Лабораторная работа №2 (1)
.docxВиниченко Анастасия, ивт-44
Вариант № 3
Лабораторная работа №2 (1).
Исследование элементарных звеньев и построение их характеристик.
clc,clear
shag = 0.01; maxznach = 10;
k = 3;
T = 0.6;
T1 = 0.6;
T2 = 0.5; e=0.03;
w=0:shag:maxznach;
W = k./(1i*T*w+1);
U=k./((T*w).^2+1);
V=-k*T*w./((T*w).^2+1);
A=sqrt(U.^2+V.^2);
fi=atan(V./U)/pi*180;
L=20*log10(k./sqrt(1+(T*w).^2));
subplot(2,1,1)
hold on, grid on
plot(U);
title('U(w)'); xlabel('w'); ylabel('U');
subplot(2,1,2)
hold on, grid on
plot(V);
title('V(w)'); xlabel('w'); ylabel('V');
figure()
subplot(3,1,1)
hold on, grid on
plot(A);
title('A(w)'); xlabel('w'); ylabel('A');
subplot(3,1,2)
hold on, grid on
plot(log10(w),L);
title('L(w)'); xlabel('log(w)'); ylabel('L');
subplot(3,1,3)
hold on, grid on
plot(fi);
title('fi(w)'); xlabel('w'); ylabel('fi');
figure()
hold on, grid on
y=20*log10(k);
line([-2 log10(1/T)], [y y]);
line([log10(1/T) log10(10*1/T)], [y y-20]);
w=0:shag:maxznach;
L=20*log10(k./sqrt(1+(T*w).^2));
plot(log10(w),L);
xlabel('log(w)'); ylabel('L'); title('L(w) точная и асимпт')
sys=tf([k],[T 1]);
ltiview('bode',sys)
Чем больше T, тем быстрее функция переходит из начального значения в конечное («сжимается» по оси абцисс).
clc, clear
shag = 0.01; maxznach = 20;
k = 3;
T = 0.6;
T1 = 0.6;
T2 = 0.5; e=0.03;
w=0:shag:maxznach;
W = k./(1i*T*w+1);
U=k*(1-(T*w).^2)./((1-(T*w).^2).^2+(2*e*T*w).^2);
V=-2*k*e*T*w./((1-(T*w).^2).^2+(2*e*T*w).^2);
A=sqrt(U.^2+V.^2);
fi=atan(V./U)/pi*180;
w=1:shag:maxznach+1;
L=20*log10(A);
w=0:shag:maxznach;
figure()
subplot(3,1,1)
plot(w,A);
title('A(w)'); xlabel('w'); ylabel('A');
subplot(3,1,2)
plot(log10(w),L);
title('L(w)'); xlabel('log(w)'); ylabel('L');
subplot(3,1,3)
hold on
w=0:shag:1/T;
fi=-atan(2*k*e*T*w./((1-(T*w).^2)))/pi*180;
plot(w,fi);
w=1/T+shag:shag:maxznach;
fi=(-pi-atan(2*k*e*T*w./((1-(T*w).^2))))/pi*180;
plot(w,fi);
title('fi(w)'); xlabel('w'); ylabel('fi');
hold on
y=20*log10(k);
line([-2 log10(1/T)], [y y]);
line([log10(1/T) log10(10*1/T)], [y y-40]);
w=0:shag:maxznach;
L=20*log10(A);
plot(log10(w),L);
xlabel('log(w)'); ylabel('L'); title('L(w) точная и асимпт');
sys = tf ([k], [T^2 (2*e*T) 1])
ltiview('bode',sys)
Чем больше T, тем быстрее функция переходит из начального значения в конечное
(«сжимается» по оси абцисс).
clc; clear; close all;
k = 3;
T = 0.6;
T1 = 0.6;
T2 = 0.5; s=0.03;
sys = @(w)( (k*1j*w) ./ ( (T1*1j*w + 1) .* (1-T2^2*w.^2+2*T2*s*1j*w) ));
W = @(d)(sys(d));
phi = @(d)(angle(W(d)));
U = @(d)(real(W(d)));
V = @(d)(imag(W(d)));
A = @(d)(sqrt(U(d).^2 + V(d).^2));
L = @(d)(20*log10(A(d)));
lw = logspace(-3, 3, 10000);
figure()
semilogx(lw, L(lw));
title('L(jw)'); xlabel('w'); grid on; hold on;
v1 = 20*log10(k*lw(1));
vend = 20*log10(k*lw(end));
u1 = -20*log10((T1*T2^2*lw(1)^2)/k);
uend = -20*log10((T1*T2^2*lw(end)^2)/k);
plot([lw(1), lw(end)], [v1 vend]);
plot([lw(1), lw(end)], [u1 uend]);
legend(["L(w)", "асимпт. 1", "асимпт. 2"]);
figure()
phi1 = @(d)(pi/2 + d*0);
phi2 = @(d)(angle(W(1 ./ (T1*1j*d+1) )));
phi3 = @(d)(angle(W(1 ./ (1-(T2*d).^2+2*T2*s*1j*d) )));
w = 0.01:0.001:10;
subplot(4, 1, 1)
plot(w, phi1(w));
grid on; ylim([-3 3]); xlabel('w'); title('k*p');
subplot(4, 1, 2)
plot(w, phi2(w));
grid on; ylim([-3 3]); xlabel('w'); title('(T1*p+1)^-^1');
subplot(4, 1, 3)
plot(w, phi3(w));
grid on; ylim([-3 3]); xlabel('w'); title('( (T2*p)^2 + 2*T2*s*p + 1)^-^1');
subplot(4, 1, 4)
plot(w, phi(w));
grid on; ylim([-3 3]); xlabel('w'); title('system');
clc,clear
k = 3;
T = 0.6;
T1 = 0.6;
T2 = 0.5; e=0.03;
sys=tf(k)
ltiview('step',sys);
sys=tf(k,[1 0])
ltiview('step',sys);
sys=tf([k 0])
ltiview('step',sys);
sys=tf(k, [T 1])
ltiview('step',sys);
sys=tf([T k])
ltiview('step',sys);
sys=tf(k,[T2^2 2*e*T2 1])
ltiview('step',sys);
sys=tf(k,[T2^2 0 1])
ltiview('step',sys);
sys=tf(k,[T1 1]) * tf(k,[T2 1])
ltiview('step',sys);
S