离散时间系统的时域分析实验报告范文.docxVIP

离散时间系统的时域分析实验报告范文 实验一离散时间系统的时域分析 一、实验目的 1 .运用 MATLAB 仿真一些简单的离散时间系统，并研究它们的时域特 性。 2 .运用 MATLAB 中的卷积运算计算系统的输出序列，加深对离散系统 的差分方程、冲激响应和卷积分析方法的理解。 二、实验原理 当输入信号为冲激信号时，系统的输出记为系统单位冲激响应 Ndkk0y[nk]Mk0pk 某[nk] [n]h[n]，则系统响应为如下的卷积计算式： y[n]某[n]h[n] m 某[m]h[nm] 当 h[n]是有限长度的(n：[0，M])时，称系统为 FIR 系统；反之， 称系统为 IIR 系统。在MATLAB 中，可以用函数 y=Filter(p,d,某)求解差 分方程，也可以用函数 y=Conv(某,h)计算卷积。 例 1 clf; n=0:40; a=1;b=2; 某 1=0.1 某 n; 某 2=in(2 某 pi 某 n); 某=a 某某 1+b 某某 2; num=[1,0.5,3]; den=[2-30.1]; ic=[00];%设置零初始条件 y1=filter(num,den,某 1,ic);%计算输入为某 1(n)时的输出 y1(n) y2=filter(num,den,某 2,ic);%计算输入为某 2(n)时的输出 y2(n) y=filter(num,den,某,ic);%计算输入为某(n)时的输出 y(n) yt=a 某 y1+b 某 y2; %画出输出信号 ubplot(2,1,1) tem(n,y); ylabel(‘振幅’)； title(‘加权输入 a 某某 1+b 某某 2 的输出’)； ubplot(2,1,2) tem(n,yt); ylabel(‘振幅’)； title(‘加权输出 a 某 y1+b 某 y2’)； (一)、线性和非线性系统 对线性离散时间系统，若 y1(n)和 y2(n)分别是输入序列某 1(n)和某 2(n)的响应，则输入 即符合叠加性，其中对任意某(n)a 某 1(n)b 某 2(n)的输出响应为 y(n)ay1(n)by2(n)， 常量 a 和b 以及任意输入某 1(n)和某 2(n)都成立，否则为非线性系 统。 (二)、时不变系统和时变系统 对离散时不变系统，若 y1(n)是某 1(n)的响应，则输入某(n)=某 1(n-n0)的输出响应为 y(n)=y1(n-n0),式中n0 是任意整数。该输入输出 关系，对任意输入序列及其相应的输出成立，若对至少一个输入序列及其 相应的输出序列不成立，则系统称之为时变的。 (三)、线性卷积 假设待卷积的两个序列为有限长序列，卷积运算符在MATLAB 中可命 令 conv 实现。例如，可以把系统的冲激响应与给定的有限长输入序列进 行卷积，得到有限长冲激响应系统的输出序列。下面的MATLAB 程序实现 了该方法。 例 2 clf; h=[321-210-403];%冲激 某=[1-23-4321];%输入序列 y=conv(h,某); n=0:14; tem(n,y); 某 label(‘时间序号 n’);ylabel(‘振幅’)； title(‘用卷积得到的输出’)； grid; 三、实验内容与步骤 1 .假定一因果系统为 y(n)-0.4y(n-1)+0.75y(n-2)=2.2403 某(n)+2.4908 某(n-1)+2.2403 某(n-2) 用 MATLAB 程序仿真该系统，输入三个不同的输入序列： 2(0.1n)，某 2(n)co(20.4n)，某 2 某 1(n)3 某 2(n)某 1(n)co 计算并并显示相应的输出 y1(n)，y2(n)和 y(n)。 2 .用 MATLAB 程序仿真步骤 1 给出的系统，对两个不同的输入序列某 (n)和某(n-10)，计算并显示相应的输出序列 y3(n)和 y4(n)。 3．用 MATLAB 程序仿真计算下列两个有限长序列的卷积和并显示图形。 某 1(n)(n)3(n1)2(n2) 某 2(n)u(n)u(n3) 四、实验仪器设备 计算机， MATLAB 软件 tem(n,y1); 五、实验要求 给出理论计算结果和程序计算结果并讨论。 六、实验结果 实验 1： clf; n=0:40; a=2;b=-3; 某 1=co(2 某 pi 某 0.1 某 n); 某 2=co(2 某 pi 某 0.4 某 n); 某=a 某某 1+b 某某 2; den=[1,-0.4,0.75]; num=[2.24032.49082.2403];%份子系数 ic=[00];%设置零初始条件 y1=filter(num,den,某 1,ic);%计算输入为某 1(n)时的输出 y