| 查看: 1218 | 回复: 6 | |||
[交流]
1stopt1.5解四元多项式方程组无稳定解，求高版本帮忙运行！！
|
|
大家好，我在算一个四个参数的非线性多项式方程组如下，得到的目标函数大概0.13，求高版本的运行结果，看是否有更准确的解。 Parameters x2[0.0001,0.00025],K1[-3.3,3.3],K2[-3.3,3.3],K3[0.1,9] ; Function -0.0002483857122054368* K1^2 - 0.00016214005342224397* K3^2 +K1* (0.11832038801446097 + 0.0004143201485611828* K3 -2.597716388604742* x2) -6791.987324581259*(-0.09254221017610412 + x2) *(0.0014465166404352989 + x2) +(-0.07341654746178977 + 2.166562300403883* x2) = 0.8738; 0.9126034755805909 - 7.408480476916429*10^(-7)* K1^2 -1.8521201192291073*10^(-7)* K2^2 - 0.0036164973644227326* K3 +0.0018405191448191696* K3^2 + 4.83607293527163*10^(-7)* K3^3 -1.3692495660106428*10^(-7)* K3^4 + K1 *(0.006461909187108292 + 7.408480476916429*10^(-7)* K2 + 1.2357726636344952*10^(-6)* K3 - 6.178863318172476*10^(-7)* K3^2 ) + K2 *(-0.003230954593554146 - 6.178863318172476*10^(-7)* K3 + 3.34494856322119*10^(-7)* K3^2 ) + 0.004893093285395663* K3^2* x2 = 0.9034; 0.9119390763671137 - 0.000014217569378577079* K1^2 - 3.5543923446442698*10^(-6)* K2^2 - 0.0158429576267302* K3 + 0.008055685207646872* K3^2 + 9.280877865753863*10^(-6)* K3^3 - 2.62771843187007*10^(-6)* K3^4 + K1* (0.02830798508143708+ 0.000014217569378577079* K2 + 0.0000237156372836732* K3 - 0.0000118578186418366* K3^2 ) + K2* (-0.01415399254071854 - 0.0000118578186418366* K3 + 6.419270242197259*10^(-6)* K3^2 ) +0.02143540054716326* K3^2 *x2= 0.8912000000000001; 0.9108915129935684 - 0.00006986738872089048*K1^2 - 0.00001746684718022262*K2^2 - 0.03512051199327993* K3 + 0.017832729459977437* K3^2 + 0.00004560770440092129* K3^3 - 0.000012913024740019884* K3^4 + K1 *(0.06275286174349187 + 0.00006986738872089048* K2 + 0.0001165423993892175* K3 - 0.00005827119969460875* K3^2 ) + K2* (-0.031376430871745936- 0.00005827119969460875* K3 + 0.0000315453111128709* K3^2 ) +0.04751778422529179* K3^2 *x2 = 0.8808; 另外，我原方程是等式两百乘以500，得到目标函数是3.05，两个情况的误差差不多，但是解却差不少，这是问什么呢？ 我得到一个解，另它作为牛顿迭代法的初值去迭代，却无解，这说明了什么呢? 急求！谢谢了 |
回复此楼» 猜你喜欢
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有2人回复
-
070200求调剂，一志愿某211，288分
已经有18人回复
-
物理学I论文润色/翻译怎么收费?
已经有112人回复
-
法国博士后职位
已经有0人回复
-
重庆交大26年硕士生招生拟调剂通知已出！欢迎加入机器视觉与3D光学成像课题组。
已经有0人回复
-
**
已经有1人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有0人回复
-
基底STO,薄膜SRO，XRD里面的振荡，是laue震荡还是kiessig振荡? 怎么判断？
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 1stopt解非线性方程组结果不稳定
已经有2人回复
- 求大神帮忙用1stopt 5.0解一个方程组，原先的版本太低，每次结果不一样。
已经有11人回复
- 用1stOpt解出非线性方程组的解作为初值用Matlab求解，解不出来
已经有4人回复
- 谁能帮忙用1Stopt运行下程序呀？谢谢，我的1.5版本的运行不了啊
已经有6人回复
- 求助1stopt拟合动力学参数
已经有4人回复
- 求助：以下非线性拟合方程可以用1stopt求解么？用什么方法
已经有16人回复
- 1stopt拟合问题求助
已经有7人回复
- 求大神用1stopt拟合数据 紧急
已经有6人回复
- 1stopt 解非线性方程组
已经有3人回复
- 求助，1stopt 解方程组为什么总是出现常数设置错误呢
已经有8人回复
- 求助，1stopt解非线性方程组。
已经有2人回复
- sos!!1stopt解非线性方程组，一直出现常数定义错误的问题！求大神帮助！
已经有12人回复
WanderingHeart
铁杆木虫 (著名写手)
- 应助: 90 (初中生)
- 金币: 13890.8
- 散金: 5
- 红花: 18
- 帖子: 1885
- 在线: 467.1小时
- 虫号: 257464
- 注册: 2006-06-04
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
2楼2015-03-30 12:59:27
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
3楼2015-03-30 13:53:41
shikang999
新虫 (小有名气)
- 计算强帖: 1
- 应助: 27 (小学生)
- 金币: 1454.7
- 红花: 7
- 帖子: 182
- 在线: 206.4小时
- 虫号: 2201981
- 注册: 2012-12-23
- 性别: GG
- 专业: 固体力学
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
小木虫: 金币+0.5, 给个红包，谢谢回帖
|
1、楼上dingd的回答可能有欠缺，4个方程4个未知数，一般情况下（不加任何服务约束）也就只有一组解，如果是线性满秩方程组的话，这是没有错的。如果是非线性方程组的话，存在多组解（均满足设定的误差规则）的情况就显得很正常了。下面给出一组在所设定的范围内找到的一组解，具体误差是否满足你的要求，可以自己测试看 x2=0.00018602746110 k1=-0.65443662426355 k2= -0.35539381460637 k3=0.10000000000929 2、用那个做初值进行牛顿迭代反而无解，只能说明你那个初值进去，采用那种牛顿迭代算法无法找到满足条件的解而已。 |
4楼2015-03-30 16:57:05
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
5楼2015-03-30 19:21:52
6楼2015-03-31 13:59:21
|
特别谢谢你，你的说法也是我在想的，这组解的误差完全满足我的需求，唉，只是我的一些参数值又变了一下，所以方程组又变了，依然用1.5解，得到的误差更大了。能否再帮我运行一下，谢谢了。 这里x2这个参数值最好在5.8*10^-5左右，目标函数误差在0.1这个order的都可以接受。若有多个相似误差的解，方便的话一块贴上来。thanks again. Parameters x2[0.000056,0.000065],K1[-3.3,3.3],K2[-3.3,3.3],K3[1,9] ; Function -0.1176618358 *K1^2 - 0.08556995644 *K3^2 + K1* (55.94823395 + 0.2071600743* K3 - 1230.554191 *x2) - 3.217406066*10^6 *(-0.09245993440 + x2)* (0.001528166642 + x2) + K3 *(-36.61948488 + 1083.281150 x2)= 436.9; 456.3017378 - 0.0003509442655 *K1^2 - 0.00008773606637* K2^2 - 1.803399598 *K3 + 0.9178216086* K3^2 + 0.0002552253695* K3^3 - 0.00007226260461* K3^4 + K1 *(3.055537706 + 0.0003509442655* K2 + 0.0006178863318* K3 - 0.0003089431659* K3^2) + K2 *(-1.527768853 - 0.0003089431659* K3 + 0.0001672474282 *K3^2) + 2.446546643 *K3^2* x2= 451.7; 455.9695382 - 0.006734949843* K1^2 - 0.001683737461* K2^2 - 7.900236206* K3 + 4.016963725* K3^2 + 0.004898014389* K3^3 - 0.001386787206* K3^4 + K1* (13.38553565 + 0.006734949843* K2 + 0.01185781864 *K3 - 0.005928909321* K3^2) + K2 *(-6.692767825 - 0.005928909321* K3 + 0.003209635121* K3^2) + 10.71770027 *K3^2* x2 = 445.6; 455.4457565 - 0.03309661069* K1^2 - 0.008274152672* K2^2 - 17.51316559 *K3 + 8.891553762* K3^2 + 0.02406961880* K3^3 - 0.006814892071 *K3^4 + K1*(29.67292323 + 0.03309661069* K2 + 0.05827119969 *K3 - 0.02913559985* K3^2) + K2* (-14.83646161 - 0.02913559985* K3 + 0.01577265556* K3^2) + 23.75889211* K3^2 *x2 = 440.4 ; |
7楼2015-03-31 14:06:27
