24小时热门版块排行榜    

查看: 1507  |  回复: 2

xtoffice

银虫 (小有名气)

[求助] istopt 拟合救助。 已有1人参与

Parameters Kc,Ks;
Variable x,y;
constant  J=0.0000245866;
ConstStr α=8/(27*β^3)+4*Ks/(3*β^3*Kc*J)-4*Ks^2*x/(3*β^3*Kc)
ConstStr β=(sprt(4/9+4*Ks/(3*Kc*J)+2*Ks^2*x/(3*Kc))
Function y=2/Ks*(β*cos(2π/3-1/3*arccos(α))+1/3) ;
data;
600        0.019145449
1200        0.032195992
1800        0.044269489
2400        0.05629646
3000        0.067416174
3600        0.078489362
4200        0.088352874
4800        0.098100071
5400        0.107800742
6000        0.117501413
6600        0.127155558
7200        0.136949281
7800        0.146696478
8400        0.156350623
拟合不了,是怎么回事。请高手指教。
回复此楼

» 猜你喜欢

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

皓小天

木虫之王 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
你自己细心点,sqrt写成了sprt肯定不对,而且beta的括号也有问题。修改如下:
Parameters Kc,Ks;
Variable x,y;
Constant J=0.0000245866;
ConstStr
b=sqrt(4/9+4*Ks/(3*Kc*J)+2*Ks^2*x/(3*Kc)),
a=8/(27*b^3)+4*Ks/(3*b^3*Kc*J)-4*Ks^2*x/(3*b^3*Kc);
Function y=2/Ks*(b*cos(2*pi/3-1/3*arccos(a))+1/3);
data;
600        0.019145449
1200        0.032195992
1800        0.044269489
2400        0.05629646
3000        0.067416174
3600        0.078489362
4200        0.088352874
4800        0.098100071
5400        0.107800742
6000        0.117501413
6600        0.127155558
7200        0.136949281
7800        0.146696478
8400        0.156350623
结果如下:
均方差(RMSE):0.00125590886248832
残差平方和(SSE):2.20822989922738E-5
相关系数(R): 0.999694045919579
相关系数之平方(R^2): 0.999388185447057
确定系数(DC): 0.999100686561502
卡方系数(Chi-Square): 0.000292076790042273
F统计(F-Statistic): 13343.5123798838

参数        最佳估算
----------        -------------
kc        1.32011205386087E-21
ks        3.44981108096092

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

No        实测值y        计算值y
1        0.019145449        0.0167668097475268
2        0.032195992        0.0303002456751706
3        0.044269489        0.0426481130554661
4        0.05629646        0.0553553257850073
5        0.067416174        0.0667882760896434
6        0.078489362        0.0785311810445057
7        0.088352874        0.0889152998412069
8        0.098100071        0.0995630840138755
9        0.107800742        0.108771270846651
10        0.117501413        0.118199568438729
11        0.127155558        0.127844491228784
12        0.136949281        0.137702694619866
13        0.146696478        0.146003126917595
14        0.156350623        0.154476897293266
请不要站内找我要书,如果需要请到书籍板块求助
2楼2015-08-31 21:25:24
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

xtoffice

银虫 (小有名气)

引用回帖:
2楼: Originally posted by 皓小天 at 2015-08-31 21:25:24
你自己细心点,sqrt写成了sprt肯定不对,而且beta的括号也有问题。修改如下:
Parameters Kc,Ks;
Variable x,y;
Constant J=0.0000245866;
ConstStr
b=sqrt(4/9+4*Ks/(3*Kc*J)+2*Ks^2*x/(3*Kc)),
a=8/(27*b^3 ...

感谢,问题已解决。
3楼2015-08-31 22:11:13
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xtoffice 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 G6APbkg8SA6w 2026-08-29 5/250 2026-08-31 01:51 by hZiFeudZyoGR
[考研] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 gy1nBQXYQJqL 2026-08-29 5/250 2026-08-31 01:38 by hZiFeudZyoGR
[考博] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 4FFAWE8HcgUD 2026-08-29 4/200 2026-08-30 20:40 by hZiFeudZyoGR
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 4FFAWE8HcgUD 2026-08-29 3/150 2026-08-30 20:19 by hZiFeudZyoGR
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 jCd0dEvKHShX 2026-08-29 5/250 2026-08-30 19:48 by hZiFeudZyoGR
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +5 jCd0dEvKHShX 2026-08-29 7/350 2026-08-30 19:45 by hZiFeudZyoGR
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 jCd0dEvKHShX 2026-08-29 6/300 2026-08-30 19:13 by hZiFeudZyoGR
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-30 17:47 by hZiFeudZyoGR
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +6 ASdOkHsho7FD 2026-08-28 9/450 2026-08-30 11:19 by l0VvVHGBGRLv
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 申请删除本帖 +6 lyz123lyz 2026-08-27 7/350 2026-08-27 17:31 by 宁静致远sy
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
信息提示
请填处理意见