【正文】
h/2+*h))*h/2。T(i+1,j+1)=4^j*T(i+1,j)/(4^j1)T(i,j)/(4^j1)。while abs(T(i+1,i+1)T(i,i))e i=i+1。h=h/2。 T(i+1,1)=T(i,1)/2+sum(feval(fname,a+h/2:h:bh/2+*h))*h/2。 for j=1:i T(i+1,j+1)=4^j*T(i+1,j)/(4^j1)T(i,j)/(4^j1)。 endend%Tt=T(i+1,j+1)。tocjacobi迭代法求方程組%P:范數(shù)的名稱,P=1,2,inf%error:近似值x的誤差%maxl迭代的最大次數(shù)function Jacobi(A,b,X0,P,error,maxl)[n n]=size(A)。X=zeros(n,1)。for k=1:maxl for j=1:n X(j)=(b(j)A(j,[1:j1,j+1:n])*X0([1:j1,j+1:n]))/A(j,j)。 end X errX=norm(XX0,P)。 X0=X。X1=A\b。 if(errXerror) disp(39。迭代數(shù)K,精確解X1和近似解X分別是:39。) K X1 X return endendif(errX=error) disp(39。請注意Jacobi迭代次數(shù)已經(jīng)超過最大迭代次數(shù)maxl 39。)endGS迭代法求方程組%P:范數(shù)的名稱,P=1,2,inf%error:近似值x的誤差%maxl迭代的最大次數(shù)function GS(A,b,X0,P,error,maxl)[n n]=size(A)。X=zeros(n,1)。for k=1:maxl for j=1:n。 XX=0。 for i=1:n if ij XX=XX+A(j,i)*X(i)。 end if ij XX=XX+A(j,i)*X0(i)。 end end X(j)=(b(j)XX)/A(j,j)。 end X errX=norm(XX0,p)。 X0=X。X1=A\b。 if(errXerror) disp(39。迭代次數(shù)k,精確解X1和近似解X分別是:39。) k X1 X return endendif(errX=error) disp(39。請注意Jacobi迭代次數(shù)已經(jīng)超過最大迭代次數(shù)maxl39。) endSOR迭代法求方程組 %P:范數(shù)的名稱,P=1,2,inf%error:近似值x的誤差%maxl迭代的最大次數(shù)%w:松弛因子function SOR(A,b,X0,P,error,maxl)[n]=size(A)。X=zeros(n,1)。for k=1:maxl for j=1:n。 XX=0。 for i=1:n if ij XX=XX+A(j,i)*X(i)。 end if ij XX=XX+A(j,i)*X0(i)。 end end X(j)=(1w)*X0(j)+w*(b(j)XX)/A(j,j)。 end errX=norm(XX0,P)。 X0=X。X1=A\b。 if(errXerror) disp(39。迭代次數(shù)k,精確解X1和近似解X分別是: 39。) k X1 X return endendif(errX=error)disp(39。請注意Jacobi迭代次數(shù)已經(jīng)超過最大迭代次數(shù)maxl 39。) end 參考文獻[1][M],中國科技大學(xué)出版社,2010[2]杜延松,[M],科學(xué)出版社,2006.[3]王沫然. [M],北京,電子工業(yè)出版社,2001,166191.[4]薛蓮,將金生. 數(shù)值計算方法[M],電子工業(yè)出版社,2007.[5]樓順天,于衛(wèi),閆華梁. Matlab程序設(shè)計語言[M],西安,西安電子科技大學(xué)出版社,1997.[6]薛定宇,[M],北京,清華大學(xué)出版社,2004.[7]李慶揚,王能超,[M],北京,清華大學(xué)出版社,1988.[8]肖筱南,趙來軍,黨林立. 現(xiàn)代數(shù)值計算方法[M], 北京,北京大學(xué)出版社2003.[9]封建湖,車剛明,[M]. 北京,科學(xué)出版社2001.[10]孫海東. EXCEL在數(shù)值計算方法方面的應(yīng)用[J],軟件世界,1999,137.致 謝畢業(yè)論文致謝詞模版本論文是在導(dǎo)老師的悉心指導(dǎo)下完成的。導(dǎo)師淵博的專業(yè)知識,嚴謹?shù)闹螌W(xué)態(tài)度,精益求精的工作作風(fēng),誨人不倦的高尚師德,嚴以律己、寬以待人的崇高風(fēng)范,樸實無華、平易近人的人格魅力對我影響深遠,他的循循善誘的教導(dǎo)和不拘一格的思路給予我無盡的啟迪。在撰寫整個畢業(yè)論文的過程當中,他為我們考慮到了每一個細節(jié),從開題報告到畢業(yè)論文的擬定修改上,鄭老師更是不厭其煩的為我們做好每一步的細心指導(dǎo)。對此,我表示衷心地感謝。沒有段老師,我的論文也不可能這么順利的完成。同時感謝這篇論文所涉及到的各位學(xué)者,本文引用了數(shù)位學(xué)者的研究文獻,如果沒有各位學(xué)者的研究成果帶給我的的幫助和啟發(fā),我將很難完成本篇論文的寫作。不僅使我樹了遠大的學(xué)術(shù)目標、掌握了基本的研究方法,還使我明白了許多待人接物與為人處世的道理。在此,謹向?qū)煴硎境绺叩木匆夂椭孕牡母兄x! 2013年5月于學(xué)院24