24小时热门版块排行榜    

查看: 1419  |  回复: 5
当前主题已经存档。
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

huhaifeng

木虫 (初入文坛)

[交流] 【求助】增广拉格朗日方法

想用增广拉格朗日方法求解一优化问题。
不过有些地方不懂,比如假设只有等式约束,
min f(x)
s.t. g(x)=0
增广拉格朗日函数P(x,lamda,r)=f(x)-lamda*g(x)+ r/2*g(x)^2;
我的问题是:
迭代出lamda后,如果用牛顿法求解x,是不是用拉格朗日函数的导数?
L(x,lamda)=f(x)-lamda*g(x)?
因为这样才满足kkt条件?
我看有的书上写的是通过设置lamda的迭代值,可以把P的导数等同于L的导数?
有点乱,看了一些文献,觉得数学性太强了,不懂。

[ Last edited by 小雨萌萌 on 2010-4-6 at 10:31 ]
回复此楼

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小雨萌萌

铜虫 (文坛精英)

优秀版主


★ ★
javeey(金币+2):谢谢解答 2010-04-06 20:07
当函数具有某些好的性质时,函数F和L的稳定点是一样的。你说的那个文献可以在网上免费下载的http://www.ebookee.net/Constrain ... Series-_212576.html
4楼2010-04-06 19:57:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答

小雨萌萌

铜虫 (文坛精英)

优秀版主


★ ★
小木虫(金币+0.5):恭喜抢沙发,给个红包
javeey(金币+1):谢谢这方面的专家提供帮助 2010-04-06 11:52
推荐一篇论文,基于增广Lagrange函数的等式约束优化算法,这个你看起来应该不难。
2楼2010-04-06 11:29:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huhaifeng

木虫 (初入文坛)

谢谢版主!
我看了一下你给的文章,他求解等式约束是Armijo搜索?这个我也不懂是什么,但是算法看了个大概。
我在袁亚湘的书里翻出来一个乘子罚函数的算法,p474-475,
里面提及的算法是:
1.设置lamda,r和x的初值
2.求解x(k+1)=arg min {P(x,lamda,r)}
3.迭代惩罚因子 r (提及的文献要避免的?为了防止r趋向于无穷大?)
4.迭代lamda
5.转 2
我在一文献里发现有人用增广lagrange算法,他用的是牛顿迭代求解第二步中的x*,但是他用到的雅克比矩阵和海森矩阵都是关于L的,而不是P,我总结了一下他的迭代是这样的:
1.设置lamda,和x的初值,r保持不变
2.迭代lamda
3.求解x(k+1)=arg min {L(x,lamda)}
4.转 2
我看袁亚湘的书,因为lamda的更新放在了第二步,所以按他的书的算法,min P变成了min L。
我只是不知道这样对吗?他本人引用的优化参考文献是 constrained optimization and lagrange multiplier methods,这个文献我查不到。
谢谢!
我不清楚这样的做法

[ Last edited by 小雨萌萌 on 2010-4-6 at 19:44 ]
3楼2010-04-06 15:05:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huhaifeng

木虫 (初入文坛)

谢谢版主!
不过可惜那个链接的文件已经被删除了
还有一个不太相关的问题:我在amazon上看到有人给这本书的评论是:
out of date... many new approaches (e.g., SQP, GRG, trust-region methods, interior point methods) have gained favor for their greater efficiency and robustness. Even so, Augmented Lagrangian methods remain relevant and interesting, in part because of their relationship to some newer approaches.
像我在一般情况下,求约束优化基本就会考虑到拉格朗日乘子法罚函数法之类的,求解非约束一般用牛顿法这样的。
但是实际上在选用某种优化算法是,是不是用评论中所说的新方法更可靠更好一些?不考虑比较简单的情况下。
5楼2010-04-07 10:28:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 让我中一个面上吧! +11 大萍1987 2026-08-20 13/650 2026-08-21 14:48 by 20120902066
[基金申请] 人气不行了 +7 fansofjerry 2026-08-21 7/350 2026-08-21 14:13 by Vivilian
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +15 医学老男孩 2026-08-20 18/900 2026-08-21 13:46 by 医学老男孩
[基金申请] 我面上完蛋了 +7 且听虎啸 2026-08-20 8/400 2026-08-21 12:31 by 酷酷墨镜
[基金申请] 只有每年这种时候来逛逛小木虫 +23 yaoyewhu2008 2026-08-20 25/1250 2026-08-21 12:19 by Poppy1104
[基金申请] 什么时候开奖? +7 CrisMessi 2026-08-18 7/350 2026-08-21 12:15 by weiyin
[基金申请] 科研孤儿太难了 +15 我4大白菜 2026-08-20 16/800 2026-08-21 11:07 by yming1319
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 时间戳今天,20号变了 +4 archvillain 2026-08-20 4/200 2026-08-21 08:10 by gatelove
[基金申请] 时间戳变了,能看出什么问题? +18 基诺咪客 2026-08-17 23/1150 2026-08-20 17:19 by Godzela
[基金申请] 重要消息,中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 2027广东省杰青 +4 奶牛小黑 2026-08-15 10/500 2026-08-19 11:02 by wanfengnew
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +6 Tide man 2026-08-14 7/350 2026-08-19 09:56 by ZJTJZ
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
信息提示
请填处理意见