[問題] 函數多次數值微分後加總已回收

看板MATLAB作者caseypie (吟遊詩人)時間15年前 (2010/07/07 09:31)推噓4(4推 0噓 5→)

留言9則, 3人參與討論串1/1

我有一個初始函數f_0(x)，有解析表達式還有一個遞迴關係： f_{n+1}(x) = x*f_n(x) + d(f_n(x))/dx 我想算 f_0 + f_1 + ...... f_100 剛剛試了直接不定數微分硬幹，結果記憶體不足...... 所以還是得用數值解法所以我先做一個x的一維陣列，然後再求得f_0(x)的一維陣列然後再用diff(f_0(x))./diff(x)做微分但是問題來了，微分後的陣列會變小，也就是從兩邊慢慢被吃掉所以那個遞迴關係因為維度不同加不起來...... 是不是只能用手動內插法把原本的f_0(x)縮小維度呢？但是我有點擔心這樣做100次以後誤差會太大...... -- http://caseypie-bard.blogspot.com/ -- ※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 59.113.76.157

→

vaakaa

07/07 10:43, , 1^F

07/07 10:43, 1^F

→

caseypie

07/07 11:29, , 2^F

07/07 11:29, 2^F

→

caseypie

07/07 11:29, , 3^F

07/07 11:29, 3^F