96 lines
1.4 KiB
Matlab
96 lines
1.4 KiB
Matlab
clear; clc; close all;
|
|
k = 2;
|
|
num1 = k;
|
|
den1 = [1 3 2];
|
|
sys1 = tf(num1, den1);
|
|
figure;
|
|
step(sys1);
|
|
grid;
|
|
[y1, t1] = step(sys1);
|
|
t = 0:0.1:6;
|
|
y2 = step(sys1, t);
|
|
|
|
T=0:0.1:50;
|
|
U(1:70)=1;
|
|
U(71:125)=2;
|
|
U(126:250)=-1;
|
|
U(251:375)=1.5;
|
|
U(376:501)=0.5;
|
|
figure;
|
|
plot (T, U, '-b');
|
|
title('sinal de entrada');
|
|
xlabel('tempo (s)');
|
|
grid;
|
|
|
|
figure;
|
|
lsim(sys1, U, T);
|
|
|
|
y3 = lsim(sys1, U, T);
|
|
u3 = U';
|
|
E3 = y3-u3;
|
|
J1 = sqrt(E3'*E3);
|
|
|
|
%% Item C
|
|
|
|
num2 = 1;
|
|
den2 = [1 3 2];
|
|
sys2 = tf(num2,den2);
|
|
K = 0:0.01:2;
|
|
figure
|
|
rlocus(sys2,K);
|
|
num3 = 2;
|
|
den3 = [1 3 4];
|
|
sys3 = tf(num3,den3);
|
|
|
|
figure;
|
|
lsim(sys3, U, T);
|
|
|
|
ka = 2;
|
|
num4 = ka * num3;
|
|
den4 = [1 3 4];
|
|
sys4 = tf(num4, den4);
|
|
figure;
|
|
lsim(sys4, U, T);
|
|
|
|
y4 = lsim(sys4, U, T);
|
|
E4 = y4-u3;
|
|
J4 = sqrt(E4'*E4);
|
|
|
|
%% Item D
|
|
kd = 1.2;
|
|
num5 = kd * 10;
|
|
den5 = [1 3 12];
|
|
sys5 = tf(num5, den5);
|
|
figure;
|
|
lsim(sys5, U, T);
|
|
|
|
num6 = [1 5];
|
|
den6 = conv([1 20], [1 3 2]);
|
|
sys6 = tf(num6, den6);
|
|
figure;
|
|
k = 0:1:500;
|
|
rlocus(sys6, k);
|
|
|
|
ka = 540/2000;
|
|
num7 = 100*ka*[1 20];
|
|
den7 = conv([1 20], [1 3 2]);
|
|
den7 = den7 + 100*[0 0 1 5];
|
|
sys7 = tf(num7, den7);
|
|
figure;
|
|
lsim(sys7, U, T);
|
|
|
|
y7 = lsim(sys7, U, T);
|
|
E7 = y7-u3;
|
|
J7 = sqrt(E7'*E7);
|
|
|
|
ka = 1040/4000;
|
|
num8 = 200*ka*[1 20];
|
|
den8 = conv([1 20], [1 3 2]);
|
|
den8 = den8 + 200*[0 0 1 5];
|
|
sys8 = tf(num8, den8);
|
|
figure;
|
|
lsim(sys8, U, T);
|
|
|
|
y8 = lsim(sys8, U, T);
|
|
E8 = y8-u3;
|
|
J8 = sqrt(E8'*E8); |