共轭梯度法.ppt

文档介绍：1
第1页/共17页
2
第2页/共17页
3
迭代系数
迭代方向
第3页/共17页
4
第4页/共17页
5
第5页/共17页
6
第6页/共17页
7
第7页/共17页
8
int main(){
double A[2][2]={{6,3}, {3,2}};
double b[2]={0,-1}, x[2]={0,0}, r[2], Ar[2];
double a;
int i,j,k;
for(k=1; k<=4; k++){
for(i=0; i<2; i++)
r[i] = b[i] - (A[i][0]*x[0] + A[i][1]*x[1]);
for(i=0; i<2; i++)
Ar[i] = A[i][0]*r[0] + A[i][1]*r[1];
a = (r[0]*r[0] + r[1]*r[1]) / (Ar[0]*r[0] + Ar[1]*r[1]);
for(i=0; i<2; i++)
x[i] = x[i] + a*r[i];
printf("k=%2d, ", k);
for(i=0; i<2; i++) printf("x%d=%, ", i, x[i]);
printf("\n");
}
}
第8页/共17页
9
第9页/共17页
10
第10页/共17页