| 查看: 1053 | 回复: 3 | ||
| 【悬赏金币】回答本帖问题，作者superl将赠送您 20 个金币 | ||
superl银虫 (小有名气)
|
[求助]
1stopt-- 求给定范围初始值，求满足方程组的初始值的最优解已有1人参与
|
|
|
我用下载的免费版本的1stopt运行，这其实就是一个反演计算，运行给出的结果很烂，基本都不符合限定相加为1的条件，第一次用这个软件，求高手指点。 Parameter a11=[2,8],a12=[1.2,2.8],a13=[0,0],a14=[0.15,5],a15=[30,70],a16=[0.2,0.5], a21=[0.35,0.85], a22=[0.1,0.3], a23=[0,0], a24=[0.01,0.5],a25=[10,28], a26=[0.15,0.25], a31=[0.3,0.7], a32=[0.09,0.23], a33=[1,1.8], a34=[0,0], a35=[0,0], a36=[0.1,0.24], a41=[1.7,3.9], a42=[0.45,0.75], a43=[4,6], a44=[1,1], a45=[0,0], a46=[0,0], a51=[10,50], a52=[1.5,2.7], a53=[0,0], a54=[0.5,5], a55=[60,120], a56=[2,4], a61=[0.1,0.9], a62=[0,0], a63=[3,5], a64=[0,0], a65=[0,0], a66=[0,0], x11=[0,1],x12=[0,1],x13=[0,1],x14=[0,1],x15=[0,1],x16=[0,1], x21=[0,1],x22=[0,1],x23=[0,1],x24=[0,1],x25=[0,1],x26=[0,1], x31=[0,1],x32=[0,1],x33=[0,1],x34=[0,1],x35=[0,1],x36=[0,1], x41=[0,1],x42=[0,1],x43=[0,1],x44=[0,1],x45=[0,1],x46=[0,1], x51=[0,1],x52=[0,1],x53=[0,1],x54=[0,1],x55=[0,1],x56=[0,1], x61=[0,1],x62=[0,1],x63=[0,1],x64=[0,1],x65=[0,1],x66=[0,1], x71=[0,1],x72=[0,1],x73=[0,1],x74=[0,1],x75=[0,1],x76=[0,1], x81=[0,1],x82=[0,1],x83=[0,1],x84=[0,1],x85=[0,1],x86=[0,1]; Constant y11=0.20435,y12=0.17197,y13=0.01914,y14=0.8352,y15=1.2722,y16=0.008722, y21=0.27297,y22=0.24468,y23=0.0215,y24=0.81946,y25=1.44834,y26=0.0131, y31=0.3553,y32=0.236,y33=0.0259,y34=0.808,y35=1.5967,y36=0.017235, y41=0.8103,y42=0.4272,y43=0.04405,y44=0.7924,y45=2.98,y46=0.04905, y51=1.5548,y52=0.61867,y53=0.0929,y54=0.7125,y55=4.654,y56=0.1598, y61=1.9693,y62=0.6804,y63=0.1079,y64=0.7732,y65=6.2797,y66=0.1936, y71=1.65822,y72=0.629438,y73=0.11123,y74=0.7224,y75=4.579118,y76=0.17304, y81=1.7544,y82=0.65046,y83=0.129776,y84=0.822139,y85=4.5582,y86=0.22376; Function a11*x11+a12*x12+a13*x13+a14*x14+a15*x15+a16*x16=y11; a21*x11+a22*x12+a23*x13+a24*x14+a25*x15+a26*x16=y12; a31*x11+a32*x12+a33*x13+a34*x14+a35*x15+a36*x16=y13; a41*x11+a42*x12+a43*x13+a44*x14+a45*x15+a46*x16=y14; a51*x11+a52*x12+a53*x13+a54*x14+a55*x15+a56*x16=y15; a61*x11+a62*x12+a63*x13+a64*x14+a65*x15+a66*x16=y16; a11*x21+a12*x22+a13*x23+a14*x24+a15*x25+a16*x26=y21; a21*x21+a22*x22+a23*x23+a24*x24+a25*x25+a26*x26=y22; a31*x21+a32*x22+a33*x23+a34*x24+a35*x25+a36*x26=y23; a41*x21+a42*x22+a43*x23+a44*x24+a45*x25+a46*x26=y24; a51*x21+a52*x22+a53*x23+a54*x24+a55*x25+a56*x26=y25; a61*x21+a62*x22+a63*x23+a64*x24+a65*x25+a66*x26=y26; a11*x31+a12*x32+a13*x33+a14*x34+a15*x35+a16*x36=y31; a21*x31+a22*x32+a23*x33+a24*x34+a25*x35+a26*x36=y32; a31*x31+a32*x32+a33*x33+a34*x34+a35*x35+a36*x36=y33; a41*x31+a42*x32+a43*x33+a44*x34+a45*x35+a46*x36=y34; a51*x31+a52*x32+a53*x33+a54*x34+a55*x35+a56*x36=y35; a61*x31+a62*x32+a63*x33+a64*x34+a65*x35+a66*x36=y36; a11*x41+a12*x42+a13*x43+a14*x44+a15*x45+a16*x46=y41; a21*x41+a22*x42+a23*x43+a24*x44+a25*x45+a26*x46=y42; a31*x41+a32*x42+a33*x43+a34*x44+a35*x45+a36*x46=y43; a41*x41+a42*x42+a43*x43+a44*x44+a45*x45+a46*x46=y44; a51*x41+a52*x42+a53*x43+a54*x44+a55*x45+a56*x46=y45; a61*x41+a62*x42+a63*x43+a64*x44+a65*x45+a66*x46=y46; a11*x51+a12*x52+a13*x53+a14*x54+a15*x55+a16*x56=y51; a21*x51+a22*x52+a23*x53+a24*x54+a25*x55+a26*x56=y52; a31*x51+a32*x52+a33*x53+a34*x54+a35*x55+a36*x56=y53; a41*x51+a42*x52+a43*x53+a44*x54+a45*x55+a46*x56=y54; a51*x51+a52*x52+a53*x53+a54*x54+a55*x55+a56*x56=y55; a61*x51+a62*x52+a63*x53+a64*x54+a65*x55+a66*x56=y56; a11*x61+a12*x62+a13*x63+a14*x64+a15*x65+a16*x66=y61; a21*x61+a22*x62+a23*x63+a24*x64+a25*x65+a26*x66=y62; a31*x61+a32*x62+a33*x63+a34*x64+a35*x65+a36*x66=y63; a41*x61+a42*x62+a43*x63+a44*x64+a45*x65+a46*x66=y64; a51*x61+a52*x62+a53*x63+a54*x64+a55*x65+a56*x66=y65; a61*x61+a62*x62+a63*x63+a64*x64+a65*x65+a66*x66=y66; a11*x71+a12*x72+a13*x73+a14*x74+a15*x75+a16*x76=y71; a21*x71+a22*x72+a23*x73+a24*x74+a25*x75+a26*x76=y72; a31*x71+a32*x72+a33*x73+a34*x74+a35*x75+a36*x76=y73; a41*x71+a42*x72+a43*x73+a44*x74+a45*x75+a46*x76=y74; a51*x71+a52*x72+a53*x73+a54*x74+a55*x75+a56*x76=y75; a61*x71+a62*x72+a63*x73+a64*x74+a65*x75+a66*x76=y76; a11*x81+a12*x82+a13*x83+a14*x84+a15*x85+a16*x86=y81; a21*x81+a22*x82+a23*x83+a24*x84+a25*x85+a26*x86=y82; a31*x81+a32*x82+a33*x83+a34*x84+a35*x85+a36*x86=y83; a41*x81+a42*x82+a43*x83+a44*x84+a45*x85+a46*x86=y84; a51*x81+a52*x82+a53*x83+a54*x84+a55*x85+a56*x86=y85; a61*x81+a62*x82+a63*x83+a64*x84+a65*x85+a66*x86=y86; x11+x12+x13+x14+x15+x16=1; x21+x22+x23+x24+x25+x26=1; x31+x32+x33+x34+x35+x36=1; x41+x42+x43+x44+x45+x46=1; x51+x52+x53+x54+x55+x56=1; x61+x62+x63+x64+x65+x66=1; x71+x72+x73+x74+x75+x76=1; x81+x82+x83+x84+x85+x86=1; [ Last edited by superl on 2014-8-6 at 20:03 ] |
回复此楼» 猜你喜欢
声子谱计算问题
已经有2人回复
-
已结束
已经有23人回复
-
物理学I论文润色/翻译怎么收费?
已经有274人回复
-
Error in routine lambda (100): wrong # or too many modes
已经有0人回复
-
祈福国自然面上
已经有57人回复
-
请问castep计算出的能带图横轴数据缩放因子怎么计算
已经有1人回复
-
Result in Physics 从见刊开始 一般多久可以被SCI收录
已经有15人回复
-
面青地会评时间？？？
已经有7人回复
-
科研求助
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 1stopt为何结果是 计算值全为0
已经有8人回复
- 求大神跑下1stOpt，一个sharemodel，谢谢！
已经有7人回复
- 1stOpt软件同一组数据每次拟合结果为什么会差距这么大呢？求高手指教~
已经有9人回复
- 用1stOpt进行复数曲线拟合时，结果出错。急求教！！！
已经有5人回复
- 1stopt求最小值
已经有37人回复
- 求助帮忙做个多元线性回归的分析spss，R均可
已经有4人回复
- 求微分方程 最大值问题
已经有4人回复
- origin关于四参数logstic曲线和曲线方程求法问题
已经有6人回复
- 怎样拟合数据求出x=f(A,B,C),y=g(A,B,C,D)中的A,B,C,D参数
已经有18人回复
- 急求一个含复数的非线性拟合。
已经有8人回复
- 请帮忙拟合一组数据，求回归系数
已经有3人回复
- 求助高版本1stopt计算非线性方程组
已经有4人回复
- 毕不了业了，发帖求助Origin8.0多元非线性拟合的问题
已经有9人回复
- 微分方程组求参数问题，求高人指教，感谢
已经有11人回复
- 1stopt中使用sharedmodel拟合曲线出现问题，求指示！！
已经有6人回复
- 求高手解答Matlab求解微分方程的初始值和边界值问题
已经有5人回复
- 求大哥大姐帮跑一下1stopt~新手啊，版本低啊
已经有9人回复
- 求助matlab拟合方程组
已经有12人回复
- 拟牛顿法求解复杂（带根号的）非线性方程组
已经有17人回复
- 大家在使用1stOpt进行数据拟合的时候有没有遇到过这种情况？
已经有8人回复
- 用1stopt或者origin非线性拟合磁导率
已经有5人回复
- 各位朋友，如何用matlab编写这个方程，并求解？
已经有15人回复
- origin8.0拟合方程中参数怎么求
已经有5人回复
- 求助 （已经-无量纲化）的非线性方程组（已经仔细检查，能算，但exitflag=0求高手
已经有4人回复
- 求高人指点用matlab求解非线性方程组，解决了追加100金币；
已经有11人回复
- 请高手谈谈，数据拟合，软件的优略（matlab，maple，mathematica，origin，1stopt等）
已经有43人回复
- 求助：origin8.0 自定义非线性拟合，参数拟合值与设定初始值一样
已经有18人回复
- 【求助】如何根据势能曲线拟合力场参数
已经有19人回复
- 【求助】求教matlab解非线性方程组
已经有9人回复
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15356.8
- 散金: 101
- 红花: 234
- 帖子: 3409
- 在线: 1222.7小时
- 虫号: 291104
- 注册: 2006-10-28
2楼2014-08-06 22:54:23
superl
银虫 (小有名气)
- 应助: 11 (小学生)
- 金币: 1293
- 散金: 214
- 帖子: 211
- 在线: 76.8小时
- 虫号: 772842
- 注册: 2009-05-17
- 性别: MM
- 专业: 同位素地球化学
|
你说的没错，你提到的几个参数都是为0，没有变化范围，只是为了查看方便，归在一起，这个问题应该不影响计算结果吧？应该是其它地方的问题，因为我运行了之后得到的结果都不满足我设定的 x11+x12+x13+x14+x15+x16=1; x21+x22+x23+x24+x25+x26=1; x31+x32+x33+x34+x35+x36=1; x41+x42+x43+x44+x45+x46=1; x51+x52+x53+x54+x55+x56=1; x61+x62+x63+x64+x65+x66=1; x71+x72+x73+x74+x75+x76=1; x81+x82+x83+x84+x85+x86=1; 相加为1的。 |
3楼2014-08-07 08:21:34
|
本帖内容被屏蔽 |
4楼2014-08-28 09:29:50
5