數(shù)值分析上機(jī)實(shí)驗(yàn)報(bào)告(參考版)
2025-08-06 23:50本頁(yè)面
【正文】 只有你自己才能把歲月描畫(huà)成一幅難以忘懷的人生畫(huà)卷。歲月是有情的，假如你奉獻(xiàn)給她的是一些色彩，它奉獻(xiàn)給你的也是一些色彩。努力過(guò)后，才知道許多事情，堅(jiān)持堅(jiān)持，就過(guò)來(lái)了。有時(shí)候覺(jué)得自己像個(gè)神經(jīng)病。在紛雜的塵世里，為自己留下一片純靜的心靈空間，不管是潮起潮落，也不管是陰晴圓缺，你都可以免去浮躁，義無(wú)反顧，勇往直前，輕松自如地走好人生路上的每一步3. 花一些時(shí)間，總會(huì)看清一些事。1. 若不給自己設(shè)限，則人生中就沒(méi)有限制你發(fā)揮的藩籬。，編程時(shí)必須注意。x=x39。b=[,]39。,。,。,。,。for k=n1:1:1 x(k)=(b(k)A(k,k+1:n)*x(k+1:n))/A(k,k)。 if flag~=0,Ab=[A,b],endendx=zeros(n,1)。 b(k+1:n)=b(k+1:n)m*b(k)。 end m=A(k+1:n,k)/A(k,k)。 if pk A([k p],:)=A([p k],:)。for k=1:(n1) [ap,p]=max(abs(A(k:n,k)))。}} 運(yùn)行結(jié)果 MATLAB上機(jī)程序unction [x]=mgauss2(A,b,flag)if nargin3,flag=0。 x[i]=(b[i]t3)/A[i][i]。j9。i=0。} x[8]=b[8]/A[8][8]。j++) A[i+1][j]=A[i+1][j]A[i][j]*t2/A[i][i]。 for(j=i。b[i+1]=t1。 t1=b[i]。A[i][j]=A[i+1][j]。ji+3。i8。 double t1=0,t2=0,t3=0。 double B[9][10]。j++) A[k][j]=u[k]*q[j]+u[j]*q[k]。k++) {for(j=0。 for(k=0。k9。 kr/=2*a。k9。} kr=0。j++) y[k]+=A[k][j]*u[j]。 for(j=0。k9。 else if(ji+1) u[j]=A[j][i]。j++) {if(j=i) u[j]=0。 for(j=0。 s=sqrt(s)。j9。i++) {s=0。 for(i=0。}void Householder(double A[9][9]) {double q[9],u[9],y[9],s,a,kr。i10。 expunction(A,b,x)。 Householder(A)。 double x[9]={}。void main(){double A[9][9]={{,},{,},{,}, {,},{,},{,}, {,},{,},{,}}。最終形式如下：(A,b）～ C語(yǔ)言程序代碼（1）比較該列的元素的絕對(duì)值的大小，將絕對(duì)值最大的元素通過(guò)行變換使其位于主對(duì)角線上；（2）進(jìn)行高斯消去法變換，把系數(shù)矩陣化成上三角形，然后回代求include include void Householder(double A[9][9])。將變換后的矩陣的第二列中第二行的元素與其下面的此列的元素逐一進(jìn)行比較，找到最大的元素，將第k行的元素與第2行的元素進(jìn)行交換，然后通過(guò)行變換，將第2列中第3到第n個(gè)元素都消成0。具體做法如下：將方陣A和向量b寫(xiě)成C=（A，b）。 理論依據(jù)：列主元素消元法是在應(yīng)用Gauss消元法的基礎(chǔ)上，憑借長(zhǎng)期經(jīng)驗(yàn)積累提出的，是線性方程組一般解法，目的是為避免在消元計(jì)算中使誤差的擴(kuò)大，甚至嚴(yán)重?fù)p失了有效數(shù)字使數(shù)據(jù)失真，而在每次初等變換前對(duì)矩陣作恰當(dāng)?shù)恼{(diào)整，以提高Gauss消元法的數(shù)字穩(wěn)定性，進(jìn)而提高計(jì)算所得數(shù)據(jù)的精確度。所以編寫(xiě)過(guò)程中一定要注意各循環(huán)的次數(shù)，以免出錯(cuò)。}} 運(yùn)行結(jié)果： 問(wèn)題討論：，精度較高，誤差小且程序相對(duì)簡(jiǎn)單，存儲(chǔ)量少。 if(j==3) printf(\n)。j=3。i5。} printf(計(jì)算結(jié)果:\n)。i=3。i++
環(huán)評(píng)公示相關(guān)推薦
鄂ICP備17016276號(hào)-1