| 查看: 1212 | 回复: 6 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
dollapym专家顾问
|
[交流]
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，两个情况的误差差不多，但是解却差不少，这是问什么呢？ 我得到一个解，另它作为牛顿迭代法的初值去迭代，却无解，这说明了什么呢? 急求！谢谢了 |
回复此楼» 猜你喜欢
天津城建大学理学院化学、环境工程、材料物理等方向接收10人
已经有0人回复
-
天津城建大学理学院化学、环境工程、材料物理等方向接收调剂，10人
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有185人回复
-
计算与实验科学工程国际会议(8月7-11日,香港) 二维材料分会S5与声学换能器分会S8
已经有0人回复
-
0702一志愿吉大B区求调剂有论文
已经有4人回复
-
欢迎调剂到赣南师范大学智能制造与未来能源学院
已经有28人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有4人回复
-
欢迎调剂到赣南师范大学智能制造与未来能源学院
已经有22人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有0人回复
-
重庆交大26年硕士生招生拟调剂通知已出！欢迎加入光子学微结构与器件课题组。
已经有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人回复
dingd
实习版主
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
5楼2015-03-30 19:21:52
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
