中心差分解兩點邊值問題-文庫吧資料
2025-03-30 06:07本頁面
【正文】 三條對角for j=2:N1 diag_1(j1)=((subs(p,{x},{(X(j+1)+X(j))/2})))/(h^2)+subs(r,{x},{X(j)})/(2*h)。 %邊界條件for j=2:N diag_0(j1)=((subs(p,{x},{(X(j+1)+X(j))/2}))+(subs(p,{x},{(X(j1)+X(j))/2})))/(h^2)+(subs(q,{x},{X(j)}))。u_a=0。%確定A的偏離對角的上對角元diag_2=zeros(N2,1)。%這是fdiag_0=zeros(N1,1)。 %將區(qū)間劃分的等分,這里控制！h=(ba)/N。%這是f函數(shù)r=0。 %這是p函數(shù)q=sin(x)+1+x。 %區(qū)間界點b=2。三．程序代碼：clcclfclfsyms x。若引進記號，設則可將（）改寫為將差分解表成 …………（）其中滿足 …………（）而滿足 …………（）先估計，由 …………（）據(jù)差分格林公式再利用柯西不等式，有常數(shù)使 …………（）將不等式（）用于（）右端，則 …………（）解差分方程（，易得）從而這樣， …
環(huán)評公示相關推薦
鄂ICP備17016276號-1