| 查看: 1608 | 回复: 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 |
回复此楼» 猜你喜欢
磺酰氟产物，毕不了业了！
已经有8人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有10人回复
-
26申博（荧光探针方向，有机合成）
已经有4人回复
-
要不要辞职读博？
已经有3人回复
-
论文终于录用啦！满足毕业条件了
已经有26人回复
-
2026年机械制造与材料应用国际会议 （ICMMMA 2026）
已经有4人回复
-
Cas 72-43-5需要30g，定制合成，能接单的留言
已经有8人回复
-
北京211副教授，35岁，想重新出发，去国外做博后，怎么样？
已经有8人回复
-
自荐读博
已经有3人回复
-
不自信的我
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 使用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
- 金币: 68123.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
- 金币: 68123.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
- 金币: 68123.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3291.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
9楼2015-01-16 09:05:26