小編給大家分享一下python中怎么生成任意頻率正弦波，相信大部分人都還不怎么了解，因此分享這篇文章給大家參考一下，希望大家閱讀完這篇文章后大有收獲，下面讓我們一起去了解一下吧！

創(chuàng)新互聯(lián)專注于企業(yè)營銷型網(wǎng)站建設(shè)、網(wǎng)站重做改版、東昌府網(wǎng)站定制設(shè)計(jì)、自適應(yīng)品牌網(wǎng)站建設(shè)、html5、商城網(wǎng)站建設(shè)、集團(tuán)公司官網(wǎng)建設(shè)、成都外貿(mào)網(wǎng)站建設(shè)公司、高端網(wǎng)站制作、響應(yīng)式網(wǎng)頁設(shè)計(jì)等建站業(yè)務(wù)，價(jià)格優(yōu)惠性價(jià)比高，為東昌府等各大城市提供網(wǎng)站開發(fā)制作服務(wù)。

如下所示：

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:為信號(hào)幅值

fi:為信號(hào)頻率

time_s:為時(shí)間長度（s）

sample:為信號(hào)采樣頻率

補(bǔ)充拓展：Python FFT合成波形實(shí)例

使用Python numpy模塊帶的FFT函數(shù)合成矩形波和方波，增加對離散傅里葉變換的理解。

導(dǎo)入模塊

import numpy as np
import matplotlib.pyplot as plt
分別是產(chǎn)生一個(gè)周期的方波和三角波程序

# 產(chǎn)生size點(diǎn)取樣的三角波，其周期為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函數(shù)第二個(gè)值是if，第三個(gè)是else

下面程序可以計(jì)算對應(yīng)的頻譜，采樣點(diǎn)數(shù)取為2的n次冪是為了便于FFT計(jì)算

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

下面對計(jì)算的頻譜進(jìn)行可視化，頻率對應(yīng)的強(qiáng)度使用工程上常用的分貝dp來表示

# 繪制三角波的FFT的前20項(xiàng)的振幅，由于不含下標(biāo)為偶數(shù)的值均為0， 因此取
# log之后無窮小，無法繪圖，用np.clip函數(shù)設(shè)置數(shù)組值的上下限，保證繪圖正確
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")

下面用正弦和余弦函數(shù)合成信號(hào)

# 取FFT計(jì)算的結(jié)果freqs中的前n項(xiàng)進(jìn)行合成，返回合成結(jié)果，計(jì)算loops個(gè)周期的波形
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 # 除去直流成分之外，其余的系數(shù)都*2
    data += np.real(p) * np.cos(k*index) # 余弦成分的系數(shù)為實(shí)數(shù)部
    data -= np.imag(p) * np.sin(k*index) # 正弦成分的系數(shù)為負(fù)的虛數(shù)部
  return index, data

其中index代表頻譜空間的采樣點(diǎn)

畫出合成信號(hào)，x坐標(biāo)使用默認(rèn)的整數(shù)表示即可

# 繪制原始的三角波和用正弦波逐級(jí)合成的結(jié)果，使用取樣點(diǎn)為x軸坐標(biāo)
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) # 計(jì)算兩個(gè)周期的合成波形
  plt.plot(data, label = "N=%s" % i)
plt.legend()
plt.title("partial Fourier series of triangle wave")
plt.show()

以上是“python中怎么生成任意頻率正弦波”這篇文章的所有內(nèi)容，感謝各位的閱讀！相信大家都有了一定的了解，希望分享的內(nèi)容對大家有所幫助，如果還想學(xué)習(xí)更多知識(shí)，歡迎關(guān)注創(chuàng)新互聯(lián)成都網(wǎng)站設(shè)計(jì)公司行業(yè)資訊頻道！

另外有需要云服務(wù)器可以了解下創(chuàng)新互聯(lián)scvps.cn，海內(nèi)外云服務(wù)器15元起步，三天無理由+7*72小時(shí)售后在線，公司持有idc許可證，提供“云服務(wù)器、裸金屬服務(wù)器、高防服務(wù)器、香港服務(wù)器、美國服務(wù)器、虛擬主機(jī)、免備案服務(wù)器”等云主機(jī)租用服務(wù)以及企業(yè)上云的綜合解決方案，具有“安全穩(wěn)定、簡單易用、服務(wù)可用性高、性價(jià)比高”等特點(diǎn)與優(yōu)勢，專為企業(yè)上云打造定制，能夠滿足用戶豐富、多元化的應(yīng)用場景需求。

文章名稱：python中怎么生成任意頻率正弦波-創(chuàng)新互聯(lián)
標(biāo)題路徑：http://www.rwnh.cn/article40/dochho.html

成都網(wǎng)站建設(shè)公司_創(chuàng)新互聯(lián)，為您提供外貿(mào)網(wǎng)站建設(shè)、手機(jī)網(wǎng)站建設(shè)、定制開發(fā)、域名注冊、網(wǎng)站內(nèi)鏈、商城網(wǎng)站

廣告

聲明：本網(wǎng)站發(fā)布的內(nèi)容（圖片、視頻和文字）以用戶投稿、用戶轉(zhuǎn)載內(nèi)容為主，如果涉及侵權(quán)請盡快告知，我們將會(huì)在第一時(shí)間刪除。文章觀點(diǎn)不代表本網(wǎng)站立場，如需處理請聯(lián)系客服。電話：028-86922220；郵箱：631063699@qq.com。內(nèi)容未經(jīng)允許不得轉(zhuǎn)載，或轉(zhuǎn)載時(shí)需注明來源： 創(chuàng)新互聯(lián)