%創(chuàng)建人：???常文靜
%?功能：????當(dāng)未知數(shù)有n個(gè)時(shí)，用罰函數(shù)法求最優(yōu)解
%???????????本程序中的目標(biāo)函數(shù)是f=（x1-0.2）^2+（x2-0.2）^2+...+（xn-0.2）^2約束條件為x1+x2+...
%???????????+xn=3
%?參數(shù)：???
%?返回值：???
%?主要思路：?
%?備注：????
%?調(diào)用方法：?
%?日期：?????2015/4/23?
%================================================================

%定義要求的未知變量x1x2....xn罰因子M
n=20;
for?i=1:n
????syms?（[‘x‘num2str（i）]）;
end
syms?M?positive;

x?=?[];
for?i=1:n
????x?=?[xsym（[‘x‘num2str（i）]）];????%將未知變量x1x2....xn存到數(shù)組x中
end


m（1）=10;????????????????????%c為遞增系數(shù)。賦初值。
c=8;
d=0;
Q=3;%需要滿足的約束條件

for?i=1:n
????d=d+（x（i）-0.2）.^2;???%目標(biāo)函數(shù)
end
f=d+M*（sum（x）-Q）.^2;????%增加罰函數(shù)的目標(biāo)函數(shù)

f0（1）=5;???%f0為所求最優(yōu)解下的f的值。賦初值5。

%求偏導(dǎo)、海森矩陣
for?i=1:n
????fx（i）=diff（fx（i））;??%f的一階偏導(dǎo)
end

for?i=1:n
????for?j=1:n
????????fxx（ij）=diff（fx（i）x（j））;??%f的二階