python实现周期方波信号频谱图
作者：落叶_小唱 发布时间：2021-04-11 00:11:46
标签：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()
输出结果如下：
由结果可以看两种不同的方波级数展开是正确的
来源：https://blog.csdn.net/ouening/article/details/53199215
0
投稿猜你喜欢
- 太多的小伙伴正在学习Python，就说自己以后要做全栈开发，大家知道这是做什么的吗？我们现在所知道的知识点，哪些是以后你要从事这个全栈所需要
- Python是一门面向对象的语言。面向对象都有三大特性：封装、继承、多态。下面分别来说说这三大特性:1、封装隐藏对象的属性和实现细节，仅对外
- 一、Python中global与nonlocal 声明如下代码a = 10 def foo(): a = 100执行foo() 结果 a
- 方法一：mylist = [1,2,2,2,2,3,3,3,4,4,4,4]myset = set(mylist)for item in m
- 引言在这篇文章中，我将介绍如何使用YOLOv5构建一个佩戴安全帽检测和识别系统。这个系统可以实时检测图像上人物是否有未佩戴安全帽，并及时进行
初学python，我们必须干点有意思的事！从微信下手吧！头像集样例如下： 大家可以发朋友圈开启辨认大赛哈哈~话不多说，直接上代码，注释我写了- 1. 语句块：{ }之间的部分即为BLOCK语句块。2. 条件语句：if ( expression ) BLOCK;if ( e
- mysql 8.0.20 winx64.zip压缩版安装教程记录如下，分享给大家1.下载MySQL官网：链接直接点击链接也可以下载：mysq
- 任务1、将指定的十进制转换为二进制、八进制、十六进制二进制转换第一种和第二种写法：num=int(input('请输入一个十进制的整
- 在正常项目开发过程中，如果MySQL版本从5.6升级到5.7版本。作为DBA在考虑数据库版本升级带来的影响时，一般会有几个注意点：sql_m
- mysql优化了所以先说说这个，其实这个，很复杂，对于不同的网站，其在线，访问量，帖子数量，网络情况，以及机器配置都有关系，优化不是一次可以
- 概述pydicom是一个常用python DICOM parser。但是，没有提供解析多帧图的示例。本文结合相关函数和DICOM知识做一个简
- 一天不小心把ROOT的权限改到最小了(只能登录,什么都做不了),这可急死我了.重装的话太麻烦，而且里面有很多的用户，一个个重新弄不知道到什么
- 下面的教程总结了Javascript在网页定位方面的相关知识。一、网页的绝对大小和相对大小首先，要明确两个基本概念。一张网页的全部面积，就是
- 这本入门手册是否合适你?我只想告诉你我非常喜欢这本书。我对Microsoft Access的经验足以让我跳过这本傻瓜系列教材，但是实际情况是
- 基于tkinter模块的GUIGUI是图形用户界面的缩写，图形化的用户界面对使用过计算机的人来说应该都不陌生，在此也无需进行赘述。Pytho
- 事务全部是关于原子性的。原子性的概念是指可以把一些事情当做一个单元来看待。从数据库的角度看，它是指应全部执行或全部都不执行的一条或多条语句的
- 一、网络请求在uni中可以调用uni.request方法进行请求网络请求需要注意的是：在小程序中网络相关的 API 在使用前需要配置域名白名
- 本文实例讲述了JS实现简单的抽奖转盘效果。分享给大家供大家参考，具体如下：<!DOCTYPE html><html lan
- preface：做着最近的任务，对数据处理，做些简单的提特征，用机器学习算法跑下程序得出结果，看看哪些特征的组合较好，这一系列流程必然要用到
