外点罚函数法.doc

《外点罚函数法.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、太原理工大学机械学院机测系课程上机实验报告课程名称：机械优化设计班级机测2班日期 成绩评定 姓名陈明鑫实验室机械馆机房 老师签名 实验名称利用外点罚函数法求解相关函数的极小值点所用软件实验目的及内容 实验目的：1.掌握并能够建立最优化基本类型问题的数学模型。2.掌握最优化方法的基本概念、基本理论和基本方法，奠定最优化的理论基础。3.能够熟练编制和调试最优化方法的程序，奠定解决实际中的优化问题的基础实验内容：理解外点罚函数法并编写相关程序求其极小值点。实验 原理步骤、实验结果及分析外点罚函数法计算结果**********外点惩罚函数法计算结果**

2、******************无约束优化方法：鲍威尔法**********++++++一维搜索方法：黄金分割法++++++**********初始惩罚因子：r0=1.00********************递增系数：c0=10.00**********初始坐标：x(0)=[0.,0.],f(0)=1.迭代轮数k=1x(1)=[0.,0.],f(1)=0.迭代精度：0.迭代轮数k=2x(2)=[0.,-0.],f(2)=0.迭代精度：0.迭代轮数k=3x(3)=[0.,0.],f(3)=0.迭代精度：0.迭代轮数k=4x(4)=[0.,

3、-0.],f(4)=0.迭代精度：0.迭代轮数k=5x(5)=[0.,0.],f(5)=0.迭代精度：0.迭代轮数k=6x(6)=[1.,-0.],f(6)=1.迭代精度：0.迭代轮数k=7x(7)=[1.,0.],f(7)=1.迭代精度：0.迭代轮数k=8x(8)=[1.,-0.],f(8)=1.迭代精度：0.迭代轮数k=9x(9)=[1.,-0.],f(9)=1.迭代精度：0.*********************外点惩罚函数法优化最优点及目标函数值为：x(*)=[1.,-0.],f(*)=1.迭代精度：0.罚函数法考核题外点：，其初始

4、迭代点为，中止误差内点：，其初始迭代点为，中止误差**********外点惩罚函数法计算结果********************无约束优化方法：鲍威尔法**********++++++一维搜索方法：黄金分割法++++++**********初始惩罚因子：r0=1.00********************递增系数：c0=10.00**********初始坐标：x(0)=[2.,2.],f(0)=18.迭代轮数k=1x(1)=[1.,0.],f(1)=0.迭代精度：1.迭代轮数k=2x(2)=[0.,0.],f(2)=1.迭代精度：0.迭代

5、轮数k=3x(3)=[0.,0.],f(3)=1.迭代精度：0.迭代轮数k=4x(4)=[0.,0.],f(4)=1.迭代精度：0.迭代轮数k=5x(5)=[0.,0.],f(5)=1.迭代精度：0.迭代轮数k=6x(6)=[0.,0.],f(6)=1.迭代精度：0.迭代轮数k=7x(7)=[0.,0.],f(7)=1.迭代精度：0.迭代轮数k=8x(8)=[0.,0.],f(8)=1.迭代精度：0.迭代轮数k=9x(9)=[0.,0.],f(9)=1.迭代精度：0.迭代轮数k=10x(10)=[0.,0.],f(10)=1.迭代精度：0.**

6、*******************外点惩罚函数法优化最优点及目标函数值为：x(*)=[0.,0.],f(*)=1.迭代精度：0.算法程序实现//*鲍威尔法*/#include#include#include#include#include#defineN2/*优化设计维数*/#defineEPSIN0./*迭代精度*/#defineH_QJ1.0/*初始区间搜索步长*/#defineY_F1/*一维搜索方法选择：1——黄金分割法*//*2——二次插值法

7、*/#defineMC10.0/*惩罚因子递增系数*/#defineNG1/*不等式约束个数*/#defineNH1/*等式约束个数*/FILE*fp;doublemr=1.0;/*外点惩罚因子*/charoutname[50]="外点惩罚函数法计算结果.txt";/*计算结果输出文件*//*给出初始点坐标*/voidcsd_x(doublex0[]){inti;for(i=0;i

8、x[0]*x[0]+x[1]*x[1];returnf;}/*不等式约束方程*/voidstrain(doublex[],doubleg[],doub