clc
data=sin(1:0.01:1000); %模拟data
t2=10; %取前10秒的资料
fs=500; %每秒纪录500个数据 即频率为500
y=data(1:fs*t2);
Y=fft(y,1024); %计算1024点的fft
f=linspace(0,fs/2,1024/2+1); %设定一个频率向量，到达Nyquist freq
figure(1),plot([1:length(y)]/fs,y)
xlabel('秒')
figure(2),semilogy(f,abs(Y(1:1024/2+1)))
网络上找到的程式 想要画出频谱图
跑程式会出现两张图 一张是sin波模拟数据图 10秒钟有接近8个全波长(频率接近0.8)
第二张图是频谱图 会发现频谱的尖峰出现在x为0.976左右的位置 而非接近0.8
请问是哪里需要修正呢？ 感谢