| 查看: 1099 | 回复: 11 | ||
[求助]
求助：帮忙做个matlab优化计算
|
| 大家帮忙给计算一下！ |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : matlab优化a.doc
2012-04-25 09:41:54, 73.5 K
» 猜你喜欢
垃圾破二本职称评审标准
已经有18人回复
-
职称评审没过，求安慰
已经有53人回复
-
毕业后当辅导员了，天天各种学生超烦
已经有5人回复
-
26申博自荐
已经有3人回复
-
A期刊撤稿
已经有4人回复
-
EST投稿状态问题
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助啊......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 (大学生)
- 金币: 18179.7
- 散金: 780
- 红花: 5
- 沙发: 2
- 帖子: 4555
- 在线: 590.4小时
- 虫号: 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 (大学生)
- 金币: 18179.7
- 散金: 780
- 红花: 5
- 沙发: 2
- 帖子: 4555
- 在线: 590.4小时
- 虫号: 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 (大学生)
- 金币: 18179.7
- 散金: 780
- 红花: 5
- 沙发: 2
- 帖子: 4555
- 在线: 590.4小时
- 虫号: 784367
- 注册: 2009-06-01
- 性别: GG
- 专业: 控制理论与方法
9楼2012-04-27 12:25:40
10楼2012-04-28 08:28:48