版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請(qǐng)進(jìn)行舉報(bào)或認(rèn)領(lǐng)
文檔簡(jiǎn)介
第二章解線性方程組的直接方法Gaus消去法簡(jiǎn)單易行,但其計(jì)算過程中,要求a稱為主元素)均不為零,因而適用范圍小,只適用于從1到n-1階順序主子式均不為零的矩陣A,計(jì)算實(shí)踐還表明,Gaus消去法的數(shù)值穩(wěn)定性差,當(dāng)出現(xiàn)小主元素時(shí),會(huì)嚴(yán)重影響計(jì)算結(jié)果的精度,甚至導(dǎo)出錯(cuò)誤的結(jié)果§2主元素法先看一個(gè)例子學(xué)魂第二章解線性方程組的直接方法例2求解方程組0.50x+1.lx+3.1x,=6.02.0x,+4.5x+0.36x=0.0205.0x,+0.96x+6.5x,=0.96式(2-10a)中所有系數(shù)均有2位有效數(shù)字[解]為減少誤差,計(jì)算過程中保留3位有效數(shù)字按Gaus消去法步驟,第一次消元得同解方程組0.50x1+1.1x2+31x3=600.100x2-12.0x3=-24.010.0x2-24.5x3=-59.0第二章解線性方程組的直接方法第二次消元得0.50x1+1.1x2+3.1x3=600.100x2-120x3=-24.01220c2460回代得解x3=202x2=240x1=-580容易驗(yàn)證,方程組(210)的準(zhǔn)確解為x1=-260,x2=1.00,x3=2.00顯然兩者相差很大但若在解方程組前,先把方程的次序?qū)W魂第二章解線性方程組的直接方法交換一下,如把(2-10a)改寫成5.0x,+0.96x,+6.5x2=0.962.0x,+4.5x+0.36x,=0.020(2-10b)0.50x1+1.1x2+3.1x3=6.0再用Gaus消去法求解,消元后得同解方程5.0x,+0.96x,+6.5x2=0.964.12x2-2.24x3=-0.3642.99x3=599學(xué)魂第二章解線性方程組的直接方法回代得解2.00,x2=1.00,x=-260與準(zhǔn)確解相同產(chǎn)生上述現(xiàn)象的原因在于舍入誤差.因?yàn)榘词?2-10)的方程順序進(jìn)行消元時(shí),主元a(=0.50,a2)=0.100都比較小,以它們?yōu)槌龜?shù)就增長(zhǎng)了舍入誤差,從而導(dǎo)致計(jì)算結(jié)果不準(zhǔn)確。為了在計(jì)算過程中,抑制舍入誤差的增長(zhǎng),應(yīng)盡量避免小主元的出現(xiàn)如例2中第二種解法,通過交換方程次序,選取絕對(duì)值大的元素作主元基于這種想法導(dǎo)出了主元素法5學(xué)魂第二章解線性方程組的直接方法2.2.1列主元素法為簡(jiǎn)便起見,我們用方程組(2-1)的增廣矩陣A,bI(2-11)b表示它,并直接在增廣矩陣上進(jìn)行運(yùn)算.學(xué)魂二章解線性方程組的直接方法列主元素法基本思想:是在每次消元前,在要消去未知數(shù)的系數(shù)中找到絕對(duì)值最大的系數(shù)作主元,通過方程對(duì)換將其換到對(duì)角線上,然后進(jìn)行消元。具體步驟如下:第一步:首先在矩陣(2-11)的第1列中選取絕對(duì)值最大的元,比如為a1,則1a;I-maxlai1將(21中第1行與第i1行互換.為方便起見記行互換后的增廣矩陣為[A①,b①],然后進(jìn)行第一次消元,得矩陣學(xué)魂第二章解線性方程組的直接方法[A2),b2]=b第二步:在矩陣[A),b2]的第2列中選主元,比如,使1a2=maxla21.將矩陣[A(,b]的第22<i<n行與第i2行互換,再進(jìn)行第二次消元,得矩陣[A0,b0第k步:在矩陣[A,b]的第k列中選主元,如,使1aik=maxla1將[A(,b的的第k行與第k行互換,進(jìn)行第k次消元8學(xué)魂第二章解線性方程組的直接方法如此經(jīng)過n-1步,增廣矩陣(2-11)被化成上三角形,最后由回代過程求解。在上述過程中,主元是按列選取的,列主元素法由此得名例2中的第二種解法就是按列主元素法進(jìn)行的2.2.2全主元素法如果不是按列選主元,而是在全體待選系數(shù)a(,j=k,k+,…,n)中選取主元,則得到全主元素法,其計(jì)算過程如下:學(xué)魂第二章解線性方程組的直接方法第一步:在全體系數(shù)an(,j=…,n)中選取絕對(duì)值最大的元作為主元,并通過行與列的互換把它換到a1的位置,然后進(jìn)行第一次消元得到矩陣[A,b1第k步:在矩陣A的右下方(n-k+1)階子矩陣的所有元素a(i,j=k,k+1,…,
溫馨提示
- 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請(qǐng)下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請(qǐng)聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
- 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁內(nèi)容里面會(huì)有圖紙預(yù)覽,若沒有圖紙預(yù)覽就沒有圖紙。
- 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
- 5. 人人文庫(kù)網(wǎng)僅提供信息存儲(chǔ)空間,僅對(duì)用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理,對(duì)用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對(duì)任何下載內(nèi)容負(fù)責(zé)。
- 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請(qǐng)與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時(shí)也不承擔(dān)用戶因使用這些下載資源對(duì)自己和他人造成任何形式的傷害或損失。
最新文檔
- 第8課 歐洲的思想解放運(yùn)動(dòng) 說課稿-2024-2025學(xué)年高一統(tǒng)編版2019必修中外歷史綱要下冊(cè)001
- 粵教版 信息技術(shù) 必修 5.1 認(rèn)識(shí)信息資源管理說課稿001
- 教學(xué)計(jì)劃(含進(jìn)度表)(說課稿)-2024-2025學(xué)年人教PEP版英語四年級(jí)上冊(cè)
- 2025版防水材料專業(yè)代理服務(wù)合同6篇
- 二零二五年度危險(xiǎn)品搬運(yùn)安全責(zé)任合同3篇
- 2024年高標(biāo)準(zhǔn)鋪裝工程承包協(xié)議版
- 2025年度電子廢棄物ROHS檢測(cè)與處理質(zhì)量保障協(xié)議
- xx省量子科技項(xiàng)目可行性研究報(bào)告
- 2024年面包磚環(huán)保生產(chǎn)技術(shù)轉(zhuǎn)讓合同
- 2025版酒店智能化改造項(xiàng)目合同范本3篇
- 中建三局住宅工程精益建造實(shí)施指南
- 分布式光伏發(fā)電項(xiàng)目并網(wǎng)驗(yàn)收意見單
- 網(wǎng)站隱私政策模板
- YY∕T 1831-2021 梅毒螺旋體抗體檢測(cè)試劑盒(免疫層析法)
- 消弧產(chǎn)品規(guī)格實(shí)用標(biāo)準(zhǔn)化規(guī)定
- 裝飾裝修工程施工合理化建議和降低成本措施提要:完整
- 第十四章35kV變電站保護(hù)整定值計(jì)算實(shí)例
- 液態(tài)模鍛工藝介紹
- 水泵水輪機(jī)結(jié)構(gòu)介紹
- 井式爐課程設(shè)計(jì)說明書
- 拼音四線三格加田字格模板(A4打印版可編輯打字)
評(píng)論
0/150
提交評(píng)論