24小时热门版块排行榜    

查看: 1762  |  回复: 7

kaiquan

金虫 (小有名气)

[求助] 1stopt 拟合隐函数的四参数

1stopt 拟合隐函数的四个参数,式子有点长,各位大侠看看,急求,在附件里 谢谢了
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : 数据拟合.doc
  • 2013-11-12 09:39:56, 17.5 K

» 猜你喜欢

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

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

feixiaolin

荣誉版主 (文坛精英)

优秀版主

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感谢参与,应助指数 +1
kaiquan: 金币+60, ★★★很有帮助 2013-11-12 10:59:50
csgt0: 金币+1, 谢谢 2013-11-12 13:51:21
相关系数(R): 0.99962478884682
相关系数之平方(R^2): 0.999249718477049
决定系数(DC): 0.999204261880148
卡方系数(Chi-Square): 1.173297317121E-5
F统计(F-Statistic): 9322.83125116381

参数                最终估算
----------                --------------
a2                 0.617537237985605
b2                 -96.0309866048836
a1                 0.452723175481844
b1                 -85.1157994040171

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

No        实测值x        计算值x
1        0.1625        0.1633055
2        0.1696        0.1698433
3        0.1788        0.1789812
4        0.1866        0.1865306
5        0.193        0.1920540
6        0.2017        0.2010975
7        0.2084        0.2078759
8        0.2162        0.2155998
9        0.2558        0.2573601

» 本帖已获得的红花(最新10朵)

2楼2013-11-12 10:04:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kaiquan

金虫 (小有名气)

送红花一朵
引用回帖:
2楼: Originally posted by feixiaolin at 2013-11-12 10:04:20
相关系数(R): 0.99962478884682
相关系数之平方(R^2): 0.999249718477049
决定系数(DC): 0.999204261880148
卡方系数(Chi-Square): 1.173297317121E-5
F统计(F-Statistic): 9322.83125116381

参数                最终估算 ...

您好,利用您拟合的四个参数,我刚才用第二组数据带入算了一下,等式左右两边相差很大啊,能不能把您的程序发一下啊,谢谢了啊
天道酬勤
3楼2013-11-12 10:42:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

feixiaolin

荣誉版主 (文坛精英)

优秀版主

【答案】应助回帖


fegg7502: 金币+1, 鼓励交流 2013-11-12 15:04:28
1stop代码
Parameters a1, b1, a2, b2;
Variable  T, x;
Function
23663.94-26.95*T-ln(x)=(1-x)^2*(0.1203*(a2+b2/T)/(x*(1-exp(-0.036*(a2+b2/T)))+exp(-0.036*(a2+b2/T)))^2+0.1203*(a1+b1/T)/ (x*(exp(-0.036*(a1+b1/T))-1)+1)^2);
data ;
297.95   0.1625
305.25   0.1696
311.25    0.1788
318.05    0.1866
328.95    0.193
333.15    0.2017
338.05    0.2084
345.55    0.2162
364.75    0.2558
4楼2013-11-12 10:53:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

feixiaolin

荣誉版主 (文坛精英)

优秀版主

【答案】应助回帖


fegg7502: 金币+1, 应助指数+1 2013-11-12 15:04:36
另外一组结果:
相关系数(R): 0.999999659568989
相关系数之平方(R^2): 0.999999319138094
决定系数(DC): 0.999277451092885
卡方系数(Chi-Square): 1.08114889392315E-5
F统计(F-Statistic): 10281079.2766421

参数        最佳估算
----------        -------------
a1                 1.26684964120538
b1                 1.8734535394872
a2                 -1.24452004120401
b2                 0.761212713683848

====== 输出结果 =====
No.         Observed x        Calculated x
1        0.1625        0.161593731597778
2        0.1696        0.168901236690022
3        0.1788        0.178358643972619
4        0.1866        0.186377729176926
5        0.193        0.192965570938433
6        0.2017        0.201897501847218
7        0.2084        0.208777207357754
8        0.2162        0.216787104934144
9        0.2558        0.25738900663584
5楼2013-11-12 10:59:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kaiquan

金虫 (小有名气)

引用回帖:
5楼: Originally posted by feixiaolin at 2013-11-12 10:59:22
另外一组结果:
相关系数(R): 0.999999659568989
相关系数之平方(R^2): 0.999999319138094
决定系数(DC): 0.999277451092885
卡方系数(Chi-Square): 1.08114889392315E-5
F统计(F-Statistic): 10281079.276642 ...

麻烦发一下代码,谢谢了
天道酬勤
6楼2013-11-12 11:08:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖


csgt0: 金币+1, 谢谢 2013-11-12 13:51:32
fegg7502: 应助指数+1, 专家考核, 鼓励交流 2013-11-12 15:04:49
看下这组:

均方差(RMSE): 0
残差平方和(SSE): 0
相关系数(R): 1
确定系数(DC): 1

参数        最佳估算
----------        -------------
a2        4.86682498921786E-8
b2        9.74676704581638E-6
a1        -4.86682497474015E-8
b1        -9.74676726037063E-6

» 本帖已获得的红花(最新10朵)

7楼2013-11-12 12:36:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kaiquan

金虫 (小有名气)

送红花一朵
引用回帖:
7楼: Originally posted by dingd at 2013-11-12 12:36:46
看下这组:

均方差(RMSE): 0
残差平方和(SSE): 0
相关系数(R): 1
确定系数(DC): 1

参数        最佳估算
----------        -------------
a2        4.86682498921786E-8
b2        9.74676704581638E-6
a1        -4.86682497474015E-8
...

不好意思啊,这个附件的式子上传时候少了点内容,新上传的式子更复杂些,谢谢您了啊
天道酬勤
8楼2013-11-12 14:17:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 kaiquan 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 jCd0dEvKHShX 2026-08-29 4/200 2026-08-30 14:05 by GkCVyODlx4mR
[基金申请] 有没有仍没收到信息的 +6 德尚中行 2026-08-27 7/350 2026-08-30 12:53 by Bxndjjd
[基金申请] 国社科又开始会评了,不知道这次命运如何 +3 雨打竹帘 2026-08-30 4/200 2026-08-30 12:35 by 白大多亮
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 ASdOkHsho7FD 2026-08-28 6/300 2026-08-30 11:18 by l0VvVHGBGRLv
[论文投稿] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 ASdOkHsho7FD 2026-08-28 6/300 2026-08-30 11:07 by l0VvVHGBGRLv
[基金申请] 29号明天会评吗 +3 笨笨唐 2026-08-28 3/150 2026-08-30 09:29 by cww8181
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +6 gdd2018 2026-08-28 11/550 2026-08-30 08:58 by 超级无敌华子
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 gy1nBQXYQJqL 2026-08-29 5/250 2026-08-30 07:05 by ZPa0EcMwuECS
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 gy1nBQXYQJqL 2026-08-29 4/200 2026-08-30 06:57 by ZPa0EcMwuECS
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +7 ASdOkHsho7FD 2026-08-28 11/550 2026-08-30 05:40 by ZPa0EcMwuECS
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 gy1nBQXYQJqL 2026-08-29 4/200 2026-08-30 01:34 by ZPa0EcMwuECS
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 基金未中,这种答复是模板吗? +5 zhaosm1982 2026-08-27 6/300 2026-08-27 16:00 by lfy8008
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
信息提示
请填处理意见