MATLAB程序 惩罚函数法 内点法
合集下载
相关主题
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
%主调函数
%fun6_2_1.m
global r
r0=1;
c=0.1;
km=10;
x0=[2;4];
for k=1:km;
r=r0*c^(k-1);
[x,fv]=fminsearch(@fun,x0);
fprintf(1,'µÚ%3.0f´Îµü´ú\n',k) fmin=fv
x0=x
end
%功能函数
function f=fun(x)
global r
f=x(1)^2+2*x(2)^2+r/(x(1)+x(2)-1);