fahanshu.m

罚函数法matlab程序.希望对大家有用.

%罚函数法 
%目标函数：fun=x1^2+2*x2^2;
%约束条件：x1+x2=1
%main

clc
clear
syms x1 x2 k
fun0=x1^2+2*x2^2
fun1=x1^2+2*x2^2+k*(x1+x2-1)^2
f0=[diff(fun1,x1);diff(fun1,x2)];
[x1,x2]=solve(f0(1),f0(2))
x1=limit(x1,k,inf);
x2=limit(x2,k,inf);
zuiyoujie=[x1 x2]


fprintf('罚函数的最优解是：\n');
zuiyoujie