24小时热门版块排行榜    

查看: 1048  |  回复: 5

hongyichen

新虫 (初入文坛)

[求助] 本人运用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版本并没有得到结果
还请各位能不能帮帮忙,帮我试着运行一下,感激不尽!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hongyichen

新虫 (初入文坛)

刚刚上传出现乱码。附件是数据

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : 程序数据.doc
  • 2017-07-03 22:14:26, 38.5 K
2楼2017-07-03 22:14:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

匿名

本帖仅楼主可见
3楼2017-07-04 00:27:51
已阅   申请计算强帖   回复此楼   编辑   查看我的主页

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
hongyichen: 金币+5, ★★★很有帮助 2017-07-04 17:16:07
把公式里的中括号“[ ]”全部换成小括号“( )”就行了
4楼2017-07-04 09:35:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hongyichen

新虫 (初入文坛)

引用回帖:
4楼: Originally posted by dingd at 2017-07-04 09:35:04
把公式里的中括号“”全部换成小括号“( )”就行了

恩恩谢谢你 可以了
5楼2017-07-04 17:13:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hongyichen

新虫 (初入文坛)

引用回帖:
3楼: Originally posted by ljk3709 at 2017-07-04 00:27:51
建议你设置悬赏金币,这样人很多看到的

恩恩谢谢你
6楼2017-07-04 17:16:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 hongyichen 的主题更新
信息提示
请填处理意见