python实现傅里叶级数展开的实现
作者：落叶_小唱 发布时间：2022-06-14 14:35:17
标签：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）
来源：https://blog.csdn.net/ouening/article/details/70339341
0
投稿猜你喜欢
- omitempty在go中的使用直接上代码：package mainimport ( "encoding/json&q
- 正在看的ORACLE教程是:在ORACLE移动数据库文件。
- 一、使用python内置commands模块执行shellcommands对Python的os.popen()进行了封装，使用SHELL命令
首先在文件夹中添加_init_.py文件，即使是空文件也可以，多层文件夹，每层文件夹中都要添加。比如我要import，a文件夹中，b文件夹下- 之前有聊过 golang 的协程，我发觉似乎还很理论，特别是在并发安全上，所以特结合网上的一些例子，来试验下go routine中 的 ch
- 前言大家好，我是辣条相信大家都能感觉到最近天气的多变，好几次出门半路天气转变。辣条也深受其扰，直接给我整感冒，就差被隔离起来了，既然天气我没
- 1. 为什么要有转义？ASCII 表中一共有 128 个字符。这里面有我们非常熟悉的字母、数字、标点符号，这些都可以从我们的键盘中输出。除此
- 导言在前面的教程里我们学习了DataList提供了一些风格样式的属性.而且我们还学习了如何定义HeadStyle, ItemStyle, A
- TCP是用来计算机之间进行通信的，通过编写客户端和服务端聊天的代码，对于服务器与客户端的工作步骤有了深刻的了解，在这里根据了不起的Node.
- hasOwnPropertyhasOwnProperty是javascript中用于检测对象是否包含某个属性的方法，返回一个布尔值。var
- python使用ctypes模块调用windows api GetVersionEx获取当前系统版本,没有使用python32#!c:/py
- pytest是python语言中一款强大的单元测试框架，用来管理和组织测试用例，可应用在单元测试、自动化测试工作中。unittest也是py
- 一个站点成功与否，不是在于所用的程序语言是PHP还是ASP，也不在于用BLOG建站或是CMS建站，内容和专业很重要。博客很好起步，CMS体系
- 插入数据insert into 表名(列名1,列名2,列名3) values(值1,值2,值3);insert into user(user
- subprocess.run()运行并等待args参数指定的指令完成，返回CompletedProcess实例。参数：(*popenargs
- 摘要1,EXEC的使用2，sp_executesql的使用MSSQL为我们提供了两种动态执行SQL语句的命令，分别是EXEC和sp_exec
- USE [DAF_DB] GO /****** Object: StoredProcedure [dbo].[PROG_WORKTASK_L
- 像这种不能创建一个.frm 文件的报错好像暗示着操作系统的文件的权限错误或者其它原因，但实际上，这些都不是的，事实上，这个mysql报错已经
- 1.在列属性中加入事件 { &
- 一、正则表达式的特殊字符介绍正则表达式^ 匹配行首 &nb
