from pylab import *

x = mgrid[-10:10:0.02]  # 这里类似于MATLAB用冒号产生步长为0.02的序列，但是语法和MATLAB不同

n = arange(1,1000)

def fourier_transform():

    a0 = (1-exp(-pi))/pi+1

    s=a0/2

    for i in range(1,100,1):

        s0 = (  (1-(-1)**i*exp(-pi))/(pi*(1+i**2))*cos(i*x)+1/pi*( (-i*(1-(-1)**i*exp(-pi)))/(1+i**2) + (1-(-1)**i)/i ) * sin(i*x) )

        s=s+s0

    plot(x,s,'orange',linewidth=0.6)

    title('fourier_transform')

    show()    

fourier_transform( )

通过以上的例子可以很好直观感受傅立叶级数的作用，在电气工程上对分析流经变压器的谐波分量、整流装置输出谐波等分析具有重要作用。

还有一点就是python的强大功能，利用for循环语句可以方便地实现级数展开的求和运算，自己感觉编程起来比MATLAB/Octave更快。