




已閱讀5頁(yè),還剩2頁(yè)未讀, 繼續(xù)免費(fèi)閱讀
版權(quán)說(shuō)明:本文檔由用戶(hù)提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請(qǐng)進(jìn)行舉報(bào)或認(rèn)領(lǐng)
文檔簡(jiǎn)介
數(shù)值分析實(shí)驗(yàn)報(bào)告一、 實(shí)驗(yàn)?zāi)康?、 了解熟悉jacobi迭代法和seidel迭代法的解法2、 將原理與matlab語(yǔ)言結(jié)合起來(lái),編程解決問(wèn)題3、 分析實(shí)驗(yàn)結(jié)果二、實(shí)驗(yàn)題目設(shè)線(xiàn)性方程組為 考察用jacobi迭代法和seidel迭代法求解該線(xiàn)性方程組的收斂情況。如果收斂,給出誤差滿(mǎn)足的解。三、實(shí)驗(yàn)原理將A 作如下分解這里Jacobi迭代矩陣為Seidel迭代矩陣為它們的迭代格式都可化為則迭代格式對(duì)任何初值都瘦臉的充要條件是迭代矩陣譜半徑Jacobi迭代矩陣的特征方程為Seidel迭代矩陣的特征方程為四、實(shí)驗(yàn)內(nèi)容用matlab編寫(xiě)計(jì)算jacobi迭代矩程序,建立m文件如下:functionM=BJ(A)D=diag(diag(A); L=tril(-A)+D; U=triu(-A)+D; M=inv(D)*(L+U); 輸入: A=20 1 -3 -1;3 18 0 7;-1 2 40 -2;1 0 -1 5; M=BJ(A)M:M = 0 -0.0500 0.1500 0.0500 -0.1667 0 0 -0.3889 0.0250 -0.0500 0 0.0500 -0.2000 0 0.2000 0則jacobi迭代矩陣為:用matlab求jacobi迭代矩陣的特征根的算法如下: A=0 -0.05 0.15 -0.05;-0.67 0 0 -0.39;0.025 -0.05 0 0.05;-0.2 0 0.2 0;V,D=eig(A)V = -0.1892 + 0.0450i -0.1892 - 0.0450i -0.3812 -0.5005 -0.9467 -0.9467 0.8867 0.5461 -0.1528 - 0.1181i -0.1528 + 0.1181i -0.2099 -0.0466 -0.1056 + 0.1325i -0.1056 - 0.1325i 0.1561 0.6701 D = -0.1774 + 0.0864i 0 0 0 0 -0.1774 - 0.0864i 0 0 0 0 0.2194 0 0 0 0 0.1355 則最大特征根為:0.2194 則,所以jacobi迭代法收斂用matlab編程jacobi迭代法求根的算法:function n,x=jacobi(A,b,X,nm,w)%用雅克比迭代法求解方程組Ax=b%輸入:A為方程組的系數(shù)矩陣,b為方程組右端的列向量,X為迭代初值構(gòu)成的列向量,nm為最大迭代次數(shù),w為誤差精度%輸出:x為求得的方程組的解構(gòu)成的列向量,n為迭代次數(shù)n=1;m=length(A);D=diag(diag(A); %令A(yù)=D-L-U,計(jì)算矩陣DL=tril(-A)+D; %令A(yù)=D-L-U,計(jì)算矩陣LU=triu(-A)+D; %令A(yù)=D-L-U,計(jì)算矩陣UM=inv(D)*(L+U); %計(jì)算迭代矩陣g=inv(D)*b; %計(jì)算迭代格式中的常數(shù)項(xiàng)%下面是迭代過(guò)程while n=nm x=M*X+g; %用迭代格式進(jìn)行迭代if norm(x-X,2) A=20 1 -3 -1;3 18 0 7;-1 2 40 -2;1 0 -1 5;b=1;2;10;-1;X=0;0;0;0; nm=100; w=1e-4; n,x=jacobi1(A,b,X,nm,w)迭代次數(shù)為n = 6方程組的解為x = 0.0687 0.1645 0.2352 -0.1667n = 6x = 0.0687 0.1645 0.2352 -0.1667所以,滿(mǎn)足精度的根為x = 0.0687 0.1645 0.2352 -0.1667 迭代次數(shù)為6次用matlab編程計(jì)算seidel迭代矩陣算法為:functionM=BS(A)D=diag(diag(A); L=tril(-A)+D; U=triu(-A)+D; M=inv(D-L)*U; 輸入: A=20 1 -3 -1;3 18 0 7;-1 2 40 -2;1 0 -1 5; M=BS(A)M = 0 -0.0500 0.1500 0.0500 0 0.0083 -0.0250 -0.3972 0 -0.0017 0.0050 0.0711 0 0.0097 -0.0290 0.0042則seidel迭代矩陣為:用matlab編程求得seidel矩陣的算法為: A=0 -0.05 0.15 0.05;0 0.0083 -0.025 -0.3972;0 -0.0017 0.005 0.0711;0 0.0097 -0.029 0.0042;V,D=eig(A)V = 1.0000 -0.1596 + 0.6750i -0.1596 - 0.6750i -0.9104 0 0.6965 0.6965 0.3924 0 -0.1250 + 0.0028i -0.1250 - 0.0028i 0.1313 0 0.0070 - 0.1348i 0.0070 + 0.1348i 0.0000 D = 0 0 0 0 0 0.0088 + 0.0768i 0 0 0 0 0.0088 - 0.0768i 0 0 0 0 -0.0001 則特征根為 0.008+0.0768i 0.0088-0.0768i -0.0001則,所以seidel迭代法收斂用seidel迭代法求根的算法為:function n,x=gaussseidel(A,b,X,nm,w)%用高斯-賽德?tīng)柕ㄇ蠼夥匠探MAx=b%輸入:A為方程組的系數(shù)矩陣,b為方程組右端的列向量,X為迭代初值構(gòu)成的列向量,nm為最大迭代次數(shù),w為誤差精度%輸出:x為求得的方程組的解構(gòu)成的列向量,n為迭代次數(shù)n=1;m=length(A);I=eye(m); %生成m*m階的單位矩陣D=diag(diag(A); %令A(yù)=D-L-U,計(jì)算矩陣DL=tril(-A)+D; %令A(yù)=D-L-U,計(jì)算矩陣LU=triu(-A)+D; %令A(yù)=D-L-U,計(jì)算矩陣UM=inv(D-L)*U; %計(jì)算迭代矩陣g=inv(I-inv(D)*L)*(inv(D)*b); %計(jì)算迭代格式中的常數(shù)項(xiàng)%下面是迭代過(guò)程while n=nm x=M*X+g; %用迭代格式進(jìn)行迭代 if norm(x-X,2) A=20 1 -3 -1;3 18 0 7;-1 2 40 -2;1 0 -1 5;b=1;2;10;-1;X=0;0;0;0;nm=100;w=1e-4;n,x=gaussseidel(A,b,X,nm,w)迭代次數(shù)為n = 5方程組的解為x = 0.0687 0.1645 0.2352 -0.1667n = 5x = 0.0687 0.1645 0.2352 -0.166滿(mǎn)足精度的根為x =0.0687 0.1645 0.2352 -0.1667 迭代次數(shù)
溫馨提示
- 1. 本站所有資源如無(wú)特殊說(shuō)明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請(qǐng)下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請(qǐng)聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶(hù)所有。
- 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁(yè)內(nèi)容里面會(huì)有圖紙預(yù)覽,若沒(méi)有圖紙預(yù)覽就沒(méi)有圖紙。
- 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
- 5. 人人文庫(kù)網(wǎng)僅提供信息存儲(chǔ)空間,僅對(duì)用戶(hù)上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理,對(duì)用戶(hù)上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對(duì)任何下載內(nèi)容負(fù)責(zé)。
- 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請(qǐng)與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時(shí)也不承擔(dān)用戶(hù)因使用這些下載資源對(duì)自己和他人造成任何形式的傷害或損失。
最新文檔
- 計(jì)算機(jī)三級(jí)軟件測(cè)試技術(shù)基礎(chǔ)與公共政策試題及答案
- 軟考網(wǎng)絡(luò)工程師技術(shù)更新與應(yīng)對(duì)試題及答案
- 西方國(guó)家醫(yī)療制度的政治分析試題及答案
- 新經(jīng)濟(jì)形態(tài)下的政策創(chuàng)新需求試題及答案
- 深入探討西方國(guó)家的公共關(guān)系與政治制度試題及答案
- 如何設(shè)計(jì)有效的公共政策干預(yù)措施試題及答案
- 項(xiàng)目管理中的利益分析與規(guī)劃試題及答案
- 機(jī)電工程跨學(xué)科整合的探索試題及答案
- 可再生能源利用試題及答案
- 公共政策編制方法探討試題及答案
- 2025-2030中國(guó)軍用無(wú)人機(jī)行業(yè)市場(chǎng)現(xiàn)狀供需分析及投資評(píng)估規(guī)劃分析研究報(bào)告
- 2025年中國(guó)工裝治具市場(chǎng)調(diào)查研究報(bào)告
- 收購(gòu)爛尾樓可行性研究報(bào)告
- 2025年國(guó)際貿(mào)易實(shí)務(wù)課程考試試題及答案
- 爆炸事故賠償協(xié)議書(shū)
- 2025華陽(yáng)新材料科技集團(tuán)有限公司招聘(500人)筆試參考題庫(kù)附帶答案詳解
- 2025年下半年度蘇州城際鐵路限公司管理崗位公開(kāi)招聘易考易錯(cuò)模擬試題(共500題)試卷后附參考答案
- 聚合物化學(xué)基礎(chǔ)試題及答案
- GB/T 12008.7-2025塑料聚氨酯生產(chǎn)用聚醚多元醇第7部分:堿性物質(zhì)含量的測(cè)定
- 蔬菜凈菜車(chē)間管理制度
- 2025年中國(guó)石油化工行業(yè)市場(chǎng)發(fā)展前景及發(fā)展趨勢(shì)與投資戰(zhàn)略研究報(bào)告
評(píng)論
0/150
提交評(píng)論