python生成任意频率正弦波方式
作者:宇宙小虾米 时间:2021-08-27 17:08:07
如下所示:
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()
来源:https://blog.csdn.net/shangxindewunai/article/details/84031218
标签:python,正弦波
![](/images/zang.png)
![](/images/jiucuo.png)
猜你喜欢
python K近邻算法的kd树实现
2022-01-09 19:05:43
Python 实现取矩阵的部分列,保存为一个新的矩阵方法
2023-07-17 23:22:30
python冒泡排序算法的实现代码
2021-09-08 03:48:27
php二分查找二种实现示例
2023-11-21 00:40:13
使用Cython中prange函数实现for循环的并行
2023-04-13 05:31:55
![](https://img.aspxhome.com/file/2023/3/76773_0s.png)
Python从ZabbixAPI获取信息及实现Zabbix-API 监控的方法
2021-01-16 18:41:38
![](https://img.aspxhome.com/file/2023/1/90711_0s.png)
ASPImage组件的实现过程[图]
2008-02-03 15:37:00
![](https://img.aspxhome.com/file/UploadPic/20082/3/200823154516685s.gif)
Python 解析库json及jsonpath pickle的实现
2021-09-10 00:27:14
![](https://img.aspxhome.com/file/2023/5/87925_0s.jpg)
Python+eval函数实现动态地计算数学表达式详解
2023-08-13 22:39:11
python实时获取外部程序输出结果的方法
2023-04-27 02:29:04
让ThinkPHP的模板引擎达到最佳效率的方法详解
2023-11-14 17:06:40
python实现列表的排序方法分享
2022-10-06 07:56:01
![](https://img.aspxhome.com/file/2023/0/95760_0s.png)
深度理解Python中Class类、Object类、Type元类
2022-07-14 03:10:13
![](https://img.aspxhome.com/file/2023/9/105059_0s.png)
基于pytorch的保存和加载模型参数的方法
2023-02-07 00:24:12
ASP代码中的容错机制
2008-10-19 17:30:00
FusionCharts图表显示双Y轴双(多)曲线
2023-08-22 17:55:38
![](https://img.aspxhome.com/file/2023/9/56199_0s.png)
python多个模块py文件的数据共享实例
2022-02-24 23:33:39
教你快速掌握SQL语言中游标的使用技巧
2009-01-08 16:24:00
python 执行终端/控制台命令的例子
2022-11-16 20:41:29
[精品]ASP中常用的22个FSO文件操作函数
2007-08-18 15:12:00
![](https://img.aspxhome.com/file/UploadPic/up/2007081816282483.jpg)