7楼: Originally posted by dingd at 2016-03-13 21:47:01
效果不错啊：
均方差(RMSE):0.077118564339768
残差平方和(SSE):0.0773145485557502
相关系数(R): 0.99955257934791
相关系数之平方(R^2): 0.999105358881059
确定系数(DC): 0.999105298809933
F统计(F-Stati ...

6楼: Originally posted by ethan1129 at 2016-03-13 21:21:52
刚才的代码有点小问题，请帮忙试一下这个。
Variable x,y;
Parameter a,b,c,d;
Constant h=6.626196*10^(-34);
Constant f=3*10^(10);
Constant k=1.3806488*10^(-23);
Constant m=h*f/k;
Function y=a*exp ...

!!!using["IMSL","luopt","math"]; //使用命名空间
Constant(::h,f,k,m)= h=6.626196*10^(-34), f=3*10^(10), k=1.3806488*10^(-23), m=h*f/k;
jf(t)= 3.27*10^(-6)+3.25*10^(-9)*t;
f(a,b,c,d:i,s,x: h,f,k,m,xArray,max)=
{
   i=0,s=0,while{i<max,
     x=xArray(i,0),
     s=s+[a*exp(-3*d*QDAGS(@jf,203.0,x))+b*(1+2/(exp(a*m/(2*x))-1))+c*(1+3/(exp(a*m/(3*x))-1)+3/(exp(a*m/(3*x))-1)^2) - xArray(i,1)]^2,
     i++
   },
   s
};
main(::xArray,max)=
{
       xArray=matrix{   //存放实验数据
"
203        768.463
253        767.922
303        767.39
353        766.649
403        765.744
428        765.285
453        764.664
478        764.044
503        763.397
528        762.797
553        762.171
603        761.024
653        759.803
"
       },
    len[xArray,0,&max], //用len函数取矩阵的行数
    ClearImslErr(),     //清空IMSL错误输出
    ERSET(0,0,0), //关闭IMSL所有警告
    Opt[@f],      //Opt函数全局优化
    ERSET(0,2,2), ERSET(0,1,0)   //恢复IMSL警告
};