24小时热门版块排行榜    

Znn3bq.jpeg
查看: 452  |  回复: 2

zzq4816

新虫 (小有名气)

[求助] 1stopt拟合方程系数,始终算不出结果,希望有大神帮忙算一遍或者看看有没有问题 已有1人参与

Parameters B[0,0.00000001],C[0,0.02],D[0,1],k[0,1];
constant wc=0.4;
constant S=350;
constant pc=3.15,pw=1,r=1;
CONSTANT r0=3/S/10/pc;
Variable t,a,a1,b;
Function a1=3*b*pw*((wc/(1+wc)-0.4*a/(1+wc))/(1/(1+wc)))^r/0.4/r0/pc/((1/(B/a^1.5+C*a^3)-r0/(D*ln(1/a)))+r0/(D*ln(1/a))*(1-a)^(-1/3)+1/k*(1-a)^(-2/3));

data;
time        a        a1        b
2.22200        0.02354        0.00423        1
2.96433        0.02718        0.00571        1
3.64655        0.03169        0.00766        1
4.24017        0.03688        0.00997        1
4.91868        0.04462        0.01293        1
5.98748        0.06109        0.01784        1
7.30983        0.08811        0.02270        1
9.00767        0.12961        0.02549        1
10.32194        0.16311        0.02523        0.983823726
11.39535        0.18973        0.02430        0.963170866
12.46181        0.21506        0.02326        0.944167732
13.59083        0.24120        0.02291        0.92594501
14.63194        0.26423        0.02119        0.910733473
16.15073        0.29432        0.01849        0.891850486
17.33270        0.31509        0.01667        0.879008325
18.91573        0.33992        0.01476        0.863137203
20.16080        0.35753        0.01356        0.851951442
21.39448        0.37363        0.01256        0.842328242
22.37030        0.38553        0.01185        0.835057445
23.06710        0.39363        0.01138        0.830070292
23.62635        0.39990        0.01103        0.826191374
24.01731        0.40417        0.01079        0.823490749
24.49386        0.40925        0.01052        0.820555243
24.96176        0.41410        0.01025        0.81756759
25.50164        0.41956        0.00994        0.814313168
25.97504        0.42421        0.00969        0.811513533
26.53211        0.42952        0.00940        0.80837147
27.39851        0.43748        0.00896        0.803816902
28.13384        0.44394        0.00861        0.800133998
37.29421        0.50917        0.00616        0.761801069
46.82561        0.56032        0.00459        0.733040157
59.37877        0.60812        0.00311        0.707264692
73.41124        0.64271        0.00189        0.689310005
88.45445        0.66483        0.00113        0.678050343
100.84466        0.67647        0.00077        0.67216066
168.23129        0.69760        0.00002        0.661701517
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
zzq4816: 金币+15, ★★★★★最佳答案 2015-11-04 20:56:19
1stOpt不区分大小写,将代码中的大些B全部改成b1:
CODE:
Parameters B1=[0,0.00000001],C=[0,0.02],D=[0,1],k=[0,1];
constant wc=0.4;
constant S=350;
constant pc=3.15,pw=1,r=1;
CONSTANT r0=3/S/10/pc;
Variable t,a,a1,b;
Function a1=3*b*pw*((wc/(1+wc)-0.4*a/(1+wc))/(1/(1+wc)))^r/0.4/r0/pc/((1/(B1/a^1.5+C*a^3)-r0/(D*ln(1/a)))+r0/(D*ln(1/a))*(1-a)^(-1/3)+1/k*(1-a)^(-2/3));

data;
time        a        a1        b
2.22200        0.02354        0.00423        1
2.96433        0.02718        0.00571        1
3.64655        0.03169        0.00766        1
。。。。。

稳定唯一:

均方差(RMSE):0.0017563846590581
残差平方和(SSE):0.000111055934540687
相关系数(R): 0.970563265117913
相关系数之平方(R^2): 0.941993051596345
确定系数(DC): 0.941185544595634
卡方系数(Chi-Square): 0.00698513257358423
F统计(F-Statistic): 181.361309789069

参数        最佳估算
----------        -------------
b1        9.872444422668E-9
c        0.02
d        1.15322313130075E-9
k        1.28829258890287E-5
1stopt拟合方程系数,始终算不出结果,希望有大神帮忙算一遍或者看看有没有问题
c194.jpg

2楼2015-11-02 19:45:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zzq4816

新虫 (小有名气)

引用回帖:
2楼: Originally posted by dingd at 2015-11-02 19:45:36
1stOpt不区分大小写,将代码中的大些B全部改成b1:

Parameters B1=,C=,D=,k=;
constant wc=0.4;
constant S=350;
constant pc=3.15,pw=1,r=1;
CONSTANT r0=3/S/10/pc;
Variable t,a,a1,b;
Function a1=3*b ...

多谢啦,解决一大问题。。。
3楼2015-11-04 20:55:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zzq4816 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 中科院总分315求调剂 +3 lallalh 2026-04-09 3/150 2026-04-10 00:51 by haironglove
[考研] 调剂 +8 卷卷卷心菜_ 2026-04-09 8/400 2026-04-09 22:16 by lbsjt
[考研] 0703调剂 +18 拾玖壹 2026-04-04 20/1000 2026-04-09 20:53 by zhouxiaoyu
[考研] 材料化工总分334求调剂 +16 Riot2025 2026-04-08 17/850 2026-04-09 20:19 by maddjdld
[考研] 0854调剂 +10 长弓傲 2026-04-09 11/550 2026-04-09 19:03 by 探123
[考研] 工科273调剂 +4 X1999 2026-04-09 5/250 2026-04-09 16:05 by wp06
[考研] 285求调剂 +9 AZMK 2026-04-07 10/500 2026-04-09 10:29 by 小物理化学
[考研] 275 求调剂 +8 Lei812514 2026-04-07 8/400 2026-04-08 12:46 by chemisry
[考研] 22408 266求调剂 +11 masss11222 2026-04-07 14/700 2026-04-08 11:06 by yulian1987
[考研] 323求调剂 +3 林zlu 2026-04-07 4/200 2026-04-07 23:21 by lbsjt
[考研] 计算机11408 287 求调剂 +3 LiLe5 2026-04-07 3/150 2026-04-07 23:15 by shanqishi
[考研] 259求调剂 +5 就爱吃土豆呀呀 2026-04-07 5/250 2026-04-07 22:40 by JourneyLucky
[考研] 305求调剂 +4 77Qi 2026-04-06 4/200 2026-04-07 20:06 by shanqishi
[考研] 338求调剂 +4 我想上岸ii 2026-04-05 4/200 2026-04-06 21:04 by 木子君1218
[考研] 材料与化工363求推荐 +11 zh096 2026-04-04 11/550 2026-04-06 19:14 by guanxin1001
[考研] 求调剂到材料 +5 程9915 2026-04-06 5/250 2026-04-06 15:21 by yulian1987
[考研] 308求调剂 +3 终不似从前 2026-04-05 3/150 2026-04-05 20:07 by 啵啵啵0119
[考研] 321求调剂 +6 认真求上学 2026-04-03 6/300 2026-04-04 19:51 by dongzh2009
[考研] 一志愿沪985,326分求调剂 +3 刘墨墨 2026-04-03 3/150 2026-04-04 11:16 by 悲伤的芋头
[考研] 350一志愿北京航空航天大学08500材料科学与工程求调剂 +5 kjnasfss 2026-04-03 5/250 2026-04-03 22:29 by 无际的草原
信息提示
请填处理意见