共軛梯度法C語言(西安交大)

1、精選優(yōu)質(zhì)文檔-傾情為你奉上精選優(yōu)質(zhì)文檔-傾情為你奉上專心-專注-專業(yè)專心-專注-專業(yè)精選優(yōu)質(zhì)文檔-傾情為你奉上專心-專注-專業(yè)#include#include#define N 10 /*定義矩陣階數(shù)*/void main()int i,j,m,ANN,BN;double XN,akvN,dkaN,rkN,dkN,pk,pkk,ak,bk;for(i=0;iN;i+) /*輸入系數(shù)矩陣A*/for(j=0;jN;j+)if(i=j)Aij=-2;else if(abs(i-j)=1) Aij=1;else Aij=0; for(i=0;iN;i+) /*輸入矩陣B*/ if(i=0|i=N-1

2、) Bi=-1; else Bi=0; printf(n); /*輸出系數(shù)矩陣A*/ printf(the Maxtr A isn); for(i=0;iN;i+) printf(n); for(j=0;jN;j+) printf(%3d,Aij); printf(n); /*輸出矩陣B*/printf(the Maxtr b isn);for(i=0;iN;i+) printf(%3d,Bi);printf(n); printf(input the Maxtr Xn); /*給X0輸入一個(gè)初始向量*/ for(i=0;iN;i+) Xi=0; printf(X chushi:n); for(

3、i=0;iN;i+) /*顯示給定的初始向量X0*/ printf(%3.0f,Xi); /*開始計(jì)算*/ for(i=0;iN;i+) /*計(jì)算rk*/ pk=0.0; for(j=0;jN;j+) pk+=Aij*Xj; akvi=pk; for(i=0;iN;i+) rki=Bi-akvi; for(i=0;iN;i+) /*給rO賦值*/ dki=rki; for(m=0;mN;m+) /*開始迭代循環(huán)*/ for(i=0;iN;i+) /*計(jì)算ak*/ pk=0.0; for(j=0;jN;j+) pk+=Aij*dkj; dkai=pk; pk=0.0;pkk=0.0; for(i

4、=0;iN;i+) pk+=dkai*dki; pkk+=rki*rki;if(pkk=0) /*如果殘差pkk=0,計(jì)算結(jié)束*/ break; ak=pkk/pk; for(i=0;iN;i+) /*計(jì)算X(k+1)*/ Xi=Xi+ak*dki; for(i=0;iN;i+) /*計(jì)算r(k+1)*/ pk=0.0; for(j=0;jN;j+) pk+=Aij*Xj; akvi=pk; for(i=0;iN;i+) rki=Bi-akvi; pk=0.0; for(i=0;iN;i+) /*計(jì)算bk*/ pk+=rki*rki; bk=pk/pkk; for(i=0;iN;i+) /*計(jì)算d(k+1)*/ dki=rki+bk*dki; printf(n); printf(X