| 查看: 1048 | 回复: 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版本并没有得到结果 还请各位能不能帮帮忙,帮我试着运行一下,感激不尽! |
» 猜你喜欢
求标准粉末衍射卡号 ICDD 01-076-1802
已经有0人回复
新西兰Robinson研究所招收全奖PhD
已经有0人回复
物理学I论文润色/翻译怎么收费?
已经有112人回复
石墨烯转移--二氧化硅衬底石墨烯
已经有0人回复
笼目材料中量子自旋液体基态的证据
已经有0人回复
数学教学论硕士可以读数学物理博士吗?
已经有0人回复
德国亥姆霍兹Hereon中心汉堡分部招镁合金腐蚀裂变SCC课题方向2026公派博士生
已经有4人回复
澳门大学 应用物理及材料工程研究院 潘晖教授课题组诚招博士后
已经有11人回复
求助NH4V4O10晶体的CIF文件
已经有0人回复
英国全奖博士招聘-深度学习与量子物理
已经有0人回复
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












回复此楼
