python实现周期方波信号频谱图
更新时间：2018年07月21日 14:33:26 作者：落叶_小唱
这篇文章主要介绍了python 周期方波信号频谱图，小编觉得挺不错的，现在分享给大家，也给大家做个参考。一起跟随小编过来看看吧
在学习傅里叶变换的时候遇到了求周期方波信号频谱图的例子，在书上和网上查阅了一些资料，发现大都是讨论的都是下图左边的周期信号的频谱，课程老师的PPT中也只列出了另一种周期信号频谱图的结论，没有在进行傅里叶变换，自己便根据定义推导了一遍，贴在这里作记录和分享之用。
关于傅立叶级数展开的另一讨论在我的另一篇文章https://www.jb51.net/article/144194.htm
2016年11月21号更新
在第二个周期方波信号的傅里叶变换里，注意是
转换为sin函数下为
之前写错了，今天更正。
对于这两种方波信号，我们也可以编程验证一下，就是用正弦函数去逼近方波信号，Python的实现代码如下：
# 分析傅里叶级数分解之后cos和sin的和项的图像输出
from numpy import mgrid,sin,cos,array,pi
from matplotlib.pyplot import plot,show,title,legend,xlabel,ylabel
x = mgrid[0:10:0.02] # 这里类似于MATLAB用冒号产生步长为0.02的序列，但是语法和MATLAB不同
# 下面的这段循环实现y=sin(x)+sin(3x)+...+sin(19x)
def cos_square():
y1 = 0;
for i in range(0,20,1):
b = (-1)**(i)*cos((2*i+1)*x)/(2*i+1)
y1=b+y1 # 这种求和的方法是从C语言移植过来的
plot(x,y1,'orange',linewidth=0.6)
title('cos_square')
xlabel('Time')
ylabel('Amplitude')
show()
def sin_square():
y2 = 0
for i in range(0,20,1):
b = sin((2*i+1)*x)/(2*i+1)
y2=b+y2 # 这种求和的方法是从C语言移植过来的
plot(x,y2,'g',linewidth=0.6)
title('sin_square')
xlabel('Time')
ylabel('Amplitude')
show()
cos_square()
sin_square()
输出结果如下：
由结果可以看两种不同的方波级数展开是正确的
以上就是本文的全部内容，希望对大家的学习有所帮助，也希望大家多多支持脚本之家。
相关文章
webbrowser模块是Python自带的标准库,无需安装,可以直接在Python中使用该模块来打开网页、PDF文件等,本文给大家详细介绍了Python webbrowser模块用法,需要的朋友可以参考下2023-08-08
深入探索Python中的反射机制和动态属性的细节,对我们来说是编写具有适应性和高可扩展性程序的关键,本篇文章旨在通过详尽的概念介绍和精心设计的代码示例,加强您对这些核心概念的把握,并助您在实践中运用自如2024-06-06
今天小编就为大家分享一篇python之文件读取一行一行的方法，具有很好的参考价值，希望对大家有所帮助。一起跟随小编过来看看吧2018-07-07
这篇文章主要介绍了selenium python浏览器多窗口处理代码示例，具有一定借鉴价值,需要的朋友可以参考下2018-01-01
这篇文章主要介绍了使用matplotlib绘制热图(heatmap)全过程，具有很好的参考价值，希望对大家有所帮助。如有错误或未考虑完全的地方，望不吝赐教2022-12-12
这篇文章主要介绍了WxPython中控件隐藏与显示的小技巧，具有很好的参考价值，希望对大家有所帮助。如有错误或未考虑完全的地方，望不吝赐教2022-11-11
这篇文章主要介绍了Python2.x版本中基本的中文编码问题解决,Python3版本中默认Unicode编码便不会遇到此类问题,需要的朋友可以参考下2015-10-10
这篇文章主要给大家介绍了如何通过一篇文章搞懂Python反斜杠的相关问题，文中通过示例代码介绍的非常详细，对大家的学习或者工作具有一定的参考学习价值，需要的朋友们下面随着小编来一起学习学习吧2021-03-03
这篇文章主要为大家介绍了Python支持异步的列表解析式示例详解，有需要的朋友可以借鉴参考下，希望能够有所帮助，祝大家多多进步，早日升职加薪2022-07-07
这篇文章主要介绍了python实时获取kafka消费队列信息,本文通过实例代码给大家介绍的非常详细，对大家的学习或工作具有一定的参考借鉴价值，需要的朋友可以参考下2023-07-07
最新评论