实验报告阻尼牛顿法.doc

太原理工大学机械学院机测系课程上机实验报告 课 程 名 称： 机械优化设计 班级 ? 日 期 ?成绩评定 ? 姓名 实验室 ?老师签名 ? 实验名称 用阻尼牛顿法程序解题 所用软件 C++ DEV 实验目的 及 内 容 ?实验目的： 1.掌握并能够建立最优化基本类型问题的数学模型。 2.掌握最优化方法的基本概念、基本理论和基本方法，奠定最优化的理论基础。 3.能够熟练编制和调试最优化方法的程序，奠定解决实际中的优化问题的基础 实验内容： 理解阻尼牛顿法并编写相关程序求其最优解。 牛顿法程序考核题 实 验 ?原 理 步 骤 、 实验原理： 实验步骤： 1，画流程图，编写程序； 2，将目标函数代入； 3，编译运行，将结果保存 实 验 结 果 及 分 析 *********阻尼牛顿法计算结果********** ++++++一维搜索方法：黄金分割法++++++ 初始坐标： x( 0)=[ 0.0000000, 0.0000000], f( 0)= 16.0000000 迭代轮数 k= 1 x( 1)=[ -1.1249988, 0.7499992], f( 1)= 9.8125000 迭代精度： 0.000009986 迭代轮数 k= 2 x( 2)=[ -1.1250000, 0.7500000], f( 2)= 9.8125000 迭代精度： 0.000000114 ********************* 阻尼牛顿法法优化最优点及目标函数值为： x( * )=[ -1.1250000, 0.7500000], f( * )= 9.8125000 迭代精度： 0.000000114 算法程序实现 /*csssqj.cpp */ #include string.h #include stdio.h #include math.h #include stdlib.h #include time.h #define N 2 /*优化设计维数*/ #define EPSIN 0.000001 /*迭代精度*/ #define H_QJ 1.0 /*初始区间搜索步长*/ #define Y_F 1 /*一维搜索方法选择：1——黄金分割法*/ /* 2——二次插值法*/ FILE *fp; char outname[50]=阻尼牛顿法计算结果.txt; /*计算结果输出文件*/ /*给出初始点坐标*/ void csd_x(double x0[]) { int i; for(i=0;iN;i++) /*初始点为坐标原点的情况*/ x0[i]=0.0; return; } /*目标函数*/ double hanshu(double x[]) { double f; f=4.0*(x[0]+1.0)*(x[0]+1.0)+2.0*(x[1]-1.0)*(x[1]-1.0) +x[0]+x[1]+10.0; /*阻尼牛顿法*/ return f; } /*函数的梯度*/ void gread(double x[],double gf[]) { /*阻尼牛顿法*/ gf[0]=8.0*(x[0]+1.0)+1.0; gf[1]=4.0*(x[1]-1.0)+1.0; return; } /*Heisen矩阵*/ void heisen(double x[],double hei[][N]) { hei[0][0]=8.0; hei[0][1]=0.0; hei[1][0]=0.0; hei[1][1]=4.0; return; } /*计算f(xk+as)*/ double xkadd(double x[],double d[],double a) { int i; double x1[N]; for(i=0;iN;i++) x1[i]=x[i]+a*d[i]; ret