内容发布更新时间 : 2024/12/26 21:18:43星期一 下面是文章的全部内容请认真阅读。

（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