《2022年信号与系统实验报告实验六离散信号与系统的Z变换分析.pdf》由会员分享，可在线阅读，更多相关《2022年信号与系统实验报告实验六离散信号与系统的Z变换分析.pdf（8页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、信号与系统实验报告实验六离散信号与系统的Z 变换分析实验六离散信号与系统的Z 变换分析学院班级姓名学号一、 实验目的1、熟悉离散信号Z变换的原理及性质2、熟悉常见信号的Z 变换3、了解正 / 反 Z 变换的 MATLAB 实现方法4、了解离散信号的Z 变换与其对应的理想抽样信号的傅氏变换与拉氏变换之间的关系5、了解利用MATLAB 实现离散系统的频率特性分析的方法二、 实验原理1、 正/反 Z 变换Z 变换分析法就是分析离散时间信号与系统的重要手段。如果以时间间隔sT对连续时间信号 f(t)进行理想抽样 ,那么 ,所得的理想抽样信号( )ft为: ( )( )*( )( )*()Tsskftf
2、 ttf ttkT理想抽样信号( )ft的双边拉普拉斯变换F(s)为: ( )( )*()()sksTstsskkFsf ttkTedtf kT e若令()( )sf kTf k,sTsze,那么( )ft的双边拉普拉斯变换F(s)为: ( )( )( )sTskz ekFsf k zF z则离散信号f(k)的 Z 变换定义为 : ( )( )kkF zf k z从上面关于Z变换的推导过程中可知,离散信号 f(k)的 Z 变换 F(z)与其对应的理想抽样信号( )ft的拉氏变换F(s)之间存在以下关系: ( )( )sTsz eFsF z精品资料 - - - 欢迎下载 - - - - - -
3、- - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 1 页，共 8 页 - - - - - - - - - - 信号与系统实验报告实验六离散信号与系统的Z 变换分析同理 ,可以推出离散信号f(k)的 Z 变换 F(z)与它对应的理想抽样信号( )ft的傅里叶变换之间的关系为()( )jTsz eFjF z如果已知信号的Z 变换 F(z),要求出所对应的原离散序列f(k),就需要进行反Z 变换 ,反 Z 变换的定义为 : 11( )( )2kf kF z zdzj?其中 ,C 为包围1( )kF z z的所有极点的闭合积分路线。在 MATLAB语言中有专门对信号进行
4、正反Z 变换的函数ztrans( ) 与 itrans( )。其调用格式分别如下 : F=ztrans( f ) 对 f(n)进行 Z 变换,其结果为F(z) F=ztrans(f,v) 对 f(n) 进行 Z 变换,其结果为F(v) F=ztrans(f,u,v) 对 f(u)进行 Z 变换 ,其结果为F(v) f=itrans ( F ) 对 F(z)进行 Z 反变换 ,其结果为f(n) f=itrans(F,u) 对 F(z)进行 Z 反变换 ,其结果为f(u) f=itrans(F,v,u )对 F(v)进行 Z 反变换 ,其结果为f(u) 注意 : 在调用函数ztran( ) 及 i
5、ztran( ) 之前,要用 syms 命令对所有需要用到的变量(如 t,u,v,w)等进行说明 ,即要将这些变量说明成符号变量。例用 MATLAB求出离散序列( )(0.5)( )kf kk的 Z 变换。MATLAB程序如下 : syms k z f=0、5k; %定义离散信号Fz=ztrans(f) %对离散信号进行Z 变换运行结果如下 : Fz = 2*z/(2*z-1) 例已知一离散信号的Z 变换式为2( )21zF zz,求出它所对应的离散信号f(k)。MATLAB程序如下 : syms k z Fz=2* z/(2*z-1); %定义 Z 变换表达式fk=iztrans(Fz,k)
6、 %求反 Z 变换运行结果如下 : fk = (1/2)k 2、离散系统的频率特性同连续系统的系统函数H(s)类似 ,离散系统的系统函数H(z)也反映了系统本身固有的特性。对于离散系统来说,如果把其系统函数H(z)中的复变量z换成j Te,那么所得的函数()j TH e就就是此离散系统的频率响应特性,即离散时间系统的频率响应为: 精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 2 页，共 8 页 - - - - - - - - - - 信号与系统实验报告实验六离散信号与系统的Z 变换分析()()()( )
7、j TjTjTjz eH eH eeH zg其中 ,()jTH e称为离散系统的幅频特性,()称为系统的相频特性。同连续系统一样,离散时间系统的幅频特性也就是频率的偶函数,相频特性也就是频率的齐函数。由于jTe就是频率的周期函数,所以离散系统的频率响应特性也就是频率的周期函数,其周期为2T,或者频率周期为2TT。 实际上 ,这就就是抽样系统的抽样频率,而其中的 T 则就是系统的抽样周期。频率响应呈现周期性就是离散系统特性区别于连续系统特性的重要特点。因此,只要分析()jTH e在2T范围内的情况 ,便可分析出系统的整个频率特性。鉴于离散系统频率响应特性的特点,为了表示方便起见,我们通常将其中的
8、T用一个变量来代替 ,即令jze代入系统函数H(z)中,用()jH e函数来表示离散系统的频率响应特性。相应地,用()jH e表示幅频特性 ,而相频特性仍用()来表示。应该特别注意 的就是 ,虽然这里的变量仍然称为频率变量,但就是它已经不就是原来意义上的角频率 概念,而实际上就是表示角度的概念。我们称之为数字角频率 。它与原来角频率的关系为:T。也就就是说 ,根据离散系统的系统函数H(z),令其中的jze,并且代入02范围内不同的频率值(实际上就是角度值),就可以逐个计算出不同频率时的响应,求出离散系统的频率响应特性。再利用离散系统频率特性的周期性特点(周期为2),求出系统的整个频率特性。离散
9、系统的幅频特性曲线与相频特性曲线能够直观地反映出系统对不同频率的输入序列的处理情况。 在函数()jH e随的变换关系中,在=0 附近 ,反映了系统对输入信号低频部分的处理情况,而在附近 ,则反映了系统对输入信号高频部分的处理情况。一般来说 ,分析离散系统频率响应特性就要绘制频率响应曲线,而这就是相当麻烦的。虽然可以通过几何矢量法来定性画出频率响应特性曲线,但一般来说这也就是很麻烦的。值得庆幸的就是 ,MATLAB为我们提供了专门用于求解离散系统频率响应的函数freqz() , 其调用格式如下 : H ,w=freqz(B,A,N)其中 ,B 与 A 分别就是表示待分析的离散系统的系统函数的分子
10、,分母多项式的向量,N 为正整数 ,返回向量 H 则包含了离散系统频率响应函数()jH e在0 范围内的 N 个频率等分点的值。向量则包含0 范围内的 N 个频率等分点。 在默认情况下 N=512。H ,w=freqqz(B,A,N,whole) 其中 ,B,A 与 N 的意义同上 ,而返回向量H 包含了频率响应函数()jH e在0 2范围内 N 个频率等分点的值。由于调用 freqz() 函数只能求出离散系统频率响应的数值,不能直接绘制曲线图,因此 ,我们可以先用freqz()函数求出系统频率响应的值,然后再利用MATLAB的 abs()与 angle()函数以及 plot()命令 ,即可绘
11、制出系统在0 或0 2范围内的幅频特性与相频特性曲线。精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 3 页，共 8 页 - - - - - - - - - - 信号与系统实验报告实验六离散信号与系统的Z 变换分析例 若离散系统的系统函数为0.5( )zH zz,请用 MATLAB计算0 频率范围内10 个等分点的频率响应()jH e的样值。MATLAB程序如下 : A=1 0; %分母多项式系数向量B=1 -0 、5; %分子多项式系数向量H,w=freqz(B,A,10) %求出对应0 范围内 10
12、个频率点的频率响应样值运行结果如下 : H = 0、5000 0、5245 + 0、1545i 0、5955 + 0、2939i 0、7061 + 0、4045i 0、8455 + 0、4755i 1、0000 + 0、5000i 1、1545 + 0、4755i 1、2939 + 0、4045i 1、4045 + 0、2939i 1、4755 + 0、1545i w = 0 0、3142 0、6283 0、9425 1、2566 1、5708 1、8850 2、1991 2、5133 2、8274 例用 MATLAB计算前面离散系统在0 2频率范围内200个频率等分点的频率响应值,并绘出相应
13、的幅频特性与相频特性曲线。MATLAB程序如下 : A=1 0; B=1 -0 、5; %H,w=freqz(B,A,200); H,w=freqz(B,A,200,whole); %求出对应0 2范围内 200 个频率点的频率响%应样值HF=abs(H); %求出幅频特性值HX=angle(H); %求出相频特性值subplot(2,1,1);plot(w,HF) %画出幅频特性曲线subplot(2,1,2);plot(w,HX) %画出相频特性曲线运行结果如下 : 精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - -
14、- - -第 4 页，共 8 页 - - - - - - - - - - 信号与系统实验报告实验六离散信号与系统的Z 变换分析三、实验内容1. 求出下列离散序列的Z 变换(1)1122( )( ) cos() ( )kkfkk解: syms k z; f = 0 、5 k * cos(k * pi / 2); Fz = ztrans(f) 结果 : Fz =4*z2/(4*z2+1) (2)223( )(1)( )( )kfkk kk解: syms k z f = k*(k - 1) * (2/3)k; Fz = ztrans(f) 结果 : Fz =24*z/(3*z-2)3 (3)3( )
15、( )(5)fkkk解: syms k z f = Heaviside(k) - Heaviside(k - 5); Fz = ztrans(f) 结果 : Fz =1+1/z+1/z2+1/z3+1/z4 (4)4( )(1)( )(5)fkk kkk解: syms k z f = k*(k-1)*(Heaviside(k) - Heaviside(k-5); Fz = ztrans(f) 精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 5 页，共 8 页 - - - - - - - - - - 信号与
16、系统实验报告实验六离散信号与系统的Z 变换分析结果 : Fz =2/z4*(z2+3*z+6) 2 已知下列单边离散序列的z 变换表达式 ,求其对应的原离散序列。(1)2121( )2zzF zzz解: syms k z Fz = (z2 + z + 1)/(z2 + z - 2); f = iztrans(Fz,k) 结果 : f =-1/2*charfcn0(k)+1/2*(-2)k+1 (2)22341111( )1Fzzzzz解: syms z k Fz = 1 + 1/z + 1/z2 + 1/z3 + 1/z4; f = iztrans(Fz,k) 结果 : f =charfcn0
17、(k)+charfcn1(k)+charfcn2(k)+charfcn3(k)+charfcn4(k) (3)2342(36)( )zzFzz解: syms z k Fz = 2*(z2 + 3*z + 6)/z4; f = iztrans(Fz) 结果 : f =12*charfcn4(k)+6*charfcn3(k)+2*charfcn2(k) (4)24(1)( )(1)(2)(3)z zzFzzzz解: syms z k Fz = z*(z2 + z + 1)/(z + 1)*(z-2)*(z+3); f = iztrans(Fz,k) 结果 : f =7/15*2k-1/6*(-1)
18、k+7/10*(-3)k 3 已知离散系统的系统函数H (z)如下 ,请绘出系统的幅频与相频特性曲线,并说明系统的作用(1)122344( )()()zH zzz解: 精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 6 页，共 8 页 - - - - - - - - - - 信号与系统实验报告实验六离散信号与系统的Z 变换分析A = 1 7/6 1/3; B = 4 4; H w = freqz(B,A,200,whole); HF = abs(H); HX = angle(H); subplot(2,1
19、,1);plot(w,HF) subplot(2,1,2);plot(w,HX) (2)221( )0.81zH zz解: A = 1 0 0 、81; B = 1 0 -1; H w = freqz(B,A,200,whole); HF = abs(H); HX = angle(H); subplot(2,1,1);plot(w,HF) subplot(2,1,2);plot(w,HX) 精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 7 页，共 8 页 - - - - - - - - - - 信号与系
20、统实验报告实验六离散信号与系统的Z 变换分析4 已知描述离散系统的差分方程为:( )1.2 (1)0.35 (2)( )0.25(1)y ky ky kf kf k请绘出系统的幅频与相频特性曲线,并说明系统的作用。解: A = 1 0 、25 0; B = 1 -1 、2 0、35; H w = freqz(B,A,200,whole); HF = abs(H); HX = angle(H); subplot(2,1,1);plot(w,HF) subplot(2,1,2);plot(w,HX) 精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 8 页，共 8 页 - - - - - - - - - -