1、精選優質文檔-傾情為你奉上#include<stdio.h>#include<math.h>#define N 10 /*定義矩陣階數*/void main()int i,j,m,ANN,BN;double XN,akvN,dkaN,rkN,dkN,pk,pkk,ak,bk;for(i=0;i<N;i+) /*輸入系數矩陣A*/for(j=0;j<N;j+)if(i=j)Aij=-2;else if(abs(i-j)=1) Aij=1;else Aij=0; for(i=0;i<N;i+) /*輸入矩陣B*/ if(i=0|i=N-1) Bi=-1;

2、else Bi=0; printf("n"); /*輸出系數矩陣A*/ printf("the Maxtr A isn"); for(i=0;i<N;i+) printf("n"); for(j=0;j<N;j+) printf("%3d",Aij); printf("n"); /*輸出矩陣B*/printf("the Maxtr b isn");for(i=0;i<N;i+) printf("%3d",Bi);printf("

3、n"); printf("input the Maxtr Xn"); /*給X0輸入一個初始向量*/ for(i=0;i<N;i+) Xi=0; printf("X chushi:n"); for(i=0;i<N;i+) /*顯示給定的初始向量X0*/ printf("%3.0f",Xi); /*開始計算*/ for(i=0;i<N;i+) /*計算rk*/ pk=0.0; for(j=0;j<N;j+) pk+=Aij*Xj; akvi=pk; for(i=0;i<N;i+) rki=Bi-a

4、kvi; for(i=0;i<N;i+) /*給rO賦值*/ dki=rki; for(m=0;m<N;m+) /*開始迭代循環*/ for(i=0;i<N;i+) /*計算ak*/ pk=0.0; for(j=0;j<N;j+) pk+=Aij*dkj; dkai=pk; pk=0.0;pkk=0.0; for(i=0;i<N;i+) pk+=dkai*dki; pkk+=rki*rki;if(pkk=0) /*如果殘差pkk=0,計算結束*/ break; ak=pkk/pk; for(i=0;i<N;i+) /*計算X(k+1)*/ Xi=Xi+ak*dki; for(i=0;i<N;i+) /*計算r(k+1)*/ pk=0.0; for(j=0;j<N;j+) pk+=Aij*Xj; akvi=pk; for(i=0;i<N;i+) rki=Bi-akvi; pk=0.0; for(i=0;i<N;i+) /*計算bk*/ pk+=rki*rki; bk=pk/pkk; for(i=0;i<N;i+) /*計算d(k+1)*/ dki=rki+bk*dki; printf("n"); printf("X cacualtate value is