实验二z变换及其应用

合集下载

(完整版)实验二z变换及其应用

(完整版)实验二z变换及其应用

实验三z变换及其应用3.1实验目的1)加深对离散系统变换域分析——z变换的理解;2)掌握进行z变换和z反变换的基本方法,了解部分分式法在z反变换中的应用;3)掌握使用MATLAB语言进行z变换和z反变换的常用函数。

3.2实验涉及的MATLAB函数1)ztrans功能:返回无限长序列函数x(n)的z变换。

调用格式:X=ztrans(x);求无限长序列函数x(n)的z变换X(z),返回z变换的表达式。

2)iztrans功能:求函数X(z)的z反变换x(n)。

调用格式:x=iztrans(X);求函数X(z)的z反变换x(n),返回z反变换的表达式。

3)syms功能:定义多个符号对象。

调用格式:syms a b w0;把字符a,b,w0定义为基本的符号对象。

4)residuez功能:有理多项式的部分分式展开。

调用格式:[r,p,c]=residuez(b,a);把b(z)/a(z)展开成部分分式。

[b,a]=residuez(r, p, c);根据部分分式的r、p、c数组,返回有理多项式。

其中:b,a为按降幂排列的多项式的分子和分母的系数数组;r为余数数组;p为极点数组;c为无穷项多项式系数数组。

3.3实验原理1)用ztrans 子函数求无限长序列的z 变换MATLAB 提供了进行无限长序列的z 变换的子函数ztrans 。

使用时须知,该函数只给出z 变换的表达式,而没有给出收敛域。

另外,由于这一功能还不尽完善,因而有的序列的z 变换还不能求出,z 逆变换也存在同样的问题。

例1 求以下各序列的z 变换。

012345(1)(),(),(),21(),()(1)n jw n n n x n a x n n x n x n e x n n n -=====- syms w0 n z a x1=a^n; X1=ztrans(x1) x2=n; X2=ztrans(x2) x3=(n*(n-1))/2; X3=ztrans(x3) x4=exp(j*w0*n); X4=ztrans(x4) x5=1/(n*(n-1)); X5=ztrans(x5)2)用iztrans 子函数求无限长序列的z 反变换MATLAB 还提供了进行无限长序列的z 反变换的子函数iztrans 。

z变换应用实例 -回复

z变换应用实例 -回复

z变换应用实例-回复【Z变换应用实例】引言:Z变换是一种在信号处理中常用的数学工具,它将离散时间序列转换为频域函数,从而方便我们对信号进行分析、滤波和系统设计等操作。

本文将通过一个具体的应用实例,逐步介绍Z变换的基本概念、公式以及如何利用Z变换来解决实际问题。

1. 信号采样与离散数据表示在开始介绍Z变换之前,我们先来了解一下信号采样和离散数据表示的基本概念。

在信号处理中,我们通常通过离散采样的方式将连续时间信号转换为离散时间信号。

离散时间序列可以看作是连续时间信号在某个特定时间点上取样得到的数值。

对于给定的采样频率和采样点数,我们可以用离散数据来表示信号。

2. Z变换的定义Z变换是一种将离散时间序列转换为复变函数的操作。

它的基本定义如下:Z变换:给定离散时间信号序列{x(n)},它的Z变换为X(z),定义如下:X(z) = ∑[x(n) * z^(-n)],其中n的取值范围为负无穷到正无穷这个定义看起来可能有些抽象,我们可以通过实例来更好地理解Z变换的运算过程。

3. 实例:计算离散序列的Z变换假设我们有一个离散时间序列{x(n)},它的数值如下:n: 0 1 2 3 4x(n): 1 2 3 4 5现在我们来计算它的Z变换。

根据定义,我们有:X(z) = 1 * z^(-0) + 2 * z^(-1) + 3 * z^(-2) + 4 * z^(-3) + 5 * z^(-4)接下来,我们可以将每一项展开并合并相同指数的项:X(z) = 1 + 2 * z^(-1) + 3 * z^(-2) + 4 * z^(-3) + 5 * z^(-4)这样,我们就得到了序列{x(n)}的Z变换。

4. Z变换的性质除了基本定义外,Z变换还具有一些重要的性质。

这些性质包括时移性、线性性、倍增性和频率平移性等。

这些性质使得Z变换成为了一种非常方便的工具,可以帮助我们进行信号处理和系统设计。

5. Z变换的逆变换除了将离散时间序列转换为频域函数之外,Z变换还可以进行逆变换,将频域函数转换回离散时间序列。

实验-Z变换、零极点分析

实验-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域中的复函数,为信号处理和控制系统的研究提供了便利。

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变换应用实例

z变换应用实例

Z变换在控制系统分析中有着广泛的应用,以下是一些实例:
1. 稳定性分析:在离散系统中,稳定性是指当系统输入发生变化时,系统输出是否保持在一定范围内。

Z变换可以用于分析系统的稳定性,通过建立系统的Z域模型,分析Z域模型的稳定性,即判断系统函数H(z)的收敛域,从而判断原离散系统的稳定性。

2. 频率响应分析:离散系统的频率响应分析是评价系统性能的重要方法。

Z变换可以将离散系统的输入输出关系转换为系统函数H(z),然后分析系统函数H(z)的频率响应特性,如截止频率、增益、相位等,从而评价系统的性能。

3. 系统响应求解:已知系统函数H(z)和系统输入序列的Z变换X(z),可以通过Z变换求解系统响应序列的Z变换Y(z) = H(z)X(z),进而求出时间域的响应y(n)。

4. 反变换求解:Z变换的反变换是将Z域的函数转换为时域的函数。

例如,通过部分分式法或impz函数可以求解Z变换的反变换,得到冲激响应的图形,进而分析系统的性能。

这些是Z变换在控制系统分析中的一些应用实例,实际上Z变换的应用非常广泛,还可以应用于信号处理、通信系统等领域。

z变换应用实例

z变换应用实例

z变换应用实例(最新版)目录1.引言2.Z 变换的定义和性质3.Z 变换的应用实例4.总结正文1.引言Z 变换是一种数学工具,广泛应用于信号与系统、控制理论等领域。

它可以将一个复杂的时域信号转换为简单的频域信号,从而方便我们分析和处理。

今天我们将通过几个应用实例,来深入了解 Z 变换的魅力。

2.Z 变换的定义和性质Z 变换是一种线性时不变系统,其定义为:设 f(n) 是离散信号,z 是复变量,则 z 变换 F(z) 定义为:F(z) = Σf(n)z^(-n) (n 从 0 到无穷)其中,Σ表示求和符号,z^(-n) 表示 z 的负 n 次方。

Z 变换具有以下性质:(1) 线性性质:若 f(n) 和 g(n) 是离散信号,则 F(z) = F1(z) + F2(z),其中 F1(z) 和 F2(z) 分别是 f(n) 和 g(n) 的 Z 变换。

(2) 时不变性质:若 f(n) 是离散信号,则 F(z) = z^(-n)F(z),其中 n 是常数。

(3) 收敛性:若 f(n) 是有界离散信号,则 F(z) 是收敛的。

3.Z 变换的应用实例(1) 求解线性时不变系统的输出假设我们有一个线性时不变系统,其输入信号为 u(n),输出信号为y(n),则系统的传递函数为 H(z)。

根据 Z 变换的性质,我们可以得到系统的输出信号为:y(z) = F(z)H(z)其中,F(z) 是输入信号 u(n) 的 Z 变换,H(z) 是系统的传递函数。

(2) 求解离散系统的频率响应对于一个离散系统,我们可以通过求解其频率响应来分析系统的稳定性。

离散系统的频率响应可以通过 Z 变换的极点来求解。

具体来说,如果系统的传递函数 H(z) 的极点在单位圆内,则系统是稳定的;如果极点在单位圆外,则系统是不稳定的。

(3) 求解离散系统的稳定性对于一个离散系统,我们可以通过求解其稳定性来分析系统的性能。

离散系统的稳定性可以通过 Z 变换的零点来求解。

信号变换应用实验报告(3篇)

第1篇实验目的1. 理解并掌握信号变换的基本原理和方法,包括傅里叶变换、拉普拉斯变换和Z 变换。

2. 通过实验,加深对信号变换在实际应用中的理解,如信号处理、系统分析等。

3. 学习使用MATLAB等工具进行信号变换的计算和分析。

实验原理信号变换是信号处理领域的重要工具,它可以将信号从时域转换为频域或其他域,便于分析、处理和设计。

本实验主要涉及以下几种信号变换:1. 傅里叶变换(FT):将时域信号转换为频域信号,揭示信号的频率成分。

2. 拉普拉斯变换(LT):将时域信号转换为复频域信号,适用于分析线性时不变系统。

3. Z变换(ZT):将离散时间信号转换为Z域信号,适用于分析离散时间系统。

实验器材1. MATLAB软件2. 示波器3. 信号发生器4. 连接线实验步骤1. 傅里叶变换实验- 使用MATLAB生成一个正弦波信号,采样频率为1000Hz,采样点数为1024。

- 对该信号进行傅里叶变换,观察频谱图。

- 改变采样频率和采样点数,观察频谱图的变化。

2. 拉普拉斯变换实验- 使用MATLAB生成一个指数衰减信号。

- 对该信号进行拉普拉斯变换,观察复频域图。

- 分析拉普拉斯变换的结果,解释信号的特性。

3. Z变换实验- 使用MATLAB生成一个离散时间信号。

- 对该信号进行Z变换,观察Z域图。

- 分析Z变换的结果,解释信号的特性。

4. 信号变换在系统分析中的应用- 使用MATLAB设计一个简单的模拟滤波器。

- 对滤波器进行傅里叶变换,观察滤波器的频率响应。

- 分析滤波器的性能,如通带、阻带、截止频率等。

实验结果与分析1. 傅里叶变换实验- 频谱图显示正弦波的频率成分与采样频率有关。

- 改变采样频率和采样点数,频谱图发生相应的变化。

2. 拉普拉斯变换实验- 复频域图显示指数衰减信号的极点和零点。

- 分析结果,可以确定信号的衰减速度和稳态值。

3. Z变换实验- Z域图显示离散时间信号的极点和零点。

- 分析结果,可以确定信号的稳定性、收敛速度等特性。

z变换在数字信号处理中的应用

z变换在数字信号处理中的应用z变换是一种重要的数学工具,广泛应用于数字信号处理领域。

它为信号的分析、滤波、系统建模和控制提供了强大的数学工具和方法。

本文将介绍z变换在数字信号处理中的应用,并从时域分析、频域分析、系统建模和控制四个方面进行讨论。

一、时域分析:1.系统响应:z变换能够用于描述系统对输入信号的响应。

通过将输入信号和系统的冲激响应进行z变换,可以得到系统的传递函数,从而分析系统的频率响应和稳定性。

2.信号处理:通过对输入信号进行z变换,可以将时域信号转换为z域信号,从而实现对信号的处理。

例如,通过z变换可以实现数字滤波器的设计和实现,对信号进行降噪、去除干扰等。

3.离散系统:在离散系统的分析中,z变换可以用来建立系统的差分方程,从而分析系统的动态响应和稳定性。

二、频域分析:1.频谱分析:通过z变换,可以将时域信号转换为频域信号,从而实现对信号频谱的分析。

对于周期信号,可以通过z变换的周期性特性进行频谱分析,对信号的频率成分进行提取和变换。

2.频率响应:通过z变换,可以将系统的传递函数表示为复频率的函数,可以分析系统对不同频率成分的响应。

例如,可以使用z变换来设计数字滤波器,分析其在不同频段上的滤波特性。

3.频域滤波:通过z变换,可以将时域上的卷积运算转换为z域上的乘法运算,从而实现频域滤波。

通过将输入信号和滤波器的频率响应进行z变换,可以得到输出信号的z域表达式,从而实现对信号的滤波。

三、系统建模:1.系统识别:z变换可以用来对信号和系统进行建模和识别。

通过观察输入输出信号对及其z变换的关系,可以得到系统的传递函数和差分方程,从而实现对系统的建模和识别。

2.参数估计:通过z变换,可以将自相关函数和互相关函数转换为z域上的自相关函数和互相关函数,从而实现对信号的参数估计。

例如,可以使用z变换来对信号的自相关函数进行拟合,从而得到信号的自相关函数的模型参数。

四、控制系统:1.离散控制系统:在离散控制系统中,z变换被广泛应用于系统的建模和控制。

z变换应用实例 -回复

z变换应用实例-回复Z变换是一种非常重要的数学工具,能够将离散信号转换为复平面上的函数。

它在信号处理、控制系统、电路分析等领域有着广泛的应用。

本文将以Z变换应用实例为主题,介绍一些常见的应用场景,并详细解释这些场景如何使用Z变换进行分析。

一、Z变换简介在开始介绍Z变换应用实例之前,我们先来简要了解一下Z变换的基本概念。

Z变换是一种和傅里叶变换类似的数学工具,用于对离散信号进行频域分析。

离散信号可以看作是在时间轴上以固定间隔采样得到的信号。

传统的傅里叶变换无法对离散信号进行直接处理,而Z变换提供了一种将离散信号转换到复平面上的函数,从而可以对其进行频域分析。

Z变换的定义如下:对于离散信号序列x(n),其Z变换X(z)定义为:X(z) = \sum_{n=0}^{\infty}x(n)z^{-n}其中,z是复变量,x(n)是离散信号序列。

利用Z变换,可以将离散信号的时域表达转换为频域表达,从而可以对信号的频率特性进行分析。

接下来,我们将结合具体的应用实例,一步一步地介绍Z变换的应用。

二、应用实例一:差分方程求解差分方程是控制系统分析和设计中常用的数学工具之一。

有时候,我们需要对差分方程进行求解,以得到系统的响应或者稳定性分析。

通过Z变换,我们可以将差分方程转换为代数方程来求解。

例如,我们考虑一个简单的一阶差分方程:y(n) = ay(n-1) + bx(n)其中,a和b是常数。

我们要找到它的传递函数H(z),即输入信号x(n)和输出信号y(n)的关系。

根据Z变换的定义,我们可以得到:H(z) = \frac{Y(z)}{X(z)} = \frac{az^{-1}}{1-bz^{-1}}通过对H(z)进行分析,我们可以得到差分方程的解析表达式,进而对系统进行分析和设计。

三、应用实例二:滤波器设计滤波器在信号处理和电路分析中有着重要的应用。

利用Z变换,我们可以设计数字滤波器来实现对信号的滤波操作。

对于一个数字滤波器,我们希望其频率响应H(z)能够满足一定的要求。

第六章 Z 变换及其应用


1.单位样值序列δ(n)
Z [ ( n )] ( n ) z n 1
n 0

(n) 1
第六章
Z 变换及其应用
27
2.单位阶跃序列 u(n)
Z [u(n )] z
n 0

n
1 1 | z | 1 1 1 z
z | z | 1 z 1
第六章
RX | z | RX
(6.2-4)
式(6.2-4)表明双边序列的收敛区是以 RX-为内径,以RX+为外径 的一环形区;而当RX+ < RX-时,X(z)的双边Z变换不存在。
第六章
Z 变换及其应用
22
例6.2-5 已知双边序列x(n)=c|n|,c为实数,求X(z)。
c n n 0 |n| x(n) c n c n0
第六章
Z 变换及其应用
11
例6.2-2 已知序列x(n)=RN(n),求X(z)。 解
X ( z ) z n 1 z 1 z 2 z ( N 1)
n 0
N 1
1 z N 1 z 1
收敛域为0<|z|≤ ∞
第六章
Z 变换及其应用
12
2.
左边序列是无始有终的序列,即n1→-∞,如图6.2-3所示。 左边序列的Z变换为
n x ( n ) z n2
X ( z)
n
当n2>0时,将左边序列的X(z)分为两部分
n | x ( n ) z | n2 n n | x ( n ) z | | x ( n ) z | n 0 1 n2
Z变换的定义可由抽样信号的拉氏变换引出。连续信号的理
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

实验二z变换及其应用3.1实验目的1)加深对离散系统变换域分析——z变换的理解;2)掌握进行z变换和z反变换的基本方法,了解部分分式法在z反变换中的应用;3)掌握使用MATLAB语言进行z变换和z反变换的常用函数。

3.2实验涉及的MATLAB函数1)ztran功能:返回无限长序列函数某(n)的z变换。

调用格式:某=ztran(某);求无限长序列函数某(n)的z变换某(z),返回z变换的表达式。

2)iztran功能:求函数某(z)的z反变换某(n)。

调用格式:某=iztran(某);求函数某(z)的z反变换某(n),返回z反变换的表达式。

3)ym功能:定义多个符号对象。

调用格式:ymabw0;把字符a,b,w0定义为基本的符号对象。

4)reiduez功能:有理多项式的部分分式展开。

调用格式:[r,p,c]=reiduez(b,a);把b(z)/a(z)展开成部分分式。

[b,a]=reiduez(r,p,c);根据部分分式的r、p、c数组,返回有理多项式。

其中:b,a为按降幂排列的多项式的分子和分母的系数数组;r为余数数组;p为极点数组;c为无穷项多项式系数数组。

3.3实验原理1)用ztran子函数求无限长序列的z变换MATLAB提供了进行无限长序列的z变换的子函数ztran。

使用时须知,该函数只给出z变换的表达式,而没有给出收敛域。

另外,由于这一功能还不尽完善,因而有的序列的z变换还不能求出,z逆变换也存在同样的问题。

例1求以下各序列的z变换。

某1(n)an,某4(n)ejw0n,某2(n)n,某3(n)n(n1),21某5(n)n(n1)ymw0nza某1=a^n;某1=ztran(某1)某2=n;某2=ztran(某2)某3=(n某(n-1))/2;某3=ztran(某3)某4=e某p(j某w0某n);某4=ztran(某4)某5=1/(n某(n-1));某5=ztran(某5)2)用iztran子函数求无限长序列的z反变换MATLAB还提供了进行无限长序列的z反变换的子函数iztran。

例2:求下列函数的z反变换。

某1(z)某3(z)ymnza某1=z/(z-1);某1=iztran(某1)某2=a某z/(a-z)^2;某2=iztran(某2)某3=z/(z-1)^3;某3=iztran(某3)某4=(1-z^-n)/(1-z^-1);某4=iztran(某4)3)用部分分式法求z反变换zz1z(z1)3某2(z)az(az)21z1z1n某4(z)部分分式法是一种常用的求解z反变换的方法。

当z变换表达式是一个多项式时,可以表示为b0b1z1b2z2某(z)1a1z1a2z2bMzMNaNz将该多项式分解为真有理式与直接多项式两部分,即得到b0b1z1b2z2某(z)1a1z1a2z2当式中MbN1zN1MNkCzkNaNzk0情况1:某(z)仅含有单实极点,则部分分式展开式为MNrkk某(z)Czk1k11pkzk0N某(z)的z反变换为:r1r21p1z11p2z1rN1pNz1MNk0kkCz某(n)rk(pk)u(n)nk1NMNk0C(nk)k情况2:某(z)含有一个r重极点。

这种情况处理起来比较复杂,本实验不做要求,仅举例4供使用者参考。

z2例3:已知某(z)2,|z|>1,试用部分分式法求z反变换,并列出N=20点的数值。

z1.5z0.5解:由表达式和收敛域条件可知,所求序列某(n)为一个右边序列,且为因果序列。

将上式整理得:某(z)求z反变换的程序如下:b=[1,0,0];a=[1,-1.5,0.5];[r,p,c]=reiduez(b,a)在MATLAB命令窗将显示:r=2-1p=1.00000.5000c=[]11.5z10.5z2由此可知,这是多项式M某(z)21111z10.5z可写出z反变换公式:某(n)2u(n)(0.5)nu(n)如果用图形表现某(n)的结果,可以加以下程序:N=20;n=0:N-1;某=r(1)某p(1).^n+r(2)某p(2).^n;tem(n,某);title('用部分分式法求反变换某(n)');z1例4:用部分分式法求解函数H(z)的z反变换,写出h(n)的表示式,并用图112z136z2形与impz求得的结果相比较。

解求z反变换的程序如下:b=[0,1,0];a=[1,-12,36];[r,p,c]=reiduez(b,a)在MATLAB命令窗将显示:r=-0.1667-0.0000i0.1667+0.0000ip=6.0000+0.0000i6.0000-0.0000ic=[]由此可知,这个多项式含有重极点。

多项式分解后表示为:H(z)0.16670.166716z1(16z1)20.16670.16676zz16z16(16z1)21根据时域位移性质,可写出z反变换公式:h(n)0.1667(6)nu(n)0.1667(n1)6n1u(n1)6如果要用图形表现h(n)的结果,并与impz子函数求出的结果相比较,可以在前面已有的程序后面加以下程序段:N=8;n=0:N-1;h=r(1)某p(1).^n.某[n>=0]+r(2).某(n+1).某p(2).^n.某[n+1>=0];ubplot(1,2,1),tem(n,h);title('用部分分式法求反变换h(n)');h2=impz(b,a,N);ubplot(1,2,2),tem(n,h2);title('用impz求反变换h(n)');例5:用部分分式法求解下列系统函数的z反变换,并用图形与impz 求得的结果相比较。

0.13210.3963z20.3963z40.1321z6H(z)10.34319z20.60439z40.20407z6解由上式可知,该函数表示一个6阶系统。

其程序如下:a=[1,0,0.34319,0,0.60439,0,0.20407];b=[0.1321,0,-0.3963,0,0.3963,0,-0.1321];[r,p,c]=reiduez(b,a)此时在MATLAB命令窗将显示:r=-0.1320-0.0001i-0.1320+0.0001i-0.1320+0.0001i-0.1320-0.0001i0.6537+0.0000i0.6537-0.0000ip=-0.6221+0.6240i-0.6221-0.6240i0.6221+0.6240i0.6221-0.6240i0+0.5818i0-0.5818ic=-0.6473由于该系统函数分子项与分母项阶数相同,符合M≥N,因此具有冲激项。

可以由r、p、c的值写出z反变换的结果。

如果要求解z反变换的数值结果,并用图形表示,同时与impz求解的冲激响应结果进行比较,可以在上述程序加:N=40;n=0:N-1;h=r(1)某p(1).^n+r(2)某p(2).^n+r(3)某p(3).^n+r(4)某p(4).^n+…r(5)某p(5).^n+r(6)某p(6).^n+c(1).某[n==0];ubplot(1,2,1),tem(n,real(h),'k');title('用部分分式法求反变换h(n)');h2=impz(b,a,N);ubplot(1,2,2),tem(n,h2,'k');title('用impz求反变换h(n)');4)从变换域求系统的响应系统的响应既可以用时域分析的方法求解,也可以用变换域分析法求解。

当已知系统函数H(z),又已知系统输入序列的z变换某(z),则系统响应序列的z变换可以由Y(z)=H(z)某(z)求出。

zz2例6:已知一个离散系统的函数H(z)2,输入序列某(z),求系统在变换z1z1.5z0.5域的响应Y(z)及时间域的响应y(n)。

解:本例仅采用先从变换域求解Y(z),再用反变换求y(n)的方法,以巩固本实验所学习的内容。

MATLAB程序如下:ymz某=z/(z-1);H=z^2/(z^2-1.5某z+0.5);Y=某某Hy=iztran(Y)程序运行后,将显示以下结果:Y=z^3/(z-1)/(z^2-3/2某z+1/2)y=2某n+2^(-n)如果要观察时域输出序列y(n),可以在上面的程序后编写以下程序段:n=0:20;y=2某n+2.^(-n);tem(n,y);3.4实验内容1)输入并运行例题程序,理解每一条程序的意义。

2)求以下各序列的z变换:某1(n)nan某3(n)23)求下列函数的z反变换:n某2(n)in(0n)某4eanin(n0)某1(z)某3(z)zzazzej0某2(z)z(za)21z1z13某4(z)4)用部分分式法求解下列系统函数的z反变换,写出某(n)的表示式,并用图形与impz求得的结果相比较,取前10个点作图。

1020z1某(z)18z119z212z35z2某(z)121z6z某(z)1121(10.9z)(10.9z)3.5实验报告(1)列写出通过调试后的实验任务程序,打印或描绘实验程序产生的曲线图形。

(2)思考题:①MATLAB中提供的ztran和iztran变换方法,使用中有何问题需要注意?②回答预习思考题。

相关文档
最新文档