| 查看: 1159 | 回复: 5 | ||
[求助]
本人运用1stopt1.5版本没有输出结果 不知道存在什么问题 可以帮忙运行一下吗 ? 已有1人参与
|
|
求助非线性函数参数估计问题,有没有大神可以用1stopt的高级版本帮忙试试? 各位大神,本人最近用1stopt进行一个非线性函数的参数估计,要估计的参数为a,b两个, 程序和数据如下: Parameter a,b; Variable x1,x2,y; Function y=[a/(2sinh(2b))]*[exp(-a*abs(x1+x2)*cosh(b))*sinh(b+a*abs(x1+x2)*sinh(b))+ exp(-a*abs(x1-x2)*cosh(b))*sinh(b+a*abs(x1-x2)*sinh(b))]; Data; x1 x2 y 0.08 0.17 0.999563647085703 0.17 0.25 0.999489815898793 0.25 0.33 0.999401949466772 0.33 0.42 0.999297491237766 0.42 0.50 0.999173634119673 0.50 0.58 0.999027456079710 0.58 0.67 0.998856178168658 0.67 0.75 0.998657595686600 0.75 0.83 0.998430729911067 0.83 0.92 0.998176719156497 0.25 0.08 0.998110258380471 0.92 1.00 0.997899895715664 0.17 0.33 0.997787633668196 0.25 0.42 0.997403928198098 0.33 0.50 0.996948442830718 0.42 0.58 0.996409743916502 0.50 0.67 0.995776444312890 0.08 0.33 0.995388492667056 0.58 0.75 0.995038572863311 0.17 0.42 0.994595154277951 0.67 0.83 0.994189729241333 0.25 0.50 0.993652858942214 0.75 0.92 0.993230158977414 0.33 0.58 0.992536907133540 0.83 1.00 0.992170684362222 0.42 0.67 0.991221891983577 0.08 0.42 0.991094253083807 0.50 0.75 0.989684128090047 0.17 0.50 0.989553121643748 0.58 0.83 0.987905614393109 0.25 0.58 0.987726579832662 0.67 0.92 0.985879890990165 0.33 0.67 0.985570817307154 0.08 0.50 0.984862896397364 0.75 1.00 0.983619850707302 0.42 0.75 0.983043217613414 0.17 0.58 0.982234083855303 0.50 0.83 0.980108088003310 0.25 0.67 0.979128444302356 0.58 0.92 0.976745379464414 0.08 0.58 0.976262467553006 1.00 1.25 0.976039408543652 1.75 2.00 0.975969455881765 0.33 0.75 0.975480433759188 0.67 1.00 0.972962704623954 0.17 0.67 0.972137693216626 0.42 0.83 0.971231626433858 1.25 1.50 0.970119603709919 1.50 1.75 0.970095924444062 0.25 0.75 0.967286767789434 0.50 0.92 0.966342074925089 0.08 0.67 0.964789651463493 0.33 0.83 0.961624903548111 0.58 1.00 0.960806373521313 0.92 1.25 0.959895117357224 0.17 0.75 0.958691077393353 0.42 0.92 0.955087318167937 0.25 0.83 0.951562642779786 0.08 0.75 0.949870998736787 0.50 1.00 0.947648920168087 0.33 0.92 0.943311398733684 0.17 0.83 0.941257805855527 0.83 1.25 0.941222933255361 0.42 1.00 0.933887601453823 0.25 0.92 0.931270539294352 0.08 0.83 0.930872948769606 0.75 1.25 0.920830067508800 0.33 1.00 0.919831636284376 0.17 0.92 0.919159818882615 0.08 0.92 0.907124967407866 0.25 1.00 0.905716352430881 0.67 1.25 0.899395312926387 1.00 1.50 0.894081844388982 1.50 2.00 0.893894595769500 0.17 1.00 0.891716845550477 1.25 1.75 0.882219469440768 0.08 1.00 0.877960196453322 0.58 1.25 0.877465510299823 0.92 1.50 0.863192093812917 0.50 1.25 0.855464853122333 0.42 1.25 0.833710800325494 0.83 1.50 0.831145567170494 0.33 1.25 0.812432190392072 0.75 1.50 0.798703457101056 0.25 1.25 0.791786854237124 0.17 1.25 0.771877322084947 0.67 1.50 0.766465395852443 1.00 1.75 0.758634340206702 1.25 2.00 0.758426234298803 0.08 1.25 0.752764081737951 0.58 1.50 0.734880031129943 0.92 1.75 0.714841869855652 0.50 1.50 0.704265162988161 0.42 1.50 0.674831290200998 0.83 1.75 0.671333534377705 0.33 1.50 0.646704660053818 0.75 1.75 0.628778946433309 0.25 1.50 0.619947748751870 1.00 2.00 0.598447990189956 0.17 1.50 0.594576362386414 0.67 1.75 0.587662668866424 0.08 1.50 0.570573305230164 0.58 1.75 0.548308272789798 0.92 2.00 0.545301962655482 0.50 1.75 0.510908106237869 0.83 2.00 0.493719235665429 0.42 1.75 0.475553068541650 0.75 2.00 0.444253113547576 0.33 1.75 0.442259218818826 0.25 1.75 0.410989899820595 0.67 2.00 0.397260183529741 0.17 1.75 0.381673210483522 0.08 1.75 0.354215266710428 0.58 2.00 0.352935252791248 0.50 2.00 0.311348248646012 0.42 2.00 0.272478068445697 0.33 2.00 0.236241028388852 0.25 2.00 0.202513336549669 0.17 2.00 0.171148009574245 0.08 2.00 0.141987092562933 用的是LM算法,可是用1.5版本并没有得到结果 还请各位能不能帮帮忙,帮我试着运行一下,感激不尽! |
» 猜你喜欢
中国科学院物理研究所谌志国研究员团队招收2027年博士研究生
已经有6人回复
2026年中德博士后交流项目 - 新型量子和磁性材料:材料制备表征和中子散射研究
已经有12人回复
物理学I论文润色/翻译怎么收费?
已经有106人回复
26申博推荐:南京航空航天大学国际前沿院光学方向招收博士生!
已经有1人回复
如何从铁电相到顺电相。
已经有1人回复
各位大佬,求一份最新Wien2k版本作为学习用途
已经有0人回复
南科大活性流体和软物质课题组诚招2027级博士、硕士研究生和博后
已经有24人回复
瑞典林雪平大学博士后招聘|PEC Water Splitting 方向
已经有22人回复
2楼2017-07-03 22:14:54
3楼2017-07-04 00:27:51
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
4楼2017-07-04 09:35:04
5楼2017-07-04 17:13:22
6楼2017-07-04 17:16:26











回复此楼
