(数字信号处理实验报告.docxVIP

实验一 信号、系统及系统响应一、 实验目的1、熟悉连续信号经理想采样前后的频谱变化关系，加深对时域采样定理的理解；2、熟悉时域离散系统的时域特性；3、利用卷积方法观察分析系统的时域特性；4、掌握序列傅立叶变换的计算机实现方法，利用序列的傅立叶变换对连续信号、离散信号及系统响应进行频域分析。二、 实验原理及方法采样是连续信号数字处理的第一个关键环节。对采样过程的研究不仅可以了解采样前后信号时域和频域特性发生变化以及信号信息不丢失的条件，而且可以加深对傅立叶变换、Z变换和序列傅立叶变换之间关系式的理解。对一个连续信号进行理想采样的过程可用下式表示： ，其中为的理想采样，p(t)为周期脉冲，即 的傅立叶变换为??上式表明为的周期延拓。其延拓周期为采样角频率（）。只有满足采样定理时，才不会发生频率混叠失真。在实验时可以用序列的傅立叶变换来计算。公式如下：离散信号和系统在时域均可用序列来表示。为了在实验中观察分析各种序列的频域特性，通常对在[0，2]上进行M点采样来观察分析。对长度为N的有限长序列x(n)，有： ?其中，，k=0,1,……M-1时域离散线性非移变系统的输入/输出关系为 上述卷积运算也可在频域实现 ?三、 实验内容及步骤1、认真复习采样理论，离散信号与系统，线性卷积，序列的傅立叶变换及性质等有关内容，阅读本实验原理与方法。2、编制实验用主程序及相应子程序。① 信号产生子程序， 用于产生实验中要用到的下列信号序列： xa(t)=Ae-at sin(Ω0t)u(t) 进行采样， 可得到采样序列 xa(n)=xa(nT)=Ae-anT sin(Ω0nT)u(n), 0≤n50其中A为幅度因子， a为衰减因子， Ω0是模拟角频率， T为采样间隔。 这些参数都要在实验过程中由键盘输入， 产生不同的xa(t)和xa(n)。 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语言中的卷积函数conv。 conv用于两个有限长度序列的卷积， 它假定两个序列都从n=0 开始。 调用格式如下： y=conv (x, h)主程序框图：3、调通并运行实验程序，完成下列实验内容：①分析采样序列的特性，产生采样信号序列，使A=444.128，a=50,。（的无失真采样频率约为1000Hz）。a.?取采样频率=1kHz，即T=1ms。观察所得采样的幅频特性||和原图中的幅频特性曲线在折叠频率附近有无明显差别。应当注意，实验中所得频谱是用序列的傅立叶变换公式求得的，所以在频率度量上存在关系：，为数字频率，为模拟频率。b.?改变采样频率，=300Hz，观察 ||的变化，并做记录（打印曲线）；进一步降低采样频率，=200Hz，观察频谱混叠是否明显存在，说明原因，并记录（打印）这时的 ||曲线。 a步实验结果如下图所示：由图形可知，当采样频率为1000Hz时，采样序列在折叠频率附近处，即=处无明显频谱混叠。b步实验结果如下图所示：?由图可知，当采样频率进一步降低时，主瓣宽度逐渐变宽，频率混叠现象也逐渐严重。存在较明显的失真现象。?②时域离散信号、系统和系统响应分析。a.? 观察信号和系统的时域和频域特性；利用线性卷积求信号通过系统的响应y(n)，比较所求响应y(n)和的时域及频域特性，注意它们之间有无差别。绘图说明，并用所学理论解释所得结果。实验结果如下图所示：?b.观察系统对信号的响应特性。利用线性卷积求系统响应y(n)，并判断y(n)图形及其非零值序列长度是否与理论结果一致，对，说出一种定性判断y(n)图形正确与否的方法。调用序列傅立叶变换数值计算子程序，求得，观察||特性曲线，定性判断结果的正确性。改变的长度，取N=5，重复该实验。注意参数变化的影响，说明变化前后的差异，并解释所得结果。实验结果如下图所示：?N=10：N=5：?欲判断结果正确与否，可以先对其进行运算，算出其卷积，再与图形对照。?当N=10时，峰值较高，且峰值很窄，变换之后图形频带主值部分比较集中；N=5时情况与之相反。?③卷积定理的验证。将实验②中的信号换为，使a=0.4，A=1，T=1，重复实验②a，打印||曲线；对主程序做简单修改，按式（10.3.9）计算，并绘出||曲线，与前面直接对y(n)进行傅立叶变换所得幅频特性曲线进行比较，验证时域卷积定理。?实验所得结果如下：四、 思考题1、在分析理想采样序列特性的实验中，采样频率不同，相应理想采样序列的傅立叶