1樓:匿名使用者
我把你的資料
來儲存到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));%下:頻譜圖
%因為頻譜圖是對稱的,所以畫一半就可以了。縱軸是幅值,橫軸是頻率
如何用matlab實現對離散資料的快速傅利葉變換
2樓:匿名使用者
1、雙來擊matlab軟體圖示,開啟matlab軟體,可以自看到matlab軟體的介面。
2、在命令列視窗中輸入:x=0:2:22;y=2*exp(x).*sin(x);建立了12個原始資料點。
3、在命令列視窗中輸入:xi=0:0.1:22;建立要進行插值的資料的橫座標。
4、在命令列視窗中輸入:yi=spline(x,y,xi);使用函式spline(x,y,xi);可以獲得對原始資料的三次樣條插值的y軸座標。
5、在命令列視窗中輸入:plot(x,y,'o',xi,yi);在影象中繪製原始資料點和三次樣條插值的資料點圖示。
6、在命令列視窗中輸入:title('三次樣條插值');xlabel('x');ylabel('y');給繪製的圖示新增標題和座標軸的標籤。
7、最後檢視繪製的原始資料點圖和三次樣條插值圖示,注意標題、橫座標、縱座標等。
3樓:匿名使用者
那是因為fft的結果數值之間差異很大,所以在圖上顯示不出來,可以只顯示部分區間的數值,就是限定座標值範圍。
4樓:匿名使用者
請問你最後解決了這個問題了嗎??我現在也遇到同樣的問題,想問下需要在matlab裡面先找到擬合的函式,在傅利葉變換嗎??還是直接用fft就行了??
matlab中如何對一組資料進行fft變換後得到頻譜圖??
5樓:用著追她
1、直接對其使用來ceil函式,源
可以向上取整:即大於等於給定數bai據的最du小整數。
2、使用zhifloor函式可以進行向下取dao整。如圖,將矩陣a中的每個元素取小於等於該元素的最大整數。
3、使用fix函式則可以去除小數部分。對於大於0的元素相當於floor,對於小於0的元素相當於ceil。
4、使用round函式可以把資料四捨五入到最近的整數。
5、round還可以帶有第二個引數,精確到位數。負數表示四捨五入到小數點左側多少位,正數表示小數點右側。0表示精確到個位,預設。
6、當round帶有第三個引數'significant'時,第二個引數表示精確位數。如圖表示每個元素均取2位有效數字。
6樓:手機使用者
將那個i輸出線分0一t支q出來輸入i到一b個kto file的模組就好了d,如果要轉到execl中2的話將他們鏈結就好e┳
怎麼用matlab做傅利葉變換
7樓:沒錯這是暱稱
1、我們使用matlab開發的傅利葉變換程式**如下:內
syms x
f = exp(-2*x^2); %our function
ezplot(f,[-2,2]) % plot of our function
ft = fourier(f) % fourier transform
將其寫容入到我們的matlab程式模組中。
2、我們執行上面的傅利葉變換程式**,將得出執行結果:ft = (2^(1/2)*pi^(1/2)*exp(-w^2/8))/2。
3、如果我們需要更高階的顯示,我們修改上述**即可,如使用ezplot(ft)作傅利葉變換折線圖。
以上就是如何用matlab做傅利葉變換的解決步驟。
8樓:匿名使用者
opencv計算機視覺14(傅利葉變換)
9樓:珞珈山石
fs=100;n=128; %取樣頻率和資料bai
點數n=0:n-1;t=n/fs; %時間序列du
x=0.5*sin(2*pi*15*t)+2*sin(2*pi*40*t); %訊號
y=fft(x,n); %對訊號進行快速
zhifourier變換
mag=abs(y); %求得daofourier變換後的振幅回
f=n*fs/n; %頻率序列
subplot(1,2,1),plot(f,mag); %繪出隨頻率變化
答的振幅
xlabel('頻率/hz');
ylabel('振幅');title('n=128');grid on;
subplot(1,2,2),plot(f(1:n/2),mag(1:n/2)); %繪出nyquist頻率之前隨頻率變化的振幅
xlabel('頻率/hz');
ylabel('振幅');title('n=128');grid on;
10樓:人文漫步者
你想利用這乙個軟體來實現福利業函式的乙個變化過程,就需要了解他的一些具體工作情況。
離散傅利葉變換,離散傅利葉變換DFT和離散時間傅利葉變換DTFT的區別
1.週期序列與有限長序列的關係 如上所述,有限長序列x n 可以看成是週期序列 只取乙個週期的結果,而週期序列 則是有限長序列x n 的週期延拓序列,即 物探數碼訊號分析與處理技術 一般稱x n 的第乙個週期從n 0到n 1的值為主值區間,所以說週期序列 是有限長序列x n 的週期開拓,而x n 是...
傅利葉變換的意義和理解
傅利葉變換的意義和理解如下 意義 傅利葉變換是數學中最深刻的見解之一,但不幸的是,它的意義深埋在一些枯燥的方程中。我們都知道傅利葉級數是一種可以把任意週期函式分解成一堆正弦波的方法。和往常一樣,這個名字來自乙個生活在很久以前的人,他叫傅利葉。在數學術語中,傅利葉變換是一種將訊號轉換成頻率的技術,即從...
為什麼要進行傅利葉變換,變換後得到的函式究竟是什麼
好問題。1.不知道你還記得傅利葉變換是怎麼來的不,至少在課本上看到的是根據週期函式的傅利葉級數的推廣 傅利葉級數告訴我們任意週期函式 這裡討論連續的情況 均可以分解為基頻及其諧波成分的疊加。而傅利葉先生當年在解決熱力學問題時將這個idea推廣了一下,就是現在的傅利葉變換。我們將週期函式的週期設為無窮...