数字信号处理课程设计实验报告

《数字信号处理课程设计实验报告》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、指导老师：***院系：***班级：***姓名：***学号：***88数字信号处理实验报告26实验一信号、系统及系统响应一、实验目的：1.熟悉理想采样的性质，了解信号采样前后的频谱变化，加深对时域采样定理的理解。2.熟悉离散信号和系统的时域特性。3.熟悉线性卷积的计算编程方法：利用卷积的方法，观察，分析系统响应的时域特性。4.掌握序列傅立叶变换的计算机实验方法，利用序列的傅里叶变换对离散信号、系统及系统响应进行频域分析。二、实验代码：l信号产生子程序：代码：clc;closeall;clearall;A=444.128;a=50*pi*sqrt(2.0);W0

2、=50*pi*sqrt(2.0);n=-50:1:50;T=1/1000;Xa=A*(exp(a*n*T)).*(sin(W0*n*T));%产生理想采样信号序列subplot(3,1,1);stem(n,Xa,'r');title('理想采样信号序列');xlabel('n');ylabel('Xa');Xb=(n==0);%产生单位脉冲序列subplot(3,1,2);stem(n,Xb,'r');title('单位脉冲序列');xlabel('n');ylabel('Xb');Xc=(n>=0)-(n>=9);%产生矩形序列subplot(3,1,3)

3、;stem(n,Xc,'r');title('矩形序列');xlabel('n');ylabel('Xc');26运行结果：26l分析理想采样信号序列的特性：代码：clc;closeall;clearall;A=444.128;a=50*pi*sqrt(2.0);W0=50*pi*sqrt(2.0);n=-50:50;T1=1/1000;Xa=A*(exp(a*n*T1)).*(sin(W0*n*T1));%T1=1/1000时得Xa序列stem(n,Xa,'r');title('Xa序列');xlabel('n');ylabel('Xa');k=-25:2

4、5;X1=Xa*(exp(-j*pi/12.5)).^(n'*k);%X1=FT[Xa]pause%暂停，观察波形subplot(3,2,1);stem(k,abs(X1),'r');%Xa的幅度谱及相位谱title('Xa的幅度谱');ylabel('幅度');subplot(3,2,2);stem(k,angle(X1),'r');title('Xa的相位谱');ylabel('相位');T2=1/300;%T2=1/300时得Xb序列Xb=A*(exp(a*n*T2)).*(sin(W0*n*T2));X2=Xb*(exp(-j*pi/12.5)).^

5、(n'*k);%X2=FT[Xb]subplot(3,2,3);stem(k,abs(X2),'r');%Xb的幅度谱及相位谱title('Xb的幅度谱');ylabel('幅度');subplot(3,2,4);stem(k,angle(X2),'r');title('Xb的相位谱');ylabel('相位');T3=1/200;%T3=1/200时得Xc序列Xc=A*(exp(a*n*T3)).*(sin(W0*n*T3));X3=Xc*(exp(-j*pi/12.5)).^(n'*k);%X3=FT[Xc]subplot(3,2,5);stem(k,a

6、bs(X3),'r');%Xc的幅度谱及相位谱title('Xc的幅度谱');ylabel('幅度');subplot(3,2,6);stem(k,angle(X3),'r');title('Xc的相位谱');ylabel('相位');26运行结果：图表1图表226l离散信号、系统和系统响应的分析：代码：1)clc;closeall;clearall;n=-50:1:50;Xb=(n==0);hb=(n==0)+2.5*(n==1)+2.5*(n==2)+(n==3);Yb=conv(Xb,hb);subplot(3,1,1);stem(n,Xb,'r');

7、%产生单位脉冲序列title('单位脉冲序列');xlabel('n');ylabel('Xb');k=-25:25;X=Xb*(exp(-j*pi/12.5)).^(n'*k);%X=FT[Xb]subplot(3,1,2);stem(k,abs(X),'r');%产生Xb的幅度谱及相位谱title('Xb的幅度谱');ylabel('幅度');subplot(3,1,3);stem(k,angle(X),'r');title('Xb的相位谱');ylabel('相位');pause%暂停观察波形subplot(3,1,1);stem(n,hb,'r');

8、title('hb序列');xlabel('n');