数字信号处理实验参考程序..docVIP

实验一 一、程序清单 %内容1：调用filter解差分方程，由系统对u(n)的响应判断系统稳定性 %========================================================= close all clear all A=[1,-0.9];B=[0.05,0.05]; x1n=[1 1 1 1 1 1 1 1 zeros(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=y1(n);tstem(y1n,y); title((c)系统对u(n)的响应y2(n)) %========================================================== %内容2：调用conv函数计算卷积 x1n=[1 1 1 1 1 1 1 1]; h1n=[ones(1,10) zeros(1,10)]; h2n=[1 2.5 2.5 1 zeros(1,10)]; y21n=conv(h1n,x1n); y22n=conv(h2n,x1n); figure(2) subplot(2,2,1);y=h1(n);tstem(h1n,y); 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); title((f)系统对单位脉冲响应h2(n)) subplot(2,2,4);y=y22(n);tstem(y22n,y); title((g)h2(n)与R8(n)的卷积y22(n)) %========================================================= %内容3：谐振器分析 un=ones(1,256); n=0:255; xsin=sin(0.014*n)+sin(0.4*n); A=[1,-1.8237,0.9801];B=[1/100.49,0,-1/100.49]; y31n=filter(B,A,un); y32n=filter(B,A,xsin); figure(3) subplot(2,1,1);y=y31(n);tstem(y31n,y); title((h)谐振器对u(n)的响应y31(n)) subplot(2,1,2);y=y32(n);tstem(y32n,y); title((i)谐振器对xsin的响应y32(n)) function tstem(xn,yn) n=0:length(xn)-1; stem(n, xn,’.’); xlabel(’n’);ylabel(’yn’); axis([0,n(end),min(xn),1.2*max(xn)]); 二、实验结果与波形 三、分析与讨论 1、综合来说，在时域求系统响应的方法有两种：一是通过解差分方程求得系统的输出，注意要合理选择初始条件；二是已知系统的单位脉冲响应，通过求输入信号和系统单位脉冲响应的线性卷积求得系统输出。用计算机求解时最好使用MATLAB语言进行。 2、实际中要检验系统的稳定性，方法是在输入端加入单位阶跃序列，观察输出波形，如果波形稳定在一个常数值上，则系统稳定，否则不稳定。上面第三个实验是稳定的。 3、谐振器具有对某个频率进行谐振的特性，本实验中的谐振器的谐振频率是0.4rad，因此稳定波形为sin(0.4n)。 4、如果输入信号为无限长序列，系统的单位脉冲响应是有限长序列，可以用分段线性卷积法求系统的响应，具体参考DFT一章内容。 如果信号经过低通滤波器，则信号的高频分量被滤掉，时域信号的变化减缓，在有阶跃处附近产生过渡带。因此，当输入矩形序列时，输出序列的开始和终了都产生了明显的过渡带，见第一个实验结果的波形。 实验二 一、程序清单 %时域采样理论验证程序 %Fs=1000Hz Tp=64/1000; Fs=1000;T=1/Fs; M=Tp*Fs;n=0:M-1; A=444.128;alph=pi*50*2^0.5;omega=pi*50*2^0.5