信号积分变换的公式推导以及几个例子

置顶 papaofdoudou

已于 2023-09-25 12:39:06 修改

阅读量4.4k

点赞数 9

分类专栏：数学工程文章标签：概率论算法微积分

于 2021-02-09 15:40:08 首次发布

本文链接：https://blog.csdn.net/tugouxp/article/details/113770573

版权

数学同时被 2 个专栏收录

147 篇文章 19 订阅

订阅专栏

工程

52 篇文章 1 订阅

订阅专栏

如果一个以为周期的函数 f(x) 在

$\bigg[-\frac{T}{2}, \frac{T}{2}\bigg]$

上满足狄利克雷条件，即:

1.除去有限个第一类间断点外，处处连续

2.分段单调，单调区间的个数有限

则 f(x) 的fourier级数表示为：

$f(t)\approx \frac{a_0}{2}+\sum_{n=1}^{\infty }\bigg[a_ncos(n\omega_0t) + b_nsin(n\omega_0t)\bigg]$

在

$\bigg[-\frac{T}{2}, \frac{T}{2}\bigg]$

上处处收敛，且在 f(x) 的连续点处收敛于 f(t) , 其中，

$\omega_0 =\frac{2\pi}{T}$

为基波频率。对上式两边求积分：

$\\ \int_{-\frac{T}{2}}^{\frac{T}{2}}f(t)dt=\int_{-\frac{T}{2}}^{\frac{T}{2}}\bigg(\frac{a_0}{2}+\sum_{n=1}^{\infty }\bigg[a_ncos(n\omega_0t) + b_nsin(n\omega_0t)\bigg]\bigg) dt=\int_{-\frac{T}{2}}^{\frac{T}{2}}\frac{a_0}{2}dt\\=\frac{a_0}{2}\cdot T$

所以:

$a_0=\frac{2}{T}\int^{\frac{T}{2}}_{-\frac{T}{2}}f(t)dt$

对于 $\\ (m=1,2,3,\cdots)$

$\\ \int_{-\frac{T}{2}}^{\frac{T}{2}}f(t)cos(m\omega_0t)dt=\int_{-\frac{T}{2}}^{\frac{T}{2}}\bigg(\frac{a_0}{2}+\sum_{n=1}^{\infty }\bigg[a_ncos(n\omega_0t) + b_nsin(n\omega_0t)\bigg]\bigg)cos(m\omega_0t) dt\\=0+a_n\int_{-\frac{T}{2}}^{\frac{T}{2}}cos^2(n\omega_0t)dt=a_n\int_{-\frac{T}{2}}^{\frac{T}{2}}\frac{1+cos(2n\omega_0t)}{2}dt=a_n\int_{-\frac{T}{2}}^{\frac{T}{2}}\frac{1}{2}dt +a_n\int_{-\frac{T}{2}}^{\frac{T}{2}} \frac{cos(2n\omega_0 t)}{2}dt\\=a_n\int_{-\frac{T}{2}}^{\frac{T}{2}}\frac{1}{2}dt +a_n\int_{-\frac{T}{2}}^{\frac{T}{2}} \frac{cos(2n\omega_0 t)}{4n\omega_0}d(2n\omega_0 t)=\frac{a_n}{2}\bigg|_{-\frac{T}{2}}^{\frac{T}{2}}+\frac{a_n}{4n\omega_0}sin(2n\omega_0t)\bigg|_{-\frac{T}{2}}^{\frac{T}{2}}\\=a_n\cdot \frac{T}{2}+0$

所以:

$\\a_n=\frac{2}{T}\int^{\frac{T}{2}}_{-\frac{T}{2}}f(t)cos(n\omega_0 t)dt \quad(n=1,2, 3,\cdots)$

$\\ \int_{-\frac{T}{2}}^{\frac{T}{2}}f(t)sin(m\omega_0t)dt=\int_{-\frac{T}{2}}^{\frac{T}{2}}\bigg(\frac{a_0}{2}+\sum_{n=1}^{\infty }\bigg[a_ncos(n\omega_0t) + b_nsin(n\omega_0t)\bigg]\bigg)sin(m\omega_0t) dt\\=0+b_n\int_{-\frac{T}{2}}^{\frac{T}{2}}sin^2(n\omega_0t)dt=b_n\int_{-\frac{T}{2}}^{\frac{T}{2}}\frac{1-cos(2n\omega_0t)}{2}dt=b_n\int_{-\frac{T}{2}}^{\frac{T}{2}}\frac{1}{2}dt -b_n\int_{-\frac{T}{2}}^{\frac{T}{2}} \frac{cos(2n\omega_0 t)}{2}dt\\=b_n\int_{-\frac{T}{2}}^{\frac{T}{2}}\frac{1}{2}dt - b_n\int_{-\frac{T}{2}}^{\frac{T}{2}} \frac{cos(2n\omega_0 t)}{4n\omega_0}d(2n\omega_0 t)=\frac{b_n}{2}\bigg|_{-\frac{T}{2}}^{\frac{T}{2}}-\frac{b_n}{4n\omega_0}sin(2n\omega_0t)\bigg|_{-\frac{T}{2}}^{\frac{T}{2}}\\=b_n\cdot \frac{T}{2}+0$

所以：

$\\b_n=\frac{2}{T}\int^{\frac{T}{2}}_{-\frac{T}{2}}f(t)sin(n\omega_0 t)dt \quad(n=1,2, 3,\cdots)$

综上：

$a_0=\frac{2}{T}\int^{\frac{T}{2}}_{-\frac{T}{2}}f(t)dt$

$\\a_n=\frac{2}{T}\int^{\frac{T}{2}}_{-\frac{T}{2}}f(t)cos(n\omega_0 t)dt \quad(n=1,2, 3,\cdots)$

$\\b_n=\frac{2}{T}\int^{\frac{T}{2}}_{-\frac{T}{2}}f(t)sin(n\omega_0 t)dt \quad(n=1,2, 3,\cdots)$

在电子通信领域，常常利用欧拉公式：

$cos(t)=\frac{e^{it}+e^{-it}}{2}$

$sin(t)=\frac{e^{it}-e^{-it}}{2i}$

所以：

$\\f(t)\approx \frac{a_0}{2}+\sum_{n=1}^{\infty }\bigg[a_ncos(n\omega_0t) + b_nsin(n\omega_0 t)\bigg]\\=\frac{a_0}{2}+\sum_{n=1}^{\infty }\bigg[ \frac{a_n}{2}(e^{in\omega_0 t}+e^{-in\omega_0 t})-i\cdot \frac{b_n}{2}(e^{in\omega_0 t}-e^{-in\omega_0 t})\bigg]\\=\frac{a_0}{2}+\sum_{n=1}^{\infty }\bigg[ \frac{a_n-ib_n}{2}e^{in\omega_0 t}+\frac{a_n+ib_n}{2}e^{-in\omega_0 t}\bigg]$

令：

$c_0=\frac{a_0}{2}$

$c_n=\frac{a_n-ib_n}{2}$

$c_{-n}=\frac{a_n+ib_n}{2}$

得到fourier级数的复指数形式：

$\\f(t)\approx\frac{a_0}{2}+\sum_{n=1}^{\infty }\bigg[ \frac{a_n-ib_n}{2}e^{in\omega_0 t}+\frac{a_n+ib_n}{2}e^{-in\omega_0 t}\bigg]\\= c_0+\sum_{n=1}^{\infty }\bigg(c_ne^{in\omega_0 t} + c_{-n}e^{-in\omega_0 t}\bigg)=\sum_{n=-\infty}^{\infty}C_ne^{in\omega_0 t}$

从公式中看出，时域信号可以看成是复指数周期基波信号

$y(t)=e^{i\omega_0t}$

以及其倍频信号：

$y_n(t)=e^{in\omega_0t}$

的线性组合。根据欧拉定理，基波复信号可以分解成实步的余弦信号和虚部的正弦信号，由于正弦信号是奇函数，在整数范围内的累加互相抵消为0，只剩下了余弦信号：

每一个高次谐波信号对T来说都是周期的，因此，一个由成谐波关系的复指数信号的现行组合形成的信号f(t)对T来说也是周期的。上式f(t)称为周期信号的傅立叶级数表述，这里面的组合系数为：

$C_0=\frac{a_0}{2}=\frac{1}{2}\cdot \frac{2}{T}\int^{\frac{T}{2}}_{-\frac{T}{2}}f(t)dt=\frac{1}{T}\int^{\frac{T}{2}}_{-\frac{T}{2}}f(t)dt$

$\\C_n=\frac{a_n-ib_n}{2}=\frac{1}{2}\cdot \frac{2}{T}\bigg[\int^{\frac{T}{2}}_{-\frac{T}{2}}f(t)cos(n\omega_0 t)dt-i\int^{\frac{T}{2}}_{-\frac{T}{2}}f(t)sin(n\omega _0 t)dt\bigg] \\=\frac{1}{T}\int^{\frac{T}{2}}_{-\frac{T}{2}}f(t)\bigg[cos(n\omega_0 t)-isin(n\omega_0 t)\bigg]dt=\frac{1}{T}\int^{\frac{T}{2}}_{-\frac{T}{2}}f(t)e^{-in\omega_0 t}dt \quad (n=1,2, 3, \cdots)$

同理:

$C_{-n}=\frac{1}{T}\int^{\frac{T}{2}}_{-\frac{T}{2}}f(t)e^{in\omega_0 t}dt \quad (n=1,2, 3, \cdots)$

上面的 $C_0, C_{-n},C_n$ 写为统一的形式为：

$C_{n}=\frac{1}{T}\int^{\frac{T}{2}}_{-\frac{T}{2}}f(t)e^{-in\omega_0 t}dt \quad (n=0,\pm1,\pm2, \pm3, \cdots)$

令

$\omega _n=n\omega_0 \quad (n=0,\pm1,\pm2, \pm3, \cdots)$

则综合上面各式，可得:

$\mathbf{f(t)\approx\frac{1}{T}\sum_{n=-\infty }^{\infty}\bigg[ \int^{\frac{T}{2}}_{-\frac{T}{2}}f(t)e^{-i\omega _nt}dt\bigg]e^{i\omega _n t}}$

拆分后得到傅里叶级数形式：

$\\\mathbf{ F(n )=\frac{1}{T}\int_{-\frac{T}{2}}^{\frac{T}{2}}f(t)e^{-in\omega_0 t}dt \quad (n \in Z)}$

$\\\mathbf{f(t)=\sum_{n=-\infty}^{\infty }F(n)\cdot e^{in\omega_0 t}\quad (n \in Z)}$

傅里叶级数推导出非周期信号的傅里叶变换：

$\\f(t)=\sum_{n=-\infty}^{\infty }F(n)\cdot e^{in\omega_0 t}=\sum_{n=-\infty}^{\infty }\bigg[\frac{1}{T}\int_{-\frac{T}{2}}^{\frac{T}{2}}f(t)e^{-in\omega_0 t}dt\bigg]\cdot e^{in\omega_0 t}\\=\frac{1}{T}\sum_{n=-\infty}^{\infty }\bigg[\int_{-\frac{T}{2}}^{\frac{T}{2}}f(t)e^{-in\omega_0 t}dt\bigg]\cdot e^{in\omega_0 t}$

当 $T\rightarrow \infty$ 时，周期信号变为非周期信号,由于 $\omega_0 =\frac{2\pi}{T}$ ，　傅里叶级数为：

$\\ \lim_{T \to \infty }f(t)=\lim_{T \to \infty }\frac{1}{T}\sum_{n=-\infty}^{\infty }\bigg[\int_{-\frac{T}{2}}^{\frac{T}{2}}f(t)e^{-in\omega_0 t}dt\bigg]\cdot e^{in\omega_0 t}\\=\lim_{T \to \infty }\frac{\omega_0}{2\pi}\sum_{n=-\infty}^{\infty }\bigg[\int_{-\frac{T}{2}}^{\frac{T}{2}}f(t)e^{-in\omega_0 t}dt\bigg]\cdot e^{in\omega_0 t}\\=\lim_{T \to \infty }\frac{1}{2\pi}\sum_{n=-\infty}^{\infty }\bigg[\int_{-\frac{T}{2}}^{\frac{T}{2}}f(t)e^{-in\omega_0 t}dt\bigg]\cdot e^{in\omega_0 t}\cdot \omega_0$

当 $T\rightarrow \infty$ 时候， $\Delta \omega =n\omega_0 - (n-1)\omega_0=\omega_0\rightarrow 0$

根据微积分的微元法，外面的累加可以看成求底边为 $\omega_0$ ,高为

$\bigg[\int_{-\frac{T}{2}}^{\frac{T}{2}}f(t)e^{-in\omega_0 t}dt\bigg]\cdot e^{in\omega_0 t}$

的图形的面积：

所以：

$F(\omega)=\int_{-\infty}^{\infty}f(t)e^{-i\omega t}dt \qquad (-\infty<\omega< \infty )$

$f(t)=\frac{1}{2\pi}\int_{-\infty}^{\infty}F(\omega)e^{i\omega t}d\omega \qquad (-\infty<\omega< \infty )$

一个例子从傅里叶级数到傅里叶变换：

此函数的解析式是：

$f(x)=\left\{\begin{matrix} 0 \qquad x>-\frac{T}{2}\ and \ x < -\tau \\ 1 \qquad \ \ x >-\tau \ and \ x < \tau\\ 0 \qquad \ \ \ \ x > \tau \ and \ x < \frac{T}{2} \end{matrix}\right.$

$\\ a_0=\frac{1}{T}\int_{-\frac{T}{2}}^{\frac{T}{2}}f(x)dx = \frac{1}{T}\int_{-\frac{T}{2}}^{\frac{T}{2}}f(x)dx =\frac{1}{T}\int_{-\tau}^{\tau}1\cdot dx = \frac{1}{T}\int_{-\tau}^{\tau}1\cdot dx\\=\frac{1}{T}\bigg|^{\tau}_{-\tau}=\frac{2\tau}{T}$

$\\ b_n =0$

$\\ a_n =\frac{2}{T}\int_{-\frac{T}{2}}^{-\frac{T}{2}}f(x)cos(n\cdot \frac{2\pi}{T}\cdot x)dx=\frac{2}{T}\int_{-\tau}^{\tau}1\cdot cos(n\cdot \frac{2\pi}{T}\cdot x)dx \\=\frac{T}{2n\pi}\cdot \frac{2}{T}\int_{-\tau}^{\tau}1\cdot cos(n\cdot \frac{2\pi}{T}\cdot x)d(n\cdot \frac{2\pi}{T}\cdot x)=\frac{1}{n\pi}sin(n\cdot \frac{2\pi}{T}\cdot x)\bigg|_{-\tau}^{\tau}=\frac{2sin(\frac{2n\pi \tau}{T})}{n\pi}$

函数图形为：

python代码：


# -*- coding: utf-8 -*-
"""
Created on Mon Feb  1 13:57:21 2021
@author: czl
"""
from pylab import *
 
x = mgrid[-20:20:0.01]
 
def fourier_wave():
    a0 = 3/16
    s=a0
    
    for n in range(1,1000,1):
        bn = 0
        an = 2*sin((2*n*pi*1.5/16))/(n*pi)
        s0 = an*cos(n*x*(2*pi/16))+bn*sin(n*x*(2*pi/16))
        s=s+s0
        
    plot(x,s,'orange',linewidth=0.6)
    title('fourier_transform')
    show()    
 
fourier_wave()

复指数形式的傅里叶变换系数如下，这里的负频率的意义是单位圆的旋转方向，并不是普通意义上“负”的概念。

$f(n\omega_0)=\frac{a_n-ib_n}{2}=\frac{a_n}{2}=\frac{sin(\frac{2n\pi \tau}{T})}{n\pi}$

$f(n\cdot \frac{2\pi}{T})=\frac{sin(\frac{2n\pi \tau}{T})}{n\pi}$

密度谱：

$T\cdot f(n\cdot \frac{2\pi}{T})=2\tau \cdot \frac{sin(\frac{2n\pi \tau}{T})}{\frac{2n\pi\tau}{T}}$

$\omega_0=\frac{2\pi}{T} \ \ \omega =n\omega_0$

当 $T->\infty$ 时:

$F(\omega)=T\cdot f(n\cdot \frac{2\pi}{T})=2\tau \cdot \frac{sin(\omega\tau)}{\omega\tau}$

$F(n)=T\cdot f(n\cdot \frac{2\pi}{T})=2\tau \cdot \frac{sin(n\omega_0\tau)}{n\omega_0\tau}=2\tau \cdot \frac{sin(n\frac{2\pi\tau}{T})}{n\frac{2\pi\tau}{T}}=2\tau \cdot \frac{sin(n\cdot \pi\cdot s)}{n\cdot \pi\cdot s}$