| 查看: 2879 | 回复: 12 | ||
| 【悬赏金币】回答本帖问题，作者缱莳霰瑜将赠送您 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值，怎么编译？谢谢啦 |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：libolin3@tal.com - 附件 1 : 数据_1.txt
2014-02-17 21:35:33, 38.07 K
» 猜你喜欢
声子谱计算问题
已经有2人回复
-
已结束
已经有23人回复
-
物理学I论文润色/翻译怎么收费?
已经有289人回复
-
Error in routine lambda (100): wrong # or too many modes
已经有0人回复
-
祈福国自然面上
已经有56人回复
-
请问castep计算出的能带图横轴数据缩放因子怎么计算
已经有1人回复
-
Result in Physics 从见刊开始 一般多久可以被SCI收录
已经有15人回复
-
面青地会评时间？？？
已经有5人回复
-
科研求助
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 请教高手们一个非线性方程如何解？急...
已经有7人回复
- 请大神求解复杂非线性方程组（有报酬）
已经有22人回复
- 1stopt非线性拟合求助
已经有6人回复
- mathematica可以求解不等式方程组吗？
已经有9人回复
- 求解方程，在线等……求大神指点！
已经有7人回复
- 1stopt 1.5计算结果不唯一，求高版本跑代码验证
已经有6人回复
- 1stopt每次运行结果都不一样
已经有15人回复
- 隐函数方程组求解
已经有15人回复
- 求助：matlab求解一个方程组，其中一个方程有三角函数
已经有10人回复
- 求解非线性方程组
已经有7人回复
- 用matlab求解非线性方程组说无解，一定是方程组本身无解，还是有可能程序有问题呢？
已经有11人回复
- （急）：用1stopt跑一段程序，解决非线性方程组的解（已有程序，需要较高版本的1stopt
已经有4人回复
- 这个方程组怎么解啊
已经有14人回复
- rmse，sse计算求助
已经有3人回复
- 求助，用1stopt解方程组，两次得到的解不唯一？什么原因，大神帮忙
已经有6人回复
- 求Matlab解非线性方程代码
已经有12人回复
- 用matlab如何解方程？
已经有6人回复
- 求大哥大姐帮跑一下1stopt~新手啊，版本低啊
已经有9人回复
- 第5章_解线性方程组的直接方法.ppt
已经有13人回复
- 1stopt问题请教(附代码)
已经有5人回复
- matlab微分方程组参数拟合，以周为单位求解，汇总后以年为单位进行数值比较
已经有4人回复
- 求助解6变量的指数方程组的方法
已经有4人回复
- 高金求助matlab解微分方程组
已经有12人回复
- 【求助】求教matlab解非线性方程组
已经有9人回复
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15356.8
- 散金: 101
- 红花: 234
- 帖子: 3409
- 在线: 1222.7小时
- 虫号: 291104
- 注册: 2006-10-28
2楼2014-02-18 11:37:35
3楼2014-02-18 15:36:34
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15356.8
- 散金: 101
- 红花: 234
- 帖子: 3409
- 在线: 1222.7小时
- 虫号: 291104
- 注册: 2006-10-28
4楼2014-02-18 16:18:10
5楼2014-02-18 16:33:07
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15356.8
- 散金: 101
- 红花: 234
- 帖子: 3409
- 在线: 1222.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 (讲师)
- 金币: 15356.8
- 散金: 101
- 红花: 234
- 帖子: 3409
- 在线: 1222.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
5