应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 1stopt拟合方程系数，始终算不出结果，希望有大神帮忙算一遍或者看看有没有问题
31/1返回列表
查看: 441  |  回复: 2
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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
回复此楼

» 猜你喜欢
1楼 2015-11-02 18:29:29
已阅   回复此楼   关注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
赞一下(1人) 回复此楼
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 的主题更新
31/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 本人考085602 化学工程 专硕 +21 不知道叫什么！ 2026-03-15 23/1150 2026-03-22 00:12 by BruceLiu320
[考研] 311求调剂 +13 冬十三 2026-03-15 14/700 2026-03-21 22:10 by peike
[考研] 266求调剂 +3 哇呼哼呼哼 2026-03-20 3/150 2026-03-21 16:46 by barlinike
[考研] 求调剂 +3 白QF 2026-03-21 3/150 2026-03-21 13:12 by zhukairuo
[考研] 085601调剂 358分 +3 zzzzggh 2026-03-20 4/200 2026-03-21 10:21 by luoyongfeng
[考研] 306求调剂 +4 chuanzhu川烛 2026-03-18 4/200 2026-03-21 08:25 by laoshidan
[考研] 316求调剂 +6 梁茜雯 2026-03-19 6/300 2026-03-21 06:32 by Ecowxq666！
[考研] 083200学硕321分一志愿暨南大学求调剂 +3 innocenceF 2026-03-17 3/150 2026-03-21 02:35 by JourneyLucky
[考研] 274求调剂 +10 S.H1 2026-03-18 10/500 2026-03-20 23:51 by JourneyLucky
[考研] 南京大学化学376求调剂 +3 hisfailed 2026-03-19 6/300 2026-03-20 23:43 by hisfailed
[考研] 295求调剂 +4 一志愿京区211 2026-03-18 6/300 2026-03-20 23:41 by JourneyLucky
[考研] 319求调剂 +3 小力气珂珂 2026-03-20 3/150 2026-03-20 19:47 by JourneyLucky
[考研] 求调剂 +3 eation27 2026-03-20 3/150 2026-03-20 19:32 by JourneyLucky
[考研] 广西大学家禽遗传育种课题组2026年硕士招生（接收计算机专业调剂） +3 123阿标 2026-03-17 3/150 2026-03-20 15:58 by 飞行琦
[考研] 085600材料与化工求调剂 +6 绪幸与子 2026-03-17 6/300 2026-03-19 13:27 by houyaoxu
[考研] 材料工程专硕调剂 +5 204818@lcx 2026-03-17 6/300 2026-03-18 22:55 by 204818@lcx
[考研] 277调剂 +5 自由煎饼果子 2026-03-16 6/300 2026-03-17 19:26 by 李leezz
[考研] 一志愿南京大学，080500材料科学与工程，调剂 +4 Jy? 2026-03-16 4/200 2026-03-17 11:02 by gaoqiong
[考研] 283求调剂 +3 听风就是雨； 2026-03-16 3/150 2026-03-17 07:41 by 热情沙漠
[考研] 机械专硕325，寻找调剂院校 +3 y9999 2026-03-15 5/250 2026-03-16 19:58 by y9999
信息提示
关闭
请填处理意见
关闭