信号与系统(第三版)第五章离散时间系统的时域分析.ppt

第五章 离散时间系统的时域分析 图5.7－5　实指数序列 6. 正弦型序列的产生 例5.7－6　正弦型序列的MATLAB程序如下：　　 clear; 　　 n=0: 10; 　　 w0=pi/5; w1=pi/4; 　　 x=sin(n*w0+w1); 　　 stem(n, x); 　　 title (′正弦型序列′); 　　 line(［0, 10］, ［0, 0］); 　　 xlabel(′n′); ylabel(′x(n)′);  正弦型序列波形如图5.7－6所示。 图5.7－6　正弦型序列 7. 复指数序列的产生 例5.7－7　复指数序列的MATLAB程序如下：　　 clear; 　　 n=0: 10; 　　 delta=-0.2; w1=0.7; 　　 x=exp((delta+j*w1)*n); 　　 subplot(2, 1, 1), stem(n, real(x)); 　　 line(［0, 10］, ［0, 0］); 　　 title (′复指数序列′); 　　 ylabel(′复指数序列的实部′); 　　 subplot(2, 1, 2), stem(n, imag(x)); 　　 line(［0, 10］, ［0, 0］); 　　 ylabel(′复指数序列的虚部′); 　　 xlabel(′n′); 图5.7－7　复指数序列 8. 任意脉冲序列扩展函数　　 function［x, n］= impseq(np, ns, nf); 　　 if nsnp|nsnf|npnf; 　　 error(′输入位置参数不满足ns=np=nf′); 　　 else n=［ns: nf］; x=［(n-np)==0］; 　　 end 9. 任意阶跃序列扩展函数　　 function［x, n］= stepseq(np, ns, nf); 　　 n=［ ns: nf］; x=［(n-np)=0］; 5.7.2　序列运算的MATLAB扩展程序 1. 序列加法扩展函数　　 function ［y, n］=seqadd(x1, n1, x2, n2); 　　 n=min(min(n1), min(n2)): max(max(n1), max(n2)); 　　 y1=zeros(1, length(n)); y2=y1; 　　 y1(find((n=min(n1))(n=max(n1))==1))=x1; 　　 y2(find((n=min(n2))(n=max(n2))==1))=x2; 　　 y=y1+y2; 2. 序列乘法扩展函数　　 function ［y, n］=seqmult(x1, n1, x2, n2); 　　 n=min(min(n1), min(n2)): max(max(n1), max(n2)); 　　 y1=zeros(1, length(n)); y2=y1; 　　 y1(find((n=min(n1))(n=max(n1))==1))=x1; 　　 y2(find((n=min(n2))(n=max(n2))==1))=x2; 　　 y=y1.*y2; 3. 序列移序扩展函数　　 function［y, ny］= seqshift(x, nx, k); 　　 y=x; ny=nx+k;  4. 序列折叠扩展函数　　 function［y, ny］= seqfold(x, nx); 　　 y=fliplr(x); ny=-fliplr(nx); 5. 序列卷积扩展函数　　 function［y, ny］=convwthn(x, nx, h, nh); 　　 nys=nx(1)+nh(1); nyf=nx(end)+ nh(end); 　　 y=conv(x, h); ny=［nys: nyf］; 5.7.3　序列运算的MATLAB程序 例5.7－8　计算例5.5－1的MATLAB程序如下：　　 clear; 　　 x=［0.5, 1.5 , 1, -0.5］; n0=-1: 2; a=2; 　 ［x1, n1］=seqshift(x, n0, 2); 　　　　%移位　　 ［ym1, n］= seqmult(x1, n1, x, n0); 　%序列相乘　　 ym=a*ym1; 　%序列倍乘　　 ［ya, n］= seqadd(ym, n, x, n0); 　　%序列相加　　 subplot(4, 1, 1); stem(n0, x); ylabel(′x′); 　　 axis(［min(n), max(n), min(x