數(shù)值分析第二次大作業(yè)-文庫吧資料
2025-06-23 06:27本頁面
【正文】 ]= x[m]。m10。 gauss ( y)。 } sum= 0 。m10。m++) { sum+=u[m]*u[m]。 for (m=0。i++) u[i]=1。 for ( i=0。 double u[10],y[10]。 } x[k]= (b[k]f)/A[k][k]。j10。k) { f= 0。 for(k=8。 b[i]=b[i]mik*b[k]。j10。i++) /*將第K列對角線以下消元為零*/ { mik = A[i][k]/A[k][k]。 } for(i=k+1。 b[k]=b[mk]。 A[mk][j]=f。j++) { f = A[k][j]。 } } if(mk!=k) /* 將第K列主元素換行到對角線上*/ { for(j=k。i++) //*選擇第K列主元素 { if(fabs(mm)fabs(A[i][k])) { mm = A[i][k]。 //記錄最大組員的行號 for(i=k。k++) /*選主元并消元*/ { mm=A[k][k]。 for(k = 0。 } }} //列主元的高斯消元法求解特征向量void gauss(double b[10]){ int i,j,k,mk。j10。i10。}//構(gòu)造函數(shù)translation 用于進行原點平移void translation (double x) { int i=0,j=0。 temp= ((v1v2)?v1:v2)。 } } } //sp 5 } }double min(double v1,double v2){ return((v2v1)?v1:v2)。 j++) { mk[i][j] = mk[i][j] ur[i] * vr[j]。 i++) { for(j = 0。 for(i = 0。 i = m。 tr = tmp_tr / hr。 i = m。 } tmp_tr = 0。 pr[i] = tmp_pr / hr。 tmp_qr = tmp_qr + A[i][j] * ur[j]。 j++) { tmp_vr = tmp_vr + mk[j][i] * ur[j]。 for(j = 0。 tmp_pr = 0。 i=m。 ur[r] = mk[r][r]cr。 i = m。 i++) ur[i] = 0。//sp 3 for(i = 0。 cr = * sgn(mk[r][r]) * dr。 i++) tmp_dr = tmp_dr + mk[i][r] * mk[i][r]。 for(i = r。 if (tmp_ir==mr) continue。i=m。r++) {//sp 1 tmp_ir = 0。for (r=0。int r, i, j。 i++) mk[i][i]+= mt。 } } for(i = 0。 k++) tmp_sum+= A[i][k] * A[k][j]。 for(k = 0。 j = m。 i = m。 double tmp_sum。 k++。 qmk(mk,m,ms,mt)。 else { ms = A[m1][m1]+ A[m][m]。 else continue。 return 0。 m。 tzi[n]= s1_i。 m。 tzi[n]= s2_i。amp。 return 0。 tzr[
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1