當前位置:成語大全網 - 新華字典 - matlab 傅裏葉變換我想要對壹組我測的數據進行傅裏葉變換,已知測量的各個時間t以及對應時間的數據值x。比

matlab 傅裏葉變換我想要對壹組我測的數據進行傅裏葉變換,已知測量的各個時間t以及對應時間的數據值x。比

我把妳的數據保存到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));%下:頻譜圖

%因為頻譜圖是對稱的,所以畫壹半就可以了。縱軸是幅值,橫軸是頻率