下載本文檔
版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請(qǐng)進(jìn)行舉報(bào)或認(rèn)領(lǐng)
文檔簡介
1、C語言矩陣求逆程序(高斯-約旦法)高斯-約旦法根據(jù)代數(shù)里面的知識(shí),可以使用伴隨矩陣也可以使用初等行變換來解求解,但是這樣如果矩陣的維數(shù)較大的時(shí)候,使用這種方法,矩陣的維數(shù)變大時(shí),計(jì)算量急劇的變大,計(jì)算時(shí)間和使用內(nèi)存也會(huì)按著指數(shù)急劇上升,這樣的算法的生命力不行。使用以下這種算法的計(jì)算量和使用內(nèi)存不會(huì)發(fā)生急劇的變化,特別是矩陣在維數(shù)大的時(shí)候。高斯-約旦法(全選主元)求逆的步驟如下:首先,對(duì)于k從0到n-1作如下幾步:從第k行、第k列開始的右下角子陣中選取絕對(duì)值最大的元素,并記住次元素所在的行號(hào)和列號(hào),在通過行交換和列交換將它交換到主元素位置上。這一步稱為全選主元。m(k,k)=1/m(k,k)m(
2、k,j)=m(k,j)*m(k,k),j=0,1,.,n-1;j!=km(i,j)=m(i,j)-m(i,k)*m(k,j),i,j=0,1,.,n-1;i,j!=km(i,k)=-m(i,k)*m(k,k),i=0,1,.,n-1;i!=k最后,根據(jù)在全選主元過程中所記錄的行、列交換的信息進(jìn)行恢復(fù),恢復(fù)的原則如下:在全選主元過程中,先交換的行(列)后進(jìn)行恢復(fù);原來的行(列)交換用列(行)交換來恢復(fù)。#includestdio.h#includemalloc.h#includemath.h/數(shù)學(xué)函數(shù)voidmain()intinv(double*p,intn);doublea44=1,2,0,
3、0,2,5,0,0,0,0,3,0,0,0,0,1,*ab;ab=a0;intn=4,i=0,j;i=inv(ab,n);/調(diào)用矩陣求逆if(i!=0)/如果返回值不是0for(i=0;in;i+)/輸出結(jié)果putchar(n);for(j=0;jn;j+)printf(%f,aij);intinv(double*p,intn)voidswap(double*a,double*b);int*is,*js,i,j,k,l;for(i=0;in;i+)putchar(n);for(j=0;jn;j+)printf(%f,*(p+i*n+j);puts(nnnn);doubletemp,fmax;i
4、s=(int*)malloc(n*sizeof(int);js=(int*)malloc(n*sizeof(int);for(k=0;kn;k+)fmax=0.0;for(i=k;in;i+)for(j=k;jfmax)fmax=temp;isk=i;jsk=j;if(fmax+1.0)=1.0)free(is);free(js);printf(noinv);return(0);if(i=isk)!=k)for(j=0;jn;j+)swap(p(k*n+j),p(i*n+j);交換指針if(j=jsk)!=k)for(i=0;in;i+)swap(p(i*n+k),p(i*n+j);/交換指針pk*n+k=1.0/pk*n+k;for(j=0;jn;j+)if(j!=k)pk*n+j*=pk*n+k;for(i=0;in;i+)if(i!=k)for(j=0;jn;j+)if(j!=k)pi*n+j=pi*n+j-pi*n+k*pk*n+j;for(i=0;i=0;k-)if(j=jsk)!=k)for(i=0;in;i+)swap(p+j*n+i),(p+k*n+i);if(i=isk)!=k)for(j=0;jn;j+)swap(p+j*n+i),(p+j*n+k);f
溫馨提示
- 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. 人人文庫網(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ì)自己和他人造成任何形式的傷害或損失。
最新文檔
- 2025年滬教新版九年級(jí)歷史上冊(cè)階段測(cè)試試卷
- 2025年統(tǒng)編版九年級(jí)歷史下冊(cè)階段測(cè)試試卷含答案
- 2025年度企業(yè)內(nèi)部培訓(xùn)項(xiàng)目合作協(xié)議合同范本4篇
- 2025年人教新起點(diǎn)必修1歷史下冊(cè)階段測(cè)試試卷含答案
- 2025年人教B版選修3歷史下冊(cè)月考試卷
- 二零二五版苗木養(yǎng)護(hù)管理服務(wù)合同2篇
- 二零二五版南京二手房買賣電子合同范本4篇
- 二零二五年度二零二五互聯(lián)網(wǎng)金融服務(wù)合同規(guī)范范本4篇
- 2025年度墻體改梁與建筑智能化升級(jí)合同4篇
- 二零二五年度農(nóng)藥產(chǎn)品生產(chǎn)安全培訓(xùn)采購合同4篇
- 環(huán)境監(jiān)測(cè)對(duì)環(huán)境保護(hù)的意義
- 2023年數(shù)學(xué)競(jìng)賽AMC8試卷(含答案)
- 神經(jīng)外科課件:神經(jīng)外科急重癥
- 2024年低壓電工證理論考試題庫及答案
- 2023年十天突破公務(wù)員面試
- 《瘋狂動(dòng)物城》中英文對(duì)照(全本臺(tái)詞)
- 醫(yī)院住院醫(yī)師規(guī)范化培訓(xùn)證明(樣本)
- 小學(xué)六年級(jí)語文閱讀理解100篇(及答案)
- 氣功修煉十奧妙
- 安徽省物業(yè)服務(wù)標(biāo)準(zhǔn)
- 勾股定理的歷史與證明課件
評(píng)論
0/150
提交評(píng)論