中心差分解兩點邊值問題(完整版)
【正文】 得越細，即步長越小，數(shù)值解與精確解的誤差就越小，數(shù)值解越逼近精確解。,X(2:N),precise_value,39。%39。A=diag(diag_0)+diag(diag_1,1)+diag(diag_2,1)。 %邊界條件u_b=2。 %這是r函數(shù).N=10。所以差分方程的解的收斂速度為。則由節(jié)點又構(gòu)成的一個網(wǎng)格剖分，稱為對偶剖分。記，稱為網(wǎng)格最大步長。我們可以用消元法或迭代法求解方程組（）（）式（）用方程組展開：寫成矩陣形式為：根據(jù)（）我們引進誤差則誤差函數(shù)滿足下列差分方程：于是收斂性及收斂速度的估計問題，就歸結(jié)到通過右端（截斷誤差）估計誤差函數(shù)的問題。 %這是p函數(shù)q=sin(x)+1+x。%確定A的偏離對角的上對角元diag_2=zeros(N2,1)。 value_of_f(j1)=subs(f,{x},{X(j)})。)。最大的誤差是%f\n39。b*39。precise_value=dx.*(dx1) %精確解%fprintf(39。value_of_f(N1)=value_of_f(N1)+u_b*(subs(p,{x},{(X(N)+X(N+1))/2}))/(h^2)。u_a=0。%這是f函數(shù)r=0。若引進記號，設(shè)則可將（）改寫為將差分解表成
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1