MATLAB之傅里叶展开（五）

傅里叶（Fourier）展开

傅里叶级数的定义如下：

设函数$f(x)$在区间$[0,2\pi ]$上绝对可积，令
$$\{\begin{array}{ll}{a}_n=\frac{1}{\pi }{\int }_0^{2\pi }f(x)\cos nx\mathrm{d}x& (n=0,1,2,\cdots )\\ b_n=\frac{1}{\pi }{\int }_0^{2\pi }f(x)\sin nx\mathrm{d}x& (n=1,2,\cdots )\end{array}$$
以$a_n,b_n$为系数作三角级数
$$\frac{a_0}{2}+\sum _{n=1}^{\infty }(a_{n}cosnx+b_{n}sinnx)$$
称为$f(x)$的傅里叶级数，$a_n,b_n$称为傅里叶系数。

具体的傅里叶级数推导过程可以参考下面这篇文章：

https://zhuanlan.zhihu.com/p/41455378

根据定义可以编写计算傅里叶系数的函数：

function [a0,an,bn] = Fourierzpi(f)
%本函数旨在计算[0,2pi]上的傅里叶系数
syms x n
a0 = int(f,0,2*pi)/pi;
an = int(f*cos(n*x),0,2*pi)/pi;
bn = int(f*sin(n*x),0,2*pi)/pi;
end

利用该函数计算$\frac{sinx+e^x}{x^2}$在[0,2pi]上的傅里叶系数

syms x
f = (sin(x)+exp(x))/x^2;
[a0,an,bn] = Fourierzpi(f)

a0 =

Inf


an =

int((cos(n*x)*(exp(x) + sin(x)))/x^2, x, 0, 2*pi)/pi


bn =

int((sin(n*x)*(exp(x) + sin(x)))/x^2, x, 0, 2*pi)/pi

如果想要计算不同区间的傅里叶变化，只需要改变函数中范围，这样就可以计算不同区间，不同函数的傅里叶系数。