| 查看: 1202 | 回复: 0 | ||
| 【悬赏金币】回答本帖问题，作者翻身仗将赠送您 5 个金币 | ||
翻身仗新虫 (正式写手)
|
[求助]
基于拉格朗日乘数法求最优化的这段代码为何是死循环
|
|
|
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); 上面这段代码是基于拉格朗日乘数法做最优化的，不知道什么地方出问题了，一运行就是死循环，请高人指教  捕获.JPG |
回复此楼» 猜你喜欢
求助
已经有3人回复
-
SCI 评审意见回复 模板
已经有4人回复
-
国家重点研发计划论文可以发MDPI里的期刊吗
已经有4人回复
-
CCS Chemistry投稿求助
已经有3人回复
-
国自然青年基金只有一篇能第一标注，其他都是第二标注能结题吗？
已经有6人回复
-
兄弟们帮我看看这两个结构怎么合成
已经有5人回复
-
一审一个审稿人，小修，会怎么样呀？
已经有6人回复
-
院士建议：对勇于提出解决卡脖子问题新解决方案的学者加大资助力度
已经有7人回复
-
没.赶.上.今.年.博.士.的.同.学.看.过.来.
已经有10人回复
-
美国教授开除中国研究生：我就不该录取你！
已经有6人回复
10
