python实现傅立叶级数展开
来源:互联网 发布:全国人口数据库 编辑:程序博客网 时间:2024/06/08 11:22
相关文章:Python实现快速傅里叶变换(FFT)
傅立叶级数的介绍我就不说了,自己也是应用为主,之前一直觉得很难懂,但最近通过自己编程实现了一些函数的傅立叶级数展开之后对傅立叶 级数展开的概念比较清楚了
(1)函数如下
函数图象如下:
代码:
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()
结果如下:
(2)
转换代码:
from pylab import *def fourier1(): s=pi/2 for i in range(1,100,1): s0 = 2/pi*(1-(-1)**i)/i**2 * cos(i*x) s=s+s0 plot(x,s,'orange',linewidth=0.6) title('fourier1') show()fourier1()
结果如下:
关于傅立叶级数展开的还有一些例子在我的另一篇博客http://blog.csdn.net/ouening/article/details/53199215 也有讨论,通过以上的例子可以很好直观感受傅立叶级数的作用,在电气工程上对分析流经变压器的谐波分量、整流装置输出谐波等分析具有重要作用。
还有一点就是python的强大功能,利用for循环语句可以方便地实现级数展开的求和运算,自己感觉编程起来比MATLAB/Octave跟快
0 0
- python实现傅立叶级数展开
- 傅立叶级数
- 连续信号(模拟信号)在有限区间上的傅立叶级数展开,离散频谱
- 完全搞懂傅里叶变换和小波(5)——傅立叶级数展开之函数项级数的概念
- 完全搞懂傅里叶变换和小波(6)——傅立叶级数展开之函数项级数的性质
- 用html+JavaScript实现傅立叶级数的动态可视化
- 傅立叶级数神贴
- 傅立叶级数动图
- 傅立叶级数和傅立叶变换是什么关系?
- 傅立叶级数与傅立叶变换总结
- 傅立叶级数/变换及其应用
- 通俗易懂的傅立叶级数理解
- 泰勒级数展开
- 傅立叶变换及应用1-傅立叶级数
- 1/n的级数展开
- ln(x) 泰勒级数展开
- MAthematica 傅里叶级数展开画图
- 傅立叶外推算法的python实现和缺点
- UVALive
- 1210
- 剑指offer(C++)——从尾到头打印链表
- typedef 枚举 联合体 条件运算符和表达式
- libevent入门
- python实现傅立叶级数展开
- 封装
- eclipse自动补全的设置
- HDU2844 Coins
- 数论
- Ubuntu下Python 3.3的安装
- 子集选择
- RSA公钥私钥的简单实现
- USACO 2.3 货币系统(背包/生成函数)