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

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

2、轉(zhuǎn)step5 l=a/u; u=b-l*c； y=f-l*y；Step3：當(dāng)u時(shí)，x=y/u，否則轉(zhuǎn)step5Step4：對(duì)于：i=n-1,n-2,.,2,1，轉(zhuǎn)step6 x=(y-c*x)/uStep5：無(wú)解信息，轉(zhuǎn)step7Step6：輸出xStep7：關(guān)機(jī)2MATLAB程序function x,L,U=thomas(a,b,c,f)n=length(b);% 對(duì)A進(jìn)行分解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)+di

3、ag(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求解實(shí)例例1方程組例2方程組例3方程組四、實(shí)驗(yàn)結(jié)論對(duì)于追趕法我最先寫(xiě)的是如下的程序：5但是出現(xiàn)了如上截圖中的錯(cuò)誤，后來(lái)與同學(xué)討論還是沒(méi)能解決我的問(wèn)題，最后借鑒了她的算法得到了正確的結(jié)果。Thomas算法在課堂上老師就已經(jīng)給我們?cè)敿?xì)地講解并指導(dǎo)了我們?nèi)绾斡肕atlab編程，但是并沒(méi)有解決a矩陣的a1如何處理，對(duì)于這個(gè)問(wèn)題，我很快解決了。我最大的問(wèn)題就是