推 Maddulin:他要問的是 exact frequency component 04/01 17:04
※ 引述《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個點，看出頻率的組成嗎?
你的目的是要看出組成比例
我猜應該是指能量的比例吧
那麼就這樣做
Ptot = sum(Pyy);
for n = 1:150
P_ratio(n) = Pyy(n)/Ptot;
end
plot(f(1:150),P_ratio));
--
可能有語法錯誤
手邊沒有MATLAB Orz
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 220.135.83.97