| 查看: 650 | 回复: 4 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
Matlab二次规划语言，高分求大神帮忙修改
|
||
|
在网上看了一遍关于配料比优化的论文，完全按照论文内容试着编了一下，可是无法运行，求大神帮忙改一下。论文就不附上了，有涉嫌侵权嫌疑，把从里面摘抄的一部分附上。 1 function f=myfun(X) 2 f=0.35*X(1)+0.002*X(2)+0.02*X(3)+0.025*X(4) +0.035*X(5)+6*X(6); 3 goal=170; 4 weight=170; 5 X0=[400;170;100;550;1200;6]; 6 b=[0;0;0;0;0;0;0;-500;600]; 7 Aeq=[0.3226 1 0.4808 0.3774 0.3704 0]; 8 beq=990; 9 lb=[300;130;50;500;1000;315]; 10 ub=[500;210;200;1 000;1400;12]; 11 A=[0.2 -1 0.2 0 0 0; -0.4 1 -0.4 0 0 0; 0 0 0 -0.65 0.35 0; 0 0 0 0.6 01.4 0; -14.6 52.71 -14.6 0 0 0; 0.01 0 0 0 0 -1; -0.02 0 0 0 0 1; -1 0 -1 0 0 0; 1 0 1 0 0 0]; 12 [X,fval,attainfactor,exitflag]=fgoalattain(xy,X0,goal,weight,A,b,Aeq,beq,lb,ub)[ Last edited by troyzz on 2013-5-18 at 12:07 ] |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : 程序语言.docx
2013-05-17 16:43:10, 31.2 K
» 猜你喜欢
时间戳又变了8-15
已经有18人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有13人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有58人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有7人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有5人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有5人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有5人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有3人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有9人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 给小老板认识的老外写邮件，请求帮忙修改投稿语言，求个模板，多谢多谢
已经有7人回复
- 帮忙用Mathematica求解一个微分方程组
已经有6人回复
- 帮忙看一下这条语句是什么意思
已经有7人回复
- linux终端运行程序出错，求大神帮忙
已经有6人回复
- 求高手帮忙修改一段c++程序
已经有4人回复
- 初学matlab 求解方程 大家帮帮忙啊。。。
已经有11人回复
- Perl 生物信息学编程求助，运行出现问题，烦请帮忙查找一下错误，谢谢！
已经有3人回复
- 求高人帮忙看一下我这个Fortran程序的错误出在哪里？
已经有12人回复
- 请帮忙用matlab求矩阵特征值
已经有8人回复
- 帮忙解释一个read语句，谢谢
已经有4人回复
- 求帮忙啊~ matlab一多项式结果太长，怎么写到txt里？
已经有14人回复
- 有会看fortran程序的，能帮忙看看这个程序吗？
已经有10人回复
- 内循环问题，求大虾帮忙
已经有3人回复
- JCA审稿意见返回，编辑要求很严厉，语言的问题，想找一个公司帮忙改，大家推荐下
已经有12人回复
- 求帮忙，关于含有bessel函数的方程的编程求解问题
已经有11人回复
- 写了一个fortran90的小程序，编译通不过，请大侠帮忙
已经有59人回复
- 【求助】50金币求大家帮忙改一段程序
已经有13人回复
- 【求助】帮忙看看哪错了error: invalid lvalue in assignment
已经有6人回复
- 【求助】f90子程序怎么编译错误？请大家帮忙，谢谢
已经有3人回复
- 【求助】用fortran怎么描述下面这样的情况，我是新手，请大家帮忙
已经有8人回复
- 【求助】帮忙看一个提示，子程序还是undefined reference to
已经有4人回复
- 【求助】帮忙改写展宽程序
已经有11人回复
- 【求助】帮忙看一下，问题程序如何修改
已经有3人回复
5楼2013-05-19 11:52:02
ausir
木虫 (正式写手)
- 应助: 18 (小学生)
- 金币: 2009.2
- 红花: 4
- 帖子: 500
- 在线: 178.9小时
- 虫号: 1939287
- 注册: 2012-08-15
- 性别: GG
- 专业: 集成电路制造与封装
【答案】应助回帖
感谢参与，应助指数 +1
|
这里面有几个问题： 1. 把 myfun 这个用户自定义函数单独存在一个.m文件中。后面在 fgoalattain中调用； 2. ub 的最后一项比 lb 对应的数值小，ub中是12， lb中是315； 3. ub 矩阵中的写法，1000 写作 1 000，这个导致识别错误； 4. 矩阵A的第4行第5列的数值到底是多少？ 01.4？ 可以按照1.4来计算，但是计算迭代到700次时依然无法收敛，系统自动停止运算，没有结果。 5. 最后的关键，第12行的函数中，根本就没有提到你自己定义的函数 myfun，而是使用了没有任何定义的 xy。根本不会有任何运算。 6. 关于fgoalattain函数的调用语法问题，在myfun 前面要加一个 @符号或使用 'myfun' 。运算结果一样。都无法取得附件中中所要求的结果。 以上看法供参考。 |
2楼2013-05-18 13:35:58
ausir
木虫 (正式写手)
- 应助: 18 (小学生)
- 金币: 2009.2
- 红花: 4
- 帖子: 500
- 在线: 178.9小时
- 虫号: 1939287
- 注册: 2012-08-15
- 性别: GG
- 专业: 集成电路制造与封装
3楼2013-05-18 13:38:06
|
原来是复制过程中出现很多错误，我又仔细的看了一下那篇论文，并参考了您的意见，自己演算了一遍，附上修改后的语言还是不行，其实我对fgoalattain不熟悉。 function f=myfun(X) f=0.35*X(1)+0.002*X(2)+0.02*X(3)+0.025*X(4) +0.035*X(5)+6*X(6); goal=170; weight=170; X0=[400;170;100;550;1200;6]; b=[0;0;0;0;0;0;0;-500;600]; Aeq=[0.3226 1 0.4808 0.3774 0.3704 0]; beq=990; lb=[300;130;50;500;1000;3.5]; ub=[500;210;200;1000;1400;12]; A=[0.2 -1 0.2 0 0 0; -0.4 1 -0.4 0 0 0; 0 0 0 -0.65 0.35 0; 0 0 0 0.6 -0.4 0; -14.6 52.708 -14.6 0 0 0; 0.01 0 0 0 0 -1; -0.02 0 0 0 0 1; -1 0 -1 0 0 0; 1 0 1 0 0 0]; [X,fval,attainfactor,exitflag]=fgoalattain(@myfun,X0,goal,weight,A,b,Aeq,beq,lb,ub,nonlcon) |
4楼2013-05-18 15:01:59
