分享我写的小程序--机械工程师理解的傅里叶级数

从零开始

见原帖，机械工程师理解的傅里叶级数，方波的构造。
http://jixietop.cn/forum.php?mod=viewthread&tid=25156&extra=

以前在老社区见id“动静之机”的网友发表过，感觉自己现在可以写出来。
于是尝试用mathematica语言试了一下。
------------
以下是动图和代码，(lisp语言的括号是可怕）。

从零开始

或者拷贝一下文件。
---------------
Animate[Graphics[{PointSize[Large], Circle[{0, 0}, 4/\[Pi]],
   Circle[{(4 Cos[t])/\[Pi], (4 Sin[t])/\[Pi]}, 4/(3 \[Pi])],
   Circle[{(4 Cos[t])/\[Pi] + (4 Cos[3 t])/(3 \[Pi]), (
      4 Sin[t])/\[Pi] + (4 Sin[3 t])/(3 \[Pi])}, 4/(5 \[Pi])],
   Circle[{(4 Cos[t])/\[Pi] + (4 Cos[3 t])/(3 \[Pi]) + (4 Cos[5 t])/(
      5 \[Pi]), (4 Sin[t])/\[Pi] + (4 Sin[3 t])/(3 \[Pi]) + (
      4 Sin[5 t])/(5 \[Pi])}, 4/(7 \[Pi])],
   Point[{{(4 Cos[t])/\[Pi], (
      4 Sin[t])/\[Pi]}, {(4 Cos[t])/\[Pi] + (4 Cos[3 t])/(3 \[Pi]), (
       4 Sin[t])/\[Pi] + (4 Sin[3 t])/(3 \[Pi])}, {(
       4 Cos[t])/\[Pi] + (4 Cos[3 t])/(3 \[Pi]) + (4 Cos[5 t])/(
       5 \[Pi]), (4 Sin[t])/\[Pi] + (4 Sin[3 t])/(3 \[Pi]) + (
       4 Sin[5 t])/(5 \[Pi])}, {(4 Cos[t])/\[Pi] + (4 Cos[3 t])/(
       3 \[Pi]) + (4 Cos[5 t])/(5 \[Pi]) + (4 Cos[7 t])/(7 \[Pi]), (
       4 Sin[t])/\[Pi] + (4 Sin[3 t])/(3 \[Pi]) + (4 Sin[5 t])/(
       5 \[Pi]) + (4 Sin[7 t])/(7 \[Pi])}}], Green,
   Point[Table[{3 + k, (4 Sin[-k + t])/\[Pi] + (4 Sin[3 (-k + t)])/(
       3 \[Pi]) + (4 Sin[5 (-k + t)])/(5 \[Pi]) + (
       4 Sin[7 ( -k + t)])/(7 \[Pi])}, {k, 0, 6, .05}]], Red,
   Arrowheads[Medium],
   Arrow[{{(4 Cos[t])/\[Pi] + (4 Cos[3 t])/(3 \[Pi]) + (4 Cos[5 t])/(
       5 \[Pi]) + (4 Cos[7 t])/(7 \[Pi]), (4 Sin[t])/\[Pi] + (
       4 Sin[3 t])/(3 \[Pi]) + (4 Sin[5 t])/(5 \[Pi]) + (4 Sin[7 t])/(
       7 \[Pi])}, {3 + 0, (4 Sin[t])/\[Pi] + (4 Sin[3 t])/(
       3 \[Pi]) + (4 Sin[5 t])/(5 \[Pi]) + (4 Sin[7 t])/(7 \[Pi])}}],
   Blue, Line[{{0, 0}, {(4 Cos[t])/\[Pi], (
      4 Sin[t])/\[Pi]}, {(4 Cos[t])/\[Pi] + (4 Cos[3 t])/(3 \[Pi]), (
       4 Sin[t])/\[Pi] + (4 Sin[3 t])/(3 \[Pi])}, {(
       4 Cos[t])/\[Pi] + (4 Cos[3 t])/(3 \[Pi]) + (4 Cos[5 t])/(
       5 \[Pi]), (4 Sin[t])/\[Pi] + (4 Sin[3 t])/(3 \[Pi]) + (
       4 Sin[5 t])/(5 \[Pi])}, {(4 Cos[t])/\[Pi] + (4 Cos[3 t])/(
       3 \[Pi]) + (4 Cos[5 t])/(5 \[Pi]) + (4 Cos[7 t])/(7 \[Pi]), (
       4 Sin[t])/\[Pi] + (4 Sin[3 t])/(3 \[Pi]) + (4 Sin[5 t])/(
       5 \[Pi]) + (4 Sin[7 t])/(7 \[Pi])}}]},
  PlotRange -> {{-2.5, 10}, {-2.5`, 2.5`}}, ImageSize -> Full], {t, 0,
   6 \[Pi]}, AnimationRate -> .1, AnimationRunning -> False,
ContentSize -> {900, 320}]

2266998

不错，应该鼓励，

nanfeng36

在YouTube上也看到一个老外使用同样的操作，他是用Mathlab做的。