[問題] (急)我畫的頻譜有什麼問題嗎???已回收

看板MATLAB作者 (厚宿共)時間17年前 (2008/12/02 01:56), 編輯推噓0(003)
留言3則, 2人參與, 最新討論串1/1
N=8000; fs=16000; dt=1/fs; t=eps:dt:(N-1)*dt; fc=2000; ERB=0.108*fc+24.7; g=t.^3.*exp(-2*pi*1.019*ERB*t).*cos(2*pi*fc*t); G=fft(g); Gn=abs(G)/max(abs(G)); Gndb=20*log10(Gn); f=fs*(0:N/2)/N; plot(f,Gndb(1:N/2+1)) ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ 如題,我畫的是函式gammatone的頻譜!!演算法為g paper上寫頻譜對稱,但我畫出來又不是那一回事!! 有人知道為什麼嗎??? 謝謝^^ -- ※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 220.134.147.167

12/02 03:17, , 1F
頻譜是對稱的...plot(Gndb)
12/02 03:17, 1F

12/02 11:01, , 2F
不是這個意思...我的意思是半邊頻譜對稱QQ
12/02 11:01, 2F

12/02 13:44, , 3F
以2000為縱軸對稱!!看起來是...但實際上不是!!
12/02 13:44, 3F
文章代碼(AID): #19D2LUms (MATLAB)
文章代碼(AID): #19D2LUms (MATLAB)