| 查看: 2255 | 回复: 8 | ||
老花木虫 (著名写手)
老花
|
[求助]
帮忙给看看这个matlab优化函数 问题
|
|
主： A=[];b=[];Aeq=[];Beq=[];Lb=[];Ub=[]; [x,f]=fmincon('myobj',[0;0],[],[],[],[],[],[],'mycon'); 目标函数：function y=myobj(x) y=(-1)*8*x(2)*(1-x(1))*x(3); 非线性约束function[c,ceq]=mycon(x) cep=[x1*(1-x(1))*x(3)-x(2)*(1+x(2));2*acos(exp(-x(4)))-x(3);1.5*2/sin(x(5))-x(4);atan((1-x(1))/(1+x(2)))]; 请帮忙看看程序的错误，谢谢！ |
回复此楼» 猜你喜欢
有多少人是今天查系统知道结果的？
已经有17人回复
-
国自然面上复盘~欢迎讨论
已经有9人回复
-
看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉……
已经有6人回复
-
导师吐槽：我怎么摊上了这么个极品研究生！
已经有6人回复
-
青基如何破局
已经有4人回复
-
基金不中，共勉
已经有8人回复
-
基金未中，这种答复是模板吗？
已经有5人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有5人回复
-
怎么看青基中了没有啊
已经有5人回复
-
梦想
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- matlab程序问题，总是出错，麻烦大家看看问题在那里？
已经有6人回复
- matlab的函数调用问题
已经有3人回复
- 帮忙看看这个平面图用matlab怎么画
已经有19人回复
- 问个关于matlab里find函数的问题！
已经有4人回复
- 求回归方程参数估计的最优化算法matlab代码
已经有9人回复
- 求助有关Matlab有约束非线性最优化问题
已经有9人回复
- 帮帮忙看看我的难题
已经有15人回复
- 【求助】matlab 求解微分方程中的未知参数
已经有20人回复
- matlab做函数图象的问题
已经有6人回复
- 【求助】matlab 二次规划的优化的问题
已经有4人回复
starlee
金虫 (正式写手)
- 应助: 6 (幼儿园)
- 金币: 1275.9
- 散金: 50
- 红花: 1
- 帖子: 310
- 在线: 627.3小时
- 虫号: 572146
- 注册: 2008-06-11
- 性别: GG
- 专业: CE
2楼2012-06-28 21:27:48
老花
木虫 (著名写手)
老花
- 应助: 0 (幼儿园)
- 贵宾: 0.402
- 金币: 2974.7
- 散金: 6205
- 红花: 41
- 沙发: 4
- 帖子: 2851
- 在线: 930.8小时
- 虫号: 1170496
- 注册: 2010-12-15
- 性别: GG
- 专业: 机械动力学
3楼2012-06-28 21:46:31
starlee
金虫 (正式写手)
- 应助: 6 (幼儿园)
- 金币: 1275.9
- 散金: 50
- 红花: 1
- 帖子: 310
- 在线: 627.3小时
- 虫号: 572146
- 注册: 2008-06-11
- 性别: GG
- 专业: CE
4楼2012-06-28 22:09:02
老花
木虫 (著名写手)
老花
- 应助: 0 (幼儿园)
- 贵宾: 0.402
- 金币: 2974.7
- 散金: 6205
- 红花: 41
- 沙发: 4
- 帖子: 2851
- 在线: 930.8小时
- 虫号: 1170496
- 注册: 2010-12-15
- 性别: GG
- 专业: 机械动力学
|
function[c,ceq]=mycon(x) cep=[x1*(1-x(1))*2/pi*acos(exp(-1.5*2/sin(x(3))))-x(2)*(1+x(2));atan((1-x(1))/(1+x(2)))-x(3)]; function y=myobj(x) y=(-1)*8*x(2)*(1-x(1))*2/pi*acos(exp(-1.5*2/sin(x(3)))) ; 修改后还是这样！ 请大家指点 Error in myobj (line 2) y=(-1)*8*x(2)*(1-x(1))*2/pi*acos(exp(-1.5*2/sin(x(3)))) ; Error in fmincon (line 601) initVals.f = feval(funfcn{3},X,varargin{:}); Caused by: Failure in initial user-supplied objective function evaluation. FMINCON cannot continue. |
5楼2012-06-28 22:33:55
老花
木虫 (著名写手)
老花
- 应助: 0 (幼儿园)
- 贵宾: 0.402
- 金币: 2974.7
- 散金: 6205
- 红花: 41
- 沙发: 4
- 帖子: 2851
- 在线: 930.8小时
- 虫号: 1170496
- 注册: 2010-12-15
- 性别: GG
- 专业: 机械动力学
6楼2012-06-28 22:46:44
starlee
金虫 (正式写手)
- 应助: 6 (幼儿园)
- 金币: 1275.9
- 散金: 50
- 红花: 1
- 帖子: 310
- 在线: 627.3小时
- 虫号: 572146
- 注册: 2008-06-11
- 性别: GG
- 专业: CE
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感谢参与，应助指数 +1
老花: 金币+100, ★★★★★最佳答案, 结果有点偏差，我再研究一下，谢谢啊！ 2012-06-28 23:12:22
老花: 回帖置顶 2012-06-28 23:24:42
感谢参与，应助指数 +1
老花: 金币+100, ★★★★★最佳答案, 结果有点偏差，我再研究一下，谢谢啊！ 2012-06-28 23:12:22
老花: 回帖置顶 2012-06-28 23:24:42
|
A=[];B=[];Aeq=[];Beq=[];Lb=[];Ub=[]; [x,fval]=fmincon('myobj',rand(1,3),A,B,Aeq,Beq,Lb,Ub,'mycon') function f=myobj(x) f=-1*8*x(2)*(1-x(1))*2/pi*acos(exp(-1.5*2/sin(x(3)))); function [c,ceq]=mycon(x) c=0; ceq(1)=x(1)*(1-x(1))*2/pi*acos(exp(-1.5*2/sin(x(3))))-x(2)*(1+x(2)); ceq(2)=atan((1-x(1))/(1+x(2)))-x(3); 结果： x = 0.3187 0.1832 0.5224 fval = -0.9971 结果合不合理？ |
7楼2012-06-28 22:59:10
老花
木虫 (著名写手)
老花
- 应助: 0 (幼儿园)
- 贵宾: 0.402
- 金币: 2974.7
- 散金: 6205
- 红花: 41
- 沙发: 4
- 帖子: 2851
- 在线: 930.8小时
- 虫号: 1170496
- 注册: 2010-12-15
- 性别: GG
- 专业: 机械动力学
8楼2012-06-28 23:24:31
starlee
金虫 (正式写手)
- 应助: 6 (幼儿园)
- 金币: 1275.9
- 散金: 50
- 红花: 1
- 帖子: 310
- 在线: 627.3小时
- 虫号: 572146
- 注册: 2008-06-11
- 性别: GG
- 专业: CE
9楼2012-06-29 12:37:21
