资源描述：
《数字信号处理方法实验报告一二.docx》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、数字信号处理课程实验报告实验名称实验一、二学生姓名××学生学号××指导教师赵××学院信息科学与工程学院专业班级××实验一离散时间信号和系统响应一.实验目的1.熟悉连续信号经理想采样前后的频谱变化关系，加深对时域采样定理的理解2.掌握时域离散系统的时域特性3.利用卷积方法观察分析系统的时域特性4.掌握序列傅里叶变换的计算机实现方法，利用序列的傅里叶变换对离散信号及系统响应进行频域分析二、实验原理1.采样是连续信号数字化处理的第一个关键环节。对采样过程的研究不仅可以了解采样前后信号时域和频域特性的变化以及信号信息不丢失的条件，而且可以加深对离散傅里叶变换、Z变换和序列傅里叶变换之间关系
2、式的理解。对连续信号以T为采样间隔进行时域等间隔理想采样，形成采样信号：式中为周期冲激脉冲，为的理想采样。的傅里叶变换为：上式表明将连续信号采样后其频谱将变为周期的，周期为Ωs=2π/T。也即采样信号的频谱是原连续信号xa(t)的频谱Xa(jΩ)在频率轴上以Ωs为周期，周期延拓而成的。因此，若对连续信号进行采样，要保证采样频率fs≥2fm，fm为信号的最高频率，才可能由采样信号无失真地恢复出原模拟信号计算机实现时，利用计算机计算上式并不方便，因此我们利用采样序列的傅里叶变换来实现，即而为采样序列的傅里叶变换2.时域中，描述系统特性的方法是差分方程和单位脉冲响应，频域中可用系统函数描
3、述系统特性。已知输入信号，可以由差分方程、单位脉冲响应或系统函数求出系统对于该输入信号的响应。本实验仅在时域求解，对于差分方程可用Matlab中的工具箱函数filter()函数求解一个时域离散线性时不变系统的输出与输入间的关系为：可用Matlab中的工具箱函数conv()函数求解三、实验内容及步骤1.时域采样定理的验证给定模拟信号：式中。其幅频特性如图所示：选择三种采样频率Fs=1kHz,300Hz,200Hz,生成采样序列分别用序列表示。编写程序计算三个序列的幅频特性曲线，并绘图显示。观察在折叠频率附近与连续信号频谱有无明显差别，分析频谱混叠现象。源代码：n=0:63;A=444
4、.128;a=50*sqrt(2)*pi;w=50*sqrt(2)*pi;fs=1000;T1=1/(fs);x=A*exp(-a*n*T1).*sin(w*n*T1);subplot(2,3,1);stem(n,x,'.');title('Fs=1000Hz时域信号波形');subplot(2,3,4);X=fft(x);X1=[X,X];plot(0:127,abs(X1));title('Fs=1000Hz信号的幅度谱');fs2=300;T2=1/(fs2);x=A*exp(-a*n*T2).*sin(w*n*T2);subplot(2,3,2);stem(n,x,'.')
5、;title('Fs=300Hz时域信号波形');subplot(2,3,5);X=fft(x);X2=[X,X];plot(0:127,abs(X2));title('Fs=300Hz信号的幅度谱');fs3=200;T3=1/(fs3);x=A*exp(-a*n*T3).*sin(w*n*T3);subplot(2,3,3);stem(n,x,'.');title('Fs=200Hz时域信号波形');subplot(2,3,6);X=fft(x);X3=[X,X];plot(0:127,abs(X3));title('Fs=200Hz信号的幅度谱');图片：2.给定一个低通滤波
6、器的差分方程为：输入序列（1）分别求出和的系统响应，并画出其波形（2）求出系统的单位脉冲响应，画出其波形源代码：A=[1,-0.9];B=[0.05,0.05];x1n=[ones(1,8),zeros(1,50)];x2n=ones(1,100);y1n=filter(B,A,x1n);subplot(3,1,1);stem(y1n,'.');title('(1)系统对R8（n）的响应y1(n)');y2n=filter(B,A,x2n);subplot(3,1,2);stem(y2n,'.');title('(2)系统对u(n)的响应y2(n)');hn=impz(B,A,50
7、);subplot(3,1,3);stem(hn,'.');title('(3)系统单位脉冲响应h(n)');图片：3.给定系统的单位脉冲响应为用线性卷积法求分别对系统和的输出响应，并画出波形源代码：x1n=ones(1,8);h2n=[1,2.5,2.5,1,zeros(1,10)];y2n=conv(h2n,x1n);subplot(2,1,1);stem(h2n,'.b');title('系统单位脉冲响应h2(n)');subplot(2,1,2);stem(y
