本帖產(chǎn)生 1 個(gè) 程序強(qiáng)帖，點(diǎn)擊這里進(jìn)行查看

當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖

whfire

金蟲 (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 180.1
散金: 34
紅花: 4
帖子: 872
在線: 143.9小時(shí)
蟲號(hào): 1147444
注冊(cè): 2010-11-15
性別: GG
專業(yè): 生化分析及生物傳感

[求助] matlab 傅里葉變換頻譜分析

我想利用matlab進(jìn)行編程，進(jìn)行傅里葉頻譜分析，但是需要將頻譜轉(zhuǎn)換成功

率譜。詳情見附件

回復(fù)此樓

1樓 2011-07-25 19:45:48

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

mallory19

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
帖子: 1
在線: 19分鐘
蟲號(hào): 2888375
注冊(cè): 2013-12-22

引用回帖:

3樓: Originally posted by sudo at 2011-07-25 20:44:28
對(duì)了，還有一些想補(bǔ)充的~

首先就是直流分量的問題，如果不先對(duì)數(shù)據(jù)都減去一個(gè)均值，那么接近0Hz部分的功率會(huì)非常非常大~這是有害的，一般做數(shù)據(jù)處理的時(shí)候如果只關(guān)注交流量的話，就要先使得整個(gè)序列的均值為0，方 ...

sudo，你好。
我在用Mathematica中的Periodogram和PeriodogramArray對(duì)同一個(gè)數(shù)據(jù)列表計(jì)算時(shí)，得到的結(jié)果卻不同。我在網(wǎng)上查了一下，也有別人在問這個(gè)問題。但沒有回答。
不知道你對(duì)Mathematica是否了解，是否可以交流一下。

贊一下

回復(fù)此樓

9樓2014-04-03 23:53:13

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

查看全部 13 個(gè)回答

sudo

木蟲 (正式寫手)

程序強(qiáng)帖: 16
應(yīng)助: 6 (幼兒園)
金幣: 1297.6
散金: 1486
紅花: 20
帖子: 588
在線: 641小時(shí)
蟲號(hào): 1211394
注冊(cè): 2011-02-24
性別: GG
專業(yè): 文藝美學(xué)

★
dubo(金幣+1): 歡迎常來程序語言版討論 2011-07-31 14:25:14

唉，可能解釋起來一言難盡~建議樓主閱讀一下信號(hào)處理的入門書嗯

首先，你師兄的圖，應(yīng)該是降采樣過后的（或者是后來截取了一段，不過我傾向于認(rèn)為是降采樣），也就是說他用來做功率譜的數(shù)據(jù)采樣率只有400Hz，你給的數(shù)據(jù)的采樣率是2000Hz

其次，你師兄給的圖是沒有取dB作為縱坐標(biāo)單位的，也就是沒有計(jì)算10*log10(P)，而是直接用“功率密度”來做縱坐標(biāo)~

最后，其實(shí)MATLAB有現(xiàn)成函數(shù)可以做功率譜：

CODE:

[P,F]=periodogram(x, window, nfft, Fs)

其中輸出P和F分別是功率譜和對(duì)應(yīng)的頻率
輸入x是信號(hào)，window是數(shù)據(jù)窗口（可以寫[]用默認(rèn)矩形窗），nfft是進(jìn)行fft變換的點(diǎn)數(shù)（可以寫[]取默認(rèn)值，一般會(huì)取到最近的2的整數(shù)次冪左右的地方），F(xiàn)s是采樣率，單位是Hz

然后我對(duì)你的數(shù)據(jù)裸進(jìn)行功率譜處理的話，結(jié)果如下圖：

贊一下(8人)

回復(fù)此樓

2樓2011-07-25 20:28:50

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

sudo

木蟲 (正式寫手)

程序強(qiáng)帖: 16
應(yīng)助: 6 (幼兒園)
金幣: 1297.6
散金: 1486
紅花: 20
帖子: 588
在線: 641小時(shí)
蟲號(hào): 1211394
注冊(cè): 2011-02-24
性別: GG
專業(yè): 文藝美學(xué)

【答案】應(yīng)助回帖

whfire(金幣+20): 哥們，謝謝。太謝謝了 2011-07-27 09:13:49

對(duì)了，還有一些想補(bǔ)充的~

首先就是直流分量的問題，如果不先對(duì)數(shù)據(jù)都減去一個(gè)均值，那么接近0Hz部分的功率會(huì)非常非常大~這是有害的，一般做數(shù)據(jù)處理的時(shí)候如果只關(guān)注交流量的話，就要先使得整個(gè)序列的均值為0，方法就是給每個(gè)數(shù)減去該序列的一個(gè)均值

然后就是降采樣的方法，其實(shí)很簡單啦，在MATLAB里面的話，如果y是你原來的數(shù)據(jù)的話（注意是測量數(shù)據(jù)，你給的數(shù)據(jù)的第2列）.....還是給出命令的全部吧：

CODE:

ty = y(1:5:end); %將原來2000Hz的數(shù)據(jù)降采樣到400Hz

ty = ty - mean(ty); %減去直流分量

[P,F]=periodogram(ty,[],[],400); %求功率譜

plot(F, P) %畫圖

結(jié)果就是

贊一下(6人)

回復(fù)此樓

» 本帖已獲得的紅花（最新10朵）

獨(dú)孤_逍遙

3樓2011-07-25 20:44:28

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

sudo

木蟲 (正式寫手)

程序強(qiáng)帖: 16
應(yīng)助: 6 (幼兒園)
金幣: 1297.6
散金: 1486
紅花: 20
帖子: 588
在線: 641小時(shí)
蟲號(hào): 1211394
注冊(cè): 2011-02-24
性別: GG
專業(yè): 文藝美學(xué)

【答案】應(yīng)助回帖

★ ★ ★ ★ ★ ★ ★
xzhdty(金幣+2): 歡迎常來 2011-07-25 20:57:22
dubo(金幣+5, 程序強(qiáng)帖+1): 感謝應(yīng)助 2011-07-31 14:26:41

另外說明：MATLAB的功率譜實(shí)現(xiàn)函數(shù)是像這樣子的：

CODE:

function [P,f]=fftpsd(x,Fs)

    m = length(x);          % 整個(gè)數(shù)據(jù)長度

    n = pow2(nextpow2(m));  % 擴(kuò)展到合適長度，2的整數(shù)次冪

    y = fft(x,n);           

    f = (0:n-1)*(Fs/n);     %頻率范圍

    P = y.*conj(y)/n;   %功率譜密度

    f=f(1:n/2);

    P=P(1:n/2);

    %plot(f,10*log10(P)); %把數(shù)據(jù)轉(zhuǎn)換為dB

    plot(f,P);

    grid;

    xlabel('Frequency (Hz)');

    ylabel('Power(dB)');

    title('{\bf Periodogram}');

為了使用FFT算法（速度快），MATLAB把數(shù)據(jù)擴(kuò)展了，本來降采樣你的數(shù)據(jù)之后，一共有4000個(gè)點(diǎn)，但是實(shí)際上使用了4096個(gè)點(diǎn)。但是由于數(shù)據(jù)是實(shí)數(shù)序列，所以，僅有采樣率一半的數(shù)據(jù)是有效的（后半部分其實(shí)是鏡像），所以在圖上，F(xiàn)和P都是2048個(gè)點(diǎn)（periodogram取了2049個(gè)，這個(gè)無所謂）

贊一下(7人)

回復(fù)此樓

4樓2011-07-25 20:53:06

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[論文投稿] 急發(fā)核心期刊論文 +3	賢達(dá)問津 2026-03-23	5/250	2026-03-23 17:13 by 妹子不好惹
[考研] 328求調(diào)劑 +4	LHHL66 2026-03-23	4/200	2026-03-23 14:55 by lbsjt
[考研] 263求調(diào)劑 +6	yqdszhdap－ 2026-03-22	9/450	2026-03-23 12:57 by yqdszhdap－
[考研] 一志愿東華大學(xué)化學(xué)070300，求調(diào)劑 +7	2117205181 2026-03-21	8/400	2026-03-22 22:55 by chixmc
[考研] 315分，誠求調(diào)劑，材料與化工085600 +3	13756423260 2026-03-22	3/150	2026-03-22 20:11 by edmund7
[考研] 085600材料與化工306 +4	z1z2z3879 2026-03-21	4/200	2026-03-21 23:44 by ms629
[考研] 一志愿南大，0703化學(xué)，分?jǐn)?shù)336，求調(diào)劑 +3	收到VS 2026-03-21	3/150	2026-03-21 18:42 by 學(xué)員8dgXkO
[考研] 求調(diào)劑 +3	13341 2026-03-20	3/150	2026-03-21 18:28 by 學(xué)員8dgXkO
[考研] 一志愿深大，0703化學(xué)，總分302，求調(diào)劑 +4	七月-七七 2026-03-21	4/200	2026-03-21 18:20 by 學(xué)員8dgXkO
[考研] 材料學(xué)碩333求調(diào)劑 +3	北道巷 2026-03-18	3/150	2026-03-21 18:17 by 學(xué)員8dgXkO
[考研] 求助 +5	夢(mèng)里的無言 2026-03-21	6/300	2026-03-21 17:51 by 學(xué)員8dgXkO
[考研] 307求調(diào)劑 +3	余意卿 2026-03-18	3/150	2026-03-21 17:31 by ColorlessPI
[考研] 求調(diào)劑 +3	Ma_xt 2026-03-17	3/150	2026-03-21 02:05 by JourneyLucky
[考研] 考研調(diào)劑求學(xué)校推薦 +3	伯樂29 2026-03-18	5/250	2026-03-20 22:59 by JourneyLucky
[考研] 317求調(diào)劑 +5	申子申申 2026-03-19	9/450	2026-03-20 22:26 by JourneyLucky
[考研] 295材料求調(diào)劑，一志愿武漢理工085601專碩 +5	Charlieyq 2026-03-19	5/250	2026-03-20 20:35 by JourneyLucky
[考研] 一志愿南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩 +5	@taotao 2026-03-20	5/250	2026-03-20 20:16 by JourneyLucky
[考研] 0856調(diào)劑，是學(xué)校就去 +8	sllhht 2026-03-19	9/450	2026-03-20 14:25 by 無懈可擊111
[考研] 085410人工智能專碩317求調(diào)劑（0854都可以） +4	xbxudjdn 2026-03-18	4/200	2026-03-20 09:07 by 不168
[考研] 材料考研調(diào)劑 +3	xwt。 2026-03-19	3/150	2026-03-19 11:22 by w沐陽w