基于拉格朗日乘数法求最优化的这段代码为何是死循环

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

clc
clear
ta = 3.1;
cv = 0.8;
syms up vp  L   % up vp是边缘分布概率值，sy则是copula函数生成元中的变量。
mgb = exp( -( (-log(up))^ta + (-log(vp))^ta )^( 1/ta ) ); % mgb为 Gumbel Copula函数表达式
mi0 = diff(mgb, up);
f = -1 * diff(mi0,  vp); % copula函数的概率密度函数
h = mgb - cv;
F = f + L .* h;
s1 = diff(F, 'up');   
s2 = diff(F, 'vp');   
s3 = diff(F, 'L');   
[L,up,vp] = solve(s1 ,s2, s3);

上面这段代码是基于拉格朗日乘数法做最优化的，不知道什么地方出问题了，一运行就是死循环，请高人指教

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：拉格朗日乘数法 拉格朗日乘数 拉格朗日 拉格朗 最优化