离散系统的冲激响应、卷积和实验报告.docVIP

实验项目八：离散系统的冲激响应、卷积和 实验报告 姓名:XX 学号:2013110903017 一、 实验室名称：信号与系统实验室 二、 实验项目名称：离散系统的冲激响应、卷积和 三、 实验原理： 在离散吋间情况下，最重要的是线性时不变（LTI）系统。线性吋不变系统的输入输fH关系 可通过冲激响应h[n]表不 y[n] = x[n]^h[n] = ^x[k]h[n- k] 又?=— 其中*表示卷积运算，MATLAB提供了求卷积函数conv,即 y=conv（x，h） 这里假设4川和蜘1都是W限长序列。如果4?1仅在么+乂.-1区间lAl为非零，而 AM仅在＜n＜nh+Nh-l上为非零，那么咖]就仅在 W /, ）＜n＜（nx-^-nhNx+Nh-2 N为非零值。同时也表明conv只需要在上述区间内计算咖]的Nx+Nh-\个样本值。需要 注意的是，conv并不产生存储在y中的y[n]样本的序号，而这个序号是有意义的，因为x 和h的区间都不是conv的输入区间，这样就应负责保持这些序号之间的联系。 filter命令计算线性常系数差分方程表征的W果LTI系统在某一给定输入时的输Hi。具 体地说，考虑一个满足下列差分方程的LTI系统： N M bm x[n - m ] k=0 zn=0 式;中x[?]是系统输入，y[n]是系统输!li。若x是含在IX间＜ n ＜ nx +7VX -1内x[n]的一 个MATLAB叫量，而肉量a和b包含系数七和么，刃P么 y=filter（b，a，x） 就会得出满妃下而差分方程的因果LTl系统的输出: N M ^a(k + \)y[n-k] = :^b(m + \)x[n-tn] k=() WJ=() 注意，6/（々+ 1）二人.和/?（w + l） = /?w，因为MATLAB要求所宵的向量序号都从1开始。例 如，为了表示差分方程y[n] + 2y[n -i] = x[n]~3x[n-l]表征的系统，就应该定义a=[l 2]和 b=[l 一3]。由filter产生的输出向量y包含了 yhl在与向量x中所在样本同一区间上的样 本，即nvnnx+；VA.-l,以使得W个向量x和y中都包含了 ；VY个样本。 四、 实验目的：加深对离散系统冲激响应、卷积和分析方法的理解 五、 实验器材（设备、元器件）：计算机、MATLAB软件。 六.实验步骤: 1、启动工具箱主界Kii,选中“实验二离散系统的冲激响应、卷积和”，点击按钮“进 入实验”，启动实验二的启动界面，如图4.2-1所示。 2、 仔细阅读劣验目的和实验内容，点击按钮“进入实验”，打开实验二主界面，如 图 4.2-2。 3、 求离散信号的卷积和。设定输入序列％ = “2 和）^*么…bm ], 以及它们相应的取位范围Atr = [xmin:xmax]和町= [ymin: ymax]，点出“确史”按 钮，可以得到信号x和y的卷积结果的图形。注意输入序列的长度和相应取值范围的长 度要相等，否则会得到警告消息 4、由离散信号的差分方程求系统输出。根据线性常系数差分方程： a{}y[n] + a}y[n — 1] + ? ? ? + 人y[z?-k] = b（）x[n] + hxh — 1] + …+ bmx[n -m] 和输入序列4z?1,求得输出序列）’h1。因此输入系数向量a和b的值，以及输入 信号x的值及其取值范围，点击按钮“确定”，将得到输出信号y的图形。如果没有输 入x的取值范围，将默认为x的起始叱标点为1。 实验所得结果: 实验内容（二）、MATLAB仿真 实验步骤： 1、考虑有限长信号 fl,0z?5 叫0,飾 （a）首先用解析方法计算y[n] = x[n]*h[n]o （b）接下来利用conv计算= 的非零样木值，并将这些样木存 入句量y中。构造一个标弓?內量ny，对座向量y样本的序号。用stem（ny，y） 画出这一结果。验证其结果与（a）是否一致。 2、对以下差分方程描述的系统 y[n] = 0.5x[n] + x[n — 1] + 2x[n — 2] y[n] = Q.Sy[n -1] + 2x[n] y[n] - ().8y[n —l] = 2x[zi — l] 分别利用filter计算出输入信号x[n] = nu[n]在1幺n幺4区间内的响应y[n]o 思考问题： 考虑函数conv和filter之间的关系，试利用filter函数来实现离散时间信号的 卷积。 1、利用conv计算＞，[n] = xh]*/i[n]的非零样本值 Matlab程序源代码： a=[ones (1,6)]; h=[0,l,2,3,4,5]; y=conv(a,h); m=length(y)-1; ny=0：l：m； stem(ny,y, 1 fill1);gr