數(shù)值分析上機題_第1頁
數(shù)值分析上機題_第2頁
數(shù)值分析上機題_第3頁
數(shù)值分析上機題_第4頁
數(shù)值分析上機題_第5頁
已閱讀5頁,還剩9頁未讀, 繼續(xù)免費閱讀

下載本文檔

版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請進行舉報或認領(lǐng)

文檔簡介

1、精選優(yōu)質(zhì)文檔-傾情為你奉上數(shù)值分析上機題目3實驗一1根據(jù)Matlab語言特點,描述Jacobi迭代法、Gauss-Seidel迭代法和SOR迭代法。2編寫Jacobi迭代法、Gauss-Seidel迭代法和SOR迭代法的M文件。3給定為五對角矩陣(1)選取不同的初始向量及右端面項向量b,給定迭代誤差要求,分別用編寫的Jacobi迭代法和Gauss-Seidel迭代法程序求解,觀察得到的序列是否收斂?若收斂,通過迭代次數(shù)分析計算結(jié)果并得出你的結(jié)論。(2)用編寫的SOR迭代法程序, 對于(1)所選取的初始向量及右端面項向量b進行求解,松馳系數(shù)取1<<2的不同值,在時停止迭代,通過迭代次

2、數(shù)分析計算結(jié)果并得出你的結(jié)論。實驗11、 根據(jù)MATLAB語言特點,描述Jacobi迭代法,Gauss-Seidel迭代法和SOR迭代法。2、 編寫Jacobi迭代法,Gauss-Seidel迭代法和SOR迭代法的M文件。Jacobi迭代法function x1,k=GS_2(A,b)n=length(A);D=diag(diag(A);L=-tril(A,-1);U=-triu(A,1);x1=zeros(n,1);x0=3*ones(n,1);k=0;while norm(x1-x0,1)>10(-7)&k<100 k=k+1; x0=x1; x1=D(L+U)*x0+

3、b);endk=kx=x1Gauss-Seidel迭代法function x1,k=GS_h(A,b)n=length(A);D=diag(diag(A);L=-tril(A,-1);U=-triu(A,1);x1=zeros(n,1);x0=3*ones(n,1);k=0;while norm(x1-x0,1)>10(-7)&k<100 k=k+1; x0=x1; x1=(D-L)U*x0-Db;endk=kx=x1SOR迭代法function x1,k=SOR_h(A,b)n=length(A);D=diag(diag(A);L=-tril(A,-1);U=-triu(

4、A,1);x1=zeros(n,1);x0=3*ones(n,1);k=0;w=0.96;while norm(x1-x0,1)>10(-7)&k<100 k=k+1; x0=x1; x1=(D-w*U)(1-w)*D+w*L)*x0+w*b);endk=kx=x13、采用Jacobi迭代法,Gauss-Seidel迭代法求解五對角矩陣clear,clcA=diag(3*ones(20,1)+diag(-0.5)*ones(19,1),-1)+diag(-0.5)*ones(19,1),1)+diag(-0.25)*ones(18,1),-2)+diag(-0.25)*on

5、es(18,1),2);b=sum(A')'x1,k1=Jacob_h(A,b)x2,k2=GS_h(A,b)運行結(jié)果:兩種方法都收斂,k1=27,k2=13。說明Gauss-Seidel迭代速度比Jacobi迭代速率快4、采用SOR迭代法程序?qū)ξ鍖蔷仃囘M行求解clear,clcA=diag(3*ones(20,1)+diag(-0.5)*ones(19,1),-1)+diag(-0.5)*ones(19,1),1)+diag(-0.25)*ones(18,1),-2)+diag(-0.25)*ones(18,1),2);b=sum(A')'x3,k3=SOR

6、_h(A,b)運行結(jié)果當w=0.5時,k3=53,當w=0.96時,k3=19,當w=1.1時,k3=14,當w=1.5時,k3=33。該結(jié)果說明,當w選取合適時,可以大大加快運算速率。實驗二題目: 多項式最小二乘法摘要:對于具體實驗時,通常不是先給出函數(shù)的解析式,再進行實驗,而是通過實驗的觀察和測量給出離散的一些點,再來求出具體的函數(shù)解析式。又因為測量誤差的存在,實際真實的解析式曲線并不一定通過測量給出的所有點。最小二乘法是求解這一問題的很好的方法,本實驗運用這一方法實現(xiàn)對給定數(shù)據(jù)的擬合。數(shù)學(xué)原理:對于給定的測量數(shù)據(jù)(xi,fi)(i=1,2,,n),設(shè)函數(shù)分布為特別的,取為多項式 (j=0

7、, 1,,m)則根據(jù)最小二乘法原理,可以構(gòu)造泛函令 (k=0, 1,,m)則可以得到法方程求該解方程組,則可以得到解,因此可得到數(shù)據(jù)的最小二乘解程序設(shè)計:編寫求解多項式擬合的Matlab函數(shù)子程序?qū)嶒炓螅河米钚《朔ㄌ幚硐旅娴膶嶒灁?shù)據(jù).xi3456789fi2.012.983.505.025.476.027.05并作出的近似分布圖。分別采用一次,二次、五次和偶數(shù)次多項式來擬合數(shù)據(jù)得到相應(yīng)的擬合多項式,并分別作出它們的曲線圖。實驗2clcclearn=input('請輸入多項式擬合次數(shù):');syms tfor i=1:n f(i)=t(i);endx=3 4 5 6 7 8

8、9'y=2.01 2.98 3.50 5.02 5.47 6.02 7.05'm=length(x);for i=1:m for j=1:nhs=inline(f(j),'t'); A(i,j)=hs(x(i); endendh=ones(m,1);A=h A;A1=A'*A;y1=A'*y;x1=A1y1;f=0;for i=1:n+1 f=f+x1(i)*t(i-1);endplot(x,y,'*');hold on x1=flipud(x1);x2=linspace(min(x),max(x);y2=polyval(x1,x

9、2);tt=poly2str(x1,'x')text(5,7,0, tt)plot(x2,y2)實驗三實驗名稱:非線性方程組數(shù)值求解的Newton類方法試驗。實驗?zāi)康模河肗ewton類方法求解線性方程組 F(x)=0,理解其解的復(fù)雜性、初始點選擇策略、減少算法工作量的方法等。實驗內(nèi)容與要求:分別用Newton法用Broyden秩1校正法求解下面非線性方程組 (1) 寫出MATLAB源代碼;(2) 給出迭代五次以上的結(jié)果;(3) 嘗試不同的初值,如可?。?;(4) 計算兩種方法的用時。實驗31、2、采用Newton法·¨clear,clcx0=0,0,0'

10、;y0=f(x0);yy0=df(x0);x1=x0-yy0y0;k=1;format longwhile norm(x1-x0,1)>10(-5) & k<100 k=k+1; x0=x1; y0=f(x0); yy0=df(x0); x1=x0-yy0y0;endx=x1k=k運行結(jié)果:x= 0.0332 -0.5047 -0.2383k =45Broyden 秩1法clear,clcx0=0,0,0'y0=f(x0);A0=df(x0);x1=x0-A0y0;y1=f(x1);k=1;format longwhile norm(x1-x0,1)>10(-2) & k< k=k+1; g=y1-y0; y=x1-x0; A1=A0+(g-A0*y)/(y'*y)*y' x0=x1; x1=x0-A1y1; A0=A1; y0=f(x0); y1=f(x1);endx=x1k=k運行結(jié)果,x = 0.3973 -0.3592 -0.1240k = 53.嘗試不同初值。x0=0.1,0.1,0.1采用Newton法·¨x = 0.2459 -0.6348

溫馨提示

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

評論

0/150

提交評論