西电-数字信号处理大作业

《西电-数字信号处理大作业》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、数字信号处理上机大作业实验一：信号、系统及系统响应(1)简述实验冃的及实验原理。1•实验目的•熟悉连续信号经理想采样前后的频谱变化关系，加深对时域采样定理的理解。•熟悉时域离散系统的时域特性。•利用卷积方法观察分析系统的吋域特性。•掌握序列傅里叶变换的计算机实现方法，利用序列的傅里叶变换对连续信号、离散信号及系统响应进行频域分析。2.实验原理与方法•时域采样。•LTI系统的输入输出关系。(2)按实验步骤附上实验过程中的信号序列、系统单位脉冲响应及系统响应序列的时域和幅频特性曲线，并对所得结果进行分析和解释。Matlab源程序如b：A二1;T1二1/1000;T2=l

2、/300;T3=1/200;a二25*pi;w0=30*pi;n=0:99;xl=A*exp(-a*n*Tl).*sin(wO*n*Tl);x2=A*exp(~a*n*T2).*sin(wO*n*T2);x3二A*exp(-a*n*T3).*sin(vO*n*T3);m=linspace(-pi,pi,10000);XI二xl*exp(-j*n'*m);%n'与m构造矩阵,xi向量与矩阵每一列相乘对应元素相加,构成DTFT后的矩阵X2=x2*exp(-j*n'*m);X3=x3*exp(-j*n,*m);figure(l)；subplot(3,2,1)plot(m

3、/pi,abs(XI));xlabel('omega/n')jylabeL('

4、H(e"j飞omega)');titleC采样频率为1000Hz时的幅度谱');subplot(3,2,3)plot(m/pi,abs(X2));xlabel('omega/n')jylabel('

5、Il(ejomega)');titleC采样频率为300Hz时的幅度谱');subplot(3,2,5)plot(m/pi,abs(X3));xlabel('omega/n')jylabel('

6、H(ejomega)');title(，采样频率为200llz时的幅度谱');subp

7、lot(3,2,2)plot(n,abs(xl));xlabel('n');ylabel('xl(t)');title(，采样频率为1000Hz时的时域波形');subplot(3,2,4)plot(n,abs(x2));xlabel('n');ylabel('x2(t)');titleC采样频率为300Hz时的时域波形');subplot(3,2,6)plot(n,abs(x3));xlabel('n，);ylabel('x3(t)');titleC采样频率为200Hz时的时域波形');波形图如下:2采样解为300Hz时純酬1.51-Q>IH0.55-0.8-0.

8、6-0.4Q200.20.40.60.81o/n0.3諏X0.1°0102030405060708090100n②时域离散信号、系统和系统响应分析。Matlab源程序如b:xb=[l];xc=ones(l,10);ha二ones(1,10);hb=[l,2.5,2.5,1];y二conv(xb,hb);nl=0:length(y)-l;n2=0:length(hb)-l;figure(1)subplot(2,1,1);stem(nl,y,'filled*);xlabel('n');ylabel('y(n)');title(，y(n)的时域响应');subplot(2

9、,1,2);stem(nl,hb,'filled，);xlabel('n');ylabel('hb(n)');title(，hb(n)的时域相应');w-linspace(-pi,pi,10000);Y二y*cxp(-j*n]，*w);Hb=hb*exp(-j*n2,*w);figure(2)subplot(2,2,1);plot(w/pi,abs(Y));xlabel('omega/n');ylabel('幅度');titleCDTFT[y(n)]的幅度');subplot(2,2,2);plot(w/pi,angle(Y));xlabel('omega/n'

10、);ylabel('相位');titleCDTFT[y(n)]的相位');subplot(2,2,3);plot(w/pi,abs(lib));xlabel('omega/n');ylabelC幅度');titleCDTFT[IIb(n)]的幅度');subplot(2,2,4);plot(w/pi,angle(Hb));xlabel('omega/n');ylabel('相位');titleCDTFT[Hb(n)]的相位');z=conv(xc,ha);n3=0:length(z)-l;Z=z*exp(-j*n3,*w);figure(3);subplot