function?[x?delta]??=?BFGS_exact（f_originClambdamu）
%?BFGS法確定多個變量的最優解
%

%?輸入：
%?x0：當前點；
%?Iter_max：最大迭代次數；

%?輸出：
%?x：最優解；
%

Iter_max?=?20;?%函數參數的個數，x1x2x3...的個數
x0?=?[11];?%初始值
?eps?=?1e-2;?%迭代容忍度
ite?=?1;?%?循環次數
f111?=?f_origin-lambda.*C+mu.*（C）^2/2;

g111=gradient（f111）;
gk?=?BFGS_Gradient（g111x0）;?%?初始梯度
H?=?eye（length（x0））;?%?初始H陣
d?=?-gk*H‘;?%?初始搜索方向
alpha?=?exact_buchang（f111x0d‘）;?%?初始搜索步長
while?norm（gk）?>?eps?&&?ite
x0?=?x0?+?alpha.*d;?%?新的循環點
gkk?=?BFGS_Gradient（g111x0）;?%?新的梯度
y?=?gkk?-?gk;?%?梯度差
s?=?alpha.*d;?%?迭代點差
H?=?H?+?（1?+（y*H*y‘）/（s*y‘））*（s‘*s）/（s*y‘）?-?（s‘*y*H?+?H*y‘*s）/（s*y‘）;?%H陣更新
d?=?-gkk*H‘;?%?新的搜索方向
alpha?=?exact_buchang（f111x0d‘）;?%?新的搜索步長
gk?=?gkk;?%?更新梯度
ite?=?ite?+?1;?%?循環次數加1
end
x?=?x0?+?alpha.*d;?%?極小值
delta=BFGS_Gradient（g111x）;
?
%?梯度矩陣
function?g?=?BFGS_Gradient（g11x0）
temp?=?sym（‘x‘[12]）;
g=（double（subs（g11tempx0）））‘;
end
?
end

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????1075??2020-04-14?22:51??第五題增廣拉格朗日法\BFGS_exact.m
?????文件?????????334??2020-04-14?22:47??第五題增廣拉格朗日法\exact_buchang.m
?????文件?????????493??2020-04-14?23:16??第五題增廣拉格朗日法\main.m