→ Maddulin:整數倍不是問題關鍵,實際上也想分析的訊號也不可能滿足 04/02 11:27
→ Maddulin:重點還是sampling window, sampling rate 04/02 11:28
i網友真的講得很不錯
原則上只要記得 當DFT/FFT這種數位的轉換都是有前提的
最大的前提都是在取樣頻率 一定是該訊號頻率的整數倍
不滿足這些條件都會有side effect
事實上oppenhiem那本DSP就介紹的很完整了:)
※ 引述《invalid (everlasting)》之銘言：
: 基本上finite sampling的FFT總是會遇到spectrum leakage的問題
: 要記得，取"300"個點進行FFT分析
: 代表的是你"已經假設"你的訊號每300個samples就會重覆一次
: 除非你的X[0]=X[300]，不然一定會有不連續點
: 不連續點會使你的頻譜失真
: 而失真的程度跟你取點的多寡是相反的
: 也跟你的的time slot有關
: 一般的解決辦法除了多抓一些sample
: 就是用windowing function
: 有很多windowing function可以使用
: 每一種有不同的特性，端看你想要知道的東西是什麼
: 不過你的信號不是周期訊號
: 能量應該是連續分布在低頻的區域
: 應該也沒辦法"完全"的列出每個頻率的比例吧..
: ※ 引述《flybow (yyyyyyyy)》之銘言：
: : 請問各位電機通訊，還有程式的高手
: : 我有一個暫態訊號.x=exp(-0.5t)-exp(-t) t=linspace(0,20,300)
: : 我做了FFT(x,300)...我要怎要才能從這300個點內，完全的分出每個頻率佔多少的
: : 比例
: : Y = fft(x,300);%無論t怎樣改動，Y的值不變
: : Pyy= Y.* conj(Y) / 300;
: : f=1000*(0:150)/300; %???
: : figure(2);plot(f,Pyy(1:151))
: : title('Frequency content of y')
: : xlabel('frequency (Hz)')
: : 我用這段程式來看，只能大略看出來，卻無法詳細的得到確切數據...
: : 我想要的答案比如說...
: : 0.1HZ 佔訊號比例 75%
: : 20 HZ 0.5%
: : 的這種數據.....
: : 請各位幫幫我，我該怎麼求?
: : 我現在正在做聲波的散射的理論分析，過去都是使用簡單的單頻波，我這次使用的
: : 則是多頻率的暫態訊號..所以有這個問題在.....
: : 請問各位了...
: : 或者可以用FFT的這300個點，看出頻率的組成嗎?
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 203.69.97.52