實驗2追趕法算法設(shè)計及MATLAB實現(xiàn)

1、數(shù)值計算方法實驗序號：實驗二 實驗名稱：追趕法算法設(shè)計及 MATLAB 實現(xiàn) 實 驗 人：專業(yè)年級：教 學(xué) 班： 學(xué) 號： 實驗時間：實驗二追趕法算法設(shè)計及MATLAB實現(xiàn)一、實驗?zāi)康?初步掌握算法設(shè)計規(guī)則；2 初步掌握MATLAB?序設(shè)計規(guī)則.二、實驗內(nèi)容1 構(gòu)造利用追趕法求解三對角線性方程組的算法；2在MATLAB環(huán)境下編寫追趕法的程序（函數(shù)）;3自由選擇若干個三對角線性方程組求解。三、實驗步驟1 追趕法算法：算法名稱：thomas輸入?yún)?shù)：向量a,b,c,f輸出參數(shù)：輸出解信息x算法的自然語言：Stepl：比力小訐力；Step2：對于 i=2,3,.n;Step2.1：當(dāng) 比廠。，否則轉(zhuǎn)

2、step5l i =ai /u i j;Ui=bi-li *Ci j ;yi=fi-li*y i4;Step3：當(dāng) un =0時，xn=yn/un，否則轉(zhuǎn) step5Step4：對于：i=n-1,n-2,.,2,1,轉(zhuǎn) step6Xi=(yi-Ci*Xi J/UiStep5:無解信息，轉(zhuǎn)step7Step6:輸出 xStep7：關(guān)機 2MATLAB 程序 function x,L,U=thomas(a,b,c,f) n=length(b);%對A進行分解 u(1)=b(1);for i=2:nif(u(i-1)=0) l(i-1)=a(i-1)/u(i-1); u(i)=b(i)-l(i-1)

3、*c(i-1);elsebreak;end end L=eye(n)+diag(l,-1); U=diag(u)+diag(c,1); x=zeros(n,1); y=x;% 求解 Ly=b y(1)=f(1); for i=2:ny(i)=f(i)-l(i-1)*y(i-1);end% 求解 Ux=y if(u(n)=0)x(n)=y(n)/u(n);end for i=n-1:-1:1x(i)=(y(i)-c(i)*x(i+1)/u(i);2end3#3.求解實例例1.方程組F 気 2J;b=X 3, 3, 37 ；Clj U U* ；Eh 0, L01 : x,LU=thomas (g

4、b, c, f )(3 1、X2 3 1|X202 3 1X31UThomas (ajbj c,f)n=length(b);滋對Aififr*解u(l)=b(l);for i=2:nif u(i-l)=0u(i)b (i)-l (i-l)*c(i-l);elsebreak;endendx=zeros (g 1)；尸；y(l)=b(l);for i=2：nendif u(n)=0x (n) =y (n) /u (nJ ;endfor i.=n-l: 1: 1x(i) = (y(i)-c(i)*x(i+1)/u(i);end s=, 2, IY ; b=3,3, 3.3; c= th 13 l?; f=l, 0, LO. xj Lj. U =thom.M (a., bj Cj )? One ar nare output argWLBnt5 not assigned during call to E:TDDOWWLOADorkVthcimas1m (thomas) 但是出現(xiàn)了如上截圖中的錯誤， 后來與同學(xué)討論還是沒能解決我的問題， 最后借 鑒了她的算法得到了正確的結(jié)果。 Thomas 算法在課堂上老師就已經(jīng)給我們詳細 地講解并指導(dǎo)了我們?nèi)绾斡?Matlab編程，但是并沒有解決a矩陣的al如何處理