| 查看: 455 | 回复: 2 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
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 |
回复此楼» 猜你喜欢
B区学生调剂-兰州交通大学材料科学与工程学院
已经有12人回复
-
山西大同大学物理学专业还有调剂名额，欢迎调剂！
已经有13人回复
-
物理学I论文润色/翻译怎么收费?
已经有248人回复
-
桂林理工大学物理学专业招收调剂
已经有18人回复
-
VASP 的一组 GPU / CPU 基准测试记录
已经有0人回复
-
津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生一名
已经有0人回复
-
【原创讨论】从电子约束到物质编辑：一套可迭代的环形磁场科技树
已经有0人回复
-
【方案分享】单环磁场+轴心控制+偏转导出电子束约束系统（可行性实验）
已经有6人回复
-
【修正版】单环用磁约束低速电子实验方案（简化版）
已经有0人回复
-
桂林理工大学物理学专业招收调剂，还有三个名额！！！
已经有9人回复
3楼2015-11-04 20:55:00
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
zzq4816: 金币+15, ★★★★★最佳答案 2015-11-04 20:56:19
感谢参与，应助指数 +1
zzq4816: 金币+15, ★★★★★最佳答案 2015-11-04 20:56:19
|
1stOpt不区分大小写，将代码中的大些B全部改成b1： 稳定唯一: 均方差(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  c194.jpg |
2楼2015-11-02 19:45:36
20