| 查看: 1413 | 回复: 11 | ||
[求助]
求助：帮忙做个matlab优化计算
|
| 大家帮忙给计算一下！ |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : matlab优化a.doc
2012-04-25 09:41:54, 73.5 K
» 猜你喜欢
今天系统多次维护，明天很可能放榜！
已经有6人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有10人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有3人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有74人回复
-
岩土工程学报什么时候才能终审完呐
已经有4人回复
-
时间戳又变了8-15
已经有25人回复
-
2027广东省杰青
已经有6人回复
-
快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。
已经有3人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有26人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有16人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助啊......matlab
已经有10人回复
- 最优化计算求助！（1stopt或者Matlab）
已经有16人回复
- 求回归方程参数估计的最优化算法matlab代码
已经有9人回复
- 求助有关Matlab有约束非线性最优化问题
已经有9人回复
- 求助高手帮忙用matlab接下面这个算式，具体见附件
已经有5人回复
- 【求助】matlab 计算一个二重积分，高手们帮帮忙吧
已经有6人回复
- 【求助】用matlab最优化方法进行参数拟合
已经有17人回复
- 【求助】催化反应动力学matlab计算各基元反应的速率常数时，该如何避免较小量被忽略？
已经有3人回复
- 【求助】请教一个matlab程序，谢谢帮忙！！
已经有15人回复
- 【求助】matlab 二次规划的优化的问题
已经有4人回复
- 【求助】精通matlab最优化计算，一书中的funval函数是什么意思
已经有10人回复
2楼2012-04-25 10:40:53
nandehutu9327
至尊木虫 (职业作家)
- 应助: 279 (大学生)
- 金币: 18965.7
- 散金: 780
- 红花: 5
- 沙发: 2
- 帖子: 4561
- 在线: 591.1小时
- 虫号: 784367
- 注册: 2009-06-01
- 性别: GG
- 专业: 控制理论与方法
3楼2012-04-25 12:06:15
|
目标函数： minf=√(〖x1〗^2+〖x2〗^2+〖x3〗^2+〖x4〗^2 )-√(〖(x1-x3)〗^2+〖(x2-x4)〗^2 ) 约束条件 sub.to X0=(177 517 86 95) %初始值 49.876 (（〖x3〗^2+〖x4〗^2)(x1*x4-x2*x3))/√((〖x3〗^2+〖x4〗^2 )[(〖x1+x3)〗^2+〖(x2+x4)〗^2])-34.482(（〖x3〗^2+〖x4〗^2)(x1*x4-x2*x3))/√((〖x3〗^2+〖x4〗^2 )[(〖x1-x3)〗^2+〖(x2-x4)〗^2])-+140≤0 X1＞0，x2>0，x3>0，x4>0 编写的m文件 function f =myfun(x) %目标函数 a=sqrt(x(1)^2+x(2)^2+x(3)^2+x(4)^2) b=sqrt((x(1)-x(3))^2+(x(2)-x(3))^2) f=a-b end function [c,ceq]=mycon(x) %约束条件 a=(x(1)*x(4)-x(2)*x(3))*sqrt(x(3)^2+x(4)2) b=(x(1)-x(3))^2+(x(2)-x(4))^2 d=(x(1)+x(3))^2+(x(2)+x(4))^2 c=49.876*a/sqrt(d) -34.482*a/sqrt(b)+140 ceq=[] end 命令行： >> x0=[177 517 86 95]; >> A=[-1 0 0 0;0 -1 0 0;0 0 -1 0;0 0 0 -1]; >> b=zeros(4,1); >> [x,fval,exitflag,output,lambda,grad,hessian] =fmincon('myfun,x0,A,b,[],[],[],[]','mycon') ??? Undefined function or variable 'x'. |
4楼2012-04-26 09:26:54
5楼2012-04-26 09:27:37
nandehutu9327
至尊木虫 (职业作家)
- 应助: 279 (大学生)
- 金币: 18965.7
- 散金: 780
- 红花: 5
- 沙发: 2
- 帖子: 4561
- 在线: 591.1小时
- 虫号: 784367
- 注册: 2009-06-01
- 性别: GG
- 专业: 控制理论与方法
6楼2012-04-26 13:07:41
|
编写的m文件 function f =myfun(x) %目标函数 a=sqrt(x(1)^2+x(2)^2+x(3)^2+x(4)^2) b=sqrt((x(1)-x(3))^2+(x(2)-x(3))^2) f=a-b end function [c,ceq]=mycon(x) %约束条件 a=(x(1)*x(4)-x(2)*x(3))*sqrt(x(3)^2+x(4)2) b=(x(1)-x(3))^2+(x(2)-x(4))^2 d=(x(1)+x(3))^2+(x(2)+x(4))^2 c=49.876*a/sqrt(d) -34.482*a/sqrt(b)+140 ceq=[] end 命令行： >> x0=[177 517 86 95]; A=[-1 0 0 0;0 -1 0 0;0 0 -1 0;0 0 0 -1]; b=zeros(4,1); >> x,fval,exitflag,output,lambda,grad,hessian]=fmincon('myfun’,x0,A,b,[],[],[],[],'mycon') 改过以后结果也不对。x = 0.0000 614.4217 3.0157 0.0000 fval = 3.0157 把x代入f是0. |
7楼2012-04-27 08:35:18
8楼2012-04-27 10:18:27
nandehutu9327
至尊木虫 (职业作家)
- 应助: 279 (大学生)
- 金币: 18965.7
- 散金: 780
- 红花: 5
- 沙发: 2
- 帖子: 4561
- 在线: 591.1小时
- 虫号: 784367
- 注册: 2009-06-01
- 性别: GG
- 专业: 控制理论与方法
9楼2012-04-27 12:25:40
10楼2012-04-28 08:28:48
25