傅里叶积分

傅里叶积分公式的推导过程

傅里叶级数公式f(x) = \frac{a_0}{2} + \sum_{n=1}^\infty[a_n\cos(n\omega x)+b_n\sin(n\omega x)],

其中a_n=\frac{2}{T}\int_{t_0}^{t_0+T}f(x)\cos(n\omega x)dx,b_n=\frac{2}{T}\int_{t_0}^{t_0+T}f(x)\sin(n\omega x)dx

已知欧拉公式e^i^x=\cos(x) + i\sin(x),所以\cos(x) = \frac{e^i^x+e^{-ix}}{2},\sin(x) = \frac{-i(e^{ix}-e^{-ix})}{2}.
带入f(x) = \frac{a_0}{2} + \sum_{n=1}^\infty[a_n\cos(n\omega x)+b_n\sin(n\omega x)]f(x) = \frac{a_0}{2}+\sum_{n=1}^{\infty}[\frac{a_n-ib_n}{2}e^{in\omega x}+\frac{a_n+ib_n}{2}e^{-in\omega x}]


C_0 = \frac{a_0}{2}=\frac{1}{T}\int_{-T/2}^{T/2}f(x)dx,C_n = \frac{a_n-ib_n}{2}[\int_{-T/2}^{T/2}f(x)\cos(n\omega x)dx-i\int_{-T/2}^{T/2}f(x)\sin(n\omega x)dx]=\frac{1}{T}\int_{-T/2}^{T/2}f(x)e^{-in\omega x}dx,n=1,2,....

C_{-n} = \frac{a_n+ib_n}{2}=\frac{1}{T}\int_{-T/2}^{T/2}f(x)e^{in\omega x}dx,n=1,2,....
C_n = \frac{1}{T}\int_{-T/2}^{T/2}f(x)e^{-in\omega x}dx,n=\pm1,\pm2,.....

\omega_n = n\omega,(n=0,\pm1,\pm2,...),得到傅里叶级数的复指数形式:f(x) = C_0+\sum_{n=1}^{\infty}(C_ne^{i\omega_nx}+C_{-n}e^{-i\omega_nx}) = \sum_{-\infty}^{\infty}C_ne^i\omega_nx,n=\pm1,\pm2,...
或将得C_n = \frac{1}{T}\int_{-T/2}^{T/2}f(x)e^{-in\omega x}dx,n=\pm1,\pm2,....带入得f(x) = \frac{1}{T}\sum_{-\infty}^{\infty}[\int_{-T/2}^{T/2}f(\tau)e^{-i\omega_nx}d\tau]e^{i\omega_nx},n=\pm1,\pm2,....(T=2π)

周期函数可以展开成傅里叶级数,那么非周期函数呢?能否用一个周期函数逼近一个非周期函数呢?
一般而言,任何一个非周期函数f(x)都可以看成是由某个周期函数f_T(x)当周期T → +∞时转化而来的.
Proof: 作周期为T的函数f_T(x),使其在[-\frac{T}{2},\frac{T}{2}]之内等于f(x),在[-\frac{T}{2},\frac{T}{2}]之外按周期为T的函数f_T(x)延拓出去。

即为f_T(x) = f(x), x\in[-\frac{T}{2},\frac{T}{2}]f_T(x) = f_T(x+T), x\notin[-\frac{T}{2},\frac{T}{2}].

则T越大,f_T(x)f(x)相等的范围越大,则当T\rightarrow \infty时,f_T(x)可以转换为f(x),即为\lim_{T\rightarrow\infty} f_T(x) = f(x)\Rightarrow f(x) = \lim_{T\rightarrow\infty}\frac{1}{T}\sum_{n=-\infty}^{+\infty}[\int_{-T/2}^{T/2}f_T(\tau)e^{-i\omega_n\tau}d\tau]e^{i\omega_nx}.

当n取一切N时,\omega_n所对应的点能均匀的分布在整个数轴上。

若取相邻的距离\triangle \omega = \omega_n-\omega_{n-1}=\frac{2\pi}{T}\Rightarrow T = \frac{2\pi}{\triangle\omega},当T\rightarrow \infty,有\triangle \omega \to 0

T = \frac{2\pi}{\triangle\omega}带入f(x) = \lim_{T\rightarrow\infty}\frac{1}{T}\sum_{n=-\infty}^{+\infty}[\int_{-T/2}^{T/2}f_T(\tau)e^{-i\omega_n\tau}d\tau]e^{i\omega_nx}得到f(x) = \lim_{T\rightarrow\infty}\frac{1}{2\pi}\sum_{n=-\infty}^{+\infty}[\int_{-T/2}^{T/2}f_T(\tau)e^{-i\omega_n\tau}d\tau]e^{i\omega_nx}\triangle \omega

当x固定时,\frac{1}{2\pi}[\int_{-T/2}^{T/2}f_T(\tau)e^{-i\omega_n\tau}d\tau]e^{i\omega_nx}\triangle \omega就为\omega的函数,即为\phi_T(\omega),即为\phi_T(\omega)=\frac{1}{2\pi}[\int_{-T/2}^{T/2}f_T(\tau)e^{-i\omega_n\tau}d\tau]e^{i\omega_nx}\triangle \omega.得到f(x) = \lim_{\triangle \omega\to0}\sum_{n=-\infty}^{+\infty}\phi_T(\omega_n)\triangle\omega

\triangle\omega\to0,T\to\infty,有\phi_T(\omega_n)\to\phi(\omega).此处\phi(\omega)=\frac{1}{2\pi}[\int_{-\infty}^{+\infty}f(\tau)e^{i\omega\tau}d\tau]

所以f(x)可以看成是\phi(\omega)(-\infty,+\infty)上的积分f(x)=\int_{-\infty}^{+\infty}\phi(\omega)d\omega,即为f(x) = \frac{1}{2\pi}\int_{-\infty}^{+\infty}[\int_{-\infty}^{+\infty}f(\tau)e^{-i\omega\tau}d\tau]e^{i\omega x}d\omega(傅里叶积分公式)。

傅里叶积分公式用欧拉公式变为三角形式的傅里叶积分
e^{ix} = \cos(x)+i\sin(x)
f(x) = \frac{1}{\pi}\int_{0}^{+\infty}[\int_{-\infty}^{+\infty}f(\tau)cos(\omega)(x-\tau)d\tau]d\omega

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
傅里叶积分是一种重要的数学工具,用于将一个时间域上的信号转换到频率域上。在MATLAB中,我们可以利用傅里叶变换函数进行傅里叶积分的计算。 MATLAB中的傅里叶变换函数主要有两种,分别是fft(x)和ifft(X),其中x是时间域上的信号,X是频率域上的信号。fft函数将信号从时间域变换到频率域,而ifft函数则将信号从频率域变换回时间域。 使用fft函数进行傅里叶积分,首先需要定义一个时间域上的信号x,然后通过调用fft(x)即可得到频率域上的信号X。X是一个包含正频率和负频率的复数数组,可以通过abs(X)求取其模值,angle(X)求取其相位值。 而使用ifft函数进行傅里叶逆变换,只需将频率域上的信号X作为参数传入即可,通过调用ifft(X)即可得到变换回的时间域上的信号。同样,变换回的时间域信号也是一个复数数组,可以通过real(ifft(X))仅提取其实部,imag(ifft(X))仅提取其虚部。 需要注意的是,通过fft和ifft进行的变换都是离散傅里叶变换(DFT和IDFT),所以信号的采样点数需要是有限的。在MATLAB中,fft函数默认使用了快速傅里叶变换(FFT)算法,能够快速、高效地进行计算。 总而言之,MATLAB中的傅里叶积分可以通过fft函数将时间域信号转换到频率域上,通过ifft函数将频率域信号转换回时间域上。这样的变换过程可以帮助我们更好地理解信号的频域特性,并进行相应的频域分析和处理。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值