泊松求和公式推导

1,泊松求和公式

设 x(t)x(t)x(t) 具有傅里叶变换。

频率表示
X(f)=∫−∞∞x(t)e−j2πft dt X(f)=\int_{-\infty}^\infty x(t)\mathrm{e}^{-\mathrm{j}2\pi ft}\mathrm{~d}t X(f)=∫−∞∞x(t)e−j2πft dt

则有
∑n=−∞∞x(n)=∑k=−∞∞X(k)=∑k=−∞∞∫−∞∞x(k)e−j2πkt dt \color{red} \sum_{n=-\infty}^\infty x(n)=\sum_{k=-\infty}^\infty X(k)=\sum_{k=-\infty}^\infty \int_{-\infty}^\infty x(k)\mathrm{e}^{-\mathrm{j}2\pi kt}\mathrm{~d}t n=−∞∑∞x(n)=k=−∞∑∞X(k)=k=−∞∑∞∫−∞∞x(k)e−j2πkt dt
角频率表示
X(jω)=∫−∞∞x(t)e−jωt dt X(\mathrm{j}\omega)=\int_{-\infty}^\infty x(t)\mathrm{e}^{-\mathrm{j}\omega t}\mathrm{~d}t X(jω)=∫−∞∞x(t)e−jωt dt

则有
∑n=−∞∞x(nT)=1T∑k=−∞∞X(j2πkT)=1T∑k=−∞∞∫−∞∞x(t)e−j(2πk/T)t dt \color{red} \sum_{n=-\infty}^\infty x(nT)=\frac{1}{T}\sum_{k=-\infty}^\infty X(\mathrm{j}\frac{2\pi k}{T})=\frac{1}{T}\sum_{k=-\infty}^\infty \int_{-\infty}^\infty x(t)\mathrm{e}^{-\mathrm{j}(2\pi k/T) t}\mathrm{~d}t n=−∞∑∞x(nT)=T1k=−∞∑∞X(jT2πk)=T1k=−∞∑∞∫−∞∞x(t)e−j(2πk/T)t dt

其中 TTT 是采样周期。

2,泊松求和公式证明

2.1,频率公式

构造函数 F(t)=∑n=−∞∞x(t+n)F(t)=\sum\limits_{n=-\infty}^{\infty} x(t+n)F(t)=n=−∞∑∞x(t+n) 。
F(t+1)=∑n=−∞∞x(t+1+n)=∑k=−∞∞x(t+k)=F(t) F(t+1)=\sum_{n=-\infty}^{\infty} x(t+1+n)=\sum_{k=-\infty}^{\infty} x(t+k)=F(t) F(t+1)=n=−∞∑∞x(t+1+n)=k=−∞∑∞x(t+k)=F(t)

因此 F(t)F(t)F(t) 是周期函数,周期 T=1T=1T=1 。

把 F(t)F(t)F(t) 展开成傅里叶级数
F(t)=∑k=−∞+∞akejk(2π/T)t=∑k=−∞+∞akej2πkt F(t)=\sum_{k=-\infty}^{+\infty} a_k \mathrm{e}^{\mathrm{j} k(2 \pi / T) t}=\sum_{k=-\infty}^{+\infty} a_k \mathrm{e}^{\mathrm{j} 2\pi kt} F(t)=k=−∞∑+∞akejk(2π/T)t=k=−∞∑+∞akej2πkt

其中
ak=1T∫TF(t)e−jk(2π/T)t dt=∫01F(t)e−j2πkt dt a_k=\frac{1}{T} \int_T F(t) \mathrm{e}^{-\mathrm{j} k(2 \pi / T) t} \mathrm{~d} t=\int_0^1 F(t) \mathrm{e}^{-\mathrm{j} 2\pi kt} \mathrm{~d} t ak=T1∫TF(t)e−jk(2π/T)t dt=∫01F(t)e−j2πkt dt

将 F(t)=∑n=−∞∞x(t+n)F(t)=\sum\limits_{n=-\infty}^{\infty} x(t+n)F(t)=n=−∞∑∞x(t+n) 代入得到
ak=∫01[∑n=−∞∞x(t+n)]e−j2πkt dt a_k=\int_0^1 [\sum_{n=-\infty}^{\infty} x(t+n)] \mathrm{e}^{-\mathrm{j} 2\pi kt} \mathrm{~d} t ak=∫01[n=−∞∑∞x(t+n)]e−j2πkt dt

交换积分和求和的次序
ak=∑n=−∞∞∫01x(t+n)e−j2πkt dt a_k=\sum_{n=-\infty}^{\infty} \int_0^1 x(t+n) \mathrm{e}^{-\mathrm{j} 2\pi kt} \mathrm{~d} t ak=n=−∞∑∞∫01x(t+n)e−j2πkt dt

令 t+n=mt+n=mt+n=m ,则
ak=∑n=−∞∞∫nn+1x(m)e−j2πk(m−n) dm=∑n=−∞∞∫nn+1x(m)e−j2πkm dm a_k=\sum_{n=-\infty}^{\infty} \int_n^{n+1} x(m) \mathrm{e}^{-\mathrm{j} 2\pi k(m-n)} \mathrm{~d} m=\sum_{n=-\infty}^{\infty} \int_n^{n+1} x(m) \mathrm{e}^{-\mathrm{j} 2\pi km} \mathrm{~d} m ak=n=−∞∑∞∫nn+1x(m)e−j2πk(m−n) dm=n=−∞∑∞∫nn+1x(m)e−j2πkm dm

注意到 aka_kak 是无穷多个区间为 1 的积分之和,可以进行合并
ak=∫−∞∞x(m)e−j2πkm dm a_k=\int_{-\infty}^\infty x(m) \mathrm{e}^{-\mathrm{j} 2\pi km} \mathrm{~d} m ak=∫−∞∞x(m)e−j2πkm dm

因为 x(t)x(t)x(t) 的傅里叶变换为
X(f)=∫−∞∞x(t)e−j2πft dt X(f)=\int_{-\infty}^\infty x(t)\mathrm{e}^{-\mathrm{j}2\pi ft}\mathrm{~d}t X(f)=∫−∞∞x(t)e−j2πft dt

对比上式, ak=X(k)a_k=X(k)ak=X(k) 。

代入傅里叶级数反变换
∑n=−∞∞x(t+n)=F(t)=∑k=−∞+∞akej2πkt=∑k=−∞+∞X(k)ej2πkt \sum_{n=-\infty}^{\infty} x(t+n)=F(t)=\sum_{k=-\infty}^{+\infty} a_k \mathrm{e}^{\mathrm{j} 2\pi kt}=\sum_{k=-\infty}^{+\infty} X(k) \mathrm{e}^{\mathrm{j} 2\pi kt} n=−∞∑∞x(t+n)=F(t)=k=−∞∑+∞akej2πkt=k=−∞∑+∞X(k)ej2πkt

取 t=0t=0t=0 ,得到
∑n=−∞∞x(n)=∑k=−∞+∞X(k) \sum_{n=-\infty}^{\infty} x(n)=\sum_{k=-\infty}^{+\infty} X(k) n=−∞∑∞x(n)=k=−∞∑+∞X(k)

得证。

2.2,角频率公式

思路和频率形式类似。

构造函数 F(t)=∑n=−∞∞x(t+nT)F(t)=\sum\limits_{n=-\infty}^{\infty} x(t+nT)F(t)=n=−∞∑∞x(t+nT) 。则
F(t+T)=∑n=−∞∞x(t+(n+1)T)=∑k=−∞∞x(t+kT)=F(t) F(t+T)=\sum_{n=-\infty}^{\infty} x(t+(n+1)T)=\sum_{k=-\infty}^{\infty} x(t+kT)=F(t) F(t+T)=n=−∞∑∞x(t+(n+1)T)=k=−∞∑∞x(t+kT)=F(t)

因此 F(t)F(t)F(t) 是周期函数,其周期为 TTT 。

F(t)F(t)F(t) 展开成傅里叶级数
F(t)=∑k=−∞+∞akejk(2π/T)t F(t)=\sum_{k=-\infty}^{+\infty} a_k \mathrm{e}^{\mathrm{j} k(2 \pi / T) t} F(t)=k=−∞∑+∞akejk(2π/T)t

其中
ak=1T∫TF(t)e−jk(2π/T)t dt=1T∫0TF(t)e−jk(2π/T)t dt a_k=\frac{1}{T} \int_T F(t) \mathrm{e}^{-\mathrm{j} k(2 \pi / T) t} \mathrm{~d} t=\frac{1}{T} \int_0^T F(t) \mathrm{e}^{-\mathrm{j} k(2 \pi / T) t} \mathrm{~d} t ak=T1∫TF(t)e−jk(2π/T)t dt=T1∫0TF(t)e−jk(2π/T)t dt

将 F(t)=∑n=−∞∞x(t+nT)F(t)=\sum\limits_{n=-\infty}^{\infty} x(t+nT)F(t)=n=−∞∑∞x(t+nT) 代入得到
ak=1T∫0T[∑n=−∞∞x(t+nT)]e−jk(2π/T)t dt a_k=\frac{1}{T} \int_0^T [\sum\limits_{n=-\infty}^{\infty} x(t+nT)] \mathrm{e}^{-\mathrm{j} k(2 \pi / T) t} \mathrm{~d} t ak=T1∫0T[n=−∞∑∞x(t+nT)]e−jk(2π/T)t dt

交换积分和求和的次序
ak=1T∑n=−∞∞∫0Tx(t+nT)e−jk(2π/T)t dt a_k=\frac{1}{T} \sum\limits_{n=-\infty}^{\infty} \int_0^T x(t+nT) \mathrm{e}^{-\mathrm{j} k(2 \pi / T) t} \mathrm{~d} t ak=T1n=−∞∑∞∫0Tx(t+nT)e−jk(2π/T)t dt

令 m=t+nTm=t+nTm=t+nT ,则
ak=1T∑n=−∞∞∫nT(n+1)Tx(m)e−jk(2π/T)(m−nT) dm=1T∑n=−∞∞∫nT(n+1)Tx(m)e−jk(2π/T)m dm a_k=\frac{1}{T} \sum\limits_{n=-\infty}^{\infty} \int_{nT}^{(n+1)T} x(m) \mathrm{e}^{-\mathrm{j} k(2 \pi / T) (m-nT)} \mathrm{~d} m=\frac{1}{T} \sum\limits_{n=-\infty}^{\infty} \int_{nT}^{(n+1)T} x(m) \mathrm{e}^{-\mathrm{j} k(2 \pi / T) m} \mathrm{~d} m ak=T1n=−∞∑∞∫nT(n+1)Tx(m)e−jk(2π/T)(m−nT) dm=T1n=−∞∑∞∫nT(n+1)Tx(m)e−jk(2π/T)m dm

注意到 aka_kak 是无穷多个区间为 TTT 的积分之和,可以进行合并
ak=1T∫−∞∞x(m)e−jk(2π/T)m dm a_k=\frac{1}{T}\int_{-\infty}^\infty x(m) \mathrm{e}^{-\mathrm{j} k(2\pi /T)m} \mathrm{~d} m ak=T1∫−∞∞x(m)e−jk(2π/T)m dm

因为 x(t)x(t)x(t) 的傅里叶变换为
X(jω)=∫−∞∞x(t)e−jωt dt X(\mathrm{j}\omega)=\int_{-\infty}^\infty x(t)\mathrm{e}^{-\mathrm{j}\omega t}\mathrm{~d}t X(jω)=∫−∞∞x(t)e−jωt dt

可见
ak=1TX(j2πkT) a_k=\frac{1}{T}X(\mathrm{j}\frac{2\pi k}{T}) ak=T1X(jT2πk)

代入傅里叶级数反变换
∑n=−∞∞x(t+nT)=F(t)=1T∑k=−∞+∞X(j2πkT)ejk(2π/T)t \sum_{n=-\infty}^{\infty} x(t+nT)=F(t)=\frac{1}{T}\sum_{k=-\infty}^{+\infty} X(\mathrm{j}\frac{2\pi k}{T}) \mathrm{e}^{\mathrm{j} k(2 \pi / T) t} n=−∞∑∞x(t+nT)=F(t)=T1k=−∞∑+∞X(jT2πk)ejk(2π/T)t

取 t=0t=0t=0 得到
∑n=−∞∞x(nT)=1T∑k=−∞+∞X(j2πkT) \sum_{n=-\infty}^{\infty} x(nT)=\frac{1}{T}\sum_{k=-\infty}^{+\infty} X(\mathrm{j}\frac{2\pi k}{T}) n=−∞∑∞x(nT)=T1k=−∞∑+∞X(jT2πk)

得证。

相关推荐
小森77674 天前
几种特殊的数字滤波器---原理及设计
数字信号处理·低通滤波器·滤波器
封奚泽优7 天前
数学七夕花礼(MATLAB版)
开发语言·数学·matlab·七夕·鲜花
花开莫与流年错_11 天前
音频相关数学支持
数学·音视频·音频·软件
MPCTHU15 天前
Advanced Math & Math Analysis |02 Limits
数学
CUC-MenG15 天前
2025杭电多校第十场 Cut Check Bit、Multiple and Factor 个人题解
数学·dp·位运算·数位dp·根号分治
Tisfy17 天前
LeetCode 837.新 21 点:动态规划+滑动窗口
数学·算法·leetcode·动态规划·dp·滑动窗口·概率
CUC-MenG20 天前
2025牛客多校第十场 K.神奇集合 F.老师和Yuuka逛商场 E.老师与好感度 I.矩阵 个人题解
数学·线段树·贪心·dp·线性dp·构造·强联通分量·树上背包·线段树二分
databook21 天前
把数学对象画出来:Manim Mobject类库速查手册
python·数学·动效
CUC-MenG22 天前
2025牛客多校第九场 G.排列 A.AVL树 F.军训 个人题解
数学·dfs·dp·笛卡尔树·组合数·曼哈顿距离·树上dp