數值計算方法實驗報告-矩陣直接LU分解法方程組和微分方程的數值解法上機實現_第1頁
數值計算方法實驗報告-矩陣直接LU分解法方程組和微分方程的數值解法上機實現_第2頁
數值計算方法實驗報告-矩陣直接LU分解法方程組和微分方程的數值解法上機實現_第3頁
全文預覽已結束

下載本文檔

版權說明：本文檔由用戶提供并上傳，收益歸屬內容提供方，若內容存在侵權，請進行舉報或認領

文檔簡介

數值計算方法實驗報告實驗報告題目：矩陣直接LU分解法方程組和微分方程的數值解法上機實現實驗要求矩陣直接LU分解法方程組和微分方程的數值解法實驗分析（包括數學原理，小組分析討論后確定實驗方案和實現思路）用matlab對矩陣三角分解法求解方程組實驗步驟（過程）（包括程序及上機的實現的結果）functionX=LU1(a,b)b=b';a=[a',b]',n=length(b');x=zeros(n,l);y=zeros(n,l);U=zeros(n);L=eye(n);fork=l:nU(l,k)=a(l,k);L(k,l)=a(k,l)/U(l,l);endfori=2:nfork=i:nlu=0;lul=0;forj=l:i-llu=lu+L(i,j)*U(j,k);lul=lul+L(k,j)*U(j.i);endU(i,k)=a(i,k)-lu;L(k,i)=(a(k,i)-lul)/U(i,i);endendfori=l:nly=0;forj=l:ily=ly+L(i,j)*y(j);endy(i)=B(i)-ly;endfori=n:-l:llyl=0;forj=i+l:nlyl=lyl+U(i,j)*X(j);endX(i)=(y(i)-lyl)/U(i,i);end輸入實驗內容：10-70x17-326x2=45-15x36輸入：A=[10,-7,0;-3,2,6;5,-1,5];B=[7,4,6];X=LU1(A,B)A=10-707-32645-156L=1.000000-0.30001.00000.5000-25.00001.0000U=10.0000-7.000000-0.10006.000000155.0000X=-0.0000-1.00001.0000總結（包括實驗過程遇到的情況等

溫馨提示

  • 1. 本站所有資源如無特殊說明，都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
  • 2. 本站的文檔不包含任何第三方提供的附件圖紙等，如果需要附件，請聯系上傳者。文件的所有權益歸上傳用戶所有。
  • 3. 本站RAR壓縮包中若帶圖紙，網頁內容里面會有圖紙預覽，若沒有圖紙預覽就沒有圖紙。
  • 4. 未經權益所有人同意不得將文件中的內容挪作商業或盈利用途。
  • 5. 人人文庫網僅提供信息存儲空間，僅對用戶上傳內容的表現方式做保護處理，對用戶上傳分享的文檔內容本身不做任何修改或編輯，并不能對任何下載內容負責。
  • 6. 下載文件中如有侵權或不適當內容，請與我們聯系，我們立即糾正。
  • 7. 本站不保證下載資源的準確性、安全性和完整性, 同時也不承擔用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。

評論

0/150

提交評論