信号系统上机实验二

一、函数Filter、Freqz和Freqs的使用

（a）

1、代码

x=[1,1,1,1,1,1,0,0,0,0,0];

h=[0,1,2,3,4,5];

y=filter(h,1,x);

ny=[0:10];

stem(ny,y);

xlabel('n');

xlabel('y');

2、图像

（b）

1、代码

a1=[1 -0.8 0];b1=[2 0 -1];N=4;

[h1 omega1]=freqz(b1,a1,N)

[h2 omega2]=freqz(b1,a1,N,'whole')

2、结果

 

（c）

1、代码

a=[1 3];b=3;figure;freqs(b,a);

w=linspace(0,3*pi);

h=freqs(b,a,w);

figure;plot(w,abs(h));

a1=[3 4 1];