北京邮电大学《数字信号处理》门爱东dsp02离散时间系统和离散信号的变换

DigitalSignalProcessing离散时间系统和离散信号的变换北京邮电大学《数字信号处理》精品课程·第二章·门爱东Contents本章内容概览离散时间系统与信号变换的核心理论框架与分析方法01取样与内插：连续信号的离散化02离散时间信号：序列表示与基本运算03离散系统基本概念与差分方程04离散时间傅里叶变换（DTFT）05Z变换与系统频域分析Chapter01取样与内插从连续信号到离散信号的桥梁——模拟世界与数字世界的转换机制SamplingTheorem取样定理与奈奎斯特条件取样是连续信号数字化的第一步。奈奎斯特-香农取样定理规定：对带限信号，取样频率fs必须满足fs≥2fmax才能无失真恢复原信号。违反此条件将导致频谱混叠，即高频分量折叠到低频段造成不可逆失真。CONDITION奈奎斯特条件取样频率fs≥2fmax（信号最高频率），称为奈奎斯特速率，是信号无失真恢复的最低要求满足条件时，可通过理想低通滤波器从取样信号中精确重建原始连续信号fs≥2fmaxPHENOMENON频谱混叠现象当fs<2fmax时，取样信号的频谱副本发生交叠，高频分量折叠至低频段，产生混叠失真电影中车轮"倒转"现象即为典型混叠：24帧/秒拍摄超过奈奎斯特极限的高速旋转物体fs<2fmaxPRACTICE工程实践意义实际取样前必须使用抗混叠滤波器，将信号带宽限制在fs/2以内，防止高频干扰混入CD音频取样率44.1kHz，恰好覆盖人耳听觉上限20kHz的奈奎斯特要求44.1kHzSamplingStrategy欠取样与过取样策略实际工程中，取样策略可灵活优化。带通信号可采用欠取样（按带宽两倍而非最高频率两倍取样），大幅降低取样率要求；过取样则以超高采样率换取抗混叠滤波器简化、量化噪声降低和信噪比提升，是高精度ADC设计的核心技术。欠取样（Undersampling）带通信号（频率范围B₁~B₂）无需按2B₂取样，只需取样率≥2(B₂−B₁)即可避免频谱混叠GSM蜂窝信号在900MHz频段占30kHz带宽，欠取样仅需略高于60kHz而非1.8GHz的取样频率≥2Δf取样率下限过取样（Oversampling）用远高于奈奎斯特率的频率取样，使量化噪声分布在更宽频带内，有效降低信号带内噪声功率结合噪声整形（NoiseShaping）技术将量化噪声搬移到信号带宽之外，显著提升信噪比SNR↑核心收益SignalReconstruction信号重建与内插方法内插是从离散序列恢复连续信号的过程，是取样的逆操作。理想内插使用sinc函数可实现数学上的完美重建，但因其无限长冲激响应在物理上不可实现。工程中采用零阶保持器（ZOH）或线性内插等近似方法，配合平滑滤波器实现高质量信号恢复。理想内插y(t)=Σx(n)·sinc(t/T−n)，数学上可精确重建带限信号，但sinc函数无限长，物理不可实现。Sinc零阶保持器将每个采样值保持至下一采样时刻，形成阶梯状重建波形，是DAC的基本工作原理。ZOH一阶保持器用直线连接相邻采样点，重建波形更平滑但高频分量更大，需配合低通滤波器使用。线性内插取样率关联取样率越高，相邻采样点间隔越小，任何内插方法的近似误差越小，重建质量越高。Fs↑SAMPLINGTHEOREMINPRACTICE取样定理的工程应用实例取样定理不仅是理论基石，更在工程实践中无处不在。从电影拍摄的帧频选择到通信系统的信号数字化，取样率的合理设定直接决定了信号保真度。电影帧频与车轮混叠54km/h混叠临界01电影以16帧/秒拍摄时，奈奎斯特极限为8Hz旋转频率，对应车速约54km/h，超过即产生"倒转"混叠。早期电影常出现马车车轮看似倒转的经典现象，正是取样率不足导致频率混叠的直观例证。02高速摄像机（如1000帧/秒）通过大幅提高取样率，可清晰记录汽车、跳水、乒乓球等高速运动细节，有效扩展奈奎斯特带宽，避免运动模糊与混叠失真。关键洞察：提高取样率=扩展可观测带宽通信系统中的取样策略44.1kHzCD取样标准01GSM系统在900MHz载频上传输30kHz带宽信号，采用欠取样技术仅需60kHz取样率，大幅降低ADC成本与功耗。通过带通取样原理，可用远低于载频的速率完成信号数字化。02音频CD采用44.1kHz取样率，满足人耳20kHz上限的奈奎斯特要求（>2×20kHz），同时留出过渡带给抗混叠滤波器，兼顾音质与存储效率。关键洞察：欠取样降低硬件复杂度，过取样提升信号质量Chapter02离散时间信号序列的数学表示、常用信号类型与基本运算规则SIGNALREPRESENTATION离散时间信号的表示方法离散时间信号是在离散时间点上定义的数值序列，用x(n)表示，n为整数。它有四种等价表示方式——公式表示、图形表示、集合枚举和数值向量表示，每种方式在不同分析场景中各有优势，是数字信号处理最基础的数据结构。公式表示x(n)=aⁿu(n)等数学表达式，便于理论推导和性质证明，是分析的核心工具aⁿu(n)图形表示用n轴上的竖线标注各时刻幅值，直观展示信号的时域波形特征脉冲图集合枚举x(n)={...0,1,3,5,3,1,0...}，以列表形式给出有限长序列的全部取值{…0,1,3,5…}数值向量表示在MATLAB中用数组存储序列值，便于计算机进行数值计算和仿真分析MATLABDISCRETESEQUENCES常用离散序列离散信号分析中有一组"基本序列"构成信号表示的基础，任何复杂序列均可分解为基本序列的组合。单位脉冲与阶跃序列针"，任何序列可表示为x(n)=Σx(k)δ(n−k)u(n)=Σδ(k)，是系统的标准阶跃激励δ(n)·u(n)矩形与指数序列号段，是窗函数法滤波的基础系统稳定性判断的核心信号RN(n)·aⁿ正弦与复指数序列为有理数，与连续正弦有本质差异信号ejΩt在离散域的映射DTFT基函数Discrete-TimeSignalOperations序列的基本运算序列运算构成离散信号处理的数学工具箱。移位和翻转改变信号在时间轴上的位置和方向，尺度变换实现信号的压缩与扩展，加法和乘法是逐点运算，而卷积运算则是线性时不变（LTI）系统分析的核心——系统输出等于输入信号与系统冲激响应的卷积。移位、翻转与尺度变换改变信号在时间轴上的位置、方向与尺度加法、乘法与卷积逐点运算与LTI系统核心分析工具移位x(n−n₀)序列整体右移n₀个单位（n₀>0），是时延操作的离散形式，在滤波器设计中广泛使用加法与乘法y(n)=x₁(n)+x₂(n)与y(n)=x₁(n)·x₂(n)，逐点运算，分别用于信号叠加和调制/加窗翻转x(−n)序列关于n=0轴镜像，常与移位组合为x(n₀−n)，在卷积计算中是必要步骤卷积y(n)=Σx(k)h(n−k)LTI系统输出等于输入x(n)与冲激响应h(n)的卷积和，是系统分析的核心运算PERIODICITYCRITERIA离散序列的周期性条件离散正弦序列的周期性取决于数字频率与2π之比是否为有理数，这一整数约束本质决定了DFT频率分辨率的设计基础。周期判据sin(ω₀n)为周期序列的充要条件是ω₀/2π=P/N（P、N为互质整数），此时序列周期为N。该判据是分析离散信号周期性的核心理论依据。P/N示例分析sin(πn/3)中ω₀/2π=1/6，周期N=6；sin(n)中ω₀/2π为无理数，不构成周期序列。通过具体算例验证周期判据的应用方法。N=6频率映射数字频率与模拟频率的关系为ω₀=Ω₀T=2πf₀/fs，取样周期T将模拟频率映射到数字域。该映射关系是连接连续与离散信号的关键桥梁。ω₀=Ω₀T工程意义DFT要求信号具有周期性或可通过加窗截断为整周期，否则将产生频谱泄漏现象。理解周期条件对数字信号处理系统设计具有重要指导价值。频谱泄漏CHAPTER03离散系统基本概念与差分方程系统分类、线性时不变系统、差分方程描述与信号流图DISCRETESYSTEMS离散系统的四个基本特性离散系统的特性可从四个维度判定：线性（叠加原理）、时不变性（参数不随时间变化）、因果性（输出不超前输入）和稳定性（有界输入产生有界输出）。同时满足线性和时不变性的LTI系统是数字信号处理的核心分析对象，其输入输出关系完全由冲激响应h(n)唯一确定。01LINEARITY线性（叠加原理）若T[x₁]=y₁、T[x₂]=y₂，则T[ax₁+bx₂]=ay₁+by₂，叠加原理成立，系统可分解为基本响应的组合。T[ax₁+bx₂]=ay₁+by₂02TIME-INVARIANCE时不变性（LSI）若T[x(n)]=y(n)，则T[x(n−n₀)]=y(n−n₀)，输入延迟n₀步则输出也延迟n₀步。x(n−n₀)→y(n−n₀)03CAUSALITY因果性y(n)仅取决于x(n)、x(n−1)等当前及过去输入，不包含未来值，实时系统必须满足因果性。y(n)←x(k),k≤n04BIBOSTABILITYBIBO稳定性对任何有界输入|x(n)|≤M，输出|y(n)|也有限；LTI系统稳定的充要条件是Σ|h(n)|<∞。Σ|h(n)|<∞DISCRETE-TIMESYSTEMS线性常系数差分方程线性常系数差分方程是LTI离散系统的核心数学模型，系数决定了系统的频率响应、稳定性和实现结构，递归型与非递归型存在根本差异。01一般形式y(n)+a₁y(n-1)+⋯+aₙy(n-N)=b₀x(n)+b₁x(n-1)+⋯+bₘx(n-M)，N阶方程描述N阶系统02递归与非递归aₖ不全为零时为递归（IIR）系统，含输出反馈；所有aₖ=0时为非递归（FIR）系统03求解方法迭代法适合数值计算；经典法求齐次解+特解得闭式表达；Z变换法可处理任意初始条件04系数与系统特性映射aₖ和bᵣ决定传输函数H(z)的零极点分布，从而决定频率响应和稳定性Section·SignalFlowGraph信号流图表示信号流图是差分方程的图形化等价表示，用节点、支路和延迟器三种基本元素构建。它将抽象的数学方程转化为直观的拓扑结构，便于分析系统的实现复杂度、数值精度和硬件资源需求。信号流图的转置定理保证了同一传输函数可用多种等价结构实现。01基本构成元素节点（加法器）多条输入支路汇聚于此，输出为各输入之和，对应差分方程中的求和操作∑求和02基本构成元素单位延迟器（z⁻¹）将信号延迟一个采样周期，对应差分方程中y(n-1)、x(n-1)等时延项z⁻¹延迟03流图分析与转置支路增益标注在连线上的系数值，对应差分方程中的aₖ和bᵣ系数，决定信号缩放比例aₖ·bᵣ系数04流图分析与转置转置定理将流图所有支路反向、输入输出互换，得到的新流图具有相同的传输函数，但数值特性不同H(z)等价CHAPTER04离散时间傅里叶变换（DTFT）离散信号的频域分析工具——定义、性质与系统频率响应DTFT的定义与存在条件离散时间傅里叶变换DTFT将离散序列映射为数字频率ω的连续函数X(e^jω)，是分析离散信号频谱的核心工具。正变换为无穷级数求和，反变换为一个周期内的积分。DTFT具有2π周期性——这是离散信号频谱最本质的特征，源于时间离散化导致的频谱周期延拓。变换对定义·正变换X(e^jω)=Σx(n)e^(−jωn)，将时域序列变换为数字频率ω的连续复函数，求和范围为n∈(−∞,+∞)。该变换建立了离散时间域与连续频域之间的数学桥梁。变换对定义·反变换x(n)=(1/2π)∫X(e^jω)e^(jωn)dω，从频域恢复时域序列，积分只需一个周期[−π,π]。正逆变换构成完整的分析合成对。存在条件充分条件：序列绝对可加Σ|x(n)|<∞，保证级数一致收敛。典型满足条件的序列包括衰减指数序列aⁿu(n)（|a|<1）、有限长序列等。2π周期性X(e^jω)=X(e^j(ω+2π))，这是DTFT最根本的特征，源于离散时间变量的整数性质。周期性意味着频谱信息完全包含在任意2π区间内，通常取主值区间[−π,π]进行分析。PROPERTIESDTFT的核心性质DTFT的性质体系将时域运算与频域运算建立对应关系。其中卷积定理是数字滤波器设计的理论基石——使系统处理从复杂卷积运算转化为频域简单乘法。01线性DTFT[ax₁+bx₂]=aX₁+bX₂叠加原理在频域同样成立，信号的线性组合对应频谱的线性叠加，是分析复杂信号的基础LINEARITY02时移x[n-n₀]↔e⁻ʲωⁿ⁰X(eʲω)时域延迟产生线性相位偏移，幅度谱保持不变，这一特性在系统延迟分析和同步处理中至关重要TIMESHIFT03频移eʲω₀ⁿx[n]↔X(eʲ⁽ω⁻ω₀⁾)时域乘以复指数实现频谱搬移，是调制解调技术的数学基础，广泛应用于通信系统的频分复用FREQSHIFT04卷积定理x[n]*h[n]↔X(eʲω)H(eʲω)时域卷积等于频域相乘，将复杂的卷积运算转化为简单的乘法运算，是LTI系统分析与设计的核心工具CONVOLUTION05帕塞瓦尔Σ|x[n]|²=(1/2π)∫|X|²dω时域总能量等于频域总能量，体现了信号在时频域的能量守恒特性，是功率谱估计的理论依据PARSEVALChapter20·DigitalSignalProcessingLTI系统的频率响应频率响应H(ejω)是LTI系统在频域中的完整描述，等于冲激响应h(n)的DTFT。它由幅度响应和相位响应两部分组成：幅度响应决定各频率分量的增益或衰减，相位响应决定各频率分量的时间延迟。滤波器设计的本质就是构造具有特定频率响应的系统。01频率响应定义H(ejω)=DTFT[h(n)]=Σh(n)e−jωn，系统冲激响应的DTFT即为频率响应频域分析的核心工具02幅度响应|H(ejω)|描述系统对各频率正弦分量的放大或衰减程度，是滤波器通阻带特性的直接体现决定滤波器选频特性03相位响应arg[H(ejω)]描述系统对各频率分量引入的相移，线性相位保证信号波形不失形影响信号时域波形保真度04群延迟τ(ω)=−d[argH]/dω表示信号包络的延迟时间，恒定群延迟意味着所有频率分量同步到达衡量系统时延特性的关键指标DISCRETE-TIMEFOURIERTRANSFORMDTFT计算实例与MATLAB可视化通过指数序列x(n)=(0.5)ⁿu(n)的DTFT计算，展示从公式推导到数值可视化的完整分析流程。解析求解过程Step01·定义代入代入DTFT定义式：X(ejω)=Σ(n=0→∞)(0.5)ⁿe−jωn，识别为等比级数形式Step02·级数求和等比级数求和：首项1，公比0.5e−jω，得X(ejω)=1/(1−0.5e−jω)，收敛条件恒成立Step03·结果验证验证收敛性：|0.5e−jω|=0.5<1，级数绝对收敛，DTFT存在且唯一确定|0.5e−jω|=0.5<1✓收敛域：全频域MATLAB可视化实现Step01·离散采样将ω∈[0,π]分为501个等间隔点，计算|X|、arg(X)、Re(X)、Im(X)四项频谱分量Step02·幅度特性幅度谱在ω=0处取最大值2，随ω增大单调递减至ω=π处的2/3，体现低通特性Step03·相位分析相位谱从0开始非线性变化，在截止频率附近变化剧烈，高频段趋于稳定值|X|max=2ω=0SYMMETRYPROPERTIESDTFT的对称性质实序列的DTFT具有共轭对称性，幅度谱为偶函数、相位谱为奇函数，使频谱分析只需考察[0,π]半周期即可。ConjugateSymmetry共轭对称性实序列x(n)的DTFT满足X(ejω)=X*(e−jω)，正负频率处频谱互为共轭。X(ejω)=X*(e−jω)EvenMagnitude幅度谱偶对称|X(ejω)|=|X(e−jω)|，幅度谱关于ω=0对称，只需分析[0,π]半周期。|X(ejω)|=|X(e−jω)|OddPhase相位谱奇对称arg[X(ejω)]=−arg[X(e−jω)]，相位谱关于原点反对称。arg[X]=−arg[X*]Decomposition序列分解x(n)=xe(n)+xo(n)，共轭对称部分xe对应DTFT实部，反对称部分xo对应虚部除以j。xe→Re,xo→Im/jCHAPTER05Z变换与系统频域分析从DTFT到Z变换——复频域分析工具与系统零极点理论Z-TransformZ变换的定义与收敛域（ROC）Z变换X(z)=Σx(n)z−n是DTFT的推广，通过复变量z将频谱分析从单位圆扩展到整个复平面。收敛域ROC是Z变换理论的核心。01Z变换定义双边Z变换X(z)=Σn=−∞∞x(n)·z−nz为复变量，级数在特定z域内收敛时变换存在单边Z变换X(z)=Σn=0∞x(n)·z−n仅对因果序列有效，适用于求解带初始条件的差分方程02收敛域ROC的分类右边序列（因果）ROC为|z|>R₁的圆外区域，包含z=∞例anu(n)→ROC:|z|>|a|左边序列（反因果）ROC为|z|<R₂的圆内区域，包含z=0例−bnu(−n−1)→ROC:|z|<|b|Chapter04·Z-TransformZ反变换的求解方法Z反变换是从X(z)恢复时域序列x(n)的过程，有三种主要求解方法。部分分式法最实用，长除法适合数值解，留数法适合理论推导。Method01部分分式法将X(z)/z分解为简单分式之和，如A/(z−p₁)+B/(z−p₂)，每项对应已知的Z变换对，逐项求反变换后叠加关键步骤：先除z再分解→确定各极点对应的留数系数→根据ROC确定每项对应因果/反因果序列最常用工程实践首选方法Method02·03长除法与留数法长除法：按z⁻¹的幂次做多项式除法，商的系数即为x(0)、x(1)、x(2)…适合求数值解的前几项留数法：x(n)=(1/2πj)∮X(z)zⁿ⁻¹dz，用留数定理计算围线积分，x(n)=ΣRes[X(z)zⁿ⁻¹,z=pₖ]数值解长除法理论推导留数法PROPERTIESZ变换的核心性质Z变换的性质体系在继承DTFT线性、时移、卷积等基本性质的同时，新增了初值定理、终值定理和Z域尺度变换等独特工具。这些性质使得我们可以在不做反变换的情况下直接获取序列特征值，并通过代数运算快速推导复杂信号的Z变换表达式。01线性：Z[ax₁(n)+bx₂(n)]=aX₁(z)+bX₂(z)，ROC至少为两个ROC的交集，系数组合可能扩大收敛域02时移：Z[x(n−n₀)]=z⁻ⁿ⁰X(z)，延迟n₀步对应乘以z⁻ⁿ⁰，在z=0或z=∞处可能增减极点03Z域尺度变换：Z[aⁿx(n)]=X(z/a)，序列乘以指数aⁿ等价于z平面上的尺度缩放04初值定理：因果序列x(0)=lim[z→∞]X(z)，无需做反变换即可直接获取序列初始值05终值定理：若(z−1)X(z)在单位圆上及圆外解析，则x(∞)=lim[z→1](z−1)X(z)，可预判序列稳态值MappingRelationsZ变换与拉氏变换、傅里叶变换的关系Z变换、拉普拉斯变换和傅里叶变换通过理想取样建立统一的数学框架。核心映射关系z=e^(sT)将s平面映射到z平面：虚轴映射为单位圆，左半s平面映射为单位圆内部。这一映射揭示了连续与离散系统稳定性的等价关系。核心关系s平面到z平面的映射公式核心关系z=esT：s=σ+jΩ映射为z=eσT·ejΩT，模|z|=eσT由σ决定，相角arg(z)=ΩT。该映射是连续系统与离散系统相互转换的数学桥梁，建立了s域与z域的一一对应关系。区域映射s平面三区域对应关系s左半平面(σ<0)→单位圆内(|z|<1)；s虚轴(σ=0)→单位圆(|z|=1)；s右半平面(σ>0)→单位圆外(|z|>1)。区域映射直接决定了连续系统稳定性与离散系统稳定性的等价判据。DTFT特例单位圆上的Z变换即为DTFTZ变换在单位圆上(|z|=1)的取值即为DTFT：X(ejω)=X(z)|z=e^jω，DTFT是Z变换的特例。当序列绝对可和时，Z变换在单位圆上收敛，此时频域特性完全由DTFT刻画。频谱延拓傅里叶变换取样后的周期延拓傅里叶变换经取样后频谱周期延拓，与DTFT的关系为X(ejω)=(1/T)ΣXa(j(ω/T−2πk/T))。该关系体现了时域离散化导致频域周期化的对偶特性，是采样定理的频域表达。Z-DOMAINANALYSIS传输函数H(z)与零极点分析传输函数H(z)=Y(z)/X(z)是LTI系统在Z域的完整描述，其零极点分布决定系统的频率选择特性。01传输函数H(z)=Σh(n)z−n=B(z)/A(z)，分子根为零点{zk}，分母根为极点{pk}，构成零极点图02频率响应提取H(ejω)=H(z)|z=e^jω，在单位圆上求值即得频率响应，极点靠近单位圆处产生谐振峰03稳定性判据因果系统BIBO稳定的充要条件是所有极点|pk|<1，即极点全部位于单位圆内04最小相位系统所有零极点均在单位圆内的系统具有最小群延迟，在相同幅度响应下相位失真最小Z-TransformMethod用单边Z变换解线性差分方程单边Z变换是求解带初始条件线性差分方程的系统化方法，通过移位性质自动纳入初始条件，获得闭式解。右移性质Z⁺[y(n-1)]=z⁻¹Y(z)+y(-1)，自动包含初始值z⁻¹·Y(z)初始条件二阶差分方程需要y(-1)和y(-2)两个初始条件y(-1),y(-2)取Z变换对差分方程两边取单边Z变换，代入移位性质得代数方程时域→Z域分离求解分离零输入响应与零状态响应，反变换得闭式解零输入·零状态Z-PLANEANALYSIS零极点与频率响应的几何解释零极点矢量法提供直观的图解方法：沿单位圆移动时，幅度等于零点矢量积/极点矢量积，极点越近单位圆谐振越尖锐，零点在圆上产生完全陷波。01幅度响应的几何意义：|H(ejω)|=K·∏|ejω−zₖ|/∏|ejω−pₖ|，即单位圆上点到零点的距离积除以到极点的距离积