% 实指数序列 n 0:35; a 1.2; K 0.2; x K*a.^n;stem n,x ; xlabel 'Time index n' ;ylabel 'Amplitude' ; % 正弦序列 n 0:40; f 0.1; phase 0; A 1.5; xA*cos 2*pi*f*n - phase ; clf; % Clear old graph stemn,x ; axis [0 40 -2 2] ; grid on; title 'SinusoidalSequence' ; xlabel 'Time index n' ; ylabel 'Amplitude' ;function [y,n] seqadd x1,n1,x2,n2 % 序列相加函数 % 实现yn x1 n +x2 n % y 在包含n1和n2的n点上求序列和, % x1在位置向量n1上的第一序列 % x2 在位置向量n2上的第二序列n2可与 n1不同 % y n 的长度 n min min n1 ,min n2 : max maxn1 ,max n2 ; y1 zeros 1,length n ; y2y1; % 初始化 % 具有y的长度的x1 y1 find nmin n1 & n max n1 x1; % 具有y的长度的x2 y2 find nmin n2 & n max n2 x2; % 序列相加y y1+y2; function [y,n] seqmult x1,n1,x2,n2 % 序列相乘函数 % 实现y n x1 n +x2 n % y 在包含n1和n2的n点上求序列和, %x1 在位置向量n1上的第一序列 % x2 在位置向量n2上的第二序列 n2可与 n1不同 % y n 的长度 n min min n1 ,min n2 :max max n1 ,max n2 ; y1 zeros 1,length n ; y2y1; % 初始化 % 具有y的长度的x1 y1 find nmin n1 & n max n1 x1; % 具有y的长度的x2 y2 find nmin n2 & n max n2 x2; % 序列相加 y y1 .* y2;function [y,ny] seqshift x,nx,n0 % 实现 y n x n-n0 %n0为平移样本数ny nx + n0; % 位置向量移位y x; % 序列的值不变nx 0:5; x 0.5.^nx; n0 3; [y,ny] seqshift x,nx,n0 ; subplot 2,1,1 ; stem nx,x ; axis [0 10 0 1.2] ; xlabel 'nx' ; ylabel 'x' ; subplot 2,1,2 ; stem ny,y ; axis [0 10 0 1.2] ; xlabel 'ny' ; ylabel 'y' ; function [y,ny] seqfold x,nx % 序列翻转（对n 0折叠）子程序 % 实现 y n x -n % 将序列数值左右翻转 y fliplr x ; % 将序列位置对零位置左右翻转，故同时改变正负号 ny -fliplr nx ; 序列能量： Ex sum x .* conj x ; Ex sum abs x .^ 2 ; 例：画出信号x1 n 1.5*? n+1 - ? n-3 的波形。
