泊松求和公式推导

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)

得证。

相关推荐
canonical_entropy2 天前
组合为什么优于继承:从工程实践到数学本质
后端·数学·设计模式
EQUINOX14 天前
切比雪夫求和不等式,离散形式,连续形式
数学
小鱼_yu5 天前
探索C#中LINQ的异步流处理使用`IAsyncEnumerable`提升数据查询效率
数学
2401_841495646 天前
【数值分析】插值法实验
python·数学·算法·可视化·数值分析·数学原理·插值法
花子の水晶植轮daisuki8 天前
数论上
数学·数论
Gohldg8 天前
C++算法·贪心例题讲解
c++·数学·算法·贪心算法
Juan_201214 天前
P1447题解
c++·数学·算法·题解
Juan_201215 天前
P3051题解
c++·数学·算法·题解
2401_8414956417 天前
【机器学习】朴素贝叶斯法
人工智能·python·数学·算法·机器学习·概率论·朴素贝叶斯法
大卫小东(Sheldon)1 个月前
写了一个BBP算法的实现库,欢迎讨论
数学·rust