基于matlab信号与线性系统分析实验四——线性连续时间系统的分析

第一题：

1.

num=[1,0];

den=[1,32,60];

p=roots(den);

z=roots(num);

plot(real(p),imag(p),'*');hold on;

plot(real(z),imag(z),'o');grid on

稳定

-30-25-20-15-10-50 2.

num=[1,0];

den=[1,32,60];

T=0:0.1:3;

y1=impulse(num,den,T);

y2=step(num,den,T);

U=sin(T);

y3=lsim(num,den,U,T);

subplot(1,1,1);plot(T,y1);title('脉冲响应');grid on;

-2-1.5-1-0.500.51 1.52 2.53

第二题：

1.

num=[1,0];

den=[1,32,60];

T=0:0.1:3;

y1=impulse(num,den,T);

y2=step(num,den,T);

U=sin(T);

y3=lsim(num,den,U,T);

subplot(1,1,1);plot(T,y1);title('脉冲响应');grid on;

00.51 1.52 2.53-0.20

0.2

0.4

0.6

0.8

1

1.2

2.

num=[1,0];

den=[1,-1,-6];

T=0:0.1:3;

y1=impulse(num,den,T);

y2=step(num,den,T);

U=sin(T);

y3=lsim(num,den,U,T);

subplot(1,1,1);plot(T,y1);title('脉冲响应');grid on;

00.51 1.52 2.530500

1000

1500

2000

2500

3000

3500

4000

4500

5000

第三题：

num=[1,0];

den=[1,32,60];

W=0:0.1:100;

h=freqs(num,den,W);

subplot(2,1,1);plot(W,abs(h));title('幅频特性');grid on;

subplot(2,1,2);plot(W,angle(h));title('相频特性');grid on;

010********

6070809010000.01

0.02

0.03

0.04

0102030405060708090100-2-1

1

2

相频特性

例题1：

num=[1,1];

den=[1,1.3,0.8];

T=0:0.1:3;

y1=impulse(num,den,T);

y2=step(num,den,T);

U=sin(T);

y3=lsim(num,den,U,T);

subplot(221);plot(T,y1);title('脉冲响应');

subplot(222);plot(T,y2);title('阶跃响应');

subplot(223);plot(T,y3);title('输入为u=sint 的响应');

012300.2

0.40.60.8

1脉冲响应01

2300.5

1

1.5

阶跃响应

0123

00.5

1

1.5

输入为u=sint 的响应

例题2：

num=[1,1];

den=[1,1.3,0.8];

W=0:0.01:100;

h=freqs(num,den,W);

subplot(1,2,1);plot(W,abs(h));title('•ùƵÌØÐÔ')

axis([0,20,0,1.5]);

set(gca,'xtick',[0,10,20]);set(gca,'ytick',[0/sqrt(2),1.25]);grid on; subplot(1,2,2);plot(W,angle(h));title('ÏàƵÌØÐÔ')

axis([0,20,-pi/2,0.2]);

set(gca,'xtick',[0,10,20]);set(gca,'ytick',[-pi/2,-pi/4,0]);grid on;

010200

1.25

　ù••••••01020-1.5708-0.7854

•à••••••

例题3：

num=[1,0,-4];

den=[1,2,-3,2,1];

p=roots(den);

z=roots(num);

subplot(2,2,1);plot(real(p),imag(p),'*');grid on;hold on;

subplot(2,2,2);plot(real(z),imag(z),'o');grid on;