function f=obj(x)
f=(x(1)-5)^2+(x(2)-6)^2+7;
function fai=objfunout(x)
%辅助函数
global r
fai=obj(x)+r*([2x1-x2-1;0])^2+r([x1+x2-2;0])^2;
function y=fun2pen(x)
%罚项函数
y=(max(0,1-x(1)))^2;
function penitentout
% 外点法罚函数方法
% fun2pen 约束函数
% objfun 目标函数
clc
global r
x0=[0;1];
r=0.3; c=5; ep=1e-5; k=1;
while rfun2pen(x0)>=ep
[x0,f]=fminsearch(@objfunout,x0)
r=cr; k=k+1;
end
disp('最优解x0='),disp(x0)
disp('k='),disp(k)
disp('目标函数f(x)='),disp(obj(x0))
MATLAB外点惩罚函数法不知道那错了
- 写回答
- 好问题 0 提建议
- 追加酬金
- 关注问题
- 邀请回答
-
0条回答 默认 最新
悬赏问题
- ¥20 有关区间dp的问题求解
- ¥15 多电路系统共用电源的串扰问题
- ¥15 slam rangenet++配置
- ¥15 有没有研究水声通信方面的帮我改俩matlab代码
- ¥15 对于相关问题的求解与代码
- ¥15 ubuntu子系统密码忘记
- ¥15 信号傅里叶变换在matlab上遇到的小问题请求帮助
- ¥15 保护模式-系统加载-段寄存器
- ¥15 电脑桌面设定一个区域禁止鼠标操作
- ¥15 求NPF226060磁芯的详细资料