python实现傅里叶级数展开的实现
更新时间：2018年07月21日 14:27:48 作者：落叶_小唱
这篇文章主要介绍了python实现傅里叶级数展开的实现，小编觉得挺不错的，现在分享给大家，也给大家做个参考。一起跟随小编过来看看吧
傅立叶级数的介绍我就不说了，自己也是应用为主，之前一直觉得很难懂，但最近通过自己编程实现了一些函数的傅立叶级数展开之后对傅立叶 级数展开的概念比较清楚了
(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()
结果如下：
关于傅里叶级数展开的还有一些例子在我的另一篇博客https://www.jb51.net/article/144196.htm也有讨论，通过以上的例子可以很好直观感受傅立叶级数的作用，在电气工程上对分析流经变压器的谐波分量、整流装置输出谐波等分析具有重要作用。
还有一点就是python的强大功能，利用for循环语句可以方便地实现级数展开的求和运算，自己感觉编程起来比MATLAB/Octave跟快
相关文章：Python实现快速傅里叶变换（FFT）
以上就是本文的全部内容，希望对大家的学习有所帮助，也希望大家多多支持脚本之家。
相关文章
这篇文章主要介绍了python 获取毫秒级时间问题的解决，具有很好的参考价值，希望对大家有所帮助。一起跟随小编过来看看吧2021-04-04
这篇文章主要给大家介绍了关于Anaconda彻底删除虚拟环境的正确方法,要在Anaconda中删除一个虚拟环境,可以按照本文以下步骤进行操作,需要的朋友可以参考下2023-10-10
Python3.7在anaconda里面使用IDLE编译器的步骤详解
这篇文章主要介绍了Python3.7在anaconda里面使用IDLE编译器的步骤，本文通过图文并茂的形式给大家介绍的非常详细，对大家的学习或工作具有一定的参考借鉴价值，需要的朋友参考下吧2020-04-04
这篇文章主要介绍了Python序列之list和tuple常用方法以及注意事项,sequence(序列)是一组有顺序的对象的集合,序列可以包含一个或多个元素,也可以没有任何元素,序列有两种：list (表) 和 tuple（元组）,需要的朋友可以参考下2015-01-01
在本篇文章里小编给大家整理的是一篇关于python常量折叠基础知识点讲解，对此有兴趣的朋友可以跟着学习下。2021-02-02
这篇文章主要介绍了Keras在训练期间可视化训练误差和测试误差实例，具有很好的参考价值，希望对大家有所帮助。一起跟随小编过来看看吧2020-06-06
OpenCV基础HSV颜色空间*args与**kwargs滑动条传参问题
这篇文章主要介绍了OpenCV基础中HSV颜色空间*args与**kwargs滑动条传参问题的详细解析，有需要的朋友可以借鉴参考下，希望能够有所帮助2021-10-10
这篇文章主要为大家详细介绍了python实现textrank关键词提取，具有一定的参考价值，感兴趣的小伙伴们可以参考一下2018-06-06
这篇文章主要介绍了Python面向对象编程基础,结合实例形式分析了Python面向对象编程类的定义、继承、特殊方法及模块相关原理与操作技巧,需要的朋友可以参考下2020-01-01
今天小编就为大家分享一篇tensorflow 变长序列存储实例，具有很好的参考价值，希望对大家有所帮助。一起跟随小编过来看看吧2020-01-01
最新评论