数字信号处理上机作业

数字信号处理上机作业

ID:6172350

大小:683.71 KB

页数:21页

时间:2018-01-05

数字信号处理上机作业_第1页
数字信号处理上机作业_第2页
数字信号处理上机作业_第3页
数字信号处理上机作业_第4页
数字信号处理上机作业_第5页
资源描述:

《数字信号处理上机作业》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库

1、数字信号处理上机作业学院:电子工程学院班级:021215组员:实验一:信号、系统及系统响应1、实验目的(1)熟悉连续信号经理想采样前后的频谱变化关系,加深对时域采样定理的理解。(2)熟悉时域离散系统的时域特性。(3)利用卷积方法观察分析系统的时域特性。(4)掌握序列傅里叶变换的计算机实现方法,利用序列的傅里叶变换对连续信号、离散信号及系统响应进行频域分析。2、实验原理与方法(1)时域采样。(2)LTI系统的输入输出关系。3、实验内容及步骤(1)认真复习采样理论、离散信号与系统、线性卷积、序列的傅里叶变换及性质等有关内容,

2、阅读本实验原理与方法。(2)编制实验用主程序及相应子程序。①信号产生子程序,用于产生实验中要用到的下列信号序列:a.xa(t)=A*e^-at*sin(Ω0t)u(t)b.单位脉冲序列:xb(n)=δ(n)c.矩形序列:xc(n)=RN(n),N=10②系统单位脉冲响应序列产生子程序。本实验要用到两种FIR系统。a.ha(n)=R10(n);b.hb(n)=δ(n)+2.5δ(n-1)+2.5δ(n-2)+δ(n-3)③有限长序列线性卷积子程序用于完成两个给定长度的序列的卷积。可以直接调用MATLAB语言中的卷积函数c

3、onv。conv用于两个有限长度序列的卷积,它假定两个序列都从n=0开始。调用格式如下:y=conv(x,h)4、实验结果分析①分析采样序列的特性。a.取采样频率fs=1kHz,,即T=1ms。b.改变采样频率,fs=300Hz,观察

4、X(e^jω)

5、的变化,并做记录(打印曲线);进一步降低采样频率,fs=200Hz,观察频谱混叠是否明显存在,说明原因,并记录(打印)这时的

6、X(e^jω)

7、曲线。程序代码如下:closeall;clearall;clc;A=50;a=50*sqrt(2)*pi;m=50*sqrt(2)*

8、pi;fs1=1000;fs2=300;fs3=200;T1=1/fs1;T2=1/fs2;T3=1/fs3;N=100;n=[0:N-1];x1=A*exp(-a*n*T1).*sin(m*n*T1);x2=A*exp(-a*n*T2).*sin(m*n*T2);x3=A*exp(-a*n*T3).*sin(m*n*T3);w=linspace(-pi,pi,10000);X1=x1*exp(-j*n'*w);X2=x2*exp(-j*n'*w);X3=x3*exp(-j*n'*w);figure(1)subplot(

9、1,3,1)plot(w/pi,abs(X1));xlabel('omega/π');ylabel('

10、H(e^j^omega)

11、')title('采样频率为1000Hz时的频谱图');subplot(1,3,2)plot(w/pi,abs(X2));xlabel('omega/π');ylabel('

12、H(e^j^omega)

13、')title('采样频率为300Hz时的频谱图');subplot(1,3,3)plot(w/pi,abs(X3));xlabel('omega/π');ylabel('

14、H(e^j

15、^omega)

16、')title('采样频率为200Hz时的频谱图');②时域离散信号、系统和系统响应分析。a.观察信号xb(n)和系统hb(n)的时域和频域特性;利用线性卷积求信号xb(n)通过系统hb(n)的响应y(n),比较所求响应y(n)和hb(n)的时域及频域特性,注意它们之间有无差别,绘图说明,并用所学理论解释所得结果。b.观察系统ha(n)对信号xc(n)的响应特性。程序代码如下:closeall;clearall;clc;xbn=[1];xcn=ones(1,10);han=ones(1,10);hbn=

17、[1,2.5,2.5,1];yn=conv(xbn,hbn);n1=0:length(xbn)-1;n2=0:length(hbn)-1;subplot(2,1,1);stem(n1,xbn,'.')xlabel('n');ylabel('xb(n)')title('xb(n)的时域特性曲线')subplot(2,1,2);stem(n2,hbn,'.')xlabel('n');ylabel('hb(n)')title('hb(n)的时域特性曲线')figure(2)subplot(2,1,1);n1=[0:length

18、(xbn)-1];w=linspace(-pi,pi,10000);Xbn=xbn*exp(-j*n1'*w);subplot(2,1,1);plot(w/pi,abs(Xbn));xlabel('omega/π');ylabel('幅度')title('DTFT[xb[n]的频谱');n2=[0:length(h

当前文档最多预览五页,下载文档查看全文

此文档下载收益归作者所有

当前文档最多预览五页,下载文档查看全文
温馨提示:
1. 部分包含数学公式或PPT动画的文件,查看预览时可能会显示错乱或异常,文件下载后无此问题,请放心下载。
2. 本文档由用户上传,版权归属用户,天天文库负责整理代发布。如果您对本文档版权有争议请及时联系客服。
3. 下载前请仔细阅读文档内容,确认文档内容符合您的需求后进行下载,若出现内容与标题不符可向本站投诉处理。
4. 下载文档时可能由于网络波动等原因无法下载或下载错误,付费完成后未能成功下载的用户请联系客服处理。