中心差分解兩點邊值問題-wenkub
2023-04-08 06:07:40 本頁面
【正文】 …………（） 其中 …………（）為差分算子的截斷誤差，舍去，便得逼近邊值問題（）（）的差分方程： …………（）i=1,2,…,N1，由方程（）（），截斷誤差可表示為 …………（）當(dāng)網(wǎng)格均勻，即時差分方程（）簡化為 …………（）這相當(dāng)于用一階中心差商，二階中心差商依次代替（）的一階微商和二階微商的結(jié)果。我們可以用消元法或迭代法求解方程組（）（）式（）用方程組展開：寫成矩陣形式為：根據(jù)（）我們引進(jìn)誤差則誤差函數(shù)滿足下列差分方程：于是收斂性及收斂速度的估計問題，就歸結(jié)到通過右端（截斷誤差）估計誤差函數(shù)的問題。三．程序代碼：clcclfclfsyms x。 %這是p函數(shù)q=sin(x)+1+x。 %將區(qū)間劃分的等分,這里控制！h=(ba)/N。%確定A的偏離對角的上對角元diag_2=zeros(N2,1)。 %邊界條件for j=2:N diag_0(j1)=((subs(p,{x},{(X(j+1)+X(j))/2
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1