基于matlab的Z变换与反Z变换.docxVIP

11n 1 1 n 《数字信号处理》  实验目的 使用 ztrans,iztrans 函数分别求出离散时间信号的 Z 变换和 Z 反 变换的结果,并用 pretty 函数进行结果美化。编写函数时养成良好的 注释习惯,有利于对函数的理解。复习 MATLAB 的基本应用，如： help,可以帮助查询相关的函数的使用方法,巩固理论知识中的离散时 间信号的传递函数与二次项式之间的转换，以及使用 zplane 函数画 出相关系统的零极点分布图,根据零极点的分布情况估计系统的滤波 特性。 （三）  程序的运行与截图 实验项目一Z变换 （1）求 x ( n ) ?[( ) 2  n  ?( ) ]u ( n ) 3  Z变换 clear all;close all;clc; syms n f=0.5^n+(1/3)^n; %定义离散信号 F=ztrans(f) %z变换 pretty(F); 运算结果 F （2） x ( n ) ?n  4  Z变换 clear all;close all;clc; syms n f=n^4; %定义离散信号 F=ztrans(f) %Z变换 pretty(F) 运算结果 （3） x ( n ) ?sin( an ?b)  Z 变换 clear all;close all;clc; syms a b n f = sin(a*n+b) %定义离散信号 F=ztrans(f) %Z 变换 pretty(F) 运算结果 2?1?1 2 ?1 ?1 实验项目二 Z 反变换 （1） X ( z ) ?  2 z ( z ?2)  2  Z 反变换 clear all;close all;clc; syms k z Fz=2*z/(z-2)^2; % 定义 Z 反变换表达式 fk=iztrans(Fz,k) %Z 反变换 pretty(fk); 运算结果 （2） X ( z ) ?  z ( z ?1) z ?2 z ?1  Z 反变换 clear all;close all;clc; syms k z Fz=z*(z-1)/(z^2+2*z+1); %定义 Z 反变换表达式 fk=iztrans(Fz,k) %Z 反变换 pretty(fk); 运算结果 f （3） X ( z ) ?  1 ?2 z  1 ?z cos ??z  ?2  Z 反变换 clear all;close all;clc; syms k z w Fz=(1+z^(-1))/(1-2*z^-1*cos(w)+z^-2); %定义Z反变换表达式 fk=iztrans(Fz,k) pretty(fk); 运算结果 %Z 反变换 3 2 3 2 实验项目三各种模型之间的变换 H ( z ) ?  10 z ( z ?1)( z ?2)  2  =  10 z z ?5z ?8z ?4 (1)clear all;close all;clc; b=[0 0 10 0];%分子的系数数组 a=[1 -5 8 -4]; %分母的系数数组 zplane(b,a)% 运算结果 使用zplane函数绘制如下系统的零极点分布图 (2)clear all;close all;clc; b=[0 0 10 0]; %分子的系数数组 a=[1 -5 8 -4]; %分母的系数数组 [r,p,c]=residuez(b,a) %使用matlab中的residuez函数，将 H ( z ) 分解成为 多个简单有理分式之和 运算结果 r = -15.0000 5.0000 10.0000 p = 2.0000 2.0000 1.0000 c = 0 (3)clear all;close all;clc; b=[0 0 10 0]; %分子的系数数组 a=[1 -5 8 -4]; %分母的系数数组 [z,p,k]=tf2zp(b,a) 运算结果 z= 0 p = 2.0000 2.0000 1.0000 k = 10 %使用tf2zp求出系统函数的零、极点和增益 (4)clear all;close all;clc; z=[1;-3];%零点,列向量 p=[2; -4];%极点,列向量 k=5; %增益 [b,a] = zp2tf(z,p,k) 原出 分子和分母的系数 运算结果 %根据求出的零、极点和增益，然后自学使用zp2tf还 (5)clear all;close all;clc; b=[0 0 10 0]; %分子的系数数组 a=[1 -5 8 -4]; %分母的系数数组 [sos,g]=tf2sos(b,a) 联形式 运算结果 sos = %使用tf2sos将系统函数分解成一系列二阶子系统的级 0 1.0000 0