2025年數(shù)值分析線性方程組求解高效直接法實踐探索報告_第1頁
2025年數(shù)值分析線性方程組求解高效直接法實踐探索報告_第2頁
2025年數(shù)值分析線性方程組求解高效直接法實踐探索報告_第3頁
2025年數(shù)值分析線性方程組求解高效直接法實踐探索報告_第4頁
2025年數(shù)值分析線性方程組求解高效直接法實踐探索報告_第5頁
已閱讀5頁，還剩12頁未讀， 繼續(xù)免費閱讀

下載本文檔

版權說明：本文檔由用戶提供并上傳，收益歸屬內容提供方，若內容存在侵權，請進行舉報或認領

文檔簡介

計算措施試驗二試驗匯報專業(yè)班級：姓名：學號：試驗成績：1.【試驗題目】2.【試驗目的】●掌握三角分解法與追趕法3.【試驗內容】。4.【試驗規(guī)定】(2)觀測小主元并分析對計算成果的影響。趕法公式。A=G?/b?β=c/(b-aβ)y=(f-ay_)/(b-aβ_?)a6.【源程序(帶注釋)】(1)一：列主元素消去法(cout<<"下面請輸入未知數(shù)的個數(shù)m=";cout<<"請按次序輸入增廣矩陣a:"<<endl;c=(fabs(a[j][i])>fabs(a[i]/*找列最大元素*/}/*將列最大數(shù)防在對角線上*/a[k][r]=a[k][r]-1[k][i]}x[m-1]=a[m-1][m]/a[m-1][x[i]=(a[i][m]-d)/a[icout<<"x["<<i<<"]="<<}j<m+1;j++)}一般消去法voidsolve(floatI[][100],floatu[][100],flfor(i=1;i<=n;i++)/*第一次回代過程開始*/{{y[i]=(b[i]+s1)/L[i]}for(i=n;i>=1;i--)/*第二次回代過程開始*/}}}{floata[100][100],I[100][100],u[100][100],x[100],b[100];for(i=1;i<=99;i++)/*將所有的數(shù)組置零，同步將L矩陣的對角值設為1*/{}printf("inputn:\n");/*輸入方程組的個數(shù)*/printf("inputarrayA:\n");/*讀取原矩陣A*/for(r=1;r<=n;r++)/*求解矩陣L和U*/{{}{}printf("arrayL:\n");/*輸出矩陣L*/for(i=1;i<=n;i++)printf("%7.3f",I[i][j]);}printf("arrayU:\n");/*輸出矩陣U*/{printf("%7.3f",u[i]lj]);printf("\n");}printf("}2)消元過程中，即使,用其作除數(shù)會導致計算中間結果數(shù)量級嚴重增長和舍入誤差的累積、擴大，最后使得計算結果不可靠。3)應避免采用絕對值很小的主元素ak);對一般的系數(shù)矩陣，最好保持乘數(shù)|≤1,因此，在高斯消去法中應引進選主元技巧，以便減少計算過程中舍入誤差對求解的影響。doublea[N],b[N],c[N-1],f[N],voidLUDecomposevoidbackSubs();//回代{{printf("由于該維數(shù)過于犀利，導致程序退出!");}printf("\n輸入下三角元素\n");printf("\n輸入主對角線元素\n");printf("\n輸入上三角元素\n");printf("x%

溫馨提示

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

評論

0/150

提交評論