DSP1.docVIP

《数字信号处理》课程研究性学习报告 姓名 学号 同组成员 指导教师 时间 基本概念和技能学习报告 【目的】 (1) 掌握离散信号和系统时域、频域和z域分析中的基本方法和概念； (2) 学会用计算机进行离散信号和系统时域、频域和z域分析。 (3) 培养学生自主学习能力，以及发现问题、分析问题和解决问题的能力。 利用MATLAB的filter函数，求出下列系统的单位脉冲响应，并判断系统是否稳定。讨论本题所获得的结果。 【题目目的】 1. 掌握LTI系统单位脉冲响应的基本概念、系统稳定性与单位脉冲响应的关系； 2. 学会filter函数的使用方法及用filter函数计算系统单位脉冲响应； 3. 体验有限字长对系统特性的影响。 【仿真结果】 【结果分析】 由图可知，，所以此系统为稳定系统。又由于h1[k]是因果序列，由其零极点分布图可以看出，其极点都在单位圆内，故此系统为稳定系统。同理，由图可知，在时域上不绝对可和，所以此系统为不稳定系统。又由于h1[k]是因果序列，由其零极点分布图可以看出，有一极点在单位圆上，故此系统为不稳定系统。 【问题探究】 已知LTI系统的系统函数，有哪些计算系统单位脉冲响应方法，比较这些方法的优缺点。 通过逆变换可以求得此方法较为常用。 可以写出原系统方程的z变换形式，据此又可写出原系统的常系数线性差分方程，令x[k]=δ[k],解方程所得y[k]即系统的单位脉冲响应，此方法较为复杂。 b=[1]; a=[1,-1.848,0.850586]; figure(1);zplane(b,a); x=[1,zeros(1,200)]; k=0:200; y=filter(b,a,x); figure(2);stem(k,y); b=[1]; a=[1,-1.85,0.85]; figure(1);zplane(b,a); x=[1,zeros(1,200)]; k=0:200; y=filter(b,a,x); figure(2); stem(k,y); (1)利用MATLAB语句 x=firls(511,[0 0.4 0.404 1],[1 1 0 0]); 产生一个长度为512的序列x[k]，并画出该序列的幅度频谱。 (2) 已知序列，分别画出时序列y[k]的幅度频谱。解释所得到的结果。 【题目目的】 1. 学会用MATLAB函数freqz计算序列频谱； 2. 掌握序列频谱的基本特性及分析方法。 【温磬提示】 只需知道MATLAB语句 x=firls(511,[0 0.4 0.404 1],[1 1 0 0] 产生一个长度为512的序列x[k]，该序列满足 不需知道其他细节。用函数freqz计算该序列的频谱，在画幅度频谱时，建议用归一化频率。 【仿真结果】 【结果分析】 因为，故在频域内，相当于x（k）的频谱进行左右搬移，因此会出现上面的结果。 【问题探究】 有部分的计算结果可能与理论分析的结果不一致，分析出现该现象的原因，给出解决问题方法并进行仿真实验。 【仿真程序】 x=firls(511,[0 0.4 0.404 1],[1 1 0 0]); b=x; a=[1]; w=linspace(-pi,pi,512); y=freqz(b,a,w); plot(w/pi,abs(y)); %归一化频率 axis([-2,2,-1,2]); x=firls(511,[0 0.4 0.404 1],[1 1 0 0]); k=0:511; w=linspace(-pi,pi,512); y=x.*cos(0.4*pi*k); g=fft(y); plot(w/pi,abs(g)); ylabel(幅度频谱) xlabel(w0=0.4*pi) x=firls(511,[0 0.4 0.404 1],[1 1 0 0]); k=0:511; w=linspace(-pi,pi,512