《数字信号处理》上机全部源代码调试通过-完整版.pdfVIP

《数字信号处理》上机全部源代码调试通过-完整版--第1页

《数字信号处理》上机全部源代码调试通过，完整版

（高西全，第四版）

实验一

%实验1：系统响应及系统稳定性

closeall;clearall

%调用fliter解差分方程，由系统对un的响应判断稳定性

%内容1：调用filter解差分方程，由系统对u(n)的响应判断稳定性

A=[1,-0.9];B=[0.05,0.05];

x1n=[11111111zeros(1,50)];

x2n=ones(1,128);

hn=impz(B,A,58);

subplot(2,2,1);y=h(n);tstem(hn,y);

title((a)系统单位脉冲响应h(n))

y1n=filter(B,A,x1n);

subplot(2,2,2);y=y1(n);tstem(y1n,y);

title((b)系统对R8(n)的响应y1(n))

y2n=filter(B,A,x2n);

subplot(2,2,4);y=y2(n);tstem(y2n,y);

title((c)系统对u(n)的响应y2(n))

y1n=filter(B,A,x1n);

subplot(2,2,2);y=y1(n);tstem(y1n,y);

title((b)系统对R8(n)的响应y1(n))

y2n=filter(B,A,x2n);

subplot(2,2,4);y=y2(n);tstem(y2n,y);

title((c)系统对u(n)的响应y2(n))

%内容2：调用conv函数计算卷积

x1n=[11111111];%产生信号x1n=R8n

h1n=[ones(1,10)zeros(1,10)];

h2n=[12.52.51zeros(1,10)]

y21n=conv(h1n,x1n);

y22n=conv(h2n,x1n);

figure(2)

subplot(2,2,1);y=h1(n);tstem(h1n,y);

%调用函数tstem绘图

title((d)系统单位脉冲响应h1(n))

subplot(2,2,2);y=y21(n);tstem(y21n,y);

title((e)h1(n)与R8(n)的卷积y21(n))

subplot(2,2,3);y=h2(n);tstem(h2n,y);%调用函数tstem绘图

title((f)系统单位脉冲响应h2(n))

精选

《数字信号处理》上机全部源代码调试通过-完整版--第1页

《数字信号处理》上机全部源代码调试通过-完整版--第2页

subplot(2,2,4);y=y22(n);tstem(y22n,y);

title((g)h2(n)与R8(n)的卷积y22(n))

%=====================================

%内容3：谐振器分析

un=ones(1,256);%产生信号un

n=0:255;

xs