傅里叶级数在三角插值多项式中的应用（含例题）

傅里叶级数在三角插值多项式中的应用

理论部分
用正弦和余弦函数表示任意函数始于18世纪50年代。下面介绍最佳平方三角逼近多项式和傅立叶级数。


例题部分
使用傅里叶变换确定函数f(x)=x^2*cos(x)在[-pi,pi]上的16次三角插值多项式。
（1）编写子程序fseries.m，求傅里叶系数及级数。程序如下：
% 功能：求解函数的傅里叶系数及傅里叶级数。 %an，bn为傅里叶系数，f为傅里叶多项式；
%fx为函数，x为自变量，n为三角多项式的次数，a,b为函数fx的定义域

function [an,bn,f]=fseries(fx,x,n,a,b)
if nargin==3
a=-pi;
b=pi;
end
l=(b-a)/2;
if a+b
fx=subs(fx,x,x+l+a);   %设置新的积分区间为[-(b-a)/2,(b-a)/2]
end
an=int(fx,x,-l,l)/l;
bn=[];
f=an/2;
for ii=1:n
ann=int(fx*cos(ii*pi*x/l),x,-l,l)/l;
bnn=int(fx*sin(ii*pi*x/l),x,-l,l)/l;
an=[an,ann];
bn=[bn,bnn];
f=f+ann*cos(ii*pi*x/l)+bnn*sin(ii*pi*x/l);
end
if a+b
f=subs(f,x,x-l-a);      %设置新的积分区间为[a,b]
end

    （2）编写命令文件。

syms x;
fx=x^2*cos(x);
[an,bn,f]=fseries(fx,x,16,-pi,pi)

    运行后，可得：

an =
[ -4, (1/3pi^3+1/2pi)/pi, -20/9, 5/8, -68/225, 13/72, -148/1225, 25/288, -260/3969, 41/800, -404/9801, 61/1800, -580/20449, 85/3528, -788/38025, 113/6272, -1028/65025]
bn=
bn =

[ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]
f =
-2+(1/3pi^3+1/2pi)/picos(x)-20/9cos(2x)+5/8cos(3x)-68/225cos(4x)+13/72cos(5x)-148/1225cos(6x)+25/288cos(7x)-260/3969cos(8x)+41/800cos(9x)-404/9801cos(10x)+61/1800cos(11x)-580/20449cos(12x)+85/3528cos(13x)-788/38025cos(14x)+113/6272cos(15x)-1028/65025cos(16*x)
可以发现，如果函数为偶函数，傅里叶技术bk值为0。
（3）结果验证

t=-pi:pi/100:pi;
y=t.^2.*cos(t);
plot(t,y,'r');           %绘制特定点对应的函数值来验证，为红色线条
hold on
ezplot(f,[-pi,pi])
%绘制16次三角插值多项式，用默认蓝色显示。

运行结果如下：

由图形可以看出，
两条曲线重叠效果好，16次三角插值多项式逼近原函数，误差小。可以在[-pi,pi]区间内使用三角多项式来研究原函数，进一步证明了傅里叶变换的合理性。