DSP實验报告（陈雪倩）.docVIP

DSP实验报告 姓名： 陈雪倩 学号： 实验一 题目一思考编程题： 1）、 两信号和波形是否相同？采样后的两序列和 是否相同？反映了什么现象？ 2）、 用help命令查看参考程序中你所不熟悉的Matlab函数的用法。 3）、 矩阵运算是Matlab中的基本操作，明确矩阵（矢量）的维数对能否正确进行计算是重要的。可通过在workspace中查看或者用size（）函数查看某矩阵（矢量）的维数。 参考程序第一段的语句：x=cos(5*pi*n1*T1)+4*sin(2*pi*n1*T1).*sin(3*pi*n1*T1) 中用计算符“*”代替“.*”，结果如何？如果进一步将参与计算的两正弦矢量sin(2*pi*n1*T1)和sin(3*pi*n1*T1)分别进行转置（提示：矢量y的转置为y’），再进行“*”，结果又如何？ 4）、 改用绿色画出。 答： 1）、和波形不相同，采样后的两序列和 相同，反映了混淆（Aliasing）现象。 2)、略 3）、若“*”代替“.*”，结果无信号图显示。如果进一步将参与计算的两正弦矢量sin(2*pi*n1*T1)和sin(3*pi*n1*T1)分别进行转置（提示：矢量y的转置为y’），再进行“*”，结果无信号图显示。 4）、如图： 二、问题二思考编程题： 1）、 该系统是否为时不变的？ 2）、 利用subplot命令，自上而下用四个小窗口分别画出，，和的波形。 答： 1）、该题根据参考程序所得图形可知：，故该系统为时变系统。 2）、图形分别如下： ： ： ： ： 总图: 三、利用卷积计算出输入信号通过FIR滤波器的输出，并观察输出信号的input-on暂态、input-off暂态和稳态部分。（计算卷积可用conv命令） 考虑下面两个滤波器，第一个的单位脉冲响应为，另一个的单位脉冲响应为；输入为周期方波，在一个周期内。 分别画出两个滤波器的输出的波形，并与书上p144例4.1.8的两幅图比较是否一致。 问题三思考编程题： 1）、 两个滤波器的DC gain分别是多少？ 2）、 响应的input-on暂态、input-off暂态和稳态段范围分别是多少？ 3）、 学习filter（）函数的使用。它是Matlab中Signal Processing Toolbox的常用函数，给定滤波器和输入，能用filter函数求出响应。改用filter代替参考程序中的conv，计算出的y1,y2结果有无变化？ 4）、 编程求教材p.138 Example 4.1.15给定的IIR系统的在输入为x(n)=u(n)-u(n-25)时输出，可以用filter或conv函数进行计算。 答: 1.两个滤波器的DC gain分别是多少？ 答：第一个为0.9866；第二个为0 2.响应的input-on暂态、input-off暂态和稳态段范围分别是多少？ 答：第一个滤波器 input-on：n=0,1,2....13;稳态：n=14,15...24;input-off:n=25,26...38 以上响应状态均以50为周期重复 第二个滤波器 input-on：n=0,1,2，3;稳态：n=4,5...24;input-off:n=25,26...28; 3.有一些变化。用filter函数计算卷积，输出结果只取input-on和稳态部分；用conv函数计算卷积，输出包含所有的暂态和稳态 4．编程求教材p.138 Example 4.1.15给定的IIR系统的在输入为x(n)=u(n)-u(n-25)时输出，可以用filter或conv函数进行计算。 答： clear; x=ones(1,25); h=(0.75).^(0:249); y=conv(x,h); n=0:199; plot(n,y(1:200),b) 实 验 二 一、观察序列频谱，观察信号通过系统后波形与频谱的变化 已知输入信号，其中，，，N可取5000点。 （1）画出的前100点波形 答：波形如下： （2）画出的DTFT频谱（） 答：波形如下： （3）某LTI系统，画出系统的幅度频响 答：波形如下： （4）求系统对的响应(可以自己编程也可利用卷积函数)。画出的波形,并与 的