Python实现正弦信号的时域波形和频谱图示例【基于matplotlib】
作者：落落图灵 发布时间：2023-03-18 02:18:43
标签：Python,正弦,频谱图,matplotlib
本文实例讲述了Python实现正弦信号的时域波形和频谱图。分享给大家供大家参考，具体如下：
# -*- coding: utf-8 -*-
# 正弦信号的时域波形与频谱图
import numpy as np
import matplotlib.pyplot as pl
import matplotlib
import math
import random
row = 4
col = 4
N = 500
fs = 5
n = [2*math.pi*fs*t/N for t in range(N)] # 生成了500个介于0.0-31.35之间的点
# print n
axis_x = np.linspace(0,3,num=N)
#频率为5Hz的正弦信号
x = [math.sin(i) for i in n]
pl.subplot(221)
pl.plot(axis_x,x)
pl.title(u'5Hz的正弦信号',fontproperties='SimHei')
pl.axis('tight')
#频率为5Hz、幅值为3的正弦+噪声
x1 = [random.gauss(0,0.5) for i in range(N)]
xx = []
#有没有直接两个列表对应项相加的方式？？
for i in range(len(x)):
xx.append(x[i]*3 + x1[i])
pl.subplot(222)
pl.plot(axis_x,xx)
pl.title(u'频率为5Hz、幅值为3的正弦+噪声',fontproperties='SimHei')
pl.axis('tight')
#频谱绘制
xf = np.fft.fft(x)
xf_abs = np.fft.fftshift(abs(xf))
axis_xf = np.linspace(-N/2,N/2-1,num=N)
pl.subplot(223)
pl.title(u'频率为5Hz的正弦频谱图',fontproperties='SimHei')
pl.plot(axis_xf,xf_abs)
pl.axis('tight')
#频谱绘制
xf = np.fft.fft(xx)
xf_abs = np.fft.fftshift(abs(xf))
pl.subplot(224)
pl.title(u'频率为5Hz的正弦频谱图',fontproperties='SimHei')
pl.plot(axis_xf,xf_abs)
pl.axis('tight')
pl.show()
运行效果：
希望本文所述对大家Python程序设计有所帮助。
来源：https://blog.csdn.net/matrix_google/article/details/60964638
0
投稿猜你喜欢
前段时间写了一个百度图片下载器，结果发现有很多人需要使用。说实话之前写的那一款百度图片下载器比较LOW，今天刚好有时间就做了一下升级。更新了- NumPyNumPy是一个用于科学计算和数据分析的Python库，也是机器学习的支柱。可以说NumPy奠定了Python在机器学习中的地位。
- function toBreakWord(intLen, id){ var obj=document.getElementById(id);
- 前几天小芳同学一直在群发起一些加速的话题，我已经把聊天记录抽出来，正打算整理出份像样的，没想到小芳同学非常速度的出了这篇。我的就省掉了，挖哈
- 本文实例讲述了Vue指令v-for遍历输出JavaScript数组及json对象的常见方式。分享给大家供大家参考，具体如下：定义数据：<
- 本文实例讲述了js点击事件的执行过程。分享给大家供大家参考，具体如下：js事件对象event包括很多事件类型，这里用onclick事件为例，
- 本文实例为大家分享了wxPython实现计算器的具体代码，供大家参考，具体内容如下＃ -*- coding: utf-8 -*-######
- 一、目的之前在博文SQL Server数据库最小宕机迁移方案中提到了使用了完全备份+差异备份的功能完成了数据库的转移，但是这个方法在遇到了7
- 前言对于我这种英语比较差的人来说，无论是敲代码还是看文档，那都是离不开翻译软件的，于是我想自己用python做一个翻译软件，花了一个小时，终
- 本教程为大家分享了Linux安装MySQL详细步骤，供大家参考，具体内容如下第一步： 下载MySQL安装包进入mysql官网，进入downl
- 如今，体育运动的热潮日益流行。同样，以不正确的方式进行运动的风险也在增加。有时可能会导致严重的伤害。考虑到这些原因，提出一种以分析运动员的关
- 为了让鼠标移到小图上显示大图,我利用鼠标事件新建了一个层来显示大图.当然之前最好得到XY坐标取得当前鼠标的X,Y坐标：function&nb
- 当我们经常在代码中使用 Python 字符串时，您可能需要以相反的顺序使用它们。Python 包含一些方便的工具和技术，可以在这些情况下为您
- 早在02年，国外关于SQL注入漏洞的技术文章已经很多，而国内在05年左右才开始的。 如今，谈SQL注入漏洞是否已是明日黄花，国内大大小小的网
- SQL Server所谓的分布式查询（Distributed Query）是能够访问存放在同一部计算机或不同计算机上的SQL Server或
- 前言前几天去我姐家里蹭饭吃，发现我姐家里的小外甥女儿已经上小学了，正在疯狂赶作业，加减乘除还有很多题目都经常出错。我姐说老师每天给他们布置了
- 前言Blender 并不是唯一一款允许你为场景编程和自动化任务的3D软件; 随着每一个新版本的推出，Blender 正逐渐成为一个可靠的 C
- 如下所示：$preg= '/xue[\s\S]*?om/i';preg_match_all($preg,"学并思网
- 外部数据导入导入excel文件pandas导入excel用read_excel()方法：import pandas as pdexcel_f
- 目录文件读写读文件写文件StringIO和BytesIOStringIOBytesIO操作文件和目录环境变量操作文件和目录文件读写读文件tr
