matlab求解惩罚函数

哪位能帮我用matlab编程外部惩罚函数法,求解下列非线性规划。
不能把我这个方程直接算个解给我看下吗?

第1个回答  2010-02-10
penalty=1e6;

A=[16,2,3,13;5,11,10,8;9,7,6,12;4,14,15,10];
b=[81;89;89;117];
x=A\b

%原始解是x=[1;2;3;4];
%下面限定x3=1
k=3; %代表x3
xk=1;%代表x3=1;
temp=zeros(1,length(b));
temp(k)=penalty;

%构造新的矩阵方程,当然也可以直接在A,b里修改
AA=A;
bb=b;
AA(k,:)=temp;
bb(k)=xk*penalty;
xx=AA\bb
相似回答