資源描述:
《數(shù)值計(jì)算方法上機(jī)實(shí)驗(yàn)報(bào)告》由會(huì)員上傳分享,免費(fèi)在線閱讀,更多相關(guān)內(nèi)容在行業(yè)資料-天天文庫。
1、........數(shù)值計(jì)算方法上機(jī)實(shí)驗(yàn)報(bào)告實(shí)驗(yàn)?zāi)康模簭?fù)習(xí)和鞏固數(shù)值計(jì)算方法的基本數(shù)學(xué)模型,全面掌握運(yùn)用計(jì)算機(jī)進(jìn)行數(shù)值計(jì)算的具體過程及相關(guān)問題。利用計(jì)算機(jī)語言獨(dú)立編寫、調(diào)試數(shù)值計(jì)算方法程序,培養(yǎng)學(xué)生利用計(jì)算機(jī)和所學(xué)理論知識(shí)分析解決實(shí)際問題的能力。上機(jī)練習(xí)任務(wù):利用計(jì)算機(jī)基本C語言編寫并調(diào)試一系列數(shù)值方法計(jì)算通用程序,并能正確計(jì)算給定題目,掌握調(diào)試技能。掌握文件使用編程技能,如文件的各類操作,數(shù)據(jù)格式設(shè)計(jì)、通用程序運(yùn)行過程中文件輸入輸出運(yùn)行方式設(shè)計(jì)等。一、各算法的算法原理及計(jì)算機(jī)程序框圖1.列主元高斯消去法l算法原理:高
2、斯消去法是利用現(xiàn)行方程組初等變換中的一種變換,即用一個(gè)不為零的數(shù)乘一個(gè)方程后加只另一個(gè)方程,使方程組變成同解的上三角方程組,然后再自下而上對上三角方程組求解。列選住院是當(dāng)高斯消元到第步時(shí),從列的以下(包括)的各元素中選出絕對值最大的,然后通過行交換將其交換到的位置上。交換系數(shù)矩陣中的兩行(包括常數(shù)項(xiàng)),只相當(dāng)于兩個(gè)方程的位置交換了,因此,列選主元不影響求解的結(jié)果。.專業(yè)學(xué)習(xí)資料.........l計(jì)算機(jī)程序框圖如上l源程序:#defineN200#include"stdio.h"#include"math.h"FI
3、LE*fp1,*fp2;voidLZ(){intn,i,j,k=0,l;doubled,t,t1;.專業(yè)學(xué)習(xí)資料.........staticdoublex[N],a[N][N];fp1=fopen("a1.txt","r");fp2=fopen("b1.txt","w");fscanf(fp1,"%d",&n);for(i=0;i4、])>fabs(d))/*選主元*/{d=a[i][k];l=i;}i++;.專業(yè)學(xué)習(xí)資料.........}while(i5、]-=a[i][k]*a[k][j];k++;}while(k=0;i--)/*回代*/{t1=0;for(j=i+1;j6、列線性方程組.專業(yè)學(xué)習(xí)資料.........輸入3輸出結(jié)果:方程組的根為x[1]=6.00000012-38方程組的根為x[2]=4.00000021322方程組的根為x[3]=2.00000032128l輸入變量、輸出變量說明:輸入變量:系數(shù)矩陣元素,常向量元素輸出變量:解向量元素1.杜里特爾分解法解線性方程l算法原理:求解線性方程組時(shí),當(dāng)對進(jìn)行杜里特爾分解,則等價(jià)于求解,這時(shí)可歸結(jié)為利用遞推計(jì)算相繼求解兩個(gè)三角形(系數(shù)矩陣為三角矩陣)方程組,用順代,由求出,再利用回帶,由求出。計(jì)算機(jī)程序框圖:源程序:#incl
7、ude"stdio.h"#include"math.h"FILE*fp1,*fp2;voidmain().專業(yè)學(xué)習(xí)資料.........{inti,j,k,N;doubles,A[200][200],B[200],x[200],y[200];staticdoubleL[200][200],U[200][200];fp1=fopen("a2.txt","r");fp2=fopen("b2.txt","w");fscanf(fp1,"%d",&N);for(i=0;i8、canf(fp1,"%lf",&A[i][j]);}for(i=0;i