我把妳的數據保存到data.txt,放到matlab工作目錄裏,然後打下面命令,或者保存成m文件
a=load('data.txt');%讀數據 t=a(:,1)';%時間 x=a(:,2)';%數據 dt=t(2)-t(1);%用diff(t)看了下,妳的采樣時間間隔壹樣,就只算了壹個差值 fs=1/dt;%傅裏葉變換後的坐標軸的最大值就是采樣頻率fs,與dt是倒數關系 ax=abs(fft(x));%fft是matlab中傅裏葉變換的命令,abs是求模,我畫的是幅值譜 f=(0:length(t)-1)*fs/length(t);%頻譜圖橫軸,總***49個點,最大值是fs,每個點就是fs/49 figure();%畫在壹起了 subplot(211);plot(t,x);%上:時域圖 subplot(212);plot(f(1:25),ax(1:25));%下:頻譜圖%因為頻譜圖是對稱的,所以畫壹半就可以了。縱軸是幅值,橫軸是頻率