《数字信号处理》学习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大厂笔试 | 大厂实习经历 ✨ 本系列打算持续跟新 春秋招笔试题 👏 感谢大家的订阅➕ 和 喜欢💗 和 手里的小花花🌸 ✨ 笔试合集传送们 -> 🧷春秋招笔试合集 本次的三题全部上线…...

51c自动驾驶~合集58
我自己的原文哦~ https://blog.51cto.com/whaosoft/13967107 #CCA-Attention 全局池化局部保留,CCA-Attention为LLM长文本建模带来突破性进展 琶洲实验室、华南理工大学联合推出关键上下文感知注意力机制(CCA-Attention),…...
DeepSeek 赋能智慧能源:微电网优化调度的智能革新路径
目录 一、智慧能源微电网优化调度概述1.1 智慧能源微电网概念1.2 优化调度的重要性1.3 目前面临的挑战 二、DeepSeek 技术探秘2.1 DeepSeek 技术原理2.2 DeepSeek 独特优势2.3 DeepSeek 在 AI 领域地位 三、DeepSeek 在微电网优化调度中的应用剖析3.1 数据处理与分析3.2 预测与…...

Vue3 + Element Plus + TypeScript中el-transfer穿梭框组件使用详解及示例
使用详解 Element Plus 的 el-transfer 组件是一个强大的穿梭框组件,常用于在两个集合之间进行数据转移,如权限分配、数据选择等场景。下面我将详细介绍其用法并提供一个完整示例。 核心特性与用法 基本属性 v-model:绑定右侧列表的值&…...

visual studio 2022更改主题为深色
visual studio 2022更改主题为深色 点击visual studio 上方的 工具-> 选项 在选项窗口中,选择 环境 -> 常规 ,将其中的颜色主题改成深色 点击确定,更改完成...
大模型多显卡多服务器并行计算方法与实践指南
一、分布式训练概述 大规模语言模型的训练通常需要分布式计算技术,以解决单机资源不足的问题。分布式训练主要分为两种模式: 数据并行:将数据分片到不同设备,每个设备拥有完整的模型副本 模型并行:将模型分割到不同设备,每个设备处理部分模型计算 现代大模型训练通常结合…...

python执行测试用例,allure报乱码且未成功生成报告
allure执行测试用例时显示乱码:‘allure’ �����ڲ����ⲿ���Ҳ���ǿ�&am…...
C++课设:简易日历程序(支持传统节假日 + 二十四节气 + 个人纪念日管理)
名人说:路漫漫其修远兮,吾将上下而求索。—— 屈原《离骚》 创作者:Code_流苏(CSDN)(一个喜欢古诗词和编程的Coder😊) 专栏介绍:《编程项目实战》 目录 一、为什么要开发一个日历程序?1. 深入理解时间算法2. 练习面向对象设计3. 学习数据结构应用二、核心算法深度解析…...

LabVIEW双光子成像系统技术
双光子成像技术的核心特性 双光子成像通过双低能量光子协同激发机制,展现出显著的技术优势: 深层组织穿透能力:适用于活体组织深度成像 高分辨率观测性能:满足微观结构的精细研究需求 低光毒性特点:减少对样本的损伤…...

手机平板能效生态设计指令EU 2023/1670标准解读
手机平板能效生态设计指令EU 2023/1670标准解读 以下是针对欧盟《手机和平板电脑生态设计法规》(EU) 2023/1670 的核心解读,综合法规核心要求、最新修正及企业合规要点: 一、法规背景与目标 生效与强制时间 发布于2023年8月31日(OJ公报&…...
es6+和css3新增的特性有哪些
一:ECMAScript 新特性(ES6) ES6 (2015) - 革命性更新 1,记住的方法,从一个方法里面用到了哪些技术 1,let /const块级作用域声明2,**默认参数**:函数参数可以设置默认值。3&#x…...