python生成任意频率正弦波方式


Posted in Python onFebruary 25, 2020

如下所示:

def signal_xHz(A, fi, time_s, sample):

return A * np.sin(np.linspace(0, fi * time_s * 2 * np.pi , sample* time_s))

A:为信号幅值

fi:为信号频率

time_s:为时间长度(s)

sample:为信号采样频率

补充拓展:Python FFT合成波形实例

使用Python numpy模块带的FFT函数合成矩形波和方波,增加对离散傅里叶变换的理解。

导入模块

import numpy as np
import matplotlib.pyplot as plt
分别是产生一个周期的方波和三角波程序

# 产生size点取样的三角波,其周期为1
def triangle_wave(size):
  x = np.arange(0, 1, 1.0/size)
  y = np.where(x<0.5, x, 0)
  y = np.where(x>=0.5, 1-x, y)
  return x, y
 
def square_wave(size):
  x = np.arange(0, 1, 1.0/size)
  y = np.where(x<0.5, 1.0, 0)
  return x, y

其中np.where函数第二个值是if,第三个是else

下面程序可以计算对应的频谱,采样点数取为2的n次幂是为了便于FFT计算

fft_size = 256
 
# 计算三角波和其FFT
x, y = triangle_wave(fft_size)
fy = np.fft.fft(y) / fft_size

下面对计算的频谱进行可视化,频率对应的强度使用工程上常用的分贝dp来表示

# 绘制三角波的FFT的前20项的振幅,由于不含下标为偶数的值均为0, 因此取
# log之后无穷小,无法绘图,用np.clip函数设置数组值的上下限,保证绘图正确
plt.figure()
plt.plot(np.clip(20*np.log10(np.abs(fy[:20])), -120, 120), "o")
plt.xlabel("frequency bin")
plt.ylabel("power(dB)")
plt.title("FFT result of triangle wave")

下面用正弦和余弦函数合成信号

# 取FFT计算的结果freqs中的前n项进行合成,返回合成结果,计算loops个周期的波形
def fft_combine(freqs, n, loops=1):
  length = len(freqs) * loops
  data = np.zeros(length)
  index = loops * np.arange(0, length, 1.0) / length * (2 * np.pi)
  for k, p in enumerate(freqs[:n]):
    if k != 0: p *= 2 # 除去直流成分之外,其余的系数都*2
    data += np.real(p) * np.cos(k*index) # 余弦成分的系数为实数部
    data -= np.imag(p) * np.sin(k*index) # 正弦成分的系数为负的虚数部
  return index, data

其中index代表频谱空间的采样点

画出合成信号,x坐标使用默认的整数表示即可

# 绘制原始的三角波和用正弦波逐级合成的结果,使用取样点为x轴坐标
plt.figure()
plt.plot(y, label="original triangle", linewidth=2)
for i in [0,1,3,5,7,9]:
  index, data = fft_combine(fy, i+1, 2) # 计算两个周期的合成波形
  plt.plot(data, label = "N=%s" % i)
plt.legend()
plt.title("partial Fourier series of triangle wave")
plt.show()

以上这篇python生成任意频率正弦波方式就是小编分享给大家的全部内容了,希望能给大家一个参考,也希望大家多多支持三水点靠木。

Python 相关文章推荐
python读取csv文件示例(python操作csv)
Mar 11 Python
关于Python元祖,列表,字典,集合的比较
Jan 06 Python
Python利用matplotlib生成图片背景及图例透明的效果
Apr 27 Python
Python 找到列表中满足某些条件的元素方法
Jun 26 Python
python 内置模块详解
Jan 01 Python
pycharm+PyQt5+python最新开发环境配置(踩坑)
Feb 11 Python
jupyter lab的目录调整及设置默认浏览器为chrome的方法
Apr 10 Python
Selenium使用Chrome模拟手机浏览器方法解析
Apr 10 Python
python修改微信和支付宝步数的示例代码
Oct 12 Python
jupyter使用自动补全和切换默认浏览器的方法
Nov 18 Python
python 基于opencv 绘制图像轮廓
Dec 11 Python
如何使用python-opencv批量生成带噪点噪线的数字验证码
Dec 21 Python
python numpy库linspace相同间隔采样的实现
Feb 25 #Python
Pandas时间序列:时期(period)及其算术运算详解
Feb 25 #Python
基于pygame实现童年掌机打砖块游戏
Feb 25 #Python
python GUI库图形界面开发之PyQt5拖放控件实例详解
Feb 25 #Python
python GUI库图形界面开发之PyQt5美化窗体与控件(异形窗体)实例
Feb 25 #Python
Python对wav文件的重采样实例
Feb 25 #Python
python实现打砖块游戏
Feb 25 #Python
You might like
我的论坛源代码(六)
2006/10/09 PHP
PHP新手上路(十四)
2006/10/09 PHP
探讨php define()函数及defined()函数使用详解
2013/06/09 PHP
PHP设置一边执行一边输出结果的代码
2013/09/30 PHP
Php-Redis安装测试笔记
2015/03/05 PHP
理解Javascript_09_Function与Object
2010/10/16 Javascript
jquery下将选择的checkbox的id组成字符串的方法
2010/11/28 Javascript
基于jQuery的message插件实现右下角弹出消息框
2011/01/11 Javascript
dojo学习第二天 ajax异步请求之绑定列表
2011/08/29 Javascript
JQuery 返回布尔值Is()条件判断方法代码
2012/05/14 Javascript
JS定时刷新页面及跳转页面的方法
2013/07/04 Javascript
js中取得变量绝对值的方法
2015/01/03 Javascript
jquery $(document).ready()和window.onload的区别浅析
2015/02/04 Javascript
js检测判断日期大于多少天的方法
2015/05/04 Javascript
JS实现模拟风力的雪花飘落效果
2015/05/13 Javascript
jQuery事件绑定on()、bind()与delegate() 方法详解
2015/06/03 Javascript
js字符串引用的两种方式(必看)
2016/09/18 Javascript
微信开发之调起摄像头、本地展示图片、上传下载图片实例
2016/12/08 Javascript
利用canvas实现的加载动画效果实例代码
2017/07/05 Javascript
JS弹窗 JS弹出DIV并使整个页面背景变暗功能的实现代码
2018/04/21 Javascript
如何使用50行javaScript代码实现简单版的call,apply,bind
2019/08/14 Javascript
Vue单页面应用中实现Markdown渲染
2021/02/14 Vue.js
在arcgis使用python脚本进行字段计算时是如何解决中文问题的
2015/10/18 Python
带你了解python装饰器
2017/06/15 Python
python实现整数的二进制循环移位
2019/03/08 Python
计算机二级python学习教程(2) python语言基本语法元素
2019/05/16 Python
python matplotlib:plt.scatter() 大小和颜色参数详解
2020/04/14 Python
Python pandas 列转行操作详解(类似hive中explode方法)
2020/05/18 Python
美国牛仔品牌:True Religion
2018/11/16 全球购物
一套SQL笔试题
2016/08/14 面试题
DTD的含义以及作用
2014/01/26 面试题
师范生自荐信
2013/10/27 职场文书
促销活动总结报告
2014/04/26 职场文书
OpenCV3.3+Python3.6实现图片高斯模糊
2021/05/18 Python
剑指Offer之Java算法习题精讲二叉树的构造和遍历
2022/03/21 Java/Android
Python实现科学占卜 让视频自动打码
2022/04/09 Python