Question

怎么对一组采集到的数据在MATLAB上做频谱分析？

数据是1000个，分别是在这个时间区间采集的 t=0：0.00004：0.03996；怎么写MATLAB程序？

fs=25000;
n=1000;
t=0:1/fs:0.03996;
base=importdata('A.TXT');
y1=fft(base(:,2));
y2=fftshift(y1);
f=(0:n-1)*fs/n-fs/2;
plot(f,abs(y1),'r');
plot(f,abs(y2),'b');
axis([0,10000,0,500000]);
这样对不对？麻烦大神看一下！谢谢了，在线等。

Answer 1

问题不大。改成如下
fs=25000;
n=1000;
t=0:1/fs:0.03996;
base=importdata('A.TXT');
y1=fft(base(:,2));
y2=fftshift(y1);
f=(0:n-1)*fs/n-fs/2;
figure(1)
plot(t,base(:,2),'r');%原始采样图
figure(2)
plot(f,abs(y2),'b');%频谱图
axis([0,10000,0,500000]);

Answer 2

　　代码：
　　A=keshiyin_moni; %这里的A=后面写文件的名字
FS=50;T=1/FS; % FS是采样率
L=length(A);
t=(0:L-1)*T;
figure(1);
plot(t,A);
figure(6);
NFFT = 2^nextpow2(L); % Next power of 2 from length of y
Y = fft(A,NFFT)/L;
f = FS/2*linspace(0,1,NFFT/2);
plot(f,2*abs(Y(1:NFFT/2)));