| 查看: 2282 | 回复: 7 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
使用拉格朗日优化，但Matlab 运行不出结果来
|
||
各位高手，我最近使用Matlab中的solve求解时，始终运行不出结果来，望高手给小弟帮忙解答，万分感谢。 Lagrange.png [ Last edited by 世纪乘风 on 2013-7-1 at 11:13 ] |
回复此楼» 猜你喜欢
欢迎发来filecode的Mz6后的代码验证其规律
已经有123人回复
-
今天基金会出结果吗？20260819
已经有13人回复
-
时间戳变了，能看出什么问题？
已经有16人回复
-
filecode，4个jtjc了
已经有13人回复
-
为什么余额宝的年化利率越来越低？主要原因有哪些？
已经有5人回复
-
今天放榜没戏了吧
已经有5人回复
-
今天放榜吗？
已经有15人回复
-
filecode=后面第一个是大写字母
已经有12人回复
-
2027广东省杰青
已经有10人回复
-
朋友圈看到的
已经有8人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 3D智能数字化与3D打印的一本新书
已经有12人回复
- LS-DYNA DYTRAN ABAQUS ADINA ANSYS NASTRAN 有限元软件比较
已经有0人回复
- 【小卒个人文集】写给初学者:小议形式电荷对量化初学者的误导
已经有27人回复
- 【资源】SIMULIA ABAQUS 6.7
已经有7人回复
3楼2013-07-01 13:42:29
aaron1988
木虫 (正式写手)
- 应助: 15 (小学生)
- 金币: 2261.1
- 红花: 3
- 帖子: 938
- 在线: 135.2小时
- 虫号: 413419
- 注册: 2007-06-26
- 专业: 化工系统工程
2楼2013-07-01 11:49:48
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
