Re: [問題] 傅立葉的問題已回收

看板MATLAB作者 (So do I)時間16年前 (2009/09/26 09:23), 編輯推噓0(002)
留言2則, 1人參與, 最新討論串2/2 (看更多)
※ 引述《eeden (簡單愛)》之銘言: : 關於傅立葉級數f(x)=a0+a*cos(2*pi*t*f)+b*sin(2*pi*t*f) : 查了定義,傅立葉級數可表示任何一個有週期性的函數 : 現在小弟目前手上有的是f(x),t的數據以及f=0.08371 : t= 0 0.07 0.14 0.21.......................71.61 : f(x)=0.2134 0.2127 0.2116..............................0.2139 就假設你的f(x)是periodic of period T=1/f 你寫的Fourier series只有DC term + frequency component at fundamental frequency 你的f(x)用Fourier series表示, 可能會只有這幾項嗎?? 怎麼可能你用一個DC+一個sinusoid就想表示出你的periodic signal unless your periodic is so simple that it can be decomposed into a DC term and a single sinusoidal component. : 想利用傅立葉級數把這條曲線的實際方程式求出來(就是要把a0,a,b積分出來) : 可是怎麼積都積不出來,想請教到底該如何寫呢?已經想好幾天了,請高手幫幫忙,感謝 : 以下是我自己寫的積分程式碼 : T=71.61; : t=linspace(0,71.61,1024); : f=0.08371; : ft=0.2134; : fq0='0*x+0.2134'; : a0=1/T*quad(fq0,0,T); : fq1='0*x+0.2134*cos(2*pi*0*0.08371)'; : a1=2/T*quad(fq1,0,T); : fq2='0*x+0.2134*sin(2*pi*0*0.08371)'; : b1=2/T*quad(fq2,0,T); : fq3=a0+a1*cos(2*pi*f*t)+b1*sin(2*pi*f*t); : plot(t,y_new, '-m',t,fq3, '-.g'); : 附上圖的部份http://img8.imageshack.us/img8/3302/12345fk.jpg
: 粉紅色是原數據畫出來的曲線,綠色是積分出來的曲線, -- ※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 122.116.181.164

09/27 04:02, , 1F
忘了說明 數據的這條曲線是單讀取f=0.08371畫出來的曲線
09/27 04:02, 1F

09/27 04:03, , 2F
已自行解決了,程式碼寫錯的關係,感謝指導
09/27 04:03, 2F
文章代碼(AID): #1AlMqadZ (MATLAB)
討論串 (同標題文章)
文章代碼(AID): #1AlMqadZ (MATLAB)