离散信号及Z变换.doc

【例1】 绘制冲激序列波形图。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.2n);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.8n);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(an*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 = z2/(z2-4)【例2】求下列函数的Z反变换。Z=sym(8*z-19)/(z2-5*z+6);x=iztrans(Z);simplify(x)ans = -19/6*charfcn0(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 onlegend(零点,极点)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 onxlabel(omega(rad/s),ylabel(幅度)titl