dt=t(2)-t(1);%%%%dt是点距;lt=length(t);%%%%是数据的长度;f0=2*pi/lt;%%%%f0是基频;fy=fft(yt);%%%%fy是fft变换得到的波谱;
for i=1:1:lt f(i)=f0*(i-1);%%%%f是频率; P(i)=abs(fy(i))^2;%%%P是功率谱;
endplot(f,P); %%%%画图;
首页 >
dt=t(2)-t(1);%%%%dt是点距;lt=length(t);%%%%是数据的长度;f0=2*pi/lt;%%%%f0是基频;fy=fft(yt);%%%%fy是fft变换得到的波谱;
for i=1:1:lt f(i)=f0*(i-1);%%%%f是频率; P(i)=abs(fy(i))^2;%%%P是功率谱;
endplot(f,P); %%%%画图;