| 查看: 3491 | 回复: 12 | |||
[求助]
1stopt解方程组
|
|
现有一个方程组 Variable b,a; Function b=x+z^2*(1-x)*(1-y)*(y+(1-x)*(1-y)*a)*(y+z^2*y^2); Data; 0.19530 1.00328 //b0 a0 0.18546 0.52781 //b1 a1 0.17207 0.03061 //b2 a2 假设c是一组点，可以通过一个点c1对应的6个值,即上面的数据来解b,a 通过解得x,y,z求得的目标值和实际值进行比较。 现在我想把c其它点对应的值一起放进来求其它点对应的x,y,z值，怎么编译？谢谢啦 |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : 数据_1.txt
2014-02-17 21:35:33, 38.07 K
» 猜你喜欢
B区学生调剂-兰州交通大学材料科学与工程学院
已经有12人回复
-
山西大同大学物理学专业还有调剂名额，欢迎调剂！
已经有13人回复
-
物理学I论文润色/翻译怎么收费?
已经有103人回复
-
桂林理工大学物理学专业招收调剂
已经有18人回复
-
VASP 的一组 GPU / CPU 基准测试记录
已经有0人回复
-
津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生一名
已经有0人回复
-
【原创讨论】从电子约束到物质编辑：一套可迭代的环形磁场科技树
已经有0人回复
-
【方案分享】单环磁场+轴心控制+偏转导出电子束约束系统（可行性实验）
已经有6人回复
-
【修正版】单环用磁约束低速电子实验方案（简化版）
已经有0人回复
-
桂林理工大学物理学专业招收调剂，还有三个名额！！！
已经有11人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求解方程，在线等……求大神指点！
已经有7人回复
- 求助：matlab求解一个方程组，其中一个方程有三角函数
已经有10人回复
- 用matlab求解非线性方程组说无解，一定是方程组本身无解，还是有可能程序有问题呢？
已经有11人回复
- 这个方程组怎么解啊
已经有14人回复
- rmse，sse计算求助
已经有3人回复
- 求助，用1stopt解方程组，两次得到的解不唯一？什么原因，大神帮忙
已经有6人回复
- 求Matlab解非线性方程代码
已经有12人回复
- 用matlab如何解方程？
已经有6人回复
- 求大哥大姐帮跑一下1stopt~新手啊，版本低啊
已经有9人回复
- 第5章_解线性方程组的直接方法.ppt
已经有13人回复
- matlab微分方程组参数拟合，以周为单位求解，汇总后以年为单位进行数值比较
已经有4人回复
- 求助解6变量的指数方程组的方法
已经有4人回复
- 高金求助matlab解微分方程组
已经有12人回复
- 【求助】求教matlab解非线性方程组
已经有9人回复
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
2楼2014-02-18 11:37:35
3楼2014-02-18 15:36:34
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
4楼2014-02-18 16:18:10
5楼2014-02-18 16:33:07
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
【答案】应助回帖
|
就是循环求解吧： 太多了，前10组： 循环常数 a0 循环常数 b0 循环常数 a1 循环常数 b1 循环常数 a2 循环常数 b2 目标函数值 x z y 95.736 400.21 48.45 15.891 3.101 15.504 23545.1733930783 -61.4105433019701 1.7338724840528 0.000356939921725784 99.225 401 49.225 15.891 1.938 13.953 23104.8556376321 -57.6032905084872 0.961829781946052 0.00126510256991255 100 401.78 50 15.556 2.222 14.444 23597.4500182412 -58.404029918814 0.0314462319595226 1.79701935010195 102.198 402.57 51.282 17.216 2.93 15.385 23247.9905035161 -59.9887251506348 0.157625120101228 0.0467799755934256 98.993 403.36 47.987 16.107 3.356 14.094 21483.3270780367 -63.5176155355752 -0.0132174322981234 2.508604718078 99.667 404.14 48 16 3 15.333 21644.6502221687 -61.0348067536678 0.0132327858744587 2.54676620819215 99.357 404.93 48.232 16.399 3.215 15.113 21869.1450662234 14289.5808987558 -86.2702757233859 0.999999981895826 100 405.72 49.211 15.773 2.524 14.826 23118.6454782423 -60.2819279159388 -0.143995346767101 0.735665078490286 97.583 406.5 48.338 16.012 3.323 15.408 23089.8608533821 -63.3054946038168 -0.000113023912612222 43.6971717119863 97.929 407.29 48.521 15.68 2.367 14.497 23685.5386771595 -60.8677794909951 -0.317834915126589 0.895069882125857 97.399 408.08 49.422 16.474 2.89 14.162 24480.8908333929 -62.7941450075183 0.336825876489558 0.00923229906277813 |
6楼2014-02-18 18:01:43
7楼2014-02-19 10:11:34
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
8楼2014-02-19 10:38:48
9楼2014-02-19 10:40:44
|
谢谢您，我用的是1.5版本的 按照您上面的计算输入时，结果显示是这样的 ====== 结果 ====== 迭代数: 16 计算用时(时:分:秒:毫秒): 00:00:00:32 计算中止原因: 达到收敛判定标准 优化算法: 麦夸特法(Levenberg-Marquardt) + 通用全局优化法 函数表达式 1: b0-(x+z^2*(1-x)*(1-y)*(y+(1-x)*(1-y)*a0)*(y+z^2*y^2)) 2: b1-(x+z^2*(1-x)*(1-y)*(y+(1-x)*(1-y)*a1)*(y+z^2*y^2)) 3: b2-(x+z^2*(1-x)*(1-y)*(y+(1-x)*(1-y)*a2)*(y+z^2*y^2)) 目标函数值: 0 x: 0 y: 0 z: 0 b0: 0 a0: 0 b1: 0 a1: 0 b2: 0 a2: 0 ====== 计算结束 ====== |
10楼2014-02-19 11:09:41
