| 查看: 2276 | 回复: 7 | ||
[求助]
使用拉格朗日优化，但Matlab 运行不出结果来
|
各位高手，我最近使用Matlab中的solve求解时，始终运行不出结果来，望高手给小弟帮忙解答，万分感谢。 Lagrange.png [ Last edited by 世纪乘风 on 2013-7-1 at 11:13 ] |
回复此楼» 猜你喜欢
filecode，4个jtjc了
已经有9人回复
-
今天基金会出结果吗？20260819
已经有8人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有115人回复
-
2027广东省杰青
已经有10人回复
-
朋友圈看到的
已经有8人回复
-
时间戳变了，能看出什么问题？
已经有14人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有30人回复
-
今天放榜吗？
已经有11人回复
-
快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。
已经有7人回复
-
什么时候开奖？
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 3D智能数字化与3D打印的一本新书
已经有12人回复
- LS-DYNA DYTRAN ABAQUS ADINA ANSYS NASTRAN 有限元软件比较
已经有0人回复
- 【小卒个人文集】写给初学者:小议形式电荷对量化初学者的误导
已经有27人回复
- 【资源】SIMULIA ABAQUS 6.7
已经有7人回复
aaron1988
木虫 (正式写手)
- 应助: 15 (小学生)
- 金币: 2261.1
- 红花: 3
- 帖子: 938
- 在线: 135.2小时
- 虫号: 413419
- 注册: 2007-06-26
- 专业: 化工系统工程
2楼2013-07-01 11:49:48
3楼2013-07-01 13:42:29
4楼2013-07-01 22:39:51
aaron1988
木虫 (正式写手)
- 应助: 15 (小学生)
- 金币: 2261.1
- 红花: 3
- 帖子: 938
- 在线: 135.2小时
- 虫号: 413419
- 注册: 2007-06-26
- 专业: 化工系统工程
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
世纪乘风: 金币+30, ★★★很有帮助 2013-07-02 12:50:00
感谢参与，应助指数 +1
世纪乘风: 金币+30, ★★★很有帮助 2013-07-02 12:50:00
|
syms P M r %Parameters W=21; C=10; a=500; b=10; Ss=40; Sr=20; hrs=3; Ir=0.13; Is=0.1; A=800; K=3; P1 = hrs+W*Ir; P2 = r*b*W*Ir/C/Is; % LM: -(a-bP)CIs+rbWIr = 0 LP=-b*(W-C)+b*Ss*sqrt(P1/72/Sr/P2)+b*hrs*(2-P2/A)*sqrt(Sr/8/P1/P2)-... hrs*b/A*sqrt(P2*Sr/2/P1)+b*M*C*Is-r*(-2*b+b^2/4/P2^2)*sqrt(2*Sr*P1*P2); Constraint = r*b*W*Ir/C/Is*2-a+b*W-b*M*W*Ir+b*sqrt(Sr*P1/2/P2); f = subs(LP,'M',solve(Constraint,M)); rvalue = solve(f); Pvalue = a/b-rvalue*W*Ir/C/Is 这题目居然有解析解。 你试下这个吧。 我是把LP的变量都统一了（用约束条件=0跟LM=0） 如果你硬要保持LP，LM，Constraint这样的格式的话（约束可能是不等式）。 那么请参考动态规划，我实在是懒得编程用个比较方便的 （这个问题是没办法做成线性规划的，因为有非线性元素的存在，不然的话，直接解线性矩阵就好了）。 （下回麻烦你把MATLAB的code也贴上来吧，码code好累的） 结果： P值是50 +/- 0.7 (这是非线性问题，并不保证存在唯一解）。 |
5楼2013-07-02 11:16:48
|
谢谢大侠的耐心解答，真的很谢谢！ 我之前也通过约束条件把M用P表达出来，然后带入目标函数，进而优化无约束的目标函数，可以得出结果来的。但对于约束条件中含有M平方或立方时，我就无能为力了，所以我就想到了拉格朗日方法，进而有LP=0,LM=0，Constraint=0, 联立求解的想法，呵呵。 给你添麻烦了，谢谢大侠。 另，不知 Matlab中 可有函数直接对带有约束条件的函数进行优化的函数 |
» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : Matlabcode.docx
2013-07-02 19:22:09, 12.24 K
6楼2013-07-02 19:22:24
aaron1988
木虫 (正式写手)
- 应助: 15 (小学生)
- 金币: 2261.1
- 红花: 3
- 帖子: 938
- 在线: 135.2小时
- 虫号: 413419
- 注册: 2007-06-26
- 专业: 化工系统工程
|
1. （一般）参考动态规划。 你找一下Dimitri P.Bertsekas的dynamic programming。这个只能自己读懂算法，然后写程序了。或者用Matlab自带的Optimization优化工具箱。你只要放进目标方程， 约束条件，还有参数,然后Matlab就能根据你的信息来给出数值结果（不是解析解！！，只是一个接近的解）。 2. 特殊情况是目标方程和约束条件都是quadratic的话（e.g x'*A*x = b ), matlab有quadprog 3. 如果你能简化成线性系统的话， matlab有linprog |
» 本帖已获得的红花（最新10朵）
世纪乘风7楼2013-07-03 09:23:25
8楼2013-07-03 20:24:14
