dt =9.3600e-09;n=length(v); %v是从excel中读取的1355*2的matrix文件名V=fft(v(:,2));freq=(0:n-1)/(n*dt);Pyy = V.* conj(V) / n;plot(freq,Pyy);