24小时热门版块排行榜    

查看: 837  |  回复: 4

rainman0147

新虫 (初入文坛)

[求助] 6参数的一个方程,带约束的,求1stopt代跑程序,不甚感激!! 已有1人参与

6参数的一个方程,然后6个要需要满足一个方程,低版本的1stopt无法拟合带约束的,高版本试用版的又只能做4参数的,哪位大神有正版的高版本的1stopt代小弟跑一下程序,不甚感激!!!
TiTle"test001";
Parameters A(1:6);
Variable x,y;
Function y=A1*((1+(1-exp(-x)))^(A2-1)-(1+(1-exp(-x)))^(-0.5*A2-1))+A3*((1+(1-exp(-x)))^(A4-1)-(1+(1-exp(-x)))^(-0.5*A4-1))+A5*((1+(1-exp(-x)))^(A6-1)-(1+(1-exp(-x)))^(-0.5*A6-1));
            A1*A2+A3*A4+A5*A6=2;
Data;
具体程序数据见附件
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : test001.txt
  • 2015-11-18 22:05:57, 16.7 K

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
数据太长,费时间,均匀取了50个点,结果供参考:

均方差(RMSE):0.111516434470976
残差平方和(SSE):0.646667588170212
相关系数(R): 0.999459849877
相关系数之平方(R^2): 0.998919991516156
确定系数(DC): 0.998762787065796
卡方系数(Chi-Square): 0.26117331567502
F统计(F-Statistic): 7436.06839668255
约束条件: a1*a2+a3*a4+a5*a6-2 = 6.766143201E-012

参数        最佳估算
----------        -------------
a1        5.08176004574493
a2        -12.5094351345345
a3        -7.99866254026298
a4        -6.88350829142222
a5        1.2084077175086
a6        8.69829577636614

====== 结果输出 ======

No.        实测值 y        计算值 y
1        0.00474        0.00672334878016531
2        0.00867        0.0106119135708815
3        0.24257        0.147811932566774
4        0.58853        0.370190041910542
5        0.90481        0.627791500420167
6        1.18267        0.893409636179607
7        1.42377        1.15113365316003
8        1.63259        1.39271533505806
9        1.81329        1.61440535348283
10        1.96963        1.81493691256894
11        2.10266        1.99506431823239
12        2.21877        2.15651003806406
13        2.3204        2.3013376926177
14        2.41456        2.43213268256514
15        2.50294        2.55156697437291
16        2.58787        2.66214979768776
17        2.67222        2.76649973828385
18        2.75896        2.86712144383152
19        2.84857        2.96625317270316
20        2.94391        3.06614058140696
21        3.04627        3.16889008891661
22        3.15712        3.27632031627744
23        3.27954        3.39025571212788
24        3.41364        3.51239730980969
25        3.5617        3.64413893330009
26        3.68484        3.75190296002894
27        3.6929        3.75885128870938
28        3.72544        3.78691901005312
29        3.90364        3.94207467909228
30        4.09857        4.11060172615032
31        4.30689        4.29360087002512
32        4.53102        4.49209198121528
33        4.76671        4.70670458714309
34        5.01752        4.93824765066839
35        5.27943        5.18746695208385
36        5.55617        5.45465770415145
37        5.84416        5.74037893851692
38        6.14764        6.04514176156981
39        6.46455        6.368939142335
40        6.7956        6.71211869361832
41        7.14388        7.07499339783819
42        7.50622        7.45728655180331
43        7.89096        7.85917026698921
44        8.29231        8.28079363396696
45        8.71462        8.72164328573571
46        9.15998        9.18174938949315
47        9.62134        9.66112915131039
48        10.1107        10.1590655380355
49        10.62284        10.6754778673183
50        11.15838        11.2102810683146
51        11.72316        11.7625868000411
52        11.99894        12.0318650815961
2楼2015-11-19 09:30:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

rainman0147

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by dingd at 2015-11-19 09:30:23
数据太长,费时间,均匀取了50个点,结果供参考:

均方差(RMSE):0.111516434470976
残差平方和(SSE):0.646667588170212
相关系数(R): 0.999459849877
相关系数之平方(R^2): 0.998919991516156
确定系数(DC): ...

谢谢您了!虽然我很想要全部数据拟合得到的结果,但是我知道这个工作量很大,也不好意思麻烦您了,对您的帮助我由衷的感谢!
3楼2015-11-19 21:40:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

rainman0147

新虫 (初入文坛)

引用回帖:
3楼: Originally posted by rainman0147 at 2015-11-19 21:40:37
谢谢您了!虽然我很想要全部数据拟合得到的结果,但是我知道这个工作量很大,也不好意思麻烦您了,对您的帮助我由衷的感谢!...

您好!刚才把您算得参数代入计算了一下,发现拟合得挺好!谢谢您了!我这边还有最后一个也需要拟合一下,不知道您能不能帮忙拟合一下,谢谢您了!

» 本帖附件资源列表

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

rainman0147

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by dingd at 2015-11-19 09:30:23
数据太长,费时间,均匀取了50个点,结果供参考:

均方差(RMSE):0.111516434470976
残差平方和(SSE):0.646667588170212
相关系数(R): 0.999459849877
相关系数之平方(R^2): 0.998919991516156
确定系数(DC): ...

您好!刚才把您算得参数代入计算了一下,发现拟合得挺好!谢谢您了!我这边还有最后一个也需要拟合一下,不知道您能不能帮忙拟合一下,谢谢您了!

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : test0001.txt
  • 2015-11-19 22:10:04, 16.6 K
5楼2015-11-19 22:10:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 rainman0147 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] filecode +7 布布和一二 2026-08-06 10/500 2026-08-06 16:22 by janehuiminli
[基金申请] 听说今天filecode变了 +3 布布和一二 2026-08-06 3/150 2026-08-06 16:22 by yq199112088
[基金申请] filecode +8 等待解的谜 2026-08-06 12/600 2026-08-06 16:11 by janehuiminli
[基金申请] 大家散了吧,后缀研究没有意义,别浪费时间了,过好目前的每一天,不要焦虑 +4 Tide man 2026-08-06 5/250 2026-08-06 16:08 by janehuiminli
[有机交流] 一个有机合成实验室都需要哪些设备? 50+3 kf2781974 2026-07-31 12/600 2026-08-06 15:11 by eddyin
[基金申请] 求各位大神看下 100+6 hpkpkpkp 2026-08-05 33/1650 2026-08-06 14:49 by zhiyanjiang
[教师之家] 咨询面上基金 +4 李长云 2026-07-31 7/350 2026-08-06 11:20 by 李长云
[基金申请] 影响面上的因素 +8 布布和一二 2026-08-05 11/550 2026-08-06 10:41 by 宝贝虫子
[基金申请] 8月时间戳变的,举个手。玩一下,释放压力 +9 archvillain 2026-08-04 11/550 2026-08-05 20:06 by wlwhappy
[基金申请] 面上再次挂了,太难了,躺也躺不了,倦也卷不过,小学校之殇! +22 低垂的野花 2026-07-31 30/1500 2026-08-05 18:03 by 低垂的野花
[基金申请] 好消息?这个有何含义??? +8 Tide man 2026-08-05 10/500 2026-08-05 16:14 by xmuxiaoyu
[考博] 【2027博士申请】纳米药物递送方向 20+3 13586093586 2026-08-03 4/200 2026-08-05 09:59 by lfy8008
[基金申请] 有没有H口的?有收到消息的吗? +3 超级海虾 2026-08-04 3/150 2026-08-04 17:26 by 学教育滴
[论文投稿] 十年后又回来了,论文投稿求助 +3 哈哈114477 2026-08-01 3/150 2026-08-04 15:40 by tegsgjy20
[基金申请] 纯娱乐,不喜欢勿喷 +7 Tide man 2026-08-04 10/500 2026-08-04 15:10 by loufangrui
[基金申请] 面上提前没消息,有中的吗 +14 archvillain 2026-08-02 18/900 2026-08-04 14:42 by archvillain
[基金申请] 娱乐 +4 Tide man 2026-08-03 4/200 2026-08-04 11:51 by wgch518
[基金申请] 什么时候能放榜呀? +3 Jacob678 2026-08-03 3/150 2026-08-03 16:14 by gltch
[基金申请] 2026年国自然面上资助率 +16 布布和一二 2026-07-30 22/1100 2026-08-03 09:20 by gy116024
[基金申请] 你们的时间戳变了吗 +3 archvillain 2026-07-30 4/200 2026-07-30 18:53 by levinzhwen
信息提示
请填处理意见