1kHzのデータを受け取り、FFTを行うプログラムを作成しています。
しかし、データのグラフがどうしても周波数がずれてしまいます。
プログラムは以下のように記述しました。
python
1from __future__ import unicode_literals, print_function 2 3import spidev 4import time 5import csv 6import numpy as np 7import matplotlib.pyplot as plt 8 9 10 11 12spi = spidev.SpiDev() 13spi.open(0,0) 14spi.max_speed_hz = 3600000 15 16 17 18def readAdc(channel): 19 adc = spi.xfer2([1,(8+channel)<<4,0]) 20 data = ((adc[1]&3) << 8) + adc[2] 21 return data 22 23def main(): 24 k = 512 # サンプル数 25 dt = 0.0004 # サンプリング間隔 26 s = np.arange(0, k*dt, dt) # 時間軸 27 freq = np.linspace(0, 1/dt, k) # 周波数軸 28 ch1 = 0 29 ch2 = 1 30 31 f = open('data3.csv','w') 32 dataWriter = csv.writer(f) 33 a1 = [] 34 a2 = [] 35 b1 = {} 36 x = [] 37 b1 = [] 38 39 data1 = {} 40 data2 = {} 41 cal1 = {} 42 cal2 = {} 43 44 45 for t in range(0,k): 46 47 data1[t] = readAdc(ch1)-200 48 49 a1.append(data1[t]) 50 time.sleep(dt) 51 52 for t in range(0,k): 53 54 data2[t] = readAdc(ch2) 55 a2.append(data2[t]) 56 time.sleep(dt) 57 58 59 for t in range(0,k): 60 x.append(0.01*t) 61 cal1[t] = data1[t] * 5.0 / 1024.0 62 cal2[t] = data2[t] * 5.0 / 1024.0 63 64 scal1 = str(cal1[t]) 65 scal2 = str(cal2[t]) 66 dat = (scal1+","+scal2+"\r\n") 67 f.write(dat) 68 time.sleep(dt) 69 70 f.close() 71 F = np.fft.fft(a1) 72 Amp = np.abs(F) 73 74 plt.figure() 75 plt.subplot(121) 76 plt.plot(s, a1) 77 plt.grid() 78 plt.subplot(122) 79 plt.plot(freq, Amp) 80 plt.grid() 81 plt.show() 82 83if __name__ == "__main__": 84 main() 85
1kHzの入力なので本来は1000のところに値が現れるはずがずれてしまいます。
![
正確にFFTをするためにはどのように記述すれば良いのでしょうか。
使用したもの
・AD変換器:MCP3008
・波形発生器:PicoScope2000
回答1件
あなたの回答
tips
プレビュー
バッドをするには、ログインかつ
こちらの条件を満たす必要があります。
2019/07/02 11:37
2019/07/02 13:05
2019/07/02 13:13
2019/07/02 13:32
2019/07/02 16:19
2019/07/03 12:35
2019/07/05 13:58