The separate pictures was made with this program (fourier.py):
#!/usr/bin/env python
# coding: utf-8
from __future__ import division
import pylab
from pylab import show, plot, linspace
from pylab import zeros
from pylab import sin, pi
def fourier(n=15):
x = pylab.linspace(-pi, pi, 1024)
y = pylab.zeros(1024)
pylab.grid(False)
pylab.hold(False)
for i in range(1,n+1,2):
y += sin(i*x)/i
pylab.plot(x,4/pi*y, "b")
pylab.title(u"Fourier series approximation of square wave, %02d addends" %n)
pylab.savefig("fourier%02d.png" % n)
#pylab.show()
if __name__ == "__main__":
for i in range(1,33,2):
fourier(i)
Then thise pictures was unified with ImageMagick into one gif as follows:
convert -delay 30 -loop 0 *.png fourier.gif
Llicència
Jo, el titular dels drets d'autor d'aquest treball, el public sota la següent llicència:
compartir – copiar, distribuir i comunicar públicament l'obra
adaptar – fer-ne obres derivades
Amb les condicions següents:
reconeixement – Heu de donar la informació adequada sobre l'autor, proporcionar un enllaç a la llicència i indicar si s'han realitzat canvis. Podeu fer-ho amb qualsevol mitjà raonable, però de cap manera no suggereixi que l'autor us dóna suport o aprova l'ús que en feu.
compartir igual – Si modifiqueu, transformeu, o generareu amb el material, haureu de distribuir les vostres contribucions sota una llicència similar o una de compatible com l'original