傅里叶变换没有给出预期的结果

Question

情况是从时间空间到频率空间的信号的傅立叶变换将具有与原始信号的频率对应的峰值。我的问题是，当我尝试使用numpy.fft.fft函数绘制傅立叶变换的结果并“适当地”缩放时，我得到一个对应于两倍信号频率的峰值。谁能明白为什么？

import numpy as np
import matplotlib.pyplot as plt

A=2
f0=4
t=np.linspace(0,2,200)
phase=0

s=[]
for i in t:
    s.append(A*np.sin(2*np.pi*f0*i+phase))

N = int(len(s))
yf = np.fft.fft(s)[1:N//2]
xf = np.linspace(1, N//2-1, len(yf))
yf1=yf*np.conj(yf)/N

plt.plot(xf, yf1)
plt.xlim([0,15])