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

///////////////////////////////////////// ///// vector.h 文件 ///// ///// 定义向量及其基本运算 ///// ///////////////////////////////////////// #includemath.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;in;++i){ cout请输入第i+1分量:; cinx.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;ia.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;ic.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;ia.dimension ;++i) sum+=a.elem[i]*a.elem [i]; sum=sqrt(sum) ; return sum ; } vector numMultiply (double n,vector a) { //数乘向量 for(int i=0;ia.dimension ;++i) a.elem[i]=n*a.elem[i]; return a; } void showPoint (vector x, double f){ //显示点，解. cout x= (; for(int i=0;ix.dimension ;++i){ coutx.elem[i]; if(i!=x.dimension-1)cout,; } cout) f(x)=fendl; coutendl; } ///////////////////////////////////////////////////////// ////////////// function.h 文件 ///////////////////////// //////// 函数及其导数的设置均在此文件完成//////////////// ///////////////////////////////////////////////////////// // * 无 约 束 问 题 * // // 目标函数--在vecFun (vector vx)中修改,x1 改成x[1] // // x2 改成x[2],依此类推... // //////////////////////////////////////////////////