※ 引述《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個點，看出頻率的組成嗎? clear t=linspace(0,20,300); xx=(t(2)-t(1))*10; L=(t(2)-t(1))*15; move=L-xx; delt=t(2)-t(1); f=0.000001:0.1:50; %由圖得知，輸入訊號為低頻，50Hz以後我就不考慮了，只取 0~50Hz的結果，也考慮各頻率所佔訊號能量的比重 %由0.000001開始是因為，bessely(n,a)在a=0時會發散 m1=1./((1+2*i*2*pi*f).*(1+i*2*pi*f)); %exp(-0.5t)-exp(-t)之FT m=abs(m1); %各頻率的訊號值大小，表示振幅響應，亦即傅立葉振幅 total=sum(m); %大小總和 weight=m./total; %比重 ============================================================================ 這是我想到的，不過還蠻爛的，也不精準.....如同板友所說， 訊號確實是連續分布 --

※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 140.123.66.51