實驗2追趕法算法設計與MATLAB實現_第1頁
實驗2追趕法算法設計與MATLAB實現_第2頁
實驗2追趕法算法設計與MATLAB實現_第3頁
實驗2追趕法算法設計與MATLAB實現_第4頁
實驗2追趕法算法設計與MATLAB實現_第5頁
已閱讀5頁，還剩2頁未讀， 繼續免費閱讀

下載本文檔

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

文檔簡介

1、數 值 計 算 方 法實驗報告實驗序號：實驗二 實驗名稱：追趕法算法設計及MATLAB實現實 驗 人： 專業年級： 教 學 班：學 號：實驗時間：實驗二追趕法算法設計及MATLAB實現一、實驗目的初步掌握算法設計規則；初步掌握MATLAB程序設計規則二、實驗容1構造利用追趕法求解三對角線性方程組的算法；2在MATLAB環境下編寫追趕法的程序（函數）；3自由選擇若干個三對角線性方程組求解。 三、實驗步驟1追趕法算法：算法名稱：thomas輸入參數：向量a,b,c,f輸出參數：輸出解信息x算法的自然語言：Step1：u=b,y=b;Step2：對于 i=2,3,.n; Step2.1：當u，否則轉

2、step5 l=a/u; u=b-l*c； y=f-l*y；Step3：當u時，x=y/u，否則轉step5Step4：對于：i=n-1,n-2,.,2,1，轉step6 x=(y-c*x)/uStep5：無解信息，轉step7Step6：輸出xStep7：關機2MATLAB程序function x,L,U=thomas(a,b,c,f)n=length(b);% 對A進行分解u(1)=b(1);for i=2:n if(u(i-1)=0) l(i-1)=a(i-1)/u(i-1); u(i)=b(i)-l(i-1)*c(i-1); else break; endendL=eye(n)+dia

3、g(l,-1);U=diag(u)+diag(c,1);x=zeros(n,1);y=x;% 求解Ly=by(1)=f(1);for i=2:n y(i)=f(i)-l(i-1)*y(i-1);end% 求解Ux=yif(u(n)=0) x(n)=y(n)/u(n);endfor i=n-1:-1:1 x(i)=(y(i)-c(i)*x(i+1)/u(i);end3求解實例例1方程組例2方程組例3方程組四、實驗結論對于追趕法我最先寫的是如下的程序：但是出現了如上截圖中的錯誤，后來與同學討論還是沒能解決我的問題，最后借鑒了她的算法得到了正確的結果。Thomas算法在課堂上老師就已經給我們詳細地講解并指導了我們如何用Matlab編程，但是并沒有解決a矩陣的a1如何處理，對于這個問題，我很快解決了。我最大的問題就是如上所示，說明我的編程能力還是比較差，需要多練習。對于如上的錯誤希望老師看過之后能夠給予指導，！您好，歡迎您閱讀我的文章，本WORD文檔可編輯修改，也

溫馨提示

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

評論

0/150

提交評論