离散信号与系统的Z变换分析

姓名：实验名称：离散信号与系统的Z变换分析学号：实验时间：2011.3.14评语：一.实验目的1．学会使用MATLAB表示信号的方法并绘制信号波形2．掌握使用MATLAB进行信号基本运算的指令二.实验内容1．求出下列离散序列的Z变换①1122()()cos()()kkfkk②223()(1)()()kfkkkk③3()()(5)fkkk④4()(1)()(5)fkkkkk2．已知下列单边离散序列的z变换表达式，求其对应的原离散序列。①2121()2zzFzzz②22341111()1Fzzzzz③2342(36)()zzFzz④24(1)()(1)(2)(3)zzzFzzzz3.已知离散系统的系统函数H(z)如下，请绘出系统的幅频和相频特性曲线，并说明系统的作用①122344()()()zHzzz②221()0.81zHzz4.已知描述离散系统的差分方程为：()1.2(1)0.35(2)()0.25(1)ykykykfkfk请绘出系统的幅频和相频特性曲线，并说明系统的作用。三.程序及仿真分析2（1）symskzFz=(z^2+z+1)/(z^2+z-2);%定义Z变换表达式fk=iztrans(Fz,k)%求反Z变换fk=-1/2*charfcn[0](k)+1/2*(-2)^k+1(2)symskzFz=1+1/z+1/z^2+1/z^3+1/z^4;%定义Z变换表达式fk=iztrans(Fz,k)%求反Z变换fk=charfcn[2](k)+charfcn[1](k)+charfcn[0](k)+charfcn[3](k)+charfcn[4](k)（3）symskzFz=(2*(z^2+3*z+6))/(z^4);%定义Z变换表达式fk=iztrans(Fz,k)%求反Z变换fk=12*charfcn[4](k)+6*charfcn[3](k)+2*charfcn[2](k)（4）symskzFz=(z*(z^2+z+1))/((z+1)*(z-2)*(z+3));%定义Z变换表达式fk=iztrans(Fz,k)%求反Z变换fk=-1/6*(-1)^k+7/15*2^k+7/10*(-3)^k3.(1)A=[17/61/3];B=[404];[H,w]=freqz(B,A,200,'whole');%求出对应范围内200个频率点的频率响%应样值HF=abs(H);%求出幅频特性值HX=angle(H);%求出相频特性值subplot(2,1,1);plot(w,HF)%画出幅频特性曲线title('幅频特性曲线')subplot(2,1,2);plot(w,HX)%画出相频特性曲线title('相频特性曲线')(2)A=[100.81];B=[10-1];[H,w]=freqz(B,A,200,'whole');%求出对应范围内200个频率点的频率响%应样值HF=abs(H);%求出幅频特性值HX=angle(H);%求出相频特性值subplot(2,1,1);plot(w,HF)%画出幅频特性曲线title('幅频特性曲线')subplot(2,1,2);plot(w,HX)%画出相频特性曲线title('相频特性曲线'4.A=[1-1.20.35];B=[10.250];[H,w]=freqz(B,A,200,'whole');%求出对应范围内200个频率点的频率响%应样值HF=abs(H);%求出幅频特性值HX=angle(H);%求出相频特性值subplot(2,1,1);plot(w,HF)%画出幅频特性曲线title('幅频特性曲线')subplot(2,1,2);plot(w,HX)%画出相频特性曲线title('相频特性曲线')四.实验总结