数字信号处理上机实验报告.doc

实验一： 按问题要求产生采样序列，改变采样频率观察频率响应特性，Matlab源程序如下： A=1; T1=1/1000; T2=1/300; T3=1/200; a=25*pi; w0=30*pi; n=0:99; x1=A*exp(-a*n*T1).*sin(w0*n*T1); x2=A*exp(-a*n*T2).*sin(w0*n*T2); x3=A*exp(-a*n*T3).*sin(w0*n*T3); m=linspace(-pi,pi,10000); X1=x1*exp(-j*n*m);%n与m构造矩阵，xi向量与矩阵每一列相乘对应元素相加，构成DTFT后的矩阵 X2=x2*exp(-j*n*m); X3=x3*exp(-j*n*m); figure(1); subplot(3,2,1) plot(m/pi,abs(X1)); xlabel(\omega/π);ylabel(|H(e^j^\omega)|); title(采样频率为1000Hz时的幅度谱); subplot(3,2,3) plot(m/pi,abs(X2)); xlabel(\omega/π);ylabel(|H(e^j^\omega)|); title(采样频率为300Hz时的幅度谱); subplot(3,2,5) plot(m/pi,abs(X3)); xlabel(\omega/π);ylabel(|H(e^j^\omega)|); title(采样频率为200Hz时的幅度谱); subplot(3,2,2) plot(n,abs(x1)); xlabel(n);ylabel(x1(t)); title(采样频率为1000Hz时的时域波形); subplot(3,2,4) plot(n,abs(x2)); xlabel(n);ylabel(x2(t)); title(采样频率为300Hz时的时域波形); subplot(3,2,6) plot(n,abs(x3)); xlabel(n);ylabel(x3(t)); title(采样频率为200Hz时的时域波形); 波形图如下： ② 时域离散信号、 系统和系统响应分析。 Matlab源程序如下: xb=[1]; xc=ones(1,10); ha=ones(1,10); hb=[1,2.5,2.5,1]; y=conv(xb,hb); n1=0:length(y)-1; n2=0:length(hb)-1; figure(1) subplot(2,1,1); stem(n1,y,filled); xlabel(n);ylabel(y(n)); title(y(n)的时域响应); subplot(2,1,2); stem(n1,hb,filled); xlabel(n);ylabel(hb(n)); title(hb(n)的时域相应); w=linspace(-pi,pi,10000); Y=y*exp(-j*n1*w); Hb=hb*exp(-j*n2*w); figure(2) subplot(2,2,1); plot(w/pi,abs(Y)); xlabel(\omega/π);ylabel(幅度); title(DTFT[y(n)]的幅度); subplot(2,2,2); plot(w/pi,angle(Y)); xlabel(\omega/π);ylabel(相位); title(DTFT[y(n)]的相位); subplot(2,2,3); plot(w/pi,abs(Hb)); xlabel(\omega/π);ylabel(幅度); title(DTFT[Hb(n)]的幅度); subplot(2,2,4); plot(w/pi,angle(Hb)); xlabel(\omega/π);ylabel(相位); title(DTFT[Hb(n)]的相位); z=conv(xc,ha); n3=0:length(z)-1; Z=z*exp(-j*n3*w); figure(3); subplot(2,1,1); plot(w/pi,abs(Z)); xlabel(\omega/π);ylabel(幅度); title(DTFT[z(n)]的幅度); subplot(2,1,2); plot(w/pi,angle(Z)); xlabel(\omega/π);ylabel(相位); title(DTFT[z(n)]的相位); 波形图如下： 系统ha(n)对信号xc(n)的响应特性： ③ 由题目②中图2即可验证卷积定理。 实验二： 1)Matlab源程序如下 N=64; n=0:999; fs=50; T=1/fs; x1=ones(1,4); x2=[1,2,3,4,4,3,2,1];