实验4离散信号的时域分析.docVIP

实验四 离散信号的时域分析 实验目的 学习MATLAB在信号分析和系统分析中的应用； 掌握MATLAB表示离散序列的方法、序列运算的MATLAB实现； 实验原理 离散信号的表示：在MATLAB中，可以用向量（数组）来表示有限长序列，不过这样的向量并没有包含样本位置的信息，因此，需要用两个向量（数组）来表示一个有限长序列，其中一个向量表示序列的值，称为样本向量；另一个向量表示样本的位置（即序列的序号），称为位置向量。比如序列在MATLAB中需要用以下两个向量来表示： n [-2,-1,0,1,2,3]; x [2,1,-1,5,2]; 其中 向量n表示序号，向量x表示序列的值。 当不需要样本位置信息或者序列从n 0 开始时，可以只用样本向量表示。另外由于内存的限制，MATLAB无法表示无限长序列。 离散信号的基本运算 信号相加：序列相加是对应样本相加，如果两序列长度不等或者位置向量不同，则不能用算数运算符“+”直接实现相加，必须对位置向量和长度统一处理后再相加。以下M函数可以实现任意两序列的相加运算： function [y, n] sigadd f1,n1,f2,n2 % [y n] sigadd f1,n1,f2,n2 ,Add two sequences. % Inputs: % f1 ---- the first sequence % n1 ---- index vector of f1 % f2 ---- the second sequence % n2 ---- index vector of f2 % Outputs: % y ---- the output sequence % n ---- index vector of y n min n1 1 ,n2 1 :max n1 end ,n2 end ; % index vector of y n y1 zeros 1,length n ; y2 y1; % initialization y1 n n1 1 n n1 end f1; y2 n n2 1 n n2 end f2; y y1 + y2; 序列的移位：序列移位后，样本向量没有变化，只是位置向量变了。任意序列的移位可以用以下的M函数实现： function [y,n] sigshift x,m,n0 % [y,n] sigshift x,m,n0 ,result of y x n-n0 ; % Inputs: % x ---- sequence to be shifted % m ---- the index vector of x % n0 --- shift amount % Outputs: % y ---- the output sequence % n ---- the index vector of y y x; n m + n0; 序列的反折：序列反折后，样本向量倒置，位置向量倒置且变号。以下函数可以实现序列的反折运算： function [y,n] sigfold x,m % [y,n] sigfold x,n ,result of y n x -n ; % Inputs: % x ---- sequence to be folded % m ---- the index vector of x % Outputs: % y ---- the output sequence % n ---- the index vector of y y fliplr x ; n -fliplr m ; 实验内容和步骤 序列产生及表示 在MATLAB中产生有限区间上的δ（n）或δ（n-n0），画其波形。 产生两个随机序列：x1序号始于-2，结束于7，样本值为整数；x2：序号始于0，结束于15，样本值介于 -2与2之间；对产生的序列画其波形图。 参考：a. %Generate the delta sequence function [x,n] deltaN n0,n1,n2 % [x,n] deltaN n0,n1,n2 ,generate delta n-n0 ,n1 n n2;draw % its stem graph if needed % Inputs: % n0 ---- sample position % n1 ---- beginning index of x % n2 ---- end index of x % Outputs: % x ---- the delta sequence % n ---- index vector of x n