下载文档原格式
实验2 离散系统的时域分析一、实验目的1、熟悉并掌握离散系统的差分方程表示法；2、加深对冲激响应和卷积分析方法的理解。
二、实验原理在时域中，离散时间系统对输入信号或者延迟信号进行运算处理，生成具有所需特性的输出信号，具体框图如下:其输入、输出关系可用以下差分方程描述：输入信号分解为冲激信号,记系统单位冲激响应，则系统响应为如下的卷积计算式：... . 当时，h[n]是有限长度的（），称系统为FIR 系统；反之，称系统为IIR 系统。
三、实验容1、用MATLAB 求系统响应1） 卷积的实现线性移不变系统可由它的单位脉冲响应来表征。
若已知了单位脉冲响应和系统激励就可通过卷积运算来求取系统响应，即)(*)()(n h n x n y 程序：x=input(‘Type in the input sequence=’); %输入xh=input(‘Type in the impulse response sequence=’); %输入h y=conv(x,h); % 对x ，h 进行卷积 N=length(y)-1; %求出N 的值n=0:1:N; %n 从0开始，间隔为1的取值取到N 为止 disp(‘output sequence=’); disp(y); %输出y stem(n,y); %画出n 为横轴，y 为纵轴的离散图xlabel(‘Time index n ’); ylable(‘Amplitude ’); % 规定x 轴y 轴的标签 输入为：x=[-2 0 1 -1 3] h=[1 2 0 -1]图形:2） 单位脉冲响应的求取线性时不变因果系统可用MA TLAB 的函数filter 来仿真 y=filter(b,a,x);其中，x和y是长度相等的两个矢量。
矢量x表示激励，矢量a，b表示系统函数形式滤波器的分子和分母系数，得到的响应为矢量y。
例如计算以下系统的单位脉冲响应y(n)+0.7y(n-1)-0.45y(y-2)-0.6y(y-3)=0.8x(n)-0.44x(n-1)+0.36x(n-2)+0.02x(n-3) 程序:N=input(‘Desired impuse response length=’);b=input(‘Type in the vector b=’);a=input(‘Type in the vector a=’);x=[1 zeros(1,N-1)];y=filter(b,a,x);k=0:1:N-1;stem(k,y);xlabel(’Time index n’); ylable(‘Amplitude’);输入：N=41b=[0.8 -0.44 0.36 0.02]a=[1 0.7 -0.45 -0.6]图形:2、以下程序中分别使用conv和filter函数计算h和x的卷积y和y1，运行程序，并分析y和y1是否有差别，为什么要使用x[n]补零后的x1来产生y1；具体分析当h[n]有i个值，x[n]有j个值，使用filter完成卷积功能，需要如何补零？程序:clf;h = [3 2 1 -2 1 0 -4 0 3]; %impulse responsex = [1 -2 3 -4 3 2 1]; %input sequencey = conv(h,x);n = 0:14;subplot(2,1,1);stem(n,y);..xlabel('Time index n'); ylabel('Amplitude');title('Output Obtained by Convolution'); grid;x1 = [x zeros(1,8)];y1 = filter(h,1,x1);subplot(2,1,2);stem(n,y1);xlabel('Time index n'); ylabel('Amplitude');title('Output Generated by Filtering'); grid;图形:因为在y=filter（b,a,x）中，利用给定矢量a和b对x中的数据进行滤波，结果放入y矢量中，y与x长度要相等，所以要使用x[n]补零后的x1来产生y1。
目 录第1章 信号与系统1.1 复习笔记1.2 课后习题详解1.3 名校考研真题详解第2章 连续系统的时域分析2.1 复习笔记2.2 课后习题详解2.3 名校考研真题详解第3章 离散系统的时域分析3.1 复习笔记3.2 课后习题详解3.3 名校考研真题详解第4章 傅里叶变换和系统的频域分析4.1 复习笔记4.2 课后习题详解4.3 名校考研真题详解第5章 连续系统的s域分析5.1 复习笔记5.2 课后习题详解5.3 名校考研真题详解第6章 离散系统的z域分析6.1 复习笔记6.2 课后习题详解6.3 名校考研真题详解第7章 系统函数7.1 复习笔记7.2 课后习题详解7.3 名校考研真题详解第8章 系统的状态变量分析8.1 复习笔记8.2 课后习题详解8.3 名校考研真题详解第1章 信号与系统1.1 复习笔记一、信号的基本概念与分类信号是载有信息的随时间变化的物理量或物理现象，其图像为信号的波形。
根据信号的不同特性，可对信号进行不同的分类：确定信号与随机信号；周期信号与非周期信号；连续时间信号与离散时间信号；实信号与复信号；能量信号与功率信号等。
二、信号的基本运算1加法和乘法f1（t）±f2（t）或f1（t）×f2（t）两信号f1（·）和f2（·）的相加、减、乘指同一时刻两信号之值对应相加、减、乘。
2．反转和平移（1）反转f（－t）f（－t）波形为f（t）波形以t＝0为轴反转。
图1-1（2）平移f（t＋t0）t0＞0，f（t＋t0）为f（t）波形在t轴上左移t0；t0＜0，f（t＋t0）为f（t）波形在t轴上右移t0。
图1-2平移的应用：在雷达系统中，雷达接收到的目标回波信号比发射信号延迟了时间t0，利用该延迟时间t0可以计算出目标与雷达之间的距离。
这里雷达接收到的目标回波信号就是延时信号。
3．尺度变换f（at）若a＞1，则f（at）波形为f（t）的波形在时间轴上压缩为原来的；若0＜a＜1，则f（at）波形为f（t）的波形在时间轴上扩展为原来的；若a＜0，则f（at）波形为f（t）的波形反转并压缩或展宽至。
第三章 离散系统的时域分析学习要求（1）会建立描述系统激励与响应关系的差分方程；（2）深刻理解系统的完全响应可分解为：零输入响应与零状态响应，自由响应与强迫响应，瞬态响应与稳态响应；（3）深刻理解系统的零输入线性与零状态线性，并根据关系求解相关的响应； （4）会根据系统差分方程和初始条件求解上述几种响应； （5）深刻理解单位冲激序列响应的意义，并会求解；（6）深刻理解系统起始状态与初始状态的区别，会根据系统差分方程和输入判断0时刻的跳变情况； （7）理解卷积和运算在信号与系统中的物理意义和运算规律，会计算信号的卷积和。
本章重点（1）系统数学模型（差分方程）的建立； （2）用时域经典法求系统的响应； （3）系统的单位冲激序列响应及其求解；（4）卷积和的定义、性质及运算，特别是()2(1)(1)y k y k f f k +-=--（k ）函数形式与其它信号的卷积； （5）利用零输入线性与零状态线性，求解系统的响应。
离散系统的描述方法：差分方程。
差分方程与微分方程求解方法在相当大的程度上一一对应。
与卷积类似，离散系统中占重要地位的卷积和（简称卷积）。
离散系统中的变换域方法包括z 变换、离散傅里叶变换以及其他多种离散正交变换（如沃尔什变换、离散余弦变换等等）。
与连续时间系统相比较，离散时间系统具有精度高，可靠性好，便于实现大规模集成的优点，借助于软件控制，可编程序控制器得到了广泛应用。
3.1 离散时间信号——序列在一些离散的瞬间才有定义的信号称为离散时间信号，简称为离散信号。
这里离散指信号的定义域——时间是离散的，它只是某些规定的值。
通常给出函数值的离散时刻的间隔是均匀的。
若此间隔为T ，以,2,1,0)(±±=n nT x 表示此离散时间信号，一般我们直接以)(n x 表示此序列。
n 表示个函数值在序列中出现的序号。
对应某序号n 的函数值称为在第n 个样点的“序列”。
一、离散信号的描述形式1、解析形式（闭合形式或闭式） 即用函数式表示。
例如：nn x )1(2)(1-= 2、图形形式即信号的波形。
线段的长短代表各序列值的大小。
有时可将他们的端点连接起来。
但是n 只有在整数值时才有意义。
二、序列的基本运算1、相加（相减）两序列)(),(n y n x 相加（减）即将两序列的对应序列相加（减）即可。
新序列)(n z 可表示为：)()()(n y n x n z ±=2、相乘（除）两序列)(),(n y n x 相乘（除）即将两序列的对应序列相乘（除）即可。
新序列)(n z 可表示为：)(/)()(),()()(n y n x n z n y n x n z ==3、延时（移位）指原序列逐项依次右移（左移）后给出的新序列)(n z ：)()(),()(m n x n z m n x n z +=-=显然，任何离散信号)(n x 都可以看成是离散冲激信号)(n δ的移位加权相加所构成，即：∑∞-∞=-=m m n m x n x )()()(δ4、反褶自变量变为原来的相反数，波形沿纵轴反转：)()(n x n z -= 5、展缩与连续信号的展缩不同，需按规律去除某些点或补足相应的零值，因此，也称为序列的“重排”。
[例3-1] 已知)(n x 波形如图3-1(a)所示，求)2/(),2(n x n x 的波形。
解：)2(n x 波形如图3-1(b)所示，这时，对应)(n x 波形中n 为奇数的各序列已不存在，波形压缩。
而)2/(n x 波形如图3-1(c)所示，图中，对应n 为奇数值各点应补入零值，n 为偶数值各点取得)(n x 波形中依 次对应的序列，因而波形扩展。
nn图3-13.2 LTI 离散系统的响应一、离散系统差分方程常系数线性差分方程的一般形式可表示为：)()1()1()(110N n y a N n y a n y a n y a N N -++-++-+- )()1()1()(110M n x b M n x b n x b n x b M M -++-++-+=-式中，b a ,是常数，未知函数)(n y 的位移阶次N 即此差分方程的阶次。
上式可简写为∑∑==-=-Mr r Nk kr n x b k n y a)()(求解常系数线性差分方程的方法一般有以下几种：（1）迭代法：包括手算逐次带入求解或利用计算机求解。
这种方法概念清楚，也比较简单，但只能得到其数值解，不能直接给出一个完整的解析式（闭式）作为解答。
（2）时域经典法：与微分方程的时域经典法类似，先分别求齐次解与特解，然后带入边界条件求待定系数。
这种方法便于从物理概念说明各响应分量之间的关系，但求解过程比较麻烦，在解决具体问题时不宜采用。
（3）分别求零输入响应与零状态响应：可以利用求齐次解的方法得到零输入响应，利用卷积和（卷积）的方法求零状态响应。
（4）变换域方法：类似与连续系统的拉氏变换，利用z 变换方法解差分方程。
本章着重介绍时域中求齐次解的方法和卷积方法。
下一章详细研究z 变换方法。
第十二章研究离散时间系统状态方程得求解。
一般差分方程的齐次方程的形式：0)(0=-∑=Nk kk n y a二、齐次解的求法（1）特征根没有重根的情况下，差分方程的齐次解为：nN N n n C C C n y ααα+++= 2211)(其中，N ααα ,,21为差分方程的特征根，即它应满足：00=∑=-Nk kn k a α。
N C C C ,,21是由边界条件决定的系数。
现举例说明求齐次方程的过程。
[例3-4] 已知差分方程为0)2()1()(=----n y n y n y ，1)2(,1)1(==y y ，求方程的解。
解：特征方程为：012=--αα求得特征解：251,25111-=+=αα 于是写出齐次解为：nn C C n y ⎪⎪⎭⎫⎝⎛-+⎪⎪⎭⎫ ⎝⎛+=251251)(21将1)2(,1)1(==y y 分别代入，得到方程组：⎪⎪⎩⎪⎪⎨⎧⎪⎪⎭⎫ ⎝⎛-+⎪⎪⎭⎫ ⎝⎛+=⎪⎪⎭⎫⎝⎛-+⎪⎪⎭⎫ ⎝⎛+=22212125125112512511C C C C可求出系数：51,5121-==C C因此：nnn y ⎪⎪⎭⎫⎝⎛--⎪⎪⎭⎫ ⎝⎛+=2515125151)( （2）有重根的情况下，齐次解的形式将略有不同，假定1α是特征方程的K 重根，那么，齐次解中相应于1α的部分将有K 项：n K K K K C n C n C n C 112211)(α++++---[例3-5] 求差分方程)()3(8y )2(12)1(6)(n x n n y n y n y =-+-+-+的齐次解。
解：特征方程为：0812623=+++ααα 即：0)2(3=+α可见，-2是此方程的三重特征根。
于是求得齐次解为：n C n C n C n y )2)(()(3221-++=（3）特征根为共轭复数时，齐次解的形式可以是等幅、增幅或衰减等形式的正（余）弦序列。
[例3-6] 求差分方程0)4()3(2)2(2)1(2)(=-+---+--n y n y n y n y n y 的齐次解。
已知边界条件1)1(=y ，0)2(=y ，1)3(=y ，1)5(=y 。
解：特征方程为：01222234=+-+-αααα0)1()1(22=+-αα特征根为：121==αα，j j -==43,αα因此：nn n j C j C C n C n y )()()1)(()(4321-+++=242321ππn jn jeC eC C n C -+++=⎪⎭⎫⎝⎛+⎪⎭⎫⎝⎛++=2sin 2cos 21ππn Q n P C n C 这里，)(,4343C C j Q C C P -=+=，利用边界条件：⎪⎪⎩⎪⎪⎨⎧++==-+==-+==++==QC C y Q C C y P C C y Q C C y 212121215)5(13)3(12)2(0)1(1 解得：0,1,1,021====Q P C C最后求得差分方程的解为：⎪⎭⎫⎝⎛+=2cos 1)(πn n y 三、特解的求法为求得特解，首先将激励函数代入方程式右端（自由项），观察自由项的函数形式来选择含有待定系数的特解函数式，将此特解函数代入方程后再求待定系数。
与微分方程的n t 和t e 形式对应，若自由项为kn 形式的函数，则特解选为：k k k D n D n D +++- 110；若自由项为n a 形式的函数，且a 不是此差分方程的特征根，则特解选为：nDa 。
[例3-7] 求差分方程)1()()1(2)(--=-+n x n x n y n y 的完全解。
其中激励函数2)(n n x =，且已知1)1(-=-y 。
解：（1）首先，求得它的齐次解为：nC )2(-；（2）将激励信号代入方程右端，得到自由项为12)1(22-=--n n n 。
特解可设为：10D n D +，代入差分方程，可得：12233])1([20101010-=-+=+-++n D D n D D n D D n D 对应项系数相等求得：91,3210==D D 因此完全解可表示为：9132)2()(++-=n C n y n（3）利用边界条件求系数：913221)1(1+--=-=-C y ; 98=C最后写出完全响应：9132)2(98)(++-=n n y n 四、系统响应的分解与连续时间系统的情况类似，线性时不变离散时间系统的完全响应也可分解为自由响应分量和强迫响应分量，或者零输入响应分量与零状态响应分量。
须指出，差分方程边界条件不一定由增序列)1(,),1(),0(-N y y y 给出，对于因果系统，常给定减序列)(,),2(),1(N y y y --- 为边界条件。
若激励信号在0=n 时接入系统，所谓零状态是指),1(-y)(,),2(N y y -- 都等于零，而不是指)(,),1(),0(N y y y 为零。
如果已知)(,),2(),1(N y y y --- ，欲求)(,),1(),0(N y y y ，可利用迭代法逐次求出。
[例3-8] 已知差分方程)(05.0)1(9.0)(n u n y n y =--，求（1）0)1(=-y ；（2）1)1(=-y 系统的全响应。
解：（1）由于激励在0=n 接入，且给定0)1(=-y ，因此，起始时刻系统处于零状态，由迭代法可求得：05.0)0(=y 。
由特征方程求得齐次解为nC )9.0(，根据自由项的形式，设特解为D ，将特解代入方程得到：(10.9)0.05;0.5D D -==因此：5.0)9.0()(+=nC n y将05.0)0(=y 代入，可求得：45.0-=C最后，写出完全响应：)(]5.0)9.0(45.0[)(n u n y n+⨯-=（2）先求零状态响应。
