




版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、精选优质文档-倾情为你奉上拉格朗日乘子法约束优化问题的标准形式为:约束优化算法的基本思想是:通过引入效用函数的方法将约束优化问题转换为无约束问题,再利用优化迭代过程不断地更新效用函数,以使得算法收敛。1. 罚函数法罚函数法(内点法)的主思想是:在可行域的边界上筑起一道很高的“围墙”,当迭代点靠近边界时,目标函数陡然增大,以示惩罚,阻止迭代点穿越边界,这样就可以将最优解“挡”在可行域之内了。它只适用于不等式约束:它的可行域为:对上述约束问题,其其可行域的内点可行集的情况下,引入效用函数:、其中或算法的具体步骤如下:给定控制误差,惩罚因子的缩小系数。步骤1:令,选定初始点,给定(一般取10)。步骤
2、2:以为初始点,求解无约束其中或,得最优解 步骤3:若,则为其近似最优解,停;否则,令,转步骤2.2. 拉格朗日乘子法(1)算法:(约数为等式的情况引入)效用函数为判断函数为当时迭代停止。步骤1:选定初始点,初始拉格朗日乘子向量,初始罚因子及其放大系数,控制误差与常数,令。步骤2:以为初始点,求解无约束问题:得到无约束问题最优解步骤3:当时,为所求的最优解,停;否则转步骤4.步骤4:当时,转步骤5;否则令,转步骤5.步骤5:令,转步骤1。(2) 算法(一般约束形式的松弛变量法和指数形式法)松弛变量法:乘子的修正公式为:判断函数为:当时迭代停止。3. 乘子法MATLAB程序及其作用3.1 函数3
3、.1.1程序(1):乘子法效用函数程序 函数功能:将约束优化问题,根据效用函数方法,将其转变成无约束问题。function f=AL_obj(x)%拉格朗日增广函数%N_equ 等式约束个数%N_inequ 不等式约束个数global r_al pena N_equ N_inequ;%全局变量h_equ=0;h_inequ=0;h,g=constrains(x);%等式约束部分for i=1:N_equ h_equ=h_equ+h(i)*r_al(i)+(pena/2)*h(i).2;end%不等式约束部分for i=1:N_inequ h_inequ=h_inequ+(0.5/pena)*(
4、max(0,(r_al(i)+pena*g(i).2-r_al(i).2);end%拉格朗日增广函数值f=obj(x)+h_equ+h_inequ;3.1.2 程序(2):判断函数函数功能:判断是否符合约束条件% the compare function is the stop conditionfunction f=compare(x)global r_al pena N_equ N_inequ;h_equ=0;h_inequ=0;h,g=constrains(x);%等式部分for i=1:N_equ h_equ=h_equ+h(i).2;end%不等式部分for i=1:N_inequ
5、h_inequ=h_inequ+(max(-g(i),r_al(i+N_equ)/pena).2;endf=sqrt(h_equ+h_inequ);3.1.3 程序(3)AL算法主程序函数功能:对无约束的效用函数利用拟牛顿算法求解其最优解,更新乘子。function X,FVAL=AL_main(x_al,r_al,N_equ,N_inequ)%本程序为拉格朗日乘子算法示例算法%函数输入:% x_al:初始迭代点% r_al:初始拉格朗日乘子% N-equ:等式约束个数% N_inequ:不等式约束个数%函数输出% X:最优函数点% FVAL:最优函数值%=程序开始=global r_al p
6、ena N_equ N_inequ; %参数(全局变量)pena=10; %惩罚系数c_scale=2; %乘法系数乘数cta=0.5; %下降标准系数e_al=0.005; %误差控制范围max_itera=25;out_itera=1; %迭代次数%=算法迭代开始=while out_itera<max_itera x_al0=x_al; r_al0=r_al; %判断函数 compareFlag=compare(x_al0); %无约束的拟牛顿法BFGS X,FVAL=fminunc(AL_obj,x_al0); x_al=X; %得到新迭代点 %判断停止条件 if compare
7、(x_al)<e_al disp('we get the opt point'); break end %c判断函数下降度 if compare(x_al)<cta*compareFlag pena=pena; %可以根据需要修改惩罚系数变量 else pena=min(1000,c_scale*pena); %乘法系数最大1000 disp('pena=2*pena'); end % 更新拉格朗日乘子 h,g=constrains(x_al); for i=1:N_equ %等式约束部分 r_al(i)=r_al(i)+pena*h(i); end
8、 for i=1:N_inequ %不等式约束部分 r_al(i+N_equ)=max(0,(r_al(i+N_equ)+pena*g(i); end out_itera=out_itera+1;end%+迭代结束+disp('!the iteration over!');disp('the value of the obj function');obj(x_al)disp('the value of constrains');compare(x_al)disp('the opt point'); X=x_al; FVAL=obj(X);3.1.4 乘子法函数使用方法(1) 定义目标函数及约束条件目标函数文件约束函数文件(2) 函数调用x_al=1,1,1; %初始迭代点r_al=1,1; %初始拉格朗日乘子N_equ=1; %等式约束个数 一个N_inequ=1; %不等式约束个数 一个X,FVAL=AL_main(x_al,r_al,N_equ,N_inequ)计算结果:we get the opt point!the iteration over!the value of t
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 甘蔗产业环境影响因素与生态效益评估
- 山东省济宁市田家炳中学2024年八上物理期末综合测试试题含解析
- 陕西省岐山县联考2025届九年级化学第一学期期末综合测试模拟试题含解析
- 草种种植与生态旅游合作合同
- 大数据产业在各行业的实际应用与前景预测
- 网络信息安全领域未来布局
- 游戏关卡深度设计与玩家体验优化研究
- 2025至2030国内抗氧化食品行业项目调研及市场前景预测评估报告
- 2025至2030中国自助点菜亭行业产业运行态势及投资规划深度研究报告
- 2025至2030中国自动沉降分析仪行业发展状况及投资潜力研究报告
- 异位妊娠的课件
- 血管内超声IVUS简介
- DL∕T 2528-2022 电力储能基本术语
- 上海2024年上海市教育评估院招聘笔试上岸历年典型考题与考点剖析附带答案详解
- 渣土清运综合项目施工组织设计
- 苏教版八年级生物下册期末试卷及答案【苏教版】
- 书面检查材料(通用6篇)
- 传感器与机器视觉 课件 第六章 机器视觉
- 2019疏浚工程预算定额
- RFC2326(中文版+英文版+可链接目录)-RTSP
- 2023八年级历史下册第六单元科技文化与社会生活第18课科技文化成就说课稿新人教版
评论
0/150
提交评论