微分方程數(shù)值解法實驗報告-文庫吧資料
2024-08-03 17:34本頁面
【正文】 uu=zeros(n+1)。 end u_3_0=u_3。while flag u_3=B1*u_3_0+g。u_3_0=zeros(size(b,1),1)。g=w*(Dw*L)\b。R2=triu(A,1)。D=diag(D)。GaussSeidel迭代差分解圖像39。figure(3)mesh(x,y,uu)。 endend%GaussSeidel迭代差分解圖像x=[0:1/n:1]。%迭代次數(shù)值enduu=zeros(n+1)。 end u_2_0=u_2。while flag u_2=R*u_2_0+g。%初始迭代值u_2_0=u_2_0(:)。m_2=0。R=(DL)\U。L=tril(A,1)。endD=diag(A)。Jacobi迭代差分解圖像39。figure(2)mesh(x,y,uu)。 endend%Jacobi迭代差分解圖像x=[0:1/n:1]。%迭代次數(shù)值enduu=zeros(n+1)。 end u_1_0=u_1。while flag u_1=B*u_1_0+g。%初始迭代值u_1_0=u_1_0(:)。m_1=0。B=D\(L+R)。L=tril(A,1)。四．程序設計(MATLAB實現(xiàn))Jacobi迭代法function[u_1,m_1]=Jacobi_Solve(A,b,n,err)D=diag(A)。二．實驗目的： 分別使用五點差分法（Jacobi迭代，Gauss_Seidel迭代，SOR迭代），PR交替隱式差分法（單參數(shù)，雙參數(shù)），共軛梯度法，預共軛梯度法分別求橢圓方程的數(shù)值解。問題三：取步長h=k=1/64,1/128,作五點差分格式，用共軛梯度法和預處理共軛梯度法解差分方程，近似到小數(shù)點后四位。問題二：取步長h=k=1/64,1/128,作五點差分格式，用單參數(shù)和雙參數(shù)PR法解差分方程，近似到小數(shù)點后四位。微分方程數(shù)值解法實驗報告姓名： 班級：學號：一：問題描述求解邊值問題：其精確解為問題一：取步長h=k=1/64,1/128,作五點差分格式，用Jacobi迭代法，Gauss_Seidel迭代法，SOR 迭代法（w=）。求解差分方程，以前后兩次重合到小數(shù)點后四位的迭代值作為解的近似值，比較三種解法的迭代次數(shù)以及差分解與精確解的精度。與SOR法比較精度和迭代步數(shù)。與SOR法與PR法比較精度和迭代步數(shù)。三．實驗原理：(1) Jacobi迭代法設線性方程組 (1)的系數(shù)矩陣A可逆且主對角元素均不為零,令 并將A分解成 (2)從而(1)可寫成 令 其中. (3)以為迭代矩陣的迭
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1