《数字信号处理》学习08-围线积分法(留数法)计算z 逆变换
目录
一,z逆变换相关概念
二,留数定理相关概念
三,习题
一,z逆变换相关概念
接下来开始学习z变换的反变换-z逆变换(z反变化)。
由象函数 求它的原序列
的过程就称为
逆变换。即
。
求z逆变换的方法通常有三种:围线积分法,部分分式展开法和长除法。
由于原序列 就是罗朗级数
,因此用围线积分法求z逆变换的积分公式如下👇
,
可以看到上式比较复杂,如果直接计算围线积分,会比较麻烦,因此可以借助复变函数的留数定理来计算出围线积分的结果。
二,留数定理相关概念
在使用留数定理之前,需要知道的基础知识点如下(也可以去看《复变函数》这本书):
- 复值函数:是指其值域为复数的函数(输入输出都是复数)。例如求z逆变换的积分公式中的
就是复值函数,输入(自变量)z 是复数,输出(因变量)
也是复数。为了简化运算过程,一般令
- 解析函数(也叫全纯函数):在某个区域内可以用幂级数展开的复值函数(如果复值函数在某一点可微,在该点的领域内也可微,则称之为解析函数)。
解析函数有一个很重要的性质:导数存在。
导数存在可以推出该函数具有可微性(在复分析中,如果一个解析函数的导数存在,那么该函数在其定义域内是可微的。),也可以知道该函数具有连续性。 - 留数:用 Res(复值函数,孤立奇点) 表示。由于积分公式中有复值函数
,因此可以假设复值函数为
,同时假设存在孤立奇点
,则留数可表示为
。
- 孤立奇(qi)点:是指一个复函数在某一点的邻域内不定义或不解析,但在该点的某个邻域外是解析的。奇点又称为奇异点。
假设复函数在
处是奇点, 那么孤立奇点可以分为以下三类:
1)本性奇点:极限不存在。
2)可去奇点:极限存在且有限。
3)极点:极限存在且为无穷。
在使用围线积分法求逆变换的计算中孤立奇点都找极点位置。即
- 留数定理:通过计算留数的结果,得到原序列
三,习题
例如给出一道题要求使用留数法求逆变换:求象函数
的原序列
,如下👇
题目1:已知,分别求:
1)收敛域对应的原序列
2)收敛域对应的原序列
解:
1)
// 先将 的分母因式分解
∵
= // 分子分母同时乘
,分式的大小(值)不变
= // 分母使用十字相乘法化简
又∵ 积分公式: ,且
∴
∵
// 符合双边序列的变换收敛域
,离散时间变量
// 如果z变换不清楚的可以查看下面的文章:
// 《数字信号处理》学习07-z变换_左边序列,右边序列、双边序列。-CSDN博客
所以的极点可以分为如下两种情况:
① 当 时,分子上存在一个极点,即
,令
,得
分母存在两个极点,即:
当时,得
当时,得
对应的z平面收敛域及围线C所包围的区域如下:

// 观察上图,可以发现,围线C所包围的圆里面有两个极点: 和
// 由于 是n阶的极点,因此,围线C所包含的极点需要反着取(即使用围线C外极点)
如下图:
//从上图可以看到围线C外的极点只有一个
∵
// 使用留数定理时,由于是围线C外积分,因此留数的值需要取负数:
∴
=
=
==
② 当 时,
分母存在两个极点,即:
当时,得
当时,得
但围线C只能包含一个极点,如下图:

// 使用留数定理,取的是围线C内的极点,因此,留数为正。
=
=
=
综上,原序列
2)求 收敛域对应的原序列
。
根据题目可得
象函数的收敛域
,符合右边序列的收敛域形式,由于收敛域的外部区域通常与因果序列相关(该右边序列是因果序列)。因此这里只讨论当
时的情况,
分母上的两个极点分别为:,
对应的围线C所包含的极点如下图所示:
// 使用留数定理,因为极点都在围线C内,所以留数前面为正,不加负号。
=
所以当收敛域时,对应的原序列为
。
题目2:用留数法求下面象函数 的原序列

解:
// 先将式子中z变量的指数变成正数,分子分母同时乘,式子大小不变,题目式子变为如下:
∵ 积分公式:
又∵
= // 因式分解
=
∴
∵ ,符合右边序列的z变换收敛域,且该右边序列为因果序列,此时
∴ 在z复平面上的收敛域及围线C的位置如下图所示:
// 观察上图,可以看到,极点位于围线C内,因此留数定理使用的是C内积分。
=
=
=
// 因为是因果序列,序列的离散时间变量n只分布在坐标轴的右边,所以需要加上n的取值范围
// 一般将序列乘上单位阶跃信号就可以表示该序列只在正半轴有取值。
所以当收敛域 时,对应的原序列为
。
题目3:用留数法求下面象函数 的原序列

解:
1)
∵ 积分公式:
∴
= // 分子分母同时乘
,大小不变
∵ ,符合右边序列z变换的收敛域,且该右边序列为因果序列,此时
∴在分母上存在两个极点:
在z复平面上的收敛域如下图所示:
∵极点都位于围线C内(c内极点),留数前面不用加负号。
// 根据留数定理,可求出原序列
=
=
=
=
=
=
=
// 因为是因果序列,序列的离散时间变量n只分布在坐标轴的右边,所以需要加上n的取值范围
// 一般将序列乘上单位阶跃信号就可以表示该序列只在正半轴有取值。
所以当收敛域 时,对应的原序列为
。
2)
由题(1)得
∵ ,符合双边序列z变换的收敛域,此时
∴需要进行分类讨论
当 时,分子上存在极点
∴在分母上存在两个极点:
// 因为是n阶极点,所以留数定理使用的是围线C外的极点
在z复平面上的收敛域如下图所示:

// 观察上图可以看到,围线C外的极点只有一个,此时的留数公式前需要加上负号
=
=
=
=
=
=
=
// 因为上面是在n<0时求出的结果,即单位阶跃信号翻褶之后再向左平移一个单位u(-n-1),所以需要加上定义域,上式结果乘上u(-n-1)
当 时,
// 接下来讨论n>=0的情况
当 时,不存在n阶极点z=0,收敛域依旧不变,此时围线C所包含的极点有一个
,如下

// 观察上图可以看到,围线C内的极点有一个,此时的留数公式前不需要加上负号
=
=
=
=
=
=
=
// 因为上面是在n>=0时求出的结果,所以需要加上定义域,上式结果乘上u(n)
当 时,
// 最后将 n<0 和 n>=0的结果合并在一起
综上,收敛域为 的原序列为
3)
由题(1)得
∵ ,符合左边序列z变换的收敛域,收敛域及围线C围在小于1/4的位置如下:

// 观察上图可以看到,围线C内无极点,而在围线C外存在两个极点,此时留数前需要加负号
=
=
=
=
=
=
=
// 因为上面是在n<0时求出的结果,即单位阶跃信号翻褶之后再向左平移一个单位u(-n-1),所以需要加上定义域,上式结果乘上u(-n-1)
当 时,原序列为
以上就是用留数法求z逆变换的相关内容,上述的计算也可以使用分部积分法和长除法,后面我会接着学习,有兴趣的关注专栏,有问题的请在评论区留言或者是私信我,回复时间不超过一天。
相关文章:
《数字信号处理》学习08-围线积分法(留数法)计算z 逆变换
目录 一,z逆变换相关概念 二,留数定理相关概念 三,习题 一,z逆变换相关概念 接下来开始学习z变换的反变换-z逆变换(z反变化)。 由象函数 求它的原序列 的过程就称为 逆变换。即 。 求z逆变换…...
vue3中的computed属性
模板界面: <template><div class"person"><h2>姓: <input type"text" v-model"person.firstName" /></h2><h2>名: <input type"text" v-model"person…...
C++学习笔记之vector容器
天上月,人间月,负笈求学肩上月,登高凭栏眼中月,竹篮打水碎又圆。 山间风,水边风,御剑远游脚下风,圣贤书斋翻书风,风吹浮萍又相逢。 STL(Standard Template Library,标准模板库 ) 从…...
LeNet-5(论文复现)
LeNet-5(论文复现) 本文所涉及所有资源均在传知代码平台可获取 文章目录 LeNet-5(论文复现)概述LeNet-5网络架构介绍训练过程测试过程使用方式说明 概述 LeNet是最早的卷积神经网络之一。1998年,Yann LeCun第一次将LeN…...
基于SpringBoot+Vue+Uniapp汽车保养系统小程序的设计与实现
详细视频演示 请联系我获取更详细的演示视频 项目运行截图 技术框架 后端采用SpringBoot框架 Spring Boot 是一个用于快速开发基于 Spring 框架的应用程序的开源框架。它采用约定大于配置的理念,提供了一套默认的配置,让开发者可以更专注于业务逻辑而…...
【问题实战】Jmeter中jtl格式转换图片后如何分开展示各个性能指标?
【问题实战】Jmeter中jtl格式转换图片后如何分开展示各个性能指标? 遇到的问题解决方法查看修改效果 遇到的问题 JMeter测试计划中只设置了一个性能监控器jpgc - PerfMon Metrics Collector;在这个监控器中设置几个性能监控指标,比如CPU、Di…...
解决 MySQL 连接数过多导致的 SQLNonTransientConnectionException 问题
这里写目录标题 解决 MySQL 连接数过多导致的 SQLNonTransientConnectionException 问题1. 概述2. 问题描述异常日志的关键部分: 3. 原因分析3.1. MySQL 连接数配置3.2. 连接池配置问题3.3. 代码中未正确关闭连接3.4. 高并发导致连接需求激增 4. 解决方案4.1. 增加 …...
猫头虎分享:什么是 ChatGPT 4o Canvas?
猫头虎是谁? 大家好,我是 猫头虎,猫头虎技术团队创始人,也被大家称为猫哥。我目前是COC北京城市开发者社区主理人、COC西安城市开发者社区主理人,以及云原生开发者社区主理人,在多个技术领域如云原生、前端…...
qiankun 主项目和子项目都是 vue2,部署在同一台服务器上,nginx 配置
1、主项目配置 1.1 micro.vue 组件 <template><div id"container-sub-app"></div> </template><script> import { loadMicroApp } from qiankun; import actions from /utils/actions.js;export default {name: microApp,mixins: [ac…...
深入浅出MongoDB(七)
深入浅出MongoDB(七) 文章目录 深入浅出MongoDB(七)查询优化创建索引以支持读取操作查询选择性覆盖查询 分析性能使用数据库分析器评估对数据库的操作使用db.currentOp()评估mongod操作使用explain评估查询性能 优化查询性能创建索…...
【华为】配置NAT访问互联网
1.AR1: int g0/0/0 ip ad 64.1.1.2 255.255.255.0 int g0/0/1 ip ad 110.242.68.1 255.255.255.02.AR2: (1)配置端口ip: int g0/0/1 ip ad 10.3.1.2 255.255.255.0 int g0/0/0 ip ad 64.1.1.1 255.255.255.0(2)配置默认路由: ip route-static 0.0.0.0 0.…...
Spring Boot项目使用多线程执行定时任务
我在一个Spring Boot项目中,采用定时器执行一些操作,比如10秒就发送一次数据。这些操作有2个,如下所示。我就想,虽然这两个操作各自指定了时间频率,但如果其中一个操作非常耗时,会不会影响其他操作呢&#…...
【安装JDK和Android SDK】
安装JDK和Android SDK 1 前言2 下载2.1 下载途径2.2 JDK下载和安装2.2.1 下载2.2.2 安装并配置环境变量2.2.3 验证 2.3 SDK下载和安装2.3.1 下载2.3.2 安装2.3.3 环境变量配置2.3.4 验证 1 前言 在软件开发中,Android应用开发通常使用Android Studio,但…...
汇总10个AI免费一键生成PPT的网站
一、前言 PPT幻灯片是现代办公和学习中的重要组成部分。它在工作、研究或培训中扮演着重要角色,并能够让观众更好地理解信息。随着当今人工智能技术的快速发展,现在有很多免费的AI PPT生成器可供选择,帮助用户更加便捷地制作出高效且具有较强…...
超材料光子晶体和禁带分析实例_CST电磁仿真教程
光子晶体是由周期性排列的不同折射率的介质制造的光学结构,可被视为广义超材料metamaterial的一种。本期我们演示设计一个基于光频能带(PBG,photonics band gap) 的二维光子晶体波导,能带分析方法也可适用于微波波段(EBG,electromagetic band…...
关于OceanBase数据库的poc测试连接经验(by liuhui)
poc客户给了OceanBase数据库实例如下 ob实例: ip:1xx.xx.xx 端口:2883 实例名:obm_xczjj_1_poc#cs_pool_1 用户名:root 密码:xxxxxx 问题出现:根据客户提供的OceanBase数据库配置报错。配置如下 查询数据…...
Docker部署如何修改本地mysql,redis连接信息
要修改数据库 MySQL 和缓存 Redis 的地址为 ruoyi-mysql 和 ruoyi-redis,通常需要在 Spring Boot 项目的配置文件中进行相应的修改。 ### 修改 MySQL 数据库地址为 ruoyi-mysql 1. **在 Spring Boot 项目中找到 application.properties 或 application.yml 文件**…...
PHP中的ReflectionClass常见用法
ReflectionClass是 PHP 中的一个类,它提供了有关类的信息的反射。 使用ReflectionClass可以在运行时获取关于类的各种信息,例如类的名称、方法、属性、注释等。 以下是一些常见的用法: 获取类的名称: $reflection new Reflec…...
processing像素画教程
前提:各位已经安装了processing 第一步:创建一个简单的网格 我们首先创建一个网格来定义我们作品的像素画布。网格将帮助您在适当的位置绘制每个像素。 int gridSize 20; // 每个像素的大小 int cols, rows; void setup() {size(400, 400); // 设置画…...
【秋招笔试】10.13字节跳动(已改编)秋招-三语言题解
🍭 大家好这里是 春秋招笔试突围,一起备战大厂笔试 💻 ACM金牌团队🏅️ | 多次AK大厂笔试 | 大厂实习经历 ✨ 本系列打算持续跟新 春秋招笔试题 👏 感谢大家的订阅➕ 和 喜欢💗 和 手里的小花花🌸 ✨ 笔试合集传送们 -> 🧷春秋招笔试合集 本次的三题全部上线…...
KubeSphere 容器平台高可用:环境搭建与可视化操作指南
Linux_k8s篇 欢迎来到Linux的世界,看笔记好好学多敲多打,每个人都是大神! 题目:KubeSphere 容器平台高可用:环境搭建与可视化操作指南 版本号: 1.0,0 作者: 老王要学习 日期: 2025.06.05 适用环境: Ubuntu22 文档说…...
XCTF-web-easyupload
试了试php,php7,pht,phtml等,都没有用 尝试.user.ini 抓包修改将.user.ini修改为jpg图片 在上传一个123.jpg 用蚁剑连接,得到flag...
OkHttp 中实现断点续传 demo
在 OkHttp 中实现断点续传主要通过以下步骤完成,核心是利用 HTTP 协议的 Range 请求头指定下载范围: 实现原理 Range 请求头:向服务器请求文件的特定字节范围(如 Range: bytes1024-) 本地文件记录:保存已…...
C++ 求圆面积的程序(Program to find area of a circle)
给定半径r,求圆的面积。圆的面积应精确到小数点后5位。 例子: 输入:r 5 输出:78.53982 解释:由于面积 PI * r * r 3.14159265358979323846 * 5 * 5 78.53982,因为我们只保留小数点后 5 位数字。 输…...
微信小程序云开发平台MySQL的连接方式
注:微信小程序云开发平台指的是腾讯云开发 先给结论:微信小程序云开发平台的MySQL,无法通过获取数据库连接信息的方式进行连接,连接只能通过云开发的SDK连接,具体要参考官方文档: 为什么? 因为…...
【论文阅读28】-CNN-BiLSTM-Attention-(2024)
本文把滑坡位移序列拆开、筛优质因子,再用 CNN-BiLSTM-Attention 来动态预测每个子序列,最后重构出总位移,预测效果超越传统模型。 文章目录 1 引言2 方法2.1 位移时间序列加性模型2.2 变分模态分解 (VMD) 具体步骤2.3.1 样本熵(S…...
【碎碎念】宝可梦 Mesh GO : 基于MESH网络的口袋妖怪 宝可梦GO游戏自组网系统
目录 游戏说明《宝可梦 Mesh GO》 —— 局域宝可梦探索Pokmon GO 类游戏核心理念应用场景Mesh 特性 宝可梦玩法融合设计游戏构想要素1. 地图探索(基于物理空间 广播范围)2. 野生宝可梦生成与广播3. 对战系统4. 道具与通信5. 延伸玩法 安全性设计 技术选…...
【VLNs篇】07:NavRL—在动态环境中学习安全飞行
项目内容论文标题NavRL: 在动态环境中学习安全飞行 (NavRL: Learning Safe Flight in Dynamic Environments)核心问题解决无人机在包含静态和动态障碍物的复杂环境中进行安全、高效自主导航的挑战,克服传统方法和现有强化学习方法的局限性。核心算法基于近端策略优化…...
【FTP】ftp文件传输会丢包吗?批量几百个文件传输,有一些文件没有传输完整,如何解决?
FTP(File Transfer Protocol)本身是一个基于 TCP 的协议,理论上不会丢包。但 FTP 文件传输过程中仍可能出现文件不完整、丢失或损坏的情况,主要原因包括: ✅ 一、FTP传输可能“丢包”或文件不完整的原因 原因描述网络…...
《Offer来了:Java面试核心知识点精讲》大纲
文章目录 一、《Offer来了:Java面试核心知识点精讲》的典型大纲框架Java基础并发编程JVM原理数据库与缓存分布式架构系统设计二、《Offer来了:Java面试核心知识点精讲(原理篇)》技术文章大纲核心主题:Java基础原理与面试高频考点Java虚拟机(JVM)原理Java并发编程原理Jav…...
