2o Trabalho Controle, Amaro e Pedro Fazza
2o Trabalho Controle, Amaro e Pedro Fazza
2o Trabalho Controle, Amaro e Pedro Fazza
DE JUIZ DE FORA
Juiz de Fora
2017
Amaro da Silva Lopes Júnior
Pedro Henrique Fazza Loures
s 0, 08 s 2 s
5
10
G p (s)
s s
5
1) Código:
clear
close all
clc
G_z = c2d(G,T,'zoh')
rlocus(G_z)
figure
rlocus(G_z);
axis([0,1,-2,2])
zgrid([0.7],[])
[K,p] = rlocfind(G_z)
raio = abs(p(1));
T_s = -4*T/log(raio)
theta = atan(imag(p(1)/real(p(1))));
T_p = T*pi/theta
Command Window:
G_z =
selected_point =
0.6576 + 0.2477i
K =
5.1590
p =
0.6603 + 0.2478i
0.6603 - 0.2478i
0.0414 + 0.0000i
T_s =
5.7278
T_p =
4.3757
V =
0.8657
clear
close all
clc
w = logspace(0,3,1000);
K = 2;
G_num = 10;
G_den = conv([1 0],[1 5]);
G = tf(G_num,G_den);
KG = K*G;
[Gm,Pm,Wgm,Wpm] = margin(KG);
phi = 45 + 5;
phim = phi - Pm;
a = (1-sind(phim))/(1+sind(phim));
T = 1/(wm*sqrt(a));
w_p = 1/(a*T);
w_z = 1/T;
Kc = K/a;
Gcmf = feedback(Gc*G,1);
figure()
Gmf = feedback(G,1);
step(Gmf);
hold on
step(Gcmf);
title('Resposta ao degrau')
legend('Não compensado','Compensado')
grid on
figure
step(Gcmf) % resposta ao degrau sistema analógico compensado
hold on
step(T_z) % resposta ao degrau sistema digital compensado
legend('Contínuo compensado','Discreto compensado')
xlabel('Tempo')
ylabel('c*(t)')
Command Window:
Pm1 =
92.7579
Gp_z =
0.001935 z + 0.001872
----------------------
z^2 - 1.905 z + 0.9048
Gc_z =
1.613 z - 1.506
---------------
z - 0.9466