实验信号与系统解.docx

实验信号与系统解.docx

《实验信号与系统解.docx》由会员分享，可在线阅读，更多相关《实验信号与系统解.docx（14页珍藏版）》请在冰豆网上搜索。

实验信号与系统解

实验四。

fsam=50;Tp=6;N=512;T=1/fsam;

t=0:

T:

Tp;

x=exp（-2*t）;

X=T*fft（x,N）;

plot（t,x）;

xlabel（'t'）;title（'时域波形'）;

w=（-N/2:

N/2-1）*（2*pi/N）*fsam;

y=1./（j*w+2）;

figure;plot（w,abs（fftshift（X））,w,abs（y）,'r-.'）;

title（'幅度谱'）;xlabel（'w'）;

legend（'计算值','理论值'）;

fsam=16;N=512;T=1/fsam;

t=-2:

T:

2;

x=[（t>=-1）&（t<=1）];

X=T*fft（x,N）;

plot（t,x）;

xlabel（'t'）;title（'时域波形'）;

w=（-N/2:

N/2-1）*（2*pi/N）*fsam;

y=2*sin（w）./w;

figure;plot（w,abs（fftshift（X））,w,abs（y）,'r-.'）;

title（'幅度谱’）;xlabel（'w'）;

legend（'计算值','理论值'）;

实验四

x=[1,2,0,-1,3,2];

kx=-2:

3;

h=[1,-1,1];

kh=0:

2;

y=conv（x,h）;

k=（kx

（1）+kh

（1））:

（kx（end）+kh（end））;

subplot（3,1,1）;stem（kx,x）;

subplot（3,1,2）;stem（kh,h）;

subplot（3,1,3）;stem（k,y）;

xlabel（'k'）;ylabel（'y'）;

k=-2:

10;

h=（k>=0）-（k>=8）;

x=（0.5.^k）.*h;

y=conv（x,h）;

ky=2*k

（1）:

2*k（end）;

subplot（3,1,1）;

stem（k,x）;

xlabel（‘k’）;

ylabel（‘x’）;

subplot（3,1,2）;

stem（k,h）;

xlabel（‘k’）;

ylabel（‘h’）;

subplot（3,1,3）;

stem（ky,y）;

xlabel（‘ky’）;

ylabel（‘y’）;

a=[1,4,3];

b=[1];

sys=tf（b,a）;

h=impulse（sys）;

g=step（sys）;

subplot（2,2,1）;

plot（h）;

title（‘系统的单位冲击响应h（t）’）;

subplot（2,2,2）;

plot（g）;

title（‘系统的单位冲击响应g（t）’）;

t=0:

0.01:

10;

h=（exp（-t）-exp（-3*t））/2;

g=（-exp（-t））/2+exp（-3*t）/6+1/3;

subplot（2,2,3）;

plot（h）;

title（‘系统的单位冲击响应的理论值h（t）’）;

subplot（2,2,4）;

plot（g）;

title（‘系统的单位冲击响应的理论值g（t）’）;

a=[1,3,2];

b=[8];

sys=tf（b,a）;

t=0:

0.01:

10;

x=exp（-2*t）-exp（-3*t）;

y=lsim（sys,x,t）;

plot（t,y）;

a=[1,3,2,0];

b=[4,1];

sys=tf（b.a）;

t=0:

0.01:

10;

x1=sin（t）;

x2=exp（-t）;

y1=lsim（sys,x1,t）;

y2=lsim（sys,x2,t）;

subplot（2,1,1）;

plot（t,y1）;

title（‘输入为sim（t）u（t）时系统的零状态响应’）;

subplot（2,1,2）;

plot（t,y2）;

title（‘输入为exp（-t）u（t）时系统的零状态响应’）;

shiyanwu

b=[1,0];

a=[1,1,25];

[H,w]=freqs（b,a）;

Hm=abs（H）;

phai=angle（H）;

Hr=real（H）;

Hi=imag（H）;

subplot（221）

plot（w,Hm）,gridon,title（'Magnituderesponse'）,xlabel（'Frequencyinrad/sec'）

subplot（223）

plot（w,phai）,gridon,title（'Phaseresponse'）,xlabel（'Frequencyinrad/sec'）

subplot（222）

plot（w,Hr）,gridon,title（'Realpartoffrequencyresponse'）,

xlabel（'Frequencyinrad/sec'）

subplot（224）

plot（w,Hi）,gridon,title（'Imaginarypartoffrequencyresponse'）,

xlabel（'Frequencyinrad/sec'）

b=[1,-1];

a=[1,1];

[H,w]=freqs（b,a）;

Hm=abs（H）;

phai=angle（H）;

Hr=real（H）;

Hi=imag（H）;

subplot（221）

plot（w,Hm）,gridon,title（'Magnituderesponse'）,xlabel（'Frequencyinrad/sec'）

subplot（223）

plot（w,phai）,gridon,title（'Phaseresponse'）,xlabel（'Frequencyinrad/sec'）

subplot（222）

plot（w,Hr）,gridon,title（'Realpartoffrequencyresponse'）,

xlabel（'Frequencyinrad/sec'）

subplot（224）

plot（w,Hi）,gridon,title（'Imaginarypartoffrequencyresponse'）,

xlabel（'Frequencyinrad/sec'）

b=[262];

a=[1,10,48,148,306,401,262];

[H,w]=freqs（b,a）;

Hm=abs（H）;

phai=angle（H）;

Hr=real（H）;

Hi=imag（H）;

subplot（221）

plot（w,Hm）,gridon,title（'Magnituderesponse'）,xlabel（'Frequencyinrad/sec'）

subplot（223）

plot（w,phai）,gridon,title（'Phaseresponse'）,xlabel（'Frequencyinrad/sec'）

subplot（222）

plot（w,Hr）,gridon,title（'Realpartoffrequencyresponse'）,

xlabel（'Frequencyinrad/sec'）

subplot（224）

plot（w,Hi）,gridon,title（'Imaginarypartoffrequencyresponse'）,

xlabel（'Frequencyinrad/sec'）

shiyanliu

b=[1,5,9,7];

a=[1,3,2];

z=roots（b）;

p=roots（a）;

sys=tf（b,a）;

pzmap（sys）

第二题第一问

b=[1,-0.5];

a1=[1,0.1+j*5];

a2=[1,0.1-j*5];

a=conv（a1,a2）;

z=roots（b）;

p=roots（a）;

sys=tf（b,a）;

pzmap（sys）

第二问

b=[1,-0.5];

a1=[1,0.1+j*5];

a2=[1,0.1-j*5];

a=conv（a1,a2）;

[H,w]=freqs（b,a）;

Hm=abs（H）;

plot（w,Hw）;

gridon;

title（‘Magnituderesponse’）;

xlabel（‘Frequencyinrad/sec’）;

第三问

b=[1,-0.5];

a1=[1,0.1+j*5];

a2=[1,0.1-j*5];

a=conv（a1,a2）;

sys=tf（b,a）;

t=-2:

0.001:

10;

h=impulse（sys）;

plot（h）;

title（‘系统的冲击响应’）;

第三题

Figure1

b=[1,5];

a=[1,6,11,6];

sys=tf（b,a）;

t=0:

0.001:

10;

x=exp（-4*t）;

y=lsim（sys,x,t）;

plot（t,y）;

title（‘系统的零状态响应’）;

figure2

[H,w]=freqs（b,a）;

Hm=abs（H）;

phai=angle（H）;

Hr=real（H）;

Hi=imag（H）;

subplot（2,2,1）;

plot（w,Hm）;

gridon;

title（‘Magnituderesponse’）;

xlabel（‘Frequencyinrad/sec’）;

subplot（2,2,2）;

plot（N,Hr）;

gridon;

title（‘Realpartoffreguencyresponse’）;

xlabel（‘Frequencyinrad/sec’）;

subplot（2,2,3）;

plot（w,phai）;

gridon;

title（‘phaseresponse’）;

xlabel（‘Frequencyinrad/sec’）;

subplot（2,2,4）;

plot（w,Hi）;

gridon;

title（‘Imaginarypartoffreguencyresponse’）;

xlabel（‘Frequencyinrad/sec’）;