Python科學計算（兩）——時域波形和正弦訊號的頻譜

標籤：

Python科學計算（兩）-- 時域和頻域波形為正弦波形訊號產生、計算和顯示

# -*- coding: utf-8 -*-import numpy as npimport matplotlib.pyplot as plimport matplotlibimport mathimport randomrow = 4col = 4N = 500fs = 5n = [2*math.pi*fs*t/N for t in range(N)]axis_x = np.linspace(0,1,num=N)myfont = matplotlib.font_manager.FontProperties(fname='c:\\windows\\fonts\\fzshjw_0.ttf')#頻率為5Hz的正弦訊號x = [math.sin(i) for i in n]pl.subplot(221)pl.plot(axis_x,x)pl.title(u'5Hz的正弦訊號', fontproperties=myfont)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=myfont)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=myfont)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=myfont)pl.plot(axis_xf,xf_abs)pl.axis('tight')

有疑問的留簡介討論

著作權聲明：本文部落格原創文章。部落格，未經同意，不得轉載。

Python科學計算（兩）——時域波形和正弦訊號的頻譜