| 查看: 1646 | 回复: 8 | ||||
matangtang木虫 (正式写手)
|
[求助]
matlab优化报错 已有1人参与
|
|
程序在附件中，其中dynamics_inverse_qiujie.m为优化主程序，biaozhun_function.m为优化目标函数，yueshu.m为约束函数。 程序里面的初值a0已经满足我的需要了（是手调得到的），但仍然存在一个错误。如图片所示 yueshu.m中计算出来的g的个数本来是10个，最后一次计算变成了11个，如下图所示。 还请高手解答 因为不能在原贴上编辑，所以又重新发了这个求助帖，原帖在。如果能解决问题，金币可以都奉送哈 http://muchong.com/bbs/viewthread.php?tid=8413221&page=2  1.png  2.png |
回复此楼» 猜你喜欢
一志愿985，本科211，0817化学工程与技术319求调剂
已经有4人回复
-
311求调剂
已经有6人回复
-
0703化学调剂 290分有科研经历，论文在投
已经有4人回复
-
化学工程321分求调剂
已经有6人回复
-
本人考085602 化学工程 专硕
已经有9人回复
-
材料与化工一志愿南昌大学327求调剂推荐
已经有7人回复
-
326求调剂
已经有3人回复
-
焦虑
已经有8人回复
-
308求调剂
已经有4人回复
-
NSFC申报书里申请人简历中代表性论著还需要在申报书最后的附件里面再上传一遍吗
已经有14人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 使用matlab最优化方法拟合获得多个动力学参数中的问题
已经有4人回复
- 菜鸟求助，老师给了一个分子让优化，可总是错误，求大神指点
已经有22人回复
- Matlab心得及学习方法（不断更新）
已经有26人回复
- matlab与comsol断开连接
已经有10人回复
- matlab求函数积分出错
已经有12人回复
- workbench计算报错
已经有4人回复
- matlab一运行就出现的问题
已经有7人回复
- 求教matlab线性优化求最小值的问题
已经有7人回复
- 请教MATLAB中的LMI求解问题，急！！！
已经有7人回复
- 求助：遗传算法的matlab编程求解资源优化问题
已经有7人回复
- 求助：资源优化问题的遗传算法matlab编程
已经有5人回复
- 求助matlab问题
已经有5人回复
- 求助有关Matlab有约束非线性最优化问题
已经有9人回复
- Matlab两个小问题
已经有10人回复
- 【求助】拜求高手解释下matlab的报错-----“out of memory”
已经有13人回复
- 【求助】matlab 二次规划的优化的问题
已经有4人回复
- 【分享】推荐一本好书《精通MATLAB最优化计算》
已经有880人回复
NicoleLas
木虫 (正式写手)
- 应助: 66 (初中生)
- 金币: 4288.7
- 红花: 26
- 帖子: 514
- 在线: 177.2小时
- 虫号: 2872321
- 注册: 2013-12-15
- 性别: GG
- 专业: 计算机软件
2楼2015-01-14 23:02:25
matangtang
木虫 (正式写手)
- 应助: 7 (幼儿园)
- 金币: 3131.5
- 散金: 856
- 红花: 6
- 帖子: 614
- 在线: 278.4小时
- 虫号: 1871554
- 注册: 2012-06-27
- 性别: GG
- 专业: 制造系统与自动化
3楼2015-01-15 00:33:35
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68127.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3291.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
matangtang: 金币+50, ★★★★★最佳答案 2015-01-15 10:41:09
感谢参与，应助指数 +1
matangtang: 金币+50, ★★★★★最佳答案 2015-01-15 10:41:09
|
您好！ 我大致看了一下，发现主要问题在于 yueshu.m 中 变量g 的维数问题。 按原来的程序： g=abs(double(z))-0.99; 那么 变量z 的维数就决定了 g 的维数 而 z 是用solve 函数解方程组得到的： [A B C Q H z]=solve('A=1-1.1/z','B=1+1.2/z','C=1+0.2/z',['Q=((',num2str(a(9)),')*(z-((',num2str(a(1)),')+(',num2str(a(2)),')*i))*(z-((',num2str(a(3)),')+(',num2str(a(4)),')*i)))/((z-((',num2str(a(5)),')+(',num2str(a(6)),')*i))*(z-((',num2str(a(7)),')+(',num2str(a(8)),')*i)))-B*C'],'H=z*Q-B*C-A*Q+z*B*C','z^2-(C*(Q+B*C))*z/H-((Q+B*C)*(B*C^2-z*Q*C+A*Q*C-z*B*C^2))/H^2','A,B,C,Q,H,z') 由于solve 函数会求解方程组所有的解，那么问题就来了，如果原方程组是多解的，而且解的个数不唯一，就会引起z维数的变化 通过跟踪 优化过程中 z 维数的变化可见（以下以index_xy表示z的维数）： index_xy = 6 1 index_xy = 6 1 index_xy = 6 1 index_xy = 6 1 index_xy = 6 1 index_xy = 6 1 index_xy = 11 1 可见第7次，z 的维数为11。 为了进一步验证，变量g和z 的维数就是调通程序的关键，尝试将变量z 改成以下形式： g=abs(double(z(1:6)))-0.99; 即强制 g和z 维数=6 最终程序可以跑通的，结果见附图1。当然结果肯定不好，不过这样确实说明，调整 g和z 维数，确实可以跑通程序。 基本上找到了原因，那么怎么解决呢？ 由于solve 函数会求解方程组所有的解，解的个数会随着方程的变化而不同，所以是否可以给 z 中的元素加个约束，只取出其中合适的6个，来代入优化过程。或者，不用 solve函数，用 fsolve函数代入不同的初值，求取不同的解。需要说明的是，以上只是我初步的想法，我也没有想到如何用代码具体实现。  附图1.png |
4楼2015-01-15 09:33:08
matangtang
木虫 (正式写手)
- 应助: 7 (幼儿园)
- 金币: 3131.5
- 散金: 856
- 红花: 6
- 帖子: 614
- 在线: 278.4小时
- 虫号: 1871554
- 注册: 2012-06-27
- 性别: GG
- 专业: 制造系统与自动化
5楼2015-01-15 10:48:28
matangtang
木虫 (正式写手)
- 应助: 7 (幼儿园)
- 金币: 3131.5
- 散金: 856
- 红花: 6
- 帖子: 614
- 在线: 278.4小时
- 虫号: 1871554
- 注册: 2012-06-27
- 性别: GG
- 专业: 制造系统与自动化
6楼2015-01-15 11:26:20
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68127.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3291.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
7楼2015-01-15 14:02:37
matangtang
木虫 (正式写手)
- 应助: 7 (幼儿园)
- 金币: 3131.5
- 散金: 856
- 红花: 6
- 帖子: 614
- 在线: 278.4小时
- 虫号: 1871554
- 注册: 2012-06-27
- 性别: GG
- 专业: 制造系统与自动化
8楼2015-01-15 22:47:01
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68127.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3291.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
9楼2015-01-16 09:05:26
20