实验二-离散时间信号与系统的Z变换分析
26利用Z变换分析信号和系统的频域特性
26利用Z变换分析信号和系统的频域特性Z变换是一种用于分析离散时间信号和离散时间系统频域特性的数学工具。
在这篇文章中，我们将介绍Z变换的定义、性质以及如何利用Z变换分析信号和系统的频域特性。
首先，我们来定义Z变换。
对于一个离散时间信号序列x[n]，它的Z 变换被定义为：X(z) = ∑(from n=0 to ∞) x[n] * z^(-n)其中，z为复变量。
Z变换将一个离散时间信号序列映射到一个复平面上的函数。
通过计算X(z)，我们可以得到信号x[n]的频域特性。
下面，我们来讨论一些Z变换的性质。
首先是线性性质。
对于两个离散时间信号序列x1[n]和x2[n]，以及它们的Z变换X1(z)和X2(z)，以及常量a和b，则有：Z(a*x1[n]+b*x2[n])=a*X1(z)+b*X2(z)也就是说，Z变换具有线性性质。
另一个重要的性质是时移性。
对于一个离散时间信号序列x[n-k]，以及它的Z变换X(z)，则有：Z(x[n-k])=z^(-k)*X(z)这意味着在时域上的延迟会导致复平面上的旋转。
接下来，我们来讨论如何利用Z变换分析信号的频域特性。
首先，我们需要确定信号的Z变换X(z)。
对于一个给定的离散时间信号x[n]，我们可以通过对它进行Z变换的计算得到X(z)。
然后，我们可以通过观察X(z)在复平面上的分布来分析信号的频域特性。
例如，我们可以通过计算X(z)的极点和零点来确定信号的稳定性。
如果X(z)的所有极点都位于单位圆内，那么信号是稳定的；否则，信号是不稳定的。
另外，我们还可以通过计算X(z)的幅度和相位特性来分析信号的频域特性。
信号的幅度特性可以通过计算，X(z)，来获得，而信号的相位特性可以通过计算arg(X(z))来获得。
除了分析信号的频域特性，Z变换还可以用于分析离散时间系统的频域特性。
对于一个离散时间系统的冲激响应h[n]和输入信号x[n]，它们的Z变换分别为H(z)和X(z)。
数字信号处理第2章
Z变换与拉氏变换的关系:
这一关系实际上是通过 到了Z平面。
若将Z平面用极坐标表示
标表示
，代入
将S平面的函数映射
，S平面用直角坐 ，得：
上述关系表明： z 的模 r 仅与 s 的实部 相对应， z 的幅角 则仅与 s 的虚部 对应。
映射关系：
Z变换与拉氏变换的关系
0 0,2 (S平面实轴映射到Z平面的正实轴)
解：
，求它的傅立叶变换。
其幅度谱和相位谱分别为：
典型例题
❖ 例2 已知序列的傅立叶变换如下，求它的反变换。
解：
显然序列 h(n)不是绝对可和的，而是平方可和 的 ，但其依然存在傅立叶变换。 Parseval定理
典型例题
❖ 例3 证明复指数序列 x(n) e j0n 的傅立叶变换为：
证：根据序列的傅立叶反变换定义，利用冲击函 数 的性质，有：
即序列绝对可和
某的有 立些序些叶既列序变不，列换满若虽依足引然然绝入不存对频满在可 域足。和的以见的冲上后条击条例件函件。也数，不但满满，足足其平平傅方方立可可叶和和变条，换件其傅
也存在。如
、某些周期序列，见后例。
序列傅立叶变换的定义
5．常用序列的傅立叶变换
序列
(n)
傅立叶变换
1
1
典型例题
❖ 例1 已知
A形k(式k=求0,X取1(…：z),N)B，(此z) A( z )
时
为了方bi 便z i通常利用
i0
N
1 ai z i
X(z)/z的
i 1
若序列为因果序列，且N≥M，当X(z)的N个极点都是单
极点时，可以展开成以下的部分分式的形式：
则其逆Z变换为：
离散时间信号与系统的复频域分析——z变换ppt
其他特 VIP专享精彩活动
权
VIP专属身份标识
开通VIP后可以享受不定期的VIP随时随地彰显尊贵身份。
专属客服
VIP专属客服，第一时间解决你的问题。专属客服Q全部权益：1.海量精选书免费读2.热门好书抢先看3.独家精品资源4.VIP专属身份标识5.全站去广告6.名
6.6.1 数字滤波器的概念
与模拟滤波器相对应，在离散系统中 广泛应用数字滤波器。它的作用是利用离 散时间系统的特性对输入信号波形或频谱 加工处理。或者说，把输入的数字信号通 过一定的运算关系变成所需要的输出数字 信号。
数字滤波器一般可以用两种方法来实 现：一种方法是用数字硬件装配成一台专 门的设备，这种设备称为数字信号处理机； 另一种方法就是将所需要的运算编制成程 序利用计算机软件来实现。
服务特 权
共享文档下载特权
VIP用户有效期内可使用共享文档下载特权下载任意下载券标价的文档（不含付费文档和VIP专享文档），每下载一篇共享文
档消耗一个共享文档下载特权。
年VIP
月VIP
连续包月VIP
享受100次共享文档下载特权，一次 发放，全年内有效
赠每的送次VI的发P类共放型的享决特文定权档。有下效载期特为权1自个V月IP，生发效放起数每量月由发您放购一买次，赠 V不 我I送 清 的P生每 零 设效月 。 置起1自 随5每动 时次月续 取共发费 消享放， 。文一前档次往下，我载持的特续账权有号，效-自
第6章 离散时间信号与系统的复 频域分析——z变换
6.1 z 变 换 的 定 义 6.2 常 用 序 列 的 z 变 换 6.3 z 变 换 的 性 质 6.4 逆 z 变 换 6.5 离散系统的z域分析 6.6 数 字 滤 波 器 6.7 用MATLAB进行z域分析
信号与系统 z变换
信号与系统 z变换信号与系统是电子信息学科中的一门重要课程，其中的z变换是信号与系统分析的一种重要工具。
本文将介绍信号与系统中的z变换原理及应用。
一、z变换原理z变换是一种离散域的数学变换，它将离散时间序列转换为复平面上的函数。
在信号与系统中，我们常常需要对信号进行分析和处理，而z变换提供了一种方便且有效的方式。
它将离散时间序列变换为z域函数，从而可以对信号进行频域分析。
z变换的定义是：X(z) = ∑[x(n)·z^(-n)]，其中x(n)为离散时间序列，z为复变量。
通过z变换，我们可以将离散时间序列的差分方程转化为代数方程，从而简化信号与系统的分析和计算。
此外，z变换还具有线性性质和时移性质，使得我们可以方便地进行信号的加权叠加和时间偏移操作。
二、z变换的应用1. 系统的频域分析：z变换将离散时间序列转换为z域函数，可以方便地进行频域分析。
通过计算系统的传递函数在z域中的值，我们可以得到系统的频率响应，从而了解系统对不同频率信号的响应特性。
2. 系统的稳定性判断：通过z变换，可以将系统的差分方程转化为代数方程。
我们可以通过分析代数方程的根的位置，判断系统的稳定性。
如果差分方程的根都在单位圆内，说明系统是稳定的。
3. 离散时间系统的滤波设计：z变换为我们提供了一种方便的方法来设计离散时间系统的滤波器。
通过在z域中对滤波器的传递函数进行分析和调整，我们可以设计出满足特定需求的滤波器。
4. 信号的采样与重构：在数字信号处理中，我们常常需要对连续时间信号进行采样和重构。
通过z变换，我们可以将连续时间信号转换为离散时间信号，并在z域中进行处理。
然后再通过z逆变换将离散时间信号重构为连续时间信号。
5. 离散时间系统的时域分析：z变换不仅可以进行频域分析，还可以进行时域分析。
通过z变换，我们可以将离散时间系统的差分方程转换为代数方程，并通过对代数方程的分析，得到系统的时域特性。
z变换是信号与系统分析中非常重要的工具。
离散时间信号、系统和Z变换
冲激信号的强度压缩到原信号的1/2。
第二章信号分析和处理基础
设时域离散系统的输入为x(n)，经过规定的运算，系统输出序 列用 y(n) 表示。设运算关系用 T ［· ］ 表示，输出与输入之间关 系用下式表示：
y(n)=T［x(n)］
其框图如图所示:
在时域离散系统中，最重要的是线性时不变系统，因为很多物 理过程可用这类系统表征。
e j(ω +2πM)n= e jω n,
0 0
M=0,〒1,〒2…
复指数序列具有以2π为周期的周期性。
指数信号
表达式：
f (t ) K e
直流(常数) 指数衰减
指数增长
t
f (t )
0
K
a0 a0 a0
0 0
O
t
重要特性：其对时间的微分和积分仍然是指数形式。
通常把 称为指数信号的时间常数，记作，代表 信号衰减速度，具有时间的量纲。
设输入为x1(n)和x2(n)时，输出分别为y1(n)和y2(n)，即： T[ax1(n)] =3ax1(n)+4;
例2 已知f(t)的波形如图所示，试画出f(-3t-2)的波形
1.5 1 0.5 0 -4 1.5 1 0.5 0 -4 1.5 1 0.5 0 -4 1.5 1 0.5 0 -4
f(t)
-3
-2
-1
0 f(t-2)
1
2
3
4
-3
-2
-1
0
1 f(3t-2)
2
3
4
-3
-2
-1
0
1 f(-3t-2)
2
列就是时域离散信号。 实际信号处理中，这些数字序列值按顺序放在存贮器中，此时 nT 代表
实验-Z变换、零极点分析
（一）离散时间信号的Z 变换1．利用MATLAB 实现z 域的部分分式展开式MATLAB 的信号处理工具箱提供了一个对F(Z)进行部分分式展开的函数residuez()，其调用形式为：[r,p,k]=residuez(num,den)式中，num 和den 分别为F(Z)的分子多项式和分母多项式的系数向量，r 为部分分式的系数向量，p 为极点向量，k 为多项式的系数向量。
【实例1】 利用MATLAB 计算321431818)(-----+zz z z F 的部分分式展开式。
解：利用MATLAB 计算部分分式展开式程序为% 部分分式展开式的实现程序num=[18];den=[18 3 -4 -1];[r,p,k]=residuez(num,den)2．Z 变换和Z 反变换MATLAB 的符号数学工具箱提供了计算Z 变换的函数ztrans()和Z 反变换的函数iztrans （），其调用形式为)()(F iztrans f f ztrans F ==上面两式中，右端的f 和F 分别为时域表示式和z 域表示式的符号表示，可应用函数sym 来实现，其调用格式为()A sym S =式中，A 为待分析的表示式的字符串，S 为符号化的数字或变量。
【实例2】求（1）指数序列()n u a n 的Z 变换;(2)()()2a z az z F -=的Z 反变换。
解 （1）Z 变换的MATLAB 程序% Z 变换的程序实现f=sym('a^n');F=ztrans(f)程序运行结果为：z/a/(z/a-1)可以用simplify( )化简得到 :-z/(-z+a)（2）Z 反变换的MATLAB 程序% Z 反变换实现程序F=sym('a*z/(z-a)^2');f=iztrans(F)程序运行结果为f =a^n*n（二）系统函数的零极点分析1. 系统函数的零极点分布离散时间系统的系统函数定义为系统零状态响应的z 变换与激励的z 变换之比，即)()()(z X z Y z H = （3-1）如果系统函数)(z H 的有理函数表示式为：11211121)(+-+-++++++++=n n n n m m m m a z a z a z a b z b z b z b z H （3-2） 那么，在MATLAB 中系统函数的零极点就可通过函数roots 得到，也可借助函数tf2zp 得到，tf2zp 的语句格式为：[Z,P,K]=tf2zp(B,A)其中，B 与A 分别表示)(z H 的分子与分母多项式的系数向量。
Z变换及其在离散系统中的应用
Z变换及其在离散系统中的应用Z变换是一种在信号处理和控制系统中广泛应用的数学工具。
它可以将离散时间信号转换为连续复平面上的函数，从而方便进行系统分析和设计。
本文将介绍Z变换的定义及其在离散系统中的应用。
一、Z变换的定义Z变换是一种将离散时间信号转换为连续复平面上的函数的数学变换方法。
它可以将离散时间信号转换为Z域中的复函数，为信号处理和控制系统的研究提供了便利。
Z变换的定义如下：X(z) = ∑[x(n) * z^(-n)]其中，X(z)是Z变换的结果，x(n)是离散时间信号，z是复平面上的复数。
在Z变换中，z的取值是复平面上的任意一点。
通过改变z的取值，可以得到不同的频域特性。
常见的选取方式有单位圆上的点、单位圆内的点以及单位圆外的点等。
二、Z变换的性质Z变换具有许多有用的性质，这些性质对于分析和设计离散系统非常有帮助。
以下是Z变换的几个重要性质：1. 线性性质：Z变换是线性的，即对于信号的和或差的Z变换等于该信号的Z变换的和或差。
2. 移位定理：对于离散时间序列，将序列向右或向左移动n个单位时，其Z变换结果乘以z的-n次方。
3. 初值定理：序列的初始值等于其Z变换在z=1处的值。
4. 终值定理：序列的最终值等于其Z变换在z=0处的值。
5. 延时定理：将序列推迟n个单位时，其Z变换结果乘以z的n次方。
三、Z变换在离散系统中的应用Z变换在离散系统中有广泛的应用。
它可以用来描述系统的传递函数，进而进行系统的分析和设计。
以下是几个常见的应用场景：1. 系统稳定性分析：通过对系统的传递函数进行Z变换，可以得到系统在Z域中的极点分布。
通过判断极点的位置，可以判断系统的稳定性。
2. 频率响应分析：通过将频域信号进行Z变换，可以得到系统在Z 域中的频率响应。
通过分析频率响应，可以了解系统对不同频率信号的特性。
3. 离散滤波器设计：Z变换可以用来分析和设计离散滤波器。
通过对滤波器的输入输出进行Z变换，可以得到滤波器的传递函数，并基于传递函数进行进一步设计和优化。
z变换实验报告
南昌大学实验报告(信号与系统)学生姓名：肖江学号：6100210030 专业班级：电子103班实验类型：□验证□综合□设计□创新实验日期：2012/6/1 实验成绩：Z变换、离散时间系统的Z域分析一、实验目的1、学会用matlab求解z变换与逆z变换。
2、学会离散系统零极点分布图的绘制，理解离散系统零极点分布图的含义。
3、求解离散系统的频率响应特性。
二、实验说明1、一离散系统的差分方程为y(n)-by(n-1)=x(n),若激励为x(n)=a n u(n)，起始值y(-1)=0,求响应y(n)。
2、当H(s)极点位于z平面中各方框附近的位置，画出对应的h(n)波形填入方框中。
3、求系统差分方程为y(n)-1.1y(n-1)+0.7y(n-2)=x(n-1),的系统的频率响应特性。
三、实验内容1、syms n a b z%定义符号n a b zx=a^n; %定义激励信号X=ztrans(x); %计算激励信号的变换H=1/(1-b*z^(-1)); %写出系统z变换式Y=H*X; %计算输出的变换式y1=iztrans(Y); %计算输出时域表达式y=simplify(y1) %化简表达式2、pos=[26,19,18,17,24,27,13,11,9,23,28,7,4,1,22];figure,id=1; %生成新图框，子图id初始化为1for r=0.8:0.2:1.2 %极点的幅度依次为0.8,1.0,1.2for theta=0:pi/4:pi %极点的弧度依次为0，Π/4，Π/2，3Π/4,Πp=r*exp(j*theta);if theta~=0&theta~=pip=[p;p']; %如果极点不在实轴上添加一个共轭极点end[b a]=zp2tf([],p,1); %由零极点得到传递函数subplot(4,7,pos(id));[h,t]=impz(b,a,20); %计算20个点的单位样值响应stem(t,h,'k-','MarkerSize',5);%绘制单位样值响应id=id+1; %子图序号加1end%退出弧角循环end%退出幅度循环3、a=[1,-1.1,0.7];b=[0,1];subplot(2,1,1),zplane(b,a); %绘制零极点分布图subplot(2,1,2),impz(b,a); %绘制单位样值响应figure,freqz(b,a) %绘制频率特性4、a=[1,-1.1,0.6];b=[0.6,-1.1,1];subplot(2,1,1),zplane(b,a); %绘制零极点分布图subplot(2,1,2),impz(b,a); %绘制单位样值响应figure,freqz(b,a); %绘制频率响应n=[0:40]'; %生成时间点x1=sin(0.1*pi*n); %生成单频信号x2=0*n; %准备方波信号x2(mod(n,10)<5)=1; %生成周期为10的方波信号y1=filter(b,a,x1); %分别对两个信号滤波y2=filter(b,a,x2);figuresubplot(2,1,1),stem(n,x1); %绘制单频信号及其输出波形subplot(2,1,2),stem(n,y1);figuresubplot(2,1,1),stem(n,x2); %绘制方波信号及其输出波形subplot(2,1,2),stem(n,y2);四、实验结果1、y =(a^(1+n)-b^(1+n))/(a-b)2、输出波形如下3、输出波形如下：4、输出波形如下：五、实验总结通过本次实验的学习，对离散系统有了更多的了解，通过用matlab画出离散系统的零极点分布图，使我对离散系统的零极点分布与其对用的频响特性有了深刻的了解；同时对全通网络的相频失真有了进一步了解，幅度没有失真，但对不同的频率信号的相移不同，因此单频信号输入时，其输出信号的波形没有失真，只是整个波形发生了移位，但对于方波信号，由于其中包含了各种频率的信号，因此不同频率的信号相频失真不同，因此输出波形不再是方波。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
实验二 离散时间信号与系统的Z 变换分析一、 实验目的1、熟悉离散信号Z 变换的原理及性质2、熟悉常见信号的Z 变换3、了解正/反Z 变换的MATLAB 实现方法4、了解离散信号的Z 变换与其对应的理想抽样信号的傅氏变换和拉氏变换之间的关系5、了解利用MATLAB 实现离散系统的频率特性分析的方法二、 实验原理1、正/反Z 变换Z 变换分析法是分析离散时间信号与系统的重要手段。
如果以时间间隔s T 对连续时间信号f (t)进行理想抽样，那么，所得的理想抽样信号()f t δ为：()()*()()*()Ts s k f t f t t f t t kT δδδ∞=-∞==-∑理想抽样信号()f t δ的双边拉普拉斯变换F (s)为：()()*()()s ksT sts s k k F s f t t kT e dt f kT e δδ∞∞∞---∞=-∞=-∞⎡⎤=-=⎢⎥⎣⎦∑∑⎰若令()()s f kT f k = ，ssT z e= ， 那么()f t δ的双边拉普拉斯变换F (s)为：()()()sT s k z e k F s f k z F z δ∞-==-∞==∑则离散信号f (k )的Z 变换定义为：()()k k F z f k z ∞-=-∞=∑从上面关于Z 变换的推导过程中可知，离散信号f (k )的Z 变换F(z)与其对应的理想抽样信号()f t δ的拉氏变换F (s)之间存在以下关系：()()sT s z e F s F z δ==同理，可以推出离散信号f (k )的Z 变换F(z)和它对应的理想抽样信号()f t δ的傅里叶变换之间的关系为 ()()j Ts z e F j F z δωΩ==如果已知信号的Z 变换F(z)，要求出所对应的原离散序列f (k )，就需要进行反Z 变换，反Z 变换的定义为：11()()2k f k F z z dz j π-=⎰Ñ其中，C 为包围1()k F z z -的所有极点的闭合积分路线。
在MATLAB 语言中有专门对信号进行正反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( )及iztran( )之前，要用syms 命令对所有需要用到的变量（如t,u,v,w ）等进行说明，即要将这些变量说明成符号变量。
例①．用MATLAB 求出离散序列()(0.5)()kf k k ε= 的Z 变换MATLAB 程序如下：syms k zf=0.5^k; %定义离散信号Fz=ztrans(f) %对离散信号进行Z 变换 运行结果如下：Fz = 2*z/(2*z-1)例②．已知一离散信号的Z 变换式为2()21zF z z =- ，求出它所对应的离散信号f (k) MATLAB 程序如下：syms k zFz=2* z/(2*z-1); %定义Z 变换表达式 fk=iztrans(Fz,k) %求反Z 变换 运行结果如下：fk = (1/2)^k例③：求序列()(1)(4)f k k t εε=---的Z 变换.clc;clear all syms nhn=sym('kroneckerDelta(n, 1) + kroneckerDelta(n, 2)+ kroneckerDelta(n, 3)') Hz=ztrans(hn) Hz=simplify(Hz)2、离散系统的频率特性同连续系统的系统函数H (s)类似，离散系统的系统函数H (z )也反映了系统本身固有的特性。
对于离散系统来说，如果把其系统函数H (z )中的复变量z 换成s j T j ee ωΩ=（其中s T ω=Ω），那么所得的函数()j H e ω就是此离散系统的频率响应特性，即离散时间系统的频率响应为：()()()()j j j j z e H e H e e H z ωωωϕω===g其中， ()j H e ω称为离散系统的幅频特性，()ϕω称为系统的相频特性。
同连续系统一样，离散时间系统的幅频特性也是频率的偶函数，相频特性也是频率的齐函数。
由于j eω是频率的周期函数，所以离散系统的频率响应特性也是频率的周期函数，其周期为2π，或者角频率周期为2T sT πΩ=。
实际上，这就是抽样系统的抽样频率，而其中的T 则是系统的抽样周期。
频率响应呈现周期性是离散系统特性区别于连续系统特性的重要特点。
因此，只要分析()j H e ω在2ωπ≤范围内的情况，便可分析出系统的整个频率特性。
()j H e ω函数来表示离散系统的频率响应特性， ()j H e ω表示幅频特性，而相频特性仍用()ϕω来表示。
应该特别注意的是，虽然这里的变量仍然称为频率变量，但是它已经不是原来意义上的角频率概念，而实际上是表示角度的概念。
我们称之为数字频率。
它与原来角频率的关系为：s T ω=Ω。
也就是说，根据离散系统的系统函数H (z )，令其中的j z eω=，并且代入0～2π范围内不同的频率值（实际上是角度值），就可以逐个计算出不同频率时的响应，求出离散系统的频率响应特性。
再利用离散系统频率特性的周期性特点（周期为2），求出系统的整个频率特性。
离散系统的幅频特性曲线和相频特性曲线能够直观地反映出系统对不同频率的输入序列的处理情况。
在函数()j H e ω随的变换关系中，在=0附近，反映了系统对输入信号低频部分的处理情况，而在=附近，则反映了系统对输入信号高频部分的处理情况。
一般来说，分析离散系统频率响应特性就要绘制频率响应曲线，而这是相当麻烦的。
虽然可以通过几何矢量法来定性画出频率响应特性曲线，但一般来说这也是很麻烦的。
值得庆幸的是，MATLAB 为我们提供了专门用于求解离散系统频率响应的函数freqz() ，其调用格式如下：[H ,w]=freqz(B,A,N) 其中，B 和A 分别是表示待分析的离散系统的系统函数的分子，分母多项式的向量，N 为正整数，返回向量H 则包含了离散系统频率响应函数()j H e ω在0~π范围内的N 个频率等分点的值。
向量则包含0~π范围内的N 个频率等分点。
在默认情况下N=512。
[H ,w]=freqz(B,A,N,'whole') 其中，B ，A 和N 的意义同上，而返回向量H 包含了频率响应函数()j H e ω在0~2π范围内N 个频率等分点的值。
由于调用freqz()函数只能求出离散系统频率响应的数值，不能直接绘制曲线图，因此，我们可以先用freqz()函数求出系统频率响应的值，然后再利用MATLAB 的abs()和angle()函数以及plot()命令，即可绘制出系统在0~π或0~2π范围内的幅频特性和相频特性曲线。
例①．若离散系统的系统函数为0.5()z H z z-=，请用MATLAB 计算0~π频率范围内10个等分点的频率响应()j H e ω的样值。
MATLAB 程序如下：A=[1 0]; %分母多项式系数向量 B=[1 -0.5]; %分子多项式系数向量[H,w]=freqz(B,A,10) %求出对应0~π范围内10个频率点的频率响应样值 运行结果如下： 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个频率等分点的频率响应值，并绘出相应的幅频特性和相频特性曲线。
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) %画出相频特性曲线运行结果如下：运行结果分析：从该系统的幅频特性曲线可以看出，该系统呈高通特性，是一阶高通滤波器。
三、 实验内容1． 求出下列离散序列的Z 变换① 1122()()cos()()k k f k k πε= ② 223()(1)()()k f k k k k ε=- ③ 3()()(5)f k k k εε=--④ []4()(1)()(5)f k k k k k εε=---2． 已知下列单边离散序列的z 变换表达式，求其对应的原离散序列。
①2121()2z z F z z z ++=+-②22341111()1F z z z z z=++++ ③2342(36)()z z F z z ++= ④ 24(1)()(1)(2)(3)z z z F z z z z ++=+-+3. 已知离散系统的系统函数H (z)如下，请绘出系统的幅频和相频特性曲线，并说明系统的作用 ① 122344()()()z H z z z +=++ ② 221()0.81z H z z -=+4. 已知描述离散系统的差分方程为：() 1.2(1)0.35(2)()0.25(1)y k y k y k e k e k --+-=+-请绘出系统的幅频和相频特性曲线，并说明系统的作用。
四、预习要求1、熟悉正反z变换的意义及用MATLAB软件实现的方法2、熟悉离散系统的频率响应特性及用MATLAB软件实现的方法3、编写MATLAB程序五、实验报告要求1、简述实验目的及实验原理2、计算相应z变换或反z变换的理论值，并与实验结果进行比较3、记录离散系统的频率响应特性曲线，分析系统作用4、写出程序清单5、收获与建议%参考程序%三1.①clc;clear allsyms k zf1=0.5^k*cos(k*pi/2); %定义离散信号Fz1=ztrans(f1) %对离散信号进行Z变换% 实验二1.②f2=k*(k-1)*(2/3)^k; %定义离散信号Fz2=ztrans(f2) %对离散信号进行Z变换% 实验二1.③f3=sym('kroneckerDelta(n, 1) + kroneckerDelta(n, 2)+ kroneckerDelta(n, 3)')Fz3=ztrans(f3)Fz3=simplify(Fz3)% 实验二1.④f4=k*(k-1)*sym('kroneckerDelta(k, 1) + kroneckerDelta(k, 2)+ kroneckerDelta(k, 3)'); %定义离散信号Fz4=ztrans(f4)Fz4=simplify(Fz4)。