fft谱

ljsr1001

clear
fs=1000
t=0:1/fs:0.6;
f1=200;
x=sin(2*pi*f1*t)
number=1024
y=fft(x,number);
n=0:length(y)-1;
f=fs*n/length(y);
plot(f,abs(y));
我想让MATLAB依据频谱图算出频谱，就是200，如何搞定？

[ 本帖最后由 zhangnan3509 于 2008-3-13 15:02 编辑 ]

zhangnan3509

clear；
fs=1000；
N=1024；
n=0:N-1；
t=n/fs;
f1=200
x=sin(2*pi*f1*t);
y=fft(x,N);
mag=abs(y);
f=(0:length(y)/2-1)*fs/length(y);
plot(f,2*mag(1:N/2)/N);

ljsr1001

哦，那为什么我做误差分析，就是让N值增加，他频率最大的点的误差会有振动呢，