24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

世纪乘风

银虫 (小有名气)

应助: 0 (幼儿园)
金币: 1628.5
帖子: 58
在线: 86.2小时
虫号: 1651164
注册: 2012-02-28
性别: GG
专业: 工业工程与管理

[求助] 使用拉格朗日优化，但Matlab 运行不出结果来

各位高手，我最近使用Matlab中的solve求解时，始终运行不出结果来，望高手给小弟帮忙解答，万分感谢。

Lagrange.png

[ Last edited by 世纪乘风 on 2013-7-1 at 11:13 ]

回复此楼

» 猜你喜欢

请问有评职称，把科研教学业绩算分排序的高校吗已经有6人回复
2025冷门绝学什么时候出结果已经有6人回复
Bioresource Technology期刊，第一次返修的时候被退回好几次了已经有7人回复
真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗已经有8人回复
寻求一种能扛住强氧化性腐蚀性的容器密封件已经有5人回复
请问哪里可以有青B申请的本子可以借鉴一下。已经有4人回复
请问下大家为什么这个铃木偶联几乎不反应呢已经有5人回复
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入已经有4人回复
康复大学泰山学者周祺惠团队招收博士研究生已经有6人回复
AI论文写作工具：是科研加速器还是学术作弊器？已经有3人回复

» 本主题相关价值贴推荐，对您同样有帮助:

3D智能数字化与3D打印的一本新书已经有12人回复
LS-DYNA DYTRAN ABAQUS ADINA ANSYS NASTRAN 有限元软件比较已经有0人回复
【小卒个人文集】写给初学者:小议形式电荷对量化初学者的误导已经有27人回复
【资源】SIMULIA ABAQUS 6.7 已经有7人回复

1楼2013-07-01 11:00:08

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

aaron1988

木虫 (正式写手)

应助: 15 (小学生)
金币: 2261.1
红花: 3
帖子: 938
在线: 135.2小时
虫号: 413419
注册: 2007-06-26
专业: 化工系统工程

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +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 (这是非线性问题，并不保证存在唯一解）。