| 查看: 1228 | 回复: 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，两个情况的误差差不多，但是解却差不少，这是问什么呢？ 我得到一个解，另它作为牛顿迭代法的初值去迭代，却无解，这说明了什么呢? 急求！谢谢了 |
回复此楼» 猜你喜欢
基底STO,薄膜SRO，XRD里面的振荡，是laue震荡还是kiessig振荡? 怎么判断？
已经有2人回复
-
B区学生调剂-兰州交通大学材料科学与工程学院
已经有12人回复
-
物理学I论文润色/翻译怎么收费?
已经有57人回复
-
山西大同大学物理学专业还有调剂名额，欢迎调剂！
已经有13人回复
-
桂林理工大学物理学专业招收调剂
已经有17人回复
-
VASP 的一组 GPU / CPU 基准测试记录
已经有0人回复
-
此贴已删除
已经有7人回复
-
津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生一名
已经有0人回复
-
【原创讨论】从电子约束到物质编辑：一套可迭代的环形磁场科技树
已经有0人回复
-
【方案分享】单环磁场+轴心控制+偏转导出电子束约束系统（可行性实验）
已经有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
