




版權說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權,請進行舉報或認領
文檔簡介
解線性方程組的迭代法演示文稿當前第1頁\共有33頁\編于星期五\13點優(yōu)選解線性方程組的迭代法當前第2頁\共有33頁\編于星期五\13點內(nèi)容提要引言簡單迭代法賽得爾迭代法迭代解法的收斂性
MATLAB的線性方程組求解函數(shù)2小結(jié)當前第3頁\共有33頁\編于星期五\13點
根據(jù)給定方程組,設計出一個迭代公式,構(gòu)造一數(shù)組的序列,代入迭代公式,計算出,再代入迭代公式,經(jīng)過k次迭代運算后得到,若收斂于某一極限數(shù)組xi,則xi就是方程組的近似解。迭代過程本質(zhì)上就是計算極限的過程,一般不能得到精確解。迭代法的優(yōu)點是程序簡單,適合于大型方程組求解,但缺點是要判斷迭代是否收斂和收斂速度問題。雅可比(Jacobi(1804-1851))迭代法(簡單迭代法)賽得爾(Seidel(1821-1896))迭代法迭代解法的基本思想1、引言當前第4頁\共有33頁\編于星期五\13點設線性代數(shù)方程組為2、簡單迭代法展開為當前第5頁\共有33頁\編于星期五\13點若對角元素逐一變量分離得方程組當前第6頁\共有33頁\編于星期五\13點即此即為迭代公式簡單迭代解法的過程如下:1
設定一組初值第i個變量第k次迭代2
第一次迭代:得到當前第7頁\共有33頁\編于星期五\13點3
第二次迭代:得到4
同樣做法,得到第k+1次迭代:當前第8頁\共有33頁\編于星期五\13點迭代次數(shù)k的取值與精度要求有關,按下式判斷:若滿足則停止迭代為了便于編程,迭代公式可改寫為:當前第9頁\共有33頁\編于星期五\13點function[x,iter,exitflag]=Jacobi_iter(A,b,x0,eps,iter_max)%線性方程組的Jacobi迭代求解(向量形式)%輸入?yún)?shù):%---A:線性方程組的系數(shù)矩陣%---b:線性方程組的右端項%---x0:初始向量,默認值為零向量%---eps:精度要求,默認值為1e-6%---iter_max:最大迭代次數(shù),默認值為100%輸出參數(shù):%---x:線性方程組的近似解%---iter:迭代次數(shù)%---exitflag:迭代成功與否的標志:exitflag=1表示迭代成功%exitflag=0表示迭代失敗n=length(b);ifnargin<5|isempty(iter_max);iter_max=100;endifnargin<4|isempty(eps);eps=1e-6;endifnargin<3|isempty(x0);x0=zeros(n,1);enditer=0;exitflag=1;D=diag(diag(A));L=tril(A,-1);U=triu(A,1);J=-inv(D)*(L+U);f=inv(D)*b;whileiter<iter_maxx=J*x0+f;ifnorm(x-x0,inf)<epsbreakendx0=x;iter=iter+1;endifiter==iter_maxexitflag=0;end當前第10頁\共有33頁\編于星期五\13點function[x,iter,exitflag]=Jacobi_iteration(A,b,x0,eps,iter_max)%線性方程組的Jacobi迭代求解(分量形式)%輸入?yún)?shù):%---A:線性方程組的系數(shù)矩陣%---b:線性方程組的右端項%---x0:初始向量,默認值為零向量%---eps:精度要求,默認值為1e-6%---iter_max:最大迭代次數(shù),默認值為100%輸出參數(shù):%---x:線性方程組的近似解%---iter:迭代次數(shù)%---exitflag:迭代成功與否的標志:exitflag=1表示迭代成功%exitflag=0表示迭代失敗n=length(b);ifnargin<5;iter_max=100;endifnargin<4;eps=1e-6;endifnargin<3;x0=zeros(n,1);endx=zeros(n,1);iter=0;exitflag=1;whileiter<iter_maxfori=1:nx(i)=(b(i)-A(i,[1:i-1,i+1:n])*x0([1:i-1,i+1:n]))/A(i,i);endifnorm(x-x0,inf)<epsbreakendx0=x;iter=iter+1;endifiter==iter_maxexitflag=0;end當前第11頁\共有33頁\編于星期五\13點MATLAB程序設計function[x,n]=richason(A,b,x0,eps,M)if(nargin==3)eps=1.0e-6;M=200;elseif(nargin==4)M=200;EndI=eye(size(A));x1=x0;x=(I-A)*x0+b;n=1;。。。。。。while(norm(x-x1)>eps)x1=x;x=(I-A)*x1+b;n=n+1;if(n>=M)disp('Warning:迭代次數(shù)太多,現(xiàn)在退出!');return;endend當前第12頁\共有33頁\編于星期五\13點例:求解方程組clearall;A=[1.0170-0.00920.0095;-0.00920.99030.0136;0.00950.01360.9898];b=[101]';x0=[000]';[x,n]=richason(A,b,x0)x=0.9739-0.00471.0010n=5當前第13頁\共有33頁\編于星期五\13點
賽得爾迭代法與簡單迭代法類似,只是迭代公式有所改進。3、賽得爾迭代法簡單迭代法賽得爾迭代法當前第14頁\共有33頁\編于星期五\13點MATLAB程序設計function[x,n]=gauseidel(A,b,x0,eps,M)ifnargin==3eps=1.0e-6;M=200;elseifnargin==4M=200;elseifnargin<3errorreturn;endD=diag(diag(A));%求A的對角矩陣L=-tril(A,-1);%求A的下三角陣U=-triu(A,1);%求A的上三角陣。。。。。。G=(D-L)\U;f=(D-L)\b;x=G*x0+f;n=1;%迭代次數(shù)whilenorm(x-x0)>=epsx0=x;x=G*x0+f;n=n+1;if(n>=M)disp('Warning:迭代次數(shù)太多,可能不收斂!');return;endend當前第15頁\共有33頁\編于星期五\13點例:線性代數(shù)方程組的迭代解法----賽得爾迭代法clearall;A=[953381;533813017;381301725317];b=[764893547]';x0=zeros(3,1);[x,n]=gauseidel(A,b,x0,1e-4,10)Warning:迭代次數(shù)太多,可能不收斂!x=-0.80373.3330-0.2450n=200當前第16頁\共有33頁\編于星期五\13點
迭代解法的前提條件是迭代解出的近似解序列必須具有收斂性。如果近似解序列是發(fā)散的,迭代法則不能獲得解。4、迭代解法的收斂性當前第17頁\共有33頁\編于星期五\13點以下列初值進行簡單迭代kX1X2X30000111-14-32-6981663-499-374-42944851-7149-2124當前第18頁\共有33頁\編于星期五\13點迭代收斂條件:嚴格對角占優(yōu)矩陣若不滿足收斂條件,適當調(diào)整方程次序或作一定的線性組合,就可能滿足收斂條件。當前第19頁\共有33頁\編于星期五\13點格式solve('eqn1','eqn2',...,'eqnN','var1,var2,...,varN')5、MATLAB的線性方程組求解函數(shù)2當前第20頁\共有33頁\編于星期五\13點當前第21頁\共有33頁\編于星期五\13點格式X=fsolve(FUN,X0)Matlab非線性方程組求解說明:求解方程形式F(X)=0X、F可以是向量或矩陣
X0初值當前第22頁\共有33頁\編于星期五\13點當前第23頁\共有33頁\編于星期五\13點實例:基于Matlab的透鏡中心偏測量光軸擬合光學中心偏測量儀作為精確測定和嚴格校正光學系統(tǒng)中心偏誤差的儀器,它可以指出透鏡組中的各鏡面相對于光軸的中心偏移數(shù)值大小和方向。它的測量結(jié)果具有兩個方面的意義:其一是通過根據(jù)被測光學件各面的中心誤差是否超出,來判定光學件是否合格;其二是根據(jù)測量的結(jié)果來指導光學系統(tǒng)的裝校。為消除被測件在測量儀器上的安裝定位過程帶來的誤差,必須對直接測量的數(shù)據(jù)進行修正。光軸擬合就是對測量數(shù)據(jù)的優(yōu)化和修正的過程。提出一種光軸擬合的數(shù)學模型,該數(shù)學模型結(jié)合了解析方法和數(shù)值分析方法,考慮了中心偏測量的實際情況,在嚴格的數(shù)學模型基礎上做了合理的簡化,使光軸的擬合問題最終轉(zhuǎn)化為對線性方程組的求解。當前第24頁\共有33頁\編于星期五\13點當前第25頁\共有33頁\編于星期五\13點3)應用最小二乘法得到關于四參數(shù)的線性方程組。得到各面球心的位置坐標后,按照一般直線擬合的方法,應使各球心對優(yōu)化軸距離的平方和最小,符合數(shù)學上的最小二乘法。N個球心到優(yōu)化軸距離的平方和:當前第26頁\共有33頁\編于星期五\13點當前第27頁\共有33頁\編于星期五\13點當前第28頁\共有33頁\編于星期五\13點當前第29頁\共有33頁\編于星期五\13點擴展:基于MATLAB的非線性方程組遺傳解法胡斐,趙治國(同濟大學汽車學院,上海201804)遺傳算法是一種基于自然選擇的用于求解有約束和無約束最優(yōu)問題的方法。遺傳算法反復修改包含若干個體的種群。遺傳算法在每一步中,隨機從當前種群中選擇若干個個體作為父輩,并用它們產(chǎn)生下一代子輩。在若干代之后,種群就朝著最優(yōu)解“進化”。我們可以利用遺傳
溫馨提示
- 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯(lián)系上傳者。文件的所有權益歸上傳用戶所有。
- 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁內(nèi)容里面會有圖紙預覽,若沒有圖紙預覽就沒有圖紙。
- 4. 未經(jīng)權益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
- 5. 人人文庫網(wǎng)僅提供信息存儲空間,僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護處理,對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對任何下載內(nèi)容負責。
- 6. 下載文件中如有侵權或不適當內(nèi)容,請與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準確性、安全性和完整性, 同時也不承擔用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。
最新文檔
- 優(yōu)化辦公流程管理規(guī)章制度
- 2025年巴彥淖爾貨運從業(yè)資格證考試卷
- 2025年江蘇貨運從業(yè)資格證模擬考試0題b2
- 人力資源企業(yè)勞動合同
- 2025年烏海貨運資格證考試有哪些項目
- 2025年宜昌貨運從業(yè)資格證模擬考試系統(tǒng)下載
- 2025年南京貨運從業(yè)資格考試技巧
- 動漫游戲創(chuàng)作與制作技術作業(yè)指導書
- 2025年南昌貨運從業(yè)資格證模擬考試題下載
- 2025年部編版語文小學三年級下冊口語與習作專項復習題
- 2025年湖南大眾傳媒職業(yè)技術學院單招職業(yè)技能測試題庫學生專用
- 一年級珍惜糧食主題班會學習教案
- 2025年高縣縣屬國企業(yè)公開招聘工作人員高頻重點提升(共500題)附帶答案詳解
- 2025年人教版數(shù)學五年級下冊教學計劃(含進度表)
- 海岸動力學英文課件Coastal Hydrodynamics-復習
- 碳足跡研究-洞察分析
- DB11-T 1191.3-2024 實驗室危險化學品安全管理要求 第3部分:科研單位
- 硬質(zhì)巖層組合切割開挖技術
- 2024解析:第二章聲現(xiàn)象-講核心(解析版)
- 2024年考研管理類綜合能力(199)真題及解析完整版
- 2024解析:第十章 浮力綜合應用-講核心(解析版)
評論
0/150
提交評論