[fname,pname]=uigetfile;
x=xlsread(strcat(pname,fname));
x=x(10:1400,3);
fs=500000;
fs =1400-10+1;
N=length(x);
t=0:1/fs:N/fs-1/fs;
FFTX=abs(fft(x));
FFTX=FFTX*2/N;
FFTX(1)=FFTX(1)/2;
n=1:1:N;
f=(n-1)*fs/N;
subplot(212)
%plot(f,log10(FFTX))
plot(f,FFTX)
xlim([0 fs/2])
作者:LiSQ