《数字信号处理》学习09-部分分式展开法计算z 逆变换

在之前的文章中,我已经学习了使用留数法(围线积分法)来计算z逆变换

《数字信号处理》学习08-围线积分法(留数法)计算z 逆变换-CSDN博客

接着学习第二种计算z变换的方法:部分分式展开法。 

目录

一,部分分式展开法的相关概念 

 二,习题

🐕


一,部分分式展开法的相关概念 

一个离散时间系统,其输入序列x(n)和输出序列y(n)之间的关系可以通过线性常系数差分方程表示:

y(n)+\sum_{i=1}^{N}a_{i}y(n-i)=\sum_{i=0}^{M}b_{i}x(n-i)

 现在对两边同时进行z变换,需要注意的是在时域上的延迟 i 单位对应的z域上乘上 z^{-i}  ,如下:

Z[y(n)]=Y(z)

Z[y(n-i)]=Y(z)z^{-i}       

Z[x(n-i)]=X(z)z^{-i}                                                                                                               //  因此,差分方程 y(n)+\sum_{i=1}^{N}a_{i}y(n-i)=\sum_{i=0}^{M}b_{i}x(n-i)的z变换如下:   

Y(z)+\sum_{i=1}^{N}a_{i}Y(z)z^{-i}=\sum_{i=0}^{M}b_{i}X(z)z^{-i}       

// 自变量为i,X(z) 和 Y(z) 可以看成常数提出来

Y(z)+Y(z)\sum_{i=1}^{N}a_{i}z^{-i}=X(z)\sum_{i=0}^{M}b_{i}z^{-i}     

//  对 Y(z)  进行合并同类项

Y(z)(1+\sum_{i=1}^{N}a_{i}z^{-i})=X(z)\sum_{i=0}^{M}b_{i}z^{-i}     

// 移项得到传输函数 H(z)

H(z)=\frac{Y(z)}{X(z)}= \frac{\sum_{i=0}^{M}b_{i}z^{-i}}{1+\sum_{i=1}^{N}a_{i}z^{-i}} 

x(z)的z变换和传输函数H(z)的公式类似,当X(z)是z的有理分式时,一般可以表示为:

X(z)=\frac{B(z)}{A(z)}=\frac{\sum_{i=0}^{M}b_{i}z^{-i}}{1+\sum_{i=0}^{N}a_{i}z^{-i}}=X_{1}(z)+X_{2}(z)+X_{3}(z)+...+X_{k}(z)

其中,A(z)B(z)是z的实数系数多项式,且没有公因式,则可以将X(z)展开成上述部分分式的形式。之后只需要对每一个部分分式求z逆变换,最后将各个z逆变换相加就是所求的原序列x(n) 

 二,习题

习题1

使用部分分式法计算z逆变换,如下题

解:
1) 由题可得

X(z)=\frac{3-\frac{5}{6}z^{-1}}{(1-\frac{1}{4}z^{-1})(1-\frac{1}{3}z^{-1})}

=\frac{z(3z-\frac{5}{6})}{(z-\frac{1}{4})(z-\frac{1}{3})}

\frac{X(z)}{z}=\frac{(3z-\frac{5}{6})}{(z-\frac{1}{4})(z-\frac{1}{3})}=\frac{A}{z-\frac{1}{3}}+\frac{B}{z-\frac{1}{4}}

// 使用之前的留数法,将待定系数AB求出来(注:这里用留数法只求系数,不求原序列x(n)

A|_{z=\frac{1}{3}}

=(z-\frac{1}{3})\frac{X(z)}{z}|_{z=\frac{1}{3}}

=(z-\frac{1}{3})\frac{3z-\frac{5}{6}}{(z-\frac{1}{4})(z-\frac{1}{3})}|_{z=\frac{1}{3}}

=\frac{3z-\frac{5}{6}}{(z-\frac{1}{4})}|_{z=\frac{1}{3}}

=\frac{\frac{12}{12}-\frac{10}{12}}{\frac{4}{12}-\frac{3}{12}}

=\frac{\frac{2}{12}}{\frac{1}{12}}

=2

B|_{z=\frac{1}{4}}

=(z-\frac{1}{4})\frac{X(z)}{z}|_{z=\frac{1}{4}}

=(z-\frac{1}{4})\frac{3z-\frac{5}{6}}{(z-\frac{1}{4})(z-\frac{1}{3})}|_{z=\frac{1}{4}}

=\frac{3z-\frac{5}{6}}{(z-\frac{1}{3})}|_{z=\frac{1}{4}}

=\frac{\frac{9}{12}-\frac{10}{12}}{\frac{3}{12}-\frac{4}{12}}

=\frac{\frac{-1}{12}}{\frac{-1}{12}}

=1

\frac{X(z)}{z}=\frac{2}{z-\frac{1}{3}}+\frac{1}{z-\frac{1}{4}}

X(z)==\frac{2z}{z-\frac{1}{3}}+\frac{z}{z-\frac{1}{4}}

∵ 题目要求z变换的收敛为 |z|>\frac{1}{3},右边序列,且该序列为因果序列

又∵参照z变换表

∴收敛域为 |z|>\frac{1}{3}时,对应的原序列  x(n)=2\cdot (\frac{1}{3})^{n}u(n)+ (\frac{1}{4})^{n}u(n)

2)

∵ 由题 1)得X(z)==\frac{2z}{z-\frac{1}{3}}+\frac{z}{z-\frac{1}{4}}

又∵ 题目要求z变换的收敛为 \frac{1}{4}<|z|<\frac{1}{3} ,双边序列

∴ 查表可知

当 \frac{1}{4}<|z| 时, x(z)=\frac{z}{z-\frac{1}{4}}对应的原序列为右边序列 x(n)=(\frac{1}{4})^{n}u(n)

当 |z|<\frac{1}{3} 时, x(z)=\frac{z}{z-\frac{1}{3}}对应的原序列为左边序列 x(n)=-2(\frac{1}{3})^{n}u(-n-1)

// 左边序列z变换的原序列查表可知

∴综上收敛域为  \frac{1}{4}<|z|<\frac{1}{3} 时,对应的原序列  x(n)=(\frac{1}{4})^{n}u(n)-2(\frac{1}{3})^{n}u(-n-1)

3)

∵ 由题 1)得X(z)=\frac{2z}{z-\frac{1}{3}}+\frac{z}{z-\frac{1}{4}}

又∵ 题目要求z变换的收敛为 |z|<\frac{1}{3} ,序列为左边序列

∴ 查表可知

∴ 收敛为 |z|<\frac{1}{3} 时,原序列 x(n)=-(\frac{1}{4})^{n}u(-n-1)-2(\frac{1}{3})^{n}u(-n-1)

习题2

已知X(z)=\frac{-3z^{-1}}{2-5z^{-1}+2z^{-2}},分别求:

1)收敛域0.5<|z|<2对应的原序列 x(n)

2)收敛域|z|>2对应的原序列 x(n)

 解:

1)

// 使用十字相乘法将分母简化成可以知道极点的形式

X(z)=\frac{-3z^{-1}}{2-5z^{-1}+2z^{-2}}   //  上下同乘 z^{2} ,分式大小不变

=\frac{-3z}{(2z-1)(z-2)}   // 分母十字相乘法

∵ \frac{X(z)}{z}=\frac{-3}{(2z-1)(z-2)}=(\frac{A}{z-0.5}+\frac{B}{z-2})

// 使用之前的留数法,将待定系数AB求出来(注:这里用留数法只求系数,不求原序列x(n)

A|_{z=\frac{1}{2}}

=(z-\frac{1}{2})\frac{X(z)}{z}|_{z=\frac{1}{2}}

=(z-\frac{1}{2})\cdot(\frac{-3}{2(z-\frac{1}{2})(z-2)})|_{z=\frac{1}{2}}

=\frac{-3}{z-2}|_{z=\frac{1}{2}}

=2

B|_{z=2}

=(z-2)\frac{X(z)}{z}|_{z=2}

=(z-2)\cdot(\frac{-3}{(2z-1)(z-2)})|_{z=2}

=\frac{-3}{2z-1}|_{z=2}

=-1

∴ X(z)=\frac{z}{z-0.5}-\frac{z}{z-2}

 ∵ 题目要求z变换的收敛为 0.5<|z|<2 ,序列为双边序列

0.5<|z| 时,对应的原序列为右边序列 x(n)=(\frac{1}{2})^{n}u(n)

|z|<2 时,对应的原序列为左边序列 x(n)=-2^{n}u(-n-1)

综上 ,收敛域为 0.5<|z|<2 时的原序列x(n)=2^{-n}u(n)+2^{n}u(-n-1)

 或写成x(n)=2^{-|n|},其中   _{n<0,x(n)=2^{n}}^{n>=0,x(n)=2^{-n}}

  2)收敛域|z|>2对应的原序列 x(n)

解:

由题 1)得X(z)=\frac{z}{z-0.5}-\frac{z}{z-2}

 ∵ 题目要求z变换的收敛为 |z|>2 ,序列为右边序列

又∵ 查表 ↓

x(n)=2^{-n}u(n)-2^{n}u(n) 

习题3:

用分部积分法求下面象函数X(z) 的原序列 x(n)

解:

// 先将式子中z变量的指数变成正数,分子分母同时乘z^{2},式子大小不变,题目式子变为如下:

X(z)=\frac{z^{2}-\frac{1}{2}z}{z^{2}-\frac{1}{4}}

∵ \frac{X(z)}{z}=\frac{z-\frac{1}{2}}{z^{2}-\frac{1}{4}}=\frac{z-\frac{1}{2}}{(z+\frac{1}{2})(z-\frac{1}{2})}=\frac{1}{z+\frac{1}{2}}

 ∴ X(z)=\frac{z}{z+\frac{1}{2}}

∵ 该z变换的收敛域为 |z|>\frac{1}{2},右边序列,且为因果序列

又∵ 查表 ↓

∴ 原序列  x(n)=(-\frac{1}{2})^{n}u(n)

以上就是用部分分式展开法计算z逆变换的相关知识。

如果有问题请在评论区留言或者是私信我,回复时间不超过一天。 

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/pingmian/56471.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

决策智能与强化学习:重放比率(replay ratio)

知乎&#xff1a;DILab决策实验室&#xff08;已授权&#xff09;链接&#xff1a;https://zhuanlan.zhihu.com/p/898641863 0. 概览 近年来&#xff0c;深度强化学习&#xff08;Deep Reinforcement Learning, DRL&#xff09;在诸多领域取得了显著的成果。然而&#xff0c;随…

域环境模拟实验搭建

1. 域环境搭建 总体来说下一步下一步即可 域&#xff1a;统一的管理计算机的集群&#xff0c;中心管理机器&#xff08;域控制器 DC&#xff09;管理整个内网&#xff08;域内成员机器&#xff09; 条件&#xff1a;1. 需要一台域控制器 -- windows server版本&#xff08;w…

高级英语1第四版教材全解pdf课后答案+课文翻译张汉熙

《高级英语1》是张汉熙教授编著的一本英语教材&#xff0c;广泛用于国内高校英语专业高年级学生的教学。这本书以提高学生的英语综合能力为目标&#xff0c;注重语言知识的系统性和实用性&#xff0c;同时强调跨文化交际能力的培养。书中选材丰富&#xff0c;涵盖了文学、历史、…

Cloudlog delete_oqrs_line 未授权SQL注入漏洞复现

0x01 产品简介 Cloudlog 是一个自托管的 PHP 应用程序,可让您在任何地方记录您的业余无线电联系人。使用PHP和MySQL构建的基于Web的业余无线电记录应用程序支持从HF到微波的一般站记录任务 0x02 漏洞概述 Cloudlog delete_oqrs_line 接口存在未授权SQL注入漏洞,未经身份验…

UE5学习笔记24-添加武器弹药

一、给角色的武器添加弹药 1.创建界面&#xff0c;根据笔记23的界面中添加 2.绑定界面控件 UPROPERTY(meta (Bindwidget))UTextBlock* WeaponAmmoAmount;UPROPERTY(meta (Bindwidget))UTextBlock* CarriedAmmoAmount; 3.添加武器类型枚举 3.1创建武器类型枚举头文件 3.2创建文…

Java 小游戏《超级马里奥》

文章目录 一、效果展示二、代码编写1. 素材准备2. 创建窗口类3. 创建常量类4. 创建动作类5. 创建关卡类6. 创建障碍物类7. 创建马里奥类8. 编写程序入口 一、效果展示 二、代码编写 1. 素材准备 首先创建一个基本的 java 项目&#xff0c;并将本游戏需要用到的图片素材 image…

PPT怎么转成PDF?5个超简单转换方法,快来学(包靠谱)

PPT怎么转成PDF&#xff1f;PPT&#xff08;PowerPoint演示文稿&#xff09;是一种重要的文件格式&#xff0c;通过它&#xff0c;我们能够直观地呈现信息、分享观点以及进行商务演示。然而在有些时候&#xff0c;我们需要将PPT文件转化为PDF格式&#xff0c;以便于分享、打印或…

基于PHP考研互助系统【附源码】

基于PHP考研互助系统 效果如下&#xff1a; 系统首页界面 用户注册界面 考研论坛页面 每日打卡页面 管理员登录主页面 管理员主界面 用户管理界面 备考经验界面 研究背景 近些年&#xff0c;随着中国经济发展&#xff0c;人民的生活质量逐渐提高&#xff0c;对网络的依赖性越…

软件Bug和缺陷的区别是什么?

前言 在任何软件生命周期中&#xff0c;软件缺陷的出现几乎是不可避免的。建立一套有效的缺陷管理流程的目的是为了减少软件缺陷出现的几率&#xff0c;并且大幅度降低由于软件缺陷带来的负面影响。对于缺陷管理流程的投资&#xff0c;可以大幅度的降低由于返工/修复缺陷导致的…

模版和STL

个人主页&#xff1a;Jason_from_China-CSDN博客 所属栏目&#xff1a;C系统性学习_Jason_from_China的博客-CSDN博客 所属栏目&#xff1a;C知识点的补充_Jason_from_China的博客-CSDN博客 前言 string是标准库中的一个类&#xff0c;它位于<string>头文件中。 这个类提…

Python基础:19、Python异常、模块与包

1&#xff09;python异常 什么是异常&#xff1a;异常就是程序运行的过程中出现了错误 bug是什么意思&#xff1a;bug就是指异常的意思&#xff0c;因为历史因为小虫子导致计算机失灵的案例&#xff0c;所以延续至今&#xff0c;bug就代表软件出现错误。 2&#xff09;异常的…

gitlab的基本用法之创建用户和组

安装好gitlab后要做一些基本的创建 创建顺序&#xff1a;先创建组——再创建项目——最后创建用户项目、组、用户的关系&#xff1a;把创建的用户和项目加入组中&#xff0c;只有同一组的用户才能操作该组的项目。 一、创建群组 点击设置就可以看到刚才创建的组 二、创建项…

SpringCloudStream使用StreamBridge实现延时队列

利用RabbitMQ实现消息的延迟队列 一、安装RabbitMQ 1、安装rabbitmq 安装可以看https://blog.csdn.net/qq_38618691/article/details/118223851,进行安装。 2、安装插件 安装完毕后,exchange是不支持延迟类型的,需要手动安装插件,需要和安装的rabbitmq版本一致 https:…

可以在桌面上用的倒计时提醒app下载

在忙碌的工作日常中&#xff0c;我们常常需要记住各种截止日期和重要事件。为了确保这些任务按时完成&#xff0c;一款桌面倒计时提醒应用变得尤为重要。想象一下&#xff0c;当你在电脑桌面上就能清晰地看到剩余时间&#xff0c;这无疑会增加你的工作效率和紧迫感。 敬业签就…

红日安全vulnstack (二)

目录 环境搭建 网卡设置 修改Kali网段 IP 分布 WEB渗透 Weblogin服务开启 漏洞扫描 CVE工具利用 MSF上线 内网渗透 域内信息收集 凭证横向移动 权限维持 黄金票据 参考文章 https://www.cnblogs.com/bktown/p/16904232.htmlhttps://blog.csdn.net/m0_75178803/ar…

好用的python相关的AI工具Bito介绍

插件名称&#xff1a;Bito 好用的python相关的AI工具Bito介绍 step 1:点插件step 2&#xff1a;搜索bito并安装step3 &#xff1a;需要登录&#xff0c;要有真实邮箱&#xff0c;按步骤走就行&#xff0c;完后就可以使用 step 1:点插件 step 2&#xff1a;搜索bito并安装 step3…

关于新国标强制电动车应内置北斗定位模块的规定有哪些?附北斗定位芯片对比参数

关于新国标要求电动自行车内置的北斗定位功能&#xff0c;需要符合以下几点&#xff1a; 支持UART或SPI接口至少支持接收处理北斗B1C和B2a信号具备定位信息的采集、存储和发送功能&#xff08;其中定位信息包括&#xff1a;经度、纬度、速度、定位时间&#xff09;具备采集、存…

1.项目初始化

目录 1.相关联的数据库表 2.使用gorm操作数据库 使用gen生成model和对数据库的操作 3.使用viper进行配置管理 读取配置文件 进行热更新 4.使用Pflag来进行命令行参数解析 5.使用日志slog 日志轮转与切割功能 6.错误码和http返回格式标准化 提供错误码 提供错误类型…

Leetcode—1115. 交替打印 FooBar【中等】(多线程)

2024每日刷题&#xff08;180&#xff09; Leetcode—1115. 交替打印 FooBar C实现代码 class FooBar { private:int n;sem_t fooSem;sem_t barSem;public:FooBar(int n) {this->n n;sem_init(&fooSem, 0, 1);sem_init(&barSem, 0, 0);}~FooBar() {sem_destroy(&…

免杀对抗—内存加载UUID标识IPV4地址MAC地址

前言 试想我们开辟一块内存,然后直接将shellcode写入到对应的内存中并且该内存是可读可写可执行的状态,那么这种方式太容易被AV所查杀,因此当我们如果是利用Windows自身提供的API来将加密或者封装好的shellcode写入到内存执行的话,将会大大增加查杀的难度。 参考文章&#xf…