フーリエ変換の図示の方法

Question

前提

ここに質問の内容を詳しく書いてください。
サンプリング周波数20kHzのデータをフーリエ変換したときに5kHzのデータが確認できない

実現したいこ

サンプリング周波数20kHzので録音された5kHzの信号を確認したい

発生している問題・エラーメッセージ

5kHzの周波数帯が確認できない

該当のソースコード

import numpy as np
import matplotlib.pyplot as plt
##########データ抽出#################

A=10
f = 5000.0    # 周波数 Hz
sec = 0.5  # 信号の長さ s
sf = 10000 # サンプリング周波数 Hz
count = sec*sf
count = int(count)
t=[]
y=[]

t = np.arange(0, sec, 1/sf) #サンプリング点の生成
y = A*np.sin(2*np.pi*f*t) # 正弦波の生

# Generate test signal
fs = 20000 # Sampling frequency
N = 10000  # Number of samples
t = np.linspace(0, 1/fs, num=N, endpoint=False)
x = 0.1 * np.sin(2 * np.pi * 1000 * t) + 0.01 * np.sin(2 * np.pi * 5000 * t)

# Perform FFT
X = np.fft.fft(y)
f = np.fft.fftfreq(N, 1/fs)
f = np.fft.fftshift(f)
X = np.fft.fftshift(X)

plt.plot(f, np.abs(X))
plt.xlim(0, 1000)
plt.xlabel('Frequency [Hz]')
plt.ylabel('Magnitude')
plt.show()

試したこと

ここに問題に対して試したことを記載してください。

補足情報（FW/ツールのバージョンなど）

ここにより詳細な情報を記載してください。

Comments

[jbpb0]

質問のコードをそのまま実行したら、
ValueError: x and y must have same first dimension, but have shapes (10000,) and (5000,)
というエラーが出ます

> サンプリング周波数20kHz

なのだから、

sf = 10000 # サンプリング周波数 Hz
↓ 修正
sf = 20000 # サンプリング周波数 Hz

とするのではないですかね
(そうしたらエラー出なくなる)

上記を直した場合ですが、

> plt.xlim(0, 1000)

として0〜1kHzだけ表示させたら、

> 5kHzの周波数帯が確認できない

となるのは当たり前だと思うのですが

その行を無効にしたら、グラフの-5kHzと5kHzにピークがあります

[Ayatakamaru]

すみませんでした，ただの勘違いでした
解決しました

ご指摘ありがとうございます

Answer 1

勘違いでした

jbpb0さんありがとうございます