| 查看: 356 | 回复: 1 | ||
[求助]
求救！！最好在一点之前，我这程序存在一点小问题，运行不了，求各位帮我改改
|
|
(1)编写目标函数M文件 function f=myfun(X) f=(pi* 110^2-pi X(1) )^2 )* 120 *0．25+0．25* Pi* X(1)^2 *1200+pi*pi *l／6* X(1)*X(2)^ 2／cos (atan(x(3)／(pi *x(1))))； (2)编写约束函数M文件 function[c，ceq]=mycon(X) e(1)=19601-410．9006* 3 ^0．7*X(1)^ (-0．69)*X(2)^ (2．1-2／3) *(X(1)-0．707*X (2)) ^1．39*(X(1)+0．707 *x(2)) ^(-1／3)； c(2)=162583-13．87* (0．03846*(1+(0．7071* x(2)／x(1))／(1-0．7071* X(2)／x(1))))^ (-0．5)* pi*X(1) *3*0．7071* X(2)； c(3)=0．03493 *pi* X(1)-X(3)； c(4)=X(3)-0．12278* pi* X(1)； c(5)=0．6* X(3)-X(2)； c(6)=X(2)-0．66* X(3)； ceq=x(4)- atan(X(3)／(pi *X(1)))* 180／pi； (3)在MATLAB命令窗口调用优化程序 xO=[30，3，6，3]； lb=[20，2，4，2]； ub=[60，15，18，7]； options=optimset(‘ display ’， ‘iter’， ‘ LargeScale’， ‘0ff’ )； [X，fval，exitflag，output]=fmincon(@myfun，x0， []，[]，[]，[]，lb，ub，@mycon，options) |
回复此楼» 猜你喜欢
关于代码变化问题，想知道的进来
已经有10人回复
-
基金中了
已经有13人回复
-
关于filecode
已经有6人回复
-
固定端突然变了，今天
已经有10人回复
-
听说今天filecode变了
已经有47人回复
-
filecode与中标关系的预测
已经有3人回复
-
filecode变化情况
已经有22人回复
-
十年后又回来了，论文投稿求助
已经有4人回复
-
大家散了吧，后缀研究没有意义，别浪费时间了，过好目前的每一天，不要焦虑
已经有7人回复
-
关于豆爷回答的JTJC与%2F数量
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 【求助】gaussian03运行问题 g03 <cu.gjf>cu.out&
已经有9人回复
- 【求助】高斯安装后多用户使用、多核运行及服务器配置问题
已经有15人回复
- 【求助】此程序也是没有问题了，但是依然运行不出结果
已经有7人回复
- 【求助】MS4.3在运行，却一直没有CASTEP这个程序在工作，计算结果也迟迟不出现。。。�
已经有15人回复
- 【求助】程序运行后出现的错误 access violation
已经有15人回复
- 【求助】安装问题，安装后程序无法执行
已经有12人回复
- 【求助】关于vector运行效率问题【已完结】
已经有5人回复
- 【求助】关于MBR运行中活性污泥颜色变黄且粘度过大的问题
已经有18人回复
- 紧急求救：我的纳米盘不能下载之国外下载问题
已经有8人回复
opensees520
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 647.2
- 红花: 1
- 帖子: 18
- 在线: 19.1小时
- 虫号: 1125976
- 注册: 2010-10-18
- 性别: GG
- 专业: 结构工程
【答案】应助回帖
★ ★ ★
臭水沟(金币+3): 谢谢应助~~ 2011-09-30 09:22:46
臭水沟(金币+3): 谢谢应助~~ 2011-09-30 09:22:46
|
1.x 大小請固定.不要一個大X一個小 2.l 與1 要區分 3.全形與半形字要區分 4.f=(pi* 110^2-pi X(1) )^2 ^^^^^^ 應該為pi*x(1) 修正上述問題後跑出?斫Y果: Max Line search Directional First-order Iter F-count f(x) constraint steplength derivative optimality Procedure 0 5 1.90425e+006 1.217e+005 Infeasible start point 1 10 3.22122e+006 1.325 1 4.29e+004 9.19e+005 2 15 2.75537e+006 0.1556 1 -8.08e+004 2.02e+005 Hessian modified twice 3 20 2.43364e+006 2.347e-009 1 -6.75e+004 1.8e+004 4 25 2.4159e+006 8.134e-010 1 -6.04e+004 112 5 30 2.41584e+006 1.414e-012 1 -6e+004 0.00494 x1=38.68237453582684 x2=9.84769440125371 x3=14.92074909280865 x4=6.999742560538413 |
2楼2011-09-30 08:32:51
100