2005年07月14日

数値を解析って見よう 第2段

親切に公開されてるMathematicaプログラムは後々の応用性がなく、ソースの見た目が汚いので嫌ですね(笑)

問題1は以下をMathematicaで実行すれば面倒な事無しに終ります。

　　fa[x_]=x+Pi;  fa1=-Pi; fa2=-Pi/3;
　　fb[x_]=2Pi/3; fb1=fa2; fb2=Pi/3;
　　fc[x_]=-x+Pi; fc1=fb2; fc2=Pi;
　　n=9;

　　a1[0]=Integrate[fa[x]/Pi,{x,fa1,fa2}];
　　Print[0," α1=",a1[0]]
　　Do[{a1[i]=Integrate[fa[x]Cos[i x]/Pi,{x,fa1,fa2}], 
　 　   a2[i]=Integrate[fa[x]Sin[i x]/Pi,{x,fa1,fa2}],
　 　   Print[i," α1=",Chop[a1[i]]," β1=",Chop[a2[i]]]},
　　   {i,1,n}]

　　b1[0]=Integrate[fb[x]/Pi,{x,fb1,fb2}];
　　Print[0," α2=",b1[0]]
　　Do[{b1[i]=Integrate[fb[x]Cos[i x]/Pi,{x,fb1,fb2}],
　　    b2[i]=Integrate[fb[x]Sin[i x]/Pi,{x,fb1,fb2}],
　　    Print[i," α2=",Chop[b1[i]]," β2=",Chop[b2[i]]]},
　　   {i,1,n}]

　　c1[0]=Integrate[fc[x]/Pi,{x,fc1,fc2}];
　　Print[0," α3=",c1[0]]
　　Do[{c1[i]=Integrate[fc[x]Cos[i x]/Pi,{x,fc1,fc2}],
　　    c2[i]=Integrate[fc[x]Sin[i x]/Pi,{x,fc1,fc2}],
　　    Print[i," α3=",Chop[c1[i]]," β3=",Chop[c2[i]]]},
　　   {i,1,n}]

　　go1=Plot[fa[x],{x,fa1,fa2}]
　　go2=Plot[fb[x],{x,fb1,fb2}]
　　go3=Plot[fc[x],{x,fc1,fc2}]

　　Do[{ga[x_]=a1[0]/2+Sum[a1[k]Cos[k x]+a2[k]Sin[k x],{k,1,i}];
　　     g1[i]=Plot[{ga[x]},{x,fa1,fa2}]},{i,1,n}]

　　Do[{gb[x_]=b1[0]/2+Sum[b1[k]Cos[k x]+b2[k]Sin[k x],{k,1,i}];
　　     g2[i]=Plot[{gb[x]},{x,fb1,fb2}]},{i,1,n}]

　　Do[{gc[x_]=c1[0]/2+Sum[c1[k]Cos[k x]+c2[k]Sin[k x],{k,1,i}];
　　     g3[i]=Plot[{gc[x]},{x,fc1,fc2}]},{i,1,n}]

　　　Do[{Show[g1[i],g2[i],g3[i],go1,go2,go3]},{i,1,n}]

←n=30のアニメーション

※上記ソースコードはHTML表示用のため直接貼り付けて実行することはできません。

投稿者 zeratinman : 2005年07月14日 00:00