（3）针对极点参数①，绘出系统的频响曲线。（必做） 自定义函数：
function y=sf4(t,w);
y=(t>=0).*[-1/14*exp(-2*t)+15/14*exp(-30*t)].*exp(-j*w*t);
程序：
w=linspace(-6*pi,6*pi,512); N=length(w); F=zeros(1,N); for k=1:N
F(k)=quadl('sf4',0,100000000000,[],[],w(k)); end figure(1); plot(w,real(F)); xlabel('\\omega'); ylabel('F(j\\omega)'); title('频响特性曲线');
2．利用filter命令求下面系统的冲激响应：
y(n)-0.7y(n-1)-0.6y(n-2)+y(n-3)=x(n)+0.5x(n-1)。（必做） A=[1,-0.7,-0.6,1]; B=[1,0.5]; n=[-20:100]; X=[n ==0];
y=filter(B,A,X); stem(n,y);
title('冲激响应'); grid;
周小社
2012.06.01