Python科学計算(二)——正弦波信号の時間領域波形とスペクトル図
Python科学計算(二)--正弦波信号の時間領域波形と周波数領域波形の生成、計算と表示
疑問のコメント検討しましょう
# -*- coding: utf-8 -*-
import numpy as np
import matplotlib.pyplot as pl
import matplotlib
import math
import random
row = 4
col = 4
N = 500
fs = 5
n = [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')
疑問のコメント検討しましょう