资源描述:
《中南大学数字信号处理实验报告详》由会员上传分享,免费在线阅读,更多相关内容在应用文档-天天文库。
1、课程名称:数字信号处理姓名:贺维佳成绩:班级:电子信息1201学号:0909120925日期:2014年5月13日地点:综合实验楼指导老师:张昊目录实验一信号、系统及系统响应1.实验目的………………………………………………………3222.实验原理与方法………………………………………………33.实验内容………………………………………………………4实验步骤…………………………………………………………4程序框图…………………………………………………………64.实验结论………………………………………………………7实验代码…………………………………………
2、………………7实验截图…………………………………………………………11实验二用FFT作谱分析1.实验目的………………………………………………………142.实验原理………………………………………………………143.实验步骤………………………………………………………164.上机实验内容…………………………………………………175.实验结果………………………………………………………17实验代码…………………………………………………………18实验截图…………………………………………………………191.实验目的(1)熟悉连续信号经理想采样前后的频谱变化关
3、系,加深对时域采样定理的理解。22(2)熟悉是与离散系统的时域特性。(3)利用卷积方法观察并分析系统的时域特性。(4)掌握序列傅里叶变换的计算机实现方法,利用序列的傅里叶变换对连续信号、离散信号及系统响应进行频域分析。2.实验原理与方法(1)采样是连续信号数字处理的第一个关键环节。对一个信号Xa(t)进行理想采样过程如下:其中为的理想采样,p(t)为周期冲激脉冲,即的傅里叶变换为将p(t)代入并进行傅里叶变换其中就是采样后得到的序列X(n),即X(n)的傅里叶变换为由上两式得22(2)在数字计算机上观察分析各种序列的频域特性,通常对在[0,2π
4、]上进行M点采样来观察分析。对长度为N的有限长序列X(n),有其中一个时域离散线性非事变系统的输入/输出关系为上述积分也可以在频域实现:3.实验内容实验步骤:(1)信号产生子程序,用于产生试验中要用到的下列信号序列:a.采样信号序列:对下面连续信号:进行采样,可得到采样序列:其中A为幅度因子,a为衰减因子,是模拟角频率,T为采样间隔,这些参数在实验过程中由键盘输入,产生不同的和。b.单位脉冲序列:c.矩形序列:(1)系统单位脉冲相应序列产生子程序。本实验要用到两种FIR系统。22a.b.(1)有限长序列线性卷积子程序,用于完成两个给定长度的序列
5、的卷积。可以直接调用MATLAB语言中的卷积函数conv。Conv用于两个有限长读序列的卷积,调用格式如下:其中参数x和h是两个已复制的行向量序列。(3)完成上述子程序后编制实验主程序。(4)调用并执行试验程序,完成下属实验内容:22实验主程序框图,如下:224.实验结论实验代码:t=0:1/10:8;A=input('请输入A');a=input('请输入a');w=input('请输入w0');xa=A*exp((-a)*t).*sin(w*t);plot(t,xa);xlabel('t');ylabel('xa');title('xt波形
6、');%xn的时域n=0:1:50;T=input('请输入T');xn=A*exp((-a)*n*T).*sin(w*n*T);subplot(2,1,1);stem(n,xn,'k');title('时域信号波形');xlabel('n');ylabel('xn');%xn的傅氏变换N=50;k=-200:200;w=k*pi/100;X=DFT(xn,N);subplot(2,1,2);plot(w/pi,abs(X));title('xn的傅氏变换');xlabel('w/pi');ylabel('
7、X(jw)
8、');%hbn的时域hb
9、=[1,2.5,2.5,1];i=0:3;subplot(2,2,1);stem(i,hb,'k');axis([0,3,0,2.5]);%这个是设置坐标轴刻度范围的,前面两个是x轴设置,后面两个是y轴设置;title('hb(n)的时域序列');xlabel('n');ylabel('hb(n)');%hbn的傅氏变换k=-200:200;w=k*pi/100;HB=DFT(hb,N);22sublpot(2,2,2);plot(w/pi,abs(HB));title('hb(n)的傅氏变换');xlabel('w/pi');ylabel('
10、
11、hb(jw)
12、');%xbn的时域xb=[1,0,0,0,0,0,0,0,0,0];i=0:9;subplot(2,2,3);stem(i,xb,'