 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 6512 | 回復(fù): 12 | ||||
| 本帖產(chǎn)生 1 個 程序強帖 ，點擊這里進行查看 | ||||
[求助]
matlab 傅里葉變換 頻譜分析
|
||||
|
我想利用matlab進行編程，進行傅里葉頻譜分析，但是需要將頻譜轉(zhuǎn)換成功 率譜。詳情見附件 |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 matlab | 頻譜 |
» 猜你喜歡
被我言中：新模板不強調(diào)格式了，假專家開始管格式了
已經(jīng)有3人回復(fù)
-
085600材料與化工求調(diào)劑
已經(jīng)有5人回復(fù)
-
材料與化工求調(diào)劑
已經(jīng)有5人回復(fù)
-
277調(diào)劑
已經(jīng)有4人回復(fù)
-
290求調(diào)劑
已經(jīng)有4人回復(fù)
-
材料工程專碩調(diào)劑
已經(jīng)有3人回復(fù)
-
26考研求調(diào)劑
已經(jīng)有6人回復(fù)
-
求調(diào)劑，總分315，考的生物醫(yī)藥，一志愿湖南師范大學(xué)。調(diào)劑到任何專業(yè)都可以
已經(jīng)有5人回復(fù)
-
311求調(diào)劑
已經(jīng)有7人回復(fù)
-
303求調(diào)劑
已經(jīng)有3人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 傅里葉變換，Laplace變換，Z變換
已經(jīng)有109人回復(fù)
- 用origin做快速傅里葉變換的方法
已經(jīng)有3人回復(fù)
- 【轉(zhuǎn)載】[MATLAB仿真在信號處理中的應(yīng)用].徐明遠,劉增力.高清文字版
已經(jīng)有94人回復(fù)
- 求窗口傅里葉變換的matlab程序
已經(jīng)有6人回復(fù)
- Matlab頻域分析 FFT問題求助
已經(jīng)有6人回復(fù)
- 關(guān)于matlab進行傅里葉變換和逆變換的程序，求大神指導(dǎo)！
已經(jīng)有9人回復(fù)
- 【求組】用matlab的fft求空間頻譜
已經(jīng)有8人回復(fù)
- 求教matlab傅里葉變換求解
已經(jīng)有5人回復(fù)
- 怎么做傅里葉變換
已經(jīng)有10人回復(fù)
- Radon變換的MATLAB實現(xiàn)
已經(jīng)有9人回復(fù)
- 請教傅里葉變化應(yīng)用于頻譜分析中時，相位譜中的范圍怎么計算出來？
已經(jīng)有6人回復(fù)
- 一個關(guān)于傅里葉變換的問題。
已經(jīng)有6人回復(fù)
- 傅里葉逆變換 頻域到時域
已經(jīng)有3人回復(fù)
- 關(guān)于貝葉斯線性判別分析的matlab程序
已經(jīng)有4人回復(fù)
- 求助：能不能幫忙寫一個用matlab求一個函數(shù)的傅里葉變換的程序
已經(jīng)有16人回復(fù)
- 利用MATLAB 分析圖片顏色之間的差別
已經(jīng)有4人回復(fù)
- 傅里葉變換紅外光譜法操作的常見問題及注意事項
已經(jīng)有39人回復(fù)
- 如何用matlab進行快速傅里葉變換？（請教）
已經(jīng)有5人回復(fù)
- 【求助】菜鳥求助 matlab 繪圖問題 繪制光譜圖
已經(jīng)有3人回復(fù)
sudo
木蟲 (正式寫手)
- 程序強帖: 16
- 應(yīng)助: 6 (幼兒園)
- 金幣: 1297.6
- 散金: 1486
- 紅花: 20
- 帖子: 588
- 在線: 641小時
- 蟲號: 1211394
- 注冊: 2011-02-24
- 性別: GG
- 專業(yè): 文藝美學(xué)
★
dubo(金幣+1): 歡迎常來程序語言版討論 2011-07-31 14:25:14
dubo(金幣+1): 歡迎常來程序語言版討論 2011-07-31 14:25:14
唉，可能解釋起來一言難盡~建議樓主閱讀一下信號處理的入門書嗯 首先，你師兄的圖，應(yīng)該是降采樣過后的（或者是后來截取了一段，不過我傾向于認為是降采樣），也就是說他用來做功率譜的數(shù)據(jù)采樣率只有400Hz，你給的數(shù)據(jù)的采樣率是2000Hz 其次，你師兄給的圖是沒有取dB作為縱坐標單位的，也就是沒有計算10*log10(P)，而是直接用“功率密度”來做縱坐標~ 最后，其實MATLAB有現(xiàn)成函數(shù)可以做功率譜： 其中輸出P和F分別是功率譜和對應(yīng)的頻率 輸入x是信號，window是數(shù)據(jù)窗口（可以寫[]用默認矩形窗），nfft是進行fft變換的點數(shù)（可以寫[]取默認值，一般會取到最近的2的整數(shù)次冪左右的地方），F(xiàn)s是采樣率，單位是Hz 然后我對你的數(shù)據(jù)裸進行功率譜處理的話，結(jié)果如下圖： |
2樓2011-07-25 20:28:50
sudo
木蟲 (正式寫手)
- 程序強帖: 16
- 應(yīng)助: 6 (幼兒園)
- 金幣: 1297.6
- 散金: 1486
- 紅花: 20
- 帖子: 588
- 在線: 641小時
- 蟲號: 1211394
- 注冊: 2011-02-24
- 性別: GG
- 專業(yè): 文藝美學(xué)
【答案】應(yīng)助回帖
★ ★ ★ ★ ★ ★ ★
xzhdty(金幣+2): 歡迎常來 2011-07-25 20:57:22
dubo(金幣+5, 程序強帖+1): 感謝應(yīng)助 2011-07-31 14:26:41
xzhdty(金幣+2): 歡迎常來 2011-07-25 20:57:22
dubo(金幣+5, 程序強帖+1): 感謝應(yīng)助 2011-07-31 14:26:41
|
另外說明：MATLAB的功率譜實現(xiàn)函數(shù)是像這樣子的： 為了使用FFT算法（速度快），MATLAB把數(shù)據(jù)擴展了，本來降采樣你的數(shù)據(jù)之后，一共有4000個點，但是實際上使用了4096個點。但是由于數(shù)據(jù)是實數(shù)序列，所以，僅有采樣率一半的數(shù)據(jù)是有效的（后半部分其實是鏡像），所以在圖上，F(xiàn)和P都是2048個點（periodogram取了2049個，這個無所謂） |
4樓2011-07-25 20:53:06
sudo
木蟲 (正式寫手)
- 程序強帖: 16
- 應(yīng)助: 6 (幼兒園)
- 金幣: 1297.6
- 散金: 1486
- 紅花: 20
- 帖子: 588
- 在線: 641小時
- 蟲號: 1211394
- 注冊: 2011-02-24
- 性別: GG
- 專業(yè): 文藝美學(xué)
【答案】應(yīng)助回帖
whfire(金幣+20): 哥們，謝謝。太謝謝了 2011-07-27 09:13:49
|
對了，還有一些想補充的~ 首先就是直流分量的問題，如果不先對數(shù)據(jù)都減去一個均值，那么接近0Hz部分的功率會非常非常大~這是有害的，一般做數(shù)據(jù)處理的時候如果只關(guān)注交流量的話，就要先使得整個序列的均值為0，方法就是給每個數(shù)減去該序列的一個均值 然后就是降采樣的方法，其實很簡單啦，在MATLAB里面的話，如果y是你原來的數(shù)據(jù)的話（注意是測量數(shù)據(jù)，你給的數(shù)據(jù)的第2列）.....還是給出命令的全部吧： 結(jié)果就是 |
» 本帖已獲得的紅花（最新10朵）
獨孤_逍遙3樓2011-07-25 20:44:28
mogul
金蟲 (正式寫手)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 965.5
- 散金: 51
- 帖子: 310
- 在線: 150.7小時
- 蟲號: 1350431
- 注冊: 2011-07-19
- 性別: GG
- 專業(yè): 金屬材料的凝固與結(jié)晶學(xué)
|
[P,F]=periodogram(x, window, nfft, Fs) [Pxx,w] = periodogram(x,window,nfft) uses the modified periodogram to estimate the PSD while specifying the length of the FFT with the integer nfft. "If you set nfft to the empty vector [], it takes the default value for N listed in the previous syntax. "(ref, matlab help) "If you set nfft to the empty vector [], nfft defaults to min(256,length(x)" (ref, 薛年喜. MATLAB在數(shù)字信號處理中的應(yīng)用(第2版). p286) |
6樓2013-03-13 22:04:06
sudo
木蟲 (正式寫手)
- 程序強帖: 16
- 應(yīng)助: 6 (幼兒園)
- 金幣: 1297.6
- 散金: 1486
- 紅花: 20
- 帖子: 588
- 在線: 641小時
- 蟲號: 1211394
- 注冊: 2011-02-24
- 性別: GG
- 專業(yè): 文藝美學(xué)
，不過其實第一列并不需要全部給出來 5樓2011-07-25 20:56:12
夏天summer
銅蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 207.5
- 散金: 48
- 帖子: 268
- 在線: 29.9小時
- 蟲號: 1982190
- 注冊: 2012-09-07
- 性別: GG
- 專業(yè): 控制理論與方法
7樓2013-07-10 10:33:17
8樓2014-02-28 16:18:25
9樓2014-04-03 23:53:13
10樓2014-04-04 08:12:47
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 whfire 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 290求調(diào)劑 +3 | p asserby. 2026-03-15 | 4/200 |
|
|---|---|---|---|---|
|
[考博] 26申博 +4 | 八旬速覽 2026-03-16 | 4/200 |
|
|
|
[考研] 267一志愿南京工業(yè)大學(xué)0817化工求調(diào)劑 +6 | SUICHILD 2026-03-12 | 6/300 |
|
|
|
[考研] 301求調(diào)劑 +3 | A_JiXing 2026-03-16 | 3/150 |
|
|
|
[考研] 070300化學(xué)學(xué)碩求調(diào)劑 +6 | 太想進步了0608 2026-03-16 | 6/300 |
|
|
|
[考研] 一志愿華中師范071000，325求調(diào)劑 +6 | RuitingC 2026-03-12 | 6/300 |
|
|
|
[考研] 0856專碩279求調(diào)劑 +5 | 加油加油！? 2026-03-15 | 5/250 |
|
|
|
[考研] 070305求調(diào)劑 +3 | mlpqaz03 2026-03-14 | 4/200 |
|
|
|
[考研] 255求調(diào)劑 +3 | 李嘉慧， 2026-03-12 | 4/200 |
|
|
|
[考研] 328求調(diào)劑 +3 | 5201314Lsy！ 2026-03-13 | 6/300 |
|
|
|
[考研] 266求調(diào)劑 +4 | 學(xué)員97LZgn 2026-03-13 | 4/200 |
|
|
|
[考研] 復(fù)試調(diào)劑 +9 | Copy267 2026-03-10 | 9/450 |
|
|
|
[考研] 一志愿中科院，化學(xué)方向，295求調(diào)劑 +4 | 一氧二氮 2026-03-11 | 4/200 |
|
|
|
[考研] 290求調(diào)劑 +9 | ADT 2026-03-11 | 9/450 |
|
|
|
[考研] 0703化學(xué)一志愿211 總分320求調(diào)劑 +5 | 瑪卡巴卡啊哈 2026-03-11 | 5/250 |
|
|
|
[考研] 329求調(diào)劑 +3 | miaodesi 2026-03-12 | 4/200 |
|
|
|
[考研] 281求調(diào)劑 +9 | Koxui 2026-03-12 | 11/550 |
|
|
|
[考研] 26調(diào)劑/材料科學(xué)與工程/總分295/求收留 +9 | 2026調(diào)劑俠 2026-03-12 | 9/450 |
|
|
|
[考研] 0856化學(xué)工程280分求調(diào)劑 +4 | shenzxsn 2026-03-11 | 4/200 |
|
|
|
[考研] 070303一志愿西北大學(xué)學(xué)碩310找調(diào)劑 +3 | d如愿上岸 2026-03-13 | 3/150 |
|