三元函数最小共轭梯度法.docVIP

下载本文档

47
0
约 5页
2017-07-14 发布于河南
举报
版权申诉

三元函数最小共轭梯度法.doc

1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
7、VIP文档为合作方或网友上传，每下载1次，网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

三元函数最小共轭梯度法

#includeiostream.h #includemath.h class Vector { public: Vector(){x1=0;x2=0;x3=0;} Vector(double xx1,double xx2,double xx3){x1=xx1;x2=xx2;x3=xx3;} friend Vector operator + (Vector v1,Vector v2); friend Vector operator - (Vector v1,Vector v2); friend Vector operator - (Vector v1); friend double operator * (Vector v1,Vector v2); friend Vector operator * (double c,Vector v); friend Vector operator * (Vector v,double c); ~Vector(){} void display(); //private: double x1,x2,x3; }; Vector operator + (Vector v1,Vector v2) { return Vector(v1.x1+v2.x1,v1.x2+v2.x2,v1.x3+v2.x3); } Vector operator - (Vector v1,Vector v2) { return Vector(v1.x1-v2.x1,v1.x2-v2.x2,v1.x3-v2.x3); } Vector operator - (Vector v1) { return Vector(-v1.x1,-v1.x2,-v1.x3); } double operator * (Vector v1,Vector v2) { return (v1.x1*v2.x1+v1.x2*v2.x2+v1.x3*v2.x3); } Vector operator * (double c,Vector v) { return Vector(c*v.x1,c*v.x2,c*v.x3); } Vector operator * (Vector v,double c) { return Vector(c*v.x1,c*v.x2,c*v.x3); } void Vector::display() { cout(x1,x2,x3)endl; } double ValueFunction(Vector p)//求某点处函数值的函数，返回值为一个double型数 { return p.x1*p.x1-2*p.x1*p.x2+2*p.x2*p.x2+p.x3*p.x3-p.x2*p.x3+2*p.x1+3*p.x2-p.x3; } Vector GradFuction(Vector p)//求梯度的函数，返回值为一个向量 { return Vector(2*(p.x1-p.x2)+2,-2*p.x1+4*p.x2-p.x3+3,2*p.x3-p.x2-1); } double GetModulus(Vector p) { return sqrt(p.x1*p.x1+p.x2*p.x2+p.x3*p.x3); } /*double GetFValue(double la) //求函数值f(x[k]+lamda*s[k]) { return ValueFunction(x[k]+la*s[k]); }*/ void main() { int n=3;//本题为三维问题 double ValueFunction(Vector p); Vector GradFuction(Vector p); Vector startpoint(0,0,0); Vector x[10000]={startpoint},g[10000],s[10000];double lamda[10000],lamdak; int k; cout选取初始点为：x[0]= ;x[0].display(); double precision=0.000001; step2:{g[0]=GradFuction(x[0]);s[0]=-g[0];k=0;} step3:{ double GetFValue(double la); double f1,f0,la1,la0=1,h=10,e=0.0000001; step12: la1=la0+h; f1=ValueFunction(x[k]+la1*s[k])