开发工具:
文件大小: 1kb
下载次数: 0
上传时间: 2011-06-19
详细说明： #include #include # define eps 0.01 double f (double coe[], double x[]) //返回待求函数的函数值 { return coe[0]*pow(x[0],2)+coe[1]*pow(x[1],2)+coe[2]*x[0]*x[1]+coe[3]*x[0]+coe[4]*x[1]+coe[5]; } void grads (double coe[], double x[],double grads_x[]) //求解函数的梯度 { grads_x[0] = 2*coe[0]*x[0]+coe[2]*x[1]+coe[3]; grads_x[1] = 2*coe[1]*x[1]+coe[2]*x[0]+coe[4]; } void H (double coe[], double h[])//求解海赛矩阵 { h[0] = 2.0*coe[0]; h[1] = coe[2]; h[2] = coe[2]; h[3] = 2.0*coe[1]; } void inverseH (do uble h[], double inverseh[])//求解海赛矩阵的逆矩阵 { double deter; deter = h[0]*h[3]-h[1]*h[2]; inverseh[0] = h[3]/deter; inverseh[1] = -h[1]/deter; inverseh[2] = -h[2]/deter; inverseh[3] = h[0]/deter; } void main()//主函数 { double coe[6] = {1,2,-2,-4,0,0}, x[2] = {1 , 1}, grads_x[2], h[4], inverseh[4], fx = 0; grads(coe, x, grads_x);//求梯度 while(sqrt(pow(grads_x[0],2)+pow(grads_x[1],2))-eps>0) { H(coe, h);//求海赛矩阵 inverseH(h, inverseh);//海赛矩阵求逆 x[0] = x[0]-(inverseh[0]*grads_x[0]+inverseh[1]*grads_x[1]); x[1] = x[1]-(inverseh[2]*grads_x[0]+inverseh[3]*grads_x[1]); grads(coe, x, grads_x); } fx = f(coe, x); printf("最优解为："); for(int i=0;i<2;i++) { printf("\n%f\n",x[i]); } printf("\n最优值为："); printf("\nfx=%f \n",fx); } ...展开收缩
(系统自动生成,下载前可以参看下载内容)
下载文件列表
相关说明
- 本站资源为会员上传分享交流与学习,如有侵犯您的权益,请联系我们删除.
- 本站是交换下载平台，提供交流渠道，下载内容来自于网络，除下载问题外，其它问题请自行百度。
- 本站已设置防盗链，请勿用迅雷、QQ旋风等多线程下载软件下载资源，下载后用WinRAR最新版进行解压.
- 如果您发现内容无法下载，请稍后再次尝试；或者到消费记录里找到下载记录反馈给我们.
- 下载后发现下载的内容跟说明不相乎，请到消费记录里找到下载记录反馈给我们，经确认后退回积分.
- 如下载前有疑问，可以通过点击"提供者"的名字，查看对方的联系方式，联系对方咨询.
