hnu实验五连续线性时不变系统分析

 根据RLC串联电路，可获知

1.(1)(2)

yzs= dsolve('D2y+5*Dy+6*y=6*0.5*exp(-2*t)','y(0)=0,Dy(0)=0');
t1 = -10;
t2 = 10;
p  = 0.01;
t  = t1:p:t2;
b = [6];
a = [1 5 6];
subplot(311);
impulse(b,a,t1:p:t2);
xlabel('t');ylabel('u');title('单位冲激响应');grid on
subplot(312);
step(b,a,t1:p:t2);
xlabel('t');ylabel('u');title('单位阶跃响应');grid on
subplot(313);
ezplot(yzs,[-10 10]);
xlabel('t');ylabel('yzs');title('零状态响应');grid on

1.(3)

b = [6];
a1= [1 5 6];
a2= [1 4 6];
a3= [1 2 6];
a4= [1 0.8 6];
a5= [1 0.4 6];
t1= -5;t2= 5;p= 0.1;
t = t1:p:t2;
subplot(221);
impulse(b,a1,t1:p:t2);hold on
impulse(b,a2,t1:p:t2);hold on
impulse(b,a3,t1:p:t2);hold on
impulse(b,a4,t1:p:t2);hold on
impulse(b,a5,t1:p:t2);
legend("R=5","R=4","R=2","R=0.8","R=0.4");
xlabel('t');ylabel('u');title('单位冲激响应');grid on

subplot(222);
step(b,a1,t1:p:t2);hold on
step(b,a2,t1:p:t2);hold on
step(b,a3,t1:p:t2);hold on
step(b,a4,t1:p:t2);hold on
step(b,a5,t1:p:t2);
legend("R=5","R=4","R=2","R=0.8","R=0.4");
xlabel('t');ylabel('u');title('单位阶跃响应');grid on

yzs1= dsolve('D2y+5*Dy+6*y=6*0.5*exp(-2*t)','y(0)=0,Dy(0)=0');
yzs2= dsolve('D2y+4*Dy+6*y=6*0.5*exp(-2*t)','y(0)=0,Dy(0)=0');
yzs3= dsolve('D2y+2*Dy+6*y=6*0.5*exp(-2*t)','y(0)=0,Dy(0)=0');
yzs4= dsolve('D2y+0.8*Dy+6*y=6*0.5*exp(-2*t)','y(0)=0,Dy(0)=0');
yzs5= dsolve('D2y+0.4*Dy+6*y=6*0.5*exp(-2*t)','y(0)=0,Dy(0)=0');
subplot(223);
ezplot(yzs1,[-5 5]);hold on
ezplot(yzs2,[-5 5]);hold on
ezplot(yzs3,[-5 5]);hold on
ezplot(yzs4,[-5 5]);hold on
ezplot(yzs5,[-5 5]);
legend("R=5","R=4","R=2","R=0.8","R=0.4");
xlabel('t');ylabel('u');title('零状态响应');grid on

 

 

2.

b = [6];
a1 = [1 5 6];
a2= [1 4 6];
a3= [1 2 6];
a4= [1 0.8 6];
a5= [1 0.4 6];
[H1,w1] =freqs(b,a1);
Hm1= abs(H1);
phai1= angle(H1);
[H2,w2] =freqs(b,a2);
Hm2= abs(H2);
[H3,w3] =freqs(b,a3);
Hm3= abs(H3);
[H4,w4] =freqs(b,a4);
Hm4= abs(H4);
[H5,w5] =freqs(b,a5);
Hm5= abs(H5);
subplot(221)
plot(w1,Hm1);grid on;title('幅频响应');xlabel('Omega(rad/s)');
subplot(222)
plot(w1,phai1);grid on;title('相频响应');xlabel('Omega(rad/s)');
subplot(223)
plot(w1,Hm1,w2,Hm2,w3,Hm3,w4,Hm4,w5,Hm5);grid on;title('幅频响应');xlabel('Omega(rad/s)');
legend("R=5","R=4","R=2","R=0.8","R=0.4");

 

3.

默认使用输入信号为us=0.5*exp(-2t)
下式为R=5时，求得H(s)=3/(s^3+7*s^2+16s+12)

下式为R=4时，求得H(s)=3/(s^3+6s^2+14s+12)

R=2时，求得H(s)=3/(s^3+4s^2+10s+12)
R=0.8时，求得H(s)=15/(5s^3+14s^2+38s+60)
R=0.4时，求得H(s)=15/(5s^3+52s^2+114s+60)

a1 = [1 7 16 12];
b1 = [3];
a2 = [1 6 14 12];
b2 = [3];
a3 = [1 4 10 12];
b3 = [3];
a4 = [5 14 38 60];
b4 = [15];
a5 = [5 52 114 60];
b5 = [15];
subplot(321);
zplane(b1,a1);title('零极点分布-R=5');grid on
subplot(322);
zplane(b2,a2);title('零极点分布-R=4');grid on
subplot(323);
zplane(b3,a3);title('零极点分布-R=2');grid on
subplot(324);
zplane(b4,a4);title('零极点分布-R=0.8');grid on
subplot(325);
zplane(b5,a5);title('零极点分布-R=0.4');grid on