xwfantasy 发表于 2011-5-13 21:57

谁有高版本1StOpt,帮帮忙

constant D=287 ,h=1.05457266*10^(-34),v=2700,kB=1.3806515*10^(-23);
ConstStr x=(h*w)/(kB*T),G=v/L+A*w^4+B*w^2*T*exp(-D/(3*T))+(C1*w^2)/((w01^2-w^2)^2)+ (C2*w^2)/((w02^2-w^2)^2);
Variable T,kl;
Parameter L,A,B,C1,C2,w01,w02;
Functionkl=(kB/(2*pi^2*v))*((kB*T)/h)^3*int((x^4*exp(x))/(G*(exp(x)-1)^2)*h/(KB*T),w=0,D*KB/h);
Data; //T   kl
30      5.16941
40      5.81634
50      6.43934
60      6.85724
70      6.97419
80      7.24853
90      7.23865
100   7.13575
110   7.16523
120   7.07047
130   7.02542
140   6.90191
150   6.76504
160   6.69633
170   6.36297
180   6.17548
190   6.07909
200   5.95913
210   5.80334
220   5.71028
230   5.59512
240   5.48558
250   5.37031
260   5.30022
270   5.20803
280   5.12605
290   5.09479
300   5.08835
哪位有高版本的1StOpt,帮我运行一下这个程序。万分感谢

20wangz 发表于 2011-6-26 07:11

回复 1 # xwfantasy 的帖子

好好检查一下公式,估计有问题,结果很不好。

dingd 发表于 2011-6-27 11:46

的确比较麻烦,下面是一组结果:
均方差(RMSE): 0.988621803497582
残差平方和(SSE): 27.3664459698227
相关系数(R): 0.899091983736487
相关系数之平方(R^2): 0.808366395219212
决定系数(DC): -0.686455406516489
F统计(F-Statistic): 2.07535875925088
参数         最佳估算
-------------------- -------------
l 1.16635690777819E-8
a 0
b 5.7136285006683E-18
c1 142.985235387459
c2 -55.8879107469693
w01 -1.33866548910737
w02 -15.9448648220859

xwfantasy 发表于 2011-6-29 21:52

谢谢啊!可能是实验数据误差太大了吧,公式是没问题的。你们的1STOPT是什么版本,装的话多少钱?
页: [1]
查看完整版本: 谁有高版本1StOpt,帮帮忙