python正弦_python生成任意频率正弦波方式
如下所示:
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生成任意频率正弦波方式
本文地址: http://www.cppcns.com/jiaoben/python/301433.html
python正弦_python生成任意频率正弦波方式相关推荐
- python正弦函数拟合_python生成任意频率正弦波方式
如下所示: def signal_xHz(A, fi, time_s, sample): return A * np.sin(np.linspace(0, fi * time_s * 2 * np.p ...
- FPGA 20个例程篇:20.USB2.0/RS232/LAN控制并行DAC输出任意频率正弦波、梯形波、三角波、方波(二)
通过上面的介绍相信大家对数字变频已经有了一个较为整体性的认识,下面笔者来对照XILINX的DDS IP核对数字变频技术展开更进一步的说明,做到了理论和实践很好地结合,这样大家再带入Modelsim进行 ...
- python使用python-barcode生成任意长度字符串的条形码(一维码)详细教程
使用python-barcode生成一维码 1.安装python包 pip install Pillow pip install python-barcode 2.实例程序 # -*- coding: ...
- python 生成图片_python生成带有表格的图片
因为工作中需要,需要生成一个带表格的图片 例如: 直接在html中写一个table标签,然后单独把表格部分保存成图片 或者是直接将excel中的内容保存成一个图片 刚开始的思路,是直接生成一个带有ta ...
- Python实现随机生成任意数量车牌号(完全符合2020年中国车牌号标准)
更改:之前写的还是有漏洞的,之前写的默认后五位全是数字,但是实际上不是这样的,作者百度了相关标准,如下: (一)前4位必须有1位且只能有1位是英文字母(但I和O不能使用),英文字母位置可以在4位中自由 ...
- python双色球_python生成双色球程序
偶尔去买双色球,每次购买时都是随机选球,自己平时就用python,当然一直在学习,所以就用python写了一个随机或指定出双色球的脚本. 代码如下:#!/usr/bin/python # -*- co ...
- python 中文名_python生成随机中文姓名和拼音
def sayhello(): print('hellopython') import pinyin import random def random_name(): # 删减部分,比较大众化姓氏 f ...
- python报表_Python生成报表
excel def export_teachers_excel(request): # 创建工作簿 wb = xlwt.Workbook() # 添加工作表 sheet = wb.add_sheet( ...
- python随机生成车牌_Python实现随机生成任意数量车牌号
之前做课设的时候舍友遇到了需要生成500w量级车牌号的问题,于是我便写了一个随机生成车牌号的程序,希望各位采纳. 注:Python实现 import random def chepaihao(len= ...
最新文章
- 电机测试matlab,基于MATLAB的电机故障诊断模型训练与测试平台
- python的核心数据类型_Python核心数据类型-集合
- android开发仿原生联系人,android 仿IOS联系人列表
- 接收list对象_PyTorch入门视频笔记从数组、列表对象中创建Tensor
- Bailian3179 最长单词【字符串】
- RAW数码照片处理器:SILKYPIX Developer 的功能介绍
- numpy的广播机制
- 怎么安装mysql5.6.19_centos6.5 rpm安装mysql5.6.19操作及步骤
- java double 移位_【原创】Java移位运算
- PHP 照片合并视频
- 10个超强英语学习资源网站,学霸们都在用!
- 39元超值!360超级充电器拆解与评测
- web项目js调用斑马打印机打印二维码
- APP测试基础--小工具介绍(1)
- 计算机网络的概念,组成和分类;计算机与网络信息安全的概念和防控.,计算机二级MS试题...
- Oracle CTAS
- 如何搜索百度云盘里的文件
- moment 获得上月开始日期和结束日期
- linux内核snat分析,(十)洞悉linux下的Netfilteriptables:网络地址转换原理之SNAT
- feign post 传递空值_别跟我说你真的知道 post 和 get 的本质区别?!!
热门文章
- 数据库事务的简单学习(一)
- AI智慧安监:打电话/玩手机智能检测,构筑安全生产新防线
- Matconvnet 构建自己的网络
- 降低开关电路振铃现象的缓冲器设计方法
- A股月份效应 | Python量化A股市场魔咒,5穷6绝7翻身准确吗?| 邢不行
- 腾讯 搅局者还是终结者
- “善行天下 玉满中华—中国珠宝玉石原石交易平台第三批交易品上市敲钟仪式”新闻发布会
- 批量计算机添加网络打印机方法,添加局域网打印机(网络打印机)的图文方法...
- 马斯克借势ChatGPT强推1000美元认证服务,《纽约时报》拒付被开除大V
- Go语言休眠sleep