离散信号及Z变换.docVIP

绘制冲激序列波形图。 function y=impDT(n) y=(n==0); n=-3:3; x=impDT(n); stem(n,x,.); xlabel(n),grid on; title(单位冲激序列); axis([-3 3 -0.1 1.1]); 【例2】绘制单位阶跃序列波形图。 function y=uDT(n) y=n=0; n=-3:5; x=uDT(n); stem(n,x,fill); xlabel(n),grid on; title(单位阶跃序列); axis([-3.1 5.1 -0.1 1.1]); 【例3】画出矩形序列 的波形图。 function y=uDT(n) y=n=0; n=-3:8; x=uDT(n)-uDT(n-5); stem(n,x,fill); xlabel(n),grid on; title(矩形序列); axis([-3.1 8.1 -0.1 1.1]); 【例4】求单边指数序列波形图。 function y=uDT(n) y=n=0; n=0:10; a1=1.2;a2=-1.2;a3=0.8;a4=-0.8; x1=a1.^n;x2=a2.^n;x3=a3.^n;x4=a4.^n; subplot(221); stem(n,x1,fill); xlabel(n),grid on; title(x(n)=1.2^{n}); subplot(222); stem(n,x2,fill); xlabel(n),grid on; title(x(n)=(-1.2)^{n}); subplot(223); stem(n,x3,fill); xlabel(n),grid on; title(x(n)=0.8^{n}); subplot(224); stem(n,x4,fill); xlabel(n),grid on; title(x(n)=(-0.8)^{n}); 【例5】绘制波形图 function y=uDT(n) y=n=0; n=0:39; x=sin(n*pi/6); stem(n,x,fill); xlabel(n),grid on; title(正弦序列); axis([0 40 -1.2 1.2]); 【例6】绘制复指数序列曲线 function y=uDT(n) y=n=0; n=0:30; A=2;a=-1/10;b=pi/6; x=A*exp((a+j*b)*n); subplot(221); stem(n,real(x),.),grid on; xlabel(n),title(实部); axis([0 31 -2.1 2.1]); subplot(222); stem(n,imag(x),.),grid on; xlabel(n),title(虚部); axis([0 31 -2.1 2.1]); subplot(223); stem(n,abs(x),.),grid on; xlabel(n),title(模); axis([0 31 -0.1 2.1]); subplot(224); stem(n,angle(x),.),grid on; xlabel(n),title(相角); axis([0 31 -4 4]); 【例1】求Z变换。 （1） x=sym(a^n*cos(pi*n)); Z=ztrans(x); simplify(Z) ans = z/(z+a) （2） x=sym(2^(n-1)-(-2)^(n-1)); Z=ztrans(x); simplify(Z) ans = z^2/(z^2-4) 【例2】求下列函数的Z反变换。 Z=sym((8*z-19)/(z^2-5*z+6)); x=iztrans(Z); simplify(x) ans = -19/6*charfcn[0](n)+5*3^(n-1)+3*2^(n-1) 【例4】已知一离散因果LTI系统的系统函数为。 B=[1,0,-0.36]; A=[1,-1.52,0.68]; zplane(B,A),grid on legend(零点,极点) title(零极点分布图) 【例6】绘制系统的频率响应曲线 b=[1,-0.96,0.9028]; a=[1,-1.56,0.8109]; [H,w]=freqz(b,a,400,whole); Hm=abs(H); Hp=angle(H); subplot(2,1,1); plot(w,Hm),grid on xlabel(\omega(rad/s)),ylabel(幅度) title(离散系统幅频特性曲线) subplot(2,1,2) plot