《2022年等式约束极值问题-外点罚函数法 .pdf》由会员分享，可在线阅读，更多相关《2022年等式约束极值问题-外点罚函数法 .pdf（4页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、重庆科技学院学生实验报告课程名称最优化方法与应用实验项目名称等式约束极值问题 -外点罚函数法开课学院及实验室实验日期学生姓名学号专业班级指导教师实验成绩一、实验目的和要求1. 熟悉外点罚函数法的思想和步骤；2.掌握外点罚函数法求解等式约束的多维极值问题。二、实验内容和原理1. 对所给等式约束的多维极值问题用matlab编程进行求解；三、主要仪器设备操作系统为 Windows 2000 及以上的电脑，并装有matlab 软件四、实验操作方法和步骤根据外点罚函数法 的步骤，读懂下列的程序，并用该程序求解所给多维函数极值。名师资料总结 - - -精品资料欢迎下载 - - - - - - - - -
2、- - - - - - - - - 名师精心整理 - - - - - - - 第 1 页，共 4 页 - - - - - - - - - 五、实验记录与处理（数据、图表、计算等） syms x y; minGeneralPF(x2+y2,1,1,y2-1,1000,2,x,y,0.0001) ans = 0 0 Minf=0 六、实验结果及分析运用所给程序能够求出最优解X=（0,0 ）最优值 0，名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 2 页，共 4 页 - - - - -
3、- - - - 附录 function x,minf = minGeneralPF(f,x0,h,c1,p,var,eps) format long; if nargin = 6 eps = 1.0e-4; end k = 0; FE = 0; for i=1:length(h) FE = FE + (h(i)2; end x1 = transpose(x0); x2 = inf; while 1 M = c1*p; FF = M*FE; SumF = f + FF; x2,minf = minNT(SumF,transpose(x1),var); if norm(x2 - x1)eps v
4、= subs(gradf,var,x0); tol = norm(v); pv = subs(jacf,var,x0); p = -inv(pv)*transpose(v); p = double(p); x1 = x0 + p; x0 = x1; end x = x1; minf = subs(f,var,x); format short; syms x y; minGeneralPF(x2+y2,1,1,y2-1,1000,10,x,y,0.0001) ans = 0 1.0000 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 4 页，共 4 页 - - - - - - - - -