共轭梯度法的C程序代码.pdf

///////////////////////////////////////// ///// vector.h 文件 ///// ///// 定义向量及其基本运算 ///// ///////////////////////////////////////// #include<math.h> #define MAXLENGTH 10 //向量定义 typedef struct { int tag; //行列向量标志。行向量为0，列向量为1。 int dimension ; //向量的维数 double elem[MAXLENGTH]; //向量的元素 }vector; vector vecCreat(int tag, int n){ //建立维数为n 的向量 vector x; x.tag=tag; x.dimension=n ; for(int i=0;i<n;++i){ cout<<"请输入第"<<i+1<<"分量:"; cin>>x.elem[i]; } return x; } double vecMultiply (vector a, vector b) { //行向量与列向量乘法运算 if((a.tag!=b.tag)&&(a.dimension==b.dimension)){//相乘条件 double c=0; for(int i=0;i<a.dimension ;++i) c+=a.elem [i]*b.elem[i]; return c; } } vector vecAdd (vector a, vector b) { //向量加法运算 if((a.tag==b.tag)&&(a.dimension==b.dimension)){//相加条件 vector c; c.dimension=a.dimension; c.tag=a.tag; for(int i=0;i<c.dimension ;++i) c.elem[i]=a.elem[i]+b.elem[i]; return c; } } vector vecConvert(vector a){ //向量转置 if(a.tag==0)a.tag=1; if(a.tag==1)a.tag=0; return a; } double vecMole (vector a){ //求向量模 double sum=0 ; for(int i=0;i<a.dimension ;++i) sum+=a.elem[i]*a.elem [i]; sum=sqrt(sum) ; return sum ; } vector numMultiply (double n,vector a) { //数乘向量 for(int i=0;i<a.dimension ;++i) a.elem[i]=n*a.elem[i]; return a; } void showPoint (vector x, double f){ //显示点，解. cout<<" x= ("; for(int i=0;i<x.dimension ;++i){ cout<<x.elem[i]; if(i!=x.dimension-1)cout<<","; } cout<<") f(x)="<<f<<endl; cout<<endl; } ///////////////////////////////////////////////////////// ////////////// function.h 文件 ///////////////////////// //////// 函数及其导数的设置均在此文件完成//////////////// ///////////////////////////////////////////////////////// // * 无 约 束 问 题 * // // 目标函数--在vecFun (vector vx)中修改,x1 改成x[1] // // x2 改成x[2],依此类推... // //////