ljsr1001 发表于 2008-3-13 13:22

fft谱

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 发表于 2008-3-13 15:13

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 发表于 2008-3-14 10:06

哦,那为什么我做误差分析,就是让N值增加,他频率最大的点的误差会有振动呢,
页: [1]
查看完整版本: fft谱