MAthematica 傅里叶级数展开画图

f[x_] := Which[-Pi <= x < Pi, x]
f[x_] := f[x - 2*Pi] /; x >= Pi
f[x_] := f[x - 2*Pi] /; x < -Pi
A = Plot[f[x], {x, -Pi, Pi}, PlotStyle -> Red];
k = 20;
a[n_] := (1/Pi)*Integrate[f[x]*Cos[n*x], {x, -Pi, Pi}];
b[n_] := (1/Pi)*Integrate[f[x]*Sin[n*x], {x, -Pi, Pi}];
a[0]/2 + Sum[a[n]*Cos[n*x] + b[n]*Sin[n*x], {n, 1, k}]
A2 = Plot[
   Evaluate[
    a[0]/2 + Sum[a[n]*Cos[n*x] + b[n]*Sin[n*x], {n, 1, 2}]], {x, -Pi,
    Pi}, PlotStyle -> Hue[0.5]];
A3 = Plot[
   Evaluate[
    a[0]/2 + Sum[a[n]*Cos[n x] + b[n]*Sin[n x], {n, 1, 3}]], {x, -Pi,
    Pi}, PlotStyle -> Hue[0.3]];
A4 = Plot[
   Evaluate[
    a[0]/2 + Sum[a[n]*Cos[n x] + b[n]*Sin[n x], {n, 1, 4}]], {x, -Pi,
    Pi}, PlotStyle -> Hue[0.4]];
A6 = Plot[
   Evaluate[
    a[0]/2 + Sum[a[n]*Cos[n x] + b[n]*Sin[n x], {n, 1, 6}]], {x, -Pi,
    Pi}, PlotStyle -> Hue[0.6]];
A8 = Plot[
   Evaluate[
    a[0]/2 + Sum[a[n]*Cos[n x] + b[n]*Sin[n x], {n, 1, 8}]], {x, -Pi,
    Pi}, PlotStyle -> Hue[0.8]];
A15 = Plot[
   Evaluate[
    a[0]/2 + Sum[a[n]*Cos[n*x] + b[n]*Sin[n*x], {n, 1, 15}]], {x, -Pi,
     Pi}, PlotStyle -> Hue[0.7]];
Show[A, A2, A3, A4, A6, A8, A15, PlotRange -> All]
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值