Parameters d, beta1, fai, det_p, det_pu, N;
ComplexStr = i;
ConstStr Fc=exp(-i*fai*d*q)*exp(-0.5*(fai*beta1*d)^2*q^2),
         Fa=exp(-i*(1-fai)*d*q)*exp(-0.5*((1-fai)*beta1*d)^2*q^2),
         R1=(2*det_p^2/q^2)*(((1-Fc)*(1-Fa))/(1-Fc*Fa)+(1/N)*Fa*((1-Fc)^2/(1-Fa*Fc)^2)*(1-(Fa*Fc)^N)),
         R2=(2*det_p^2)/(q^2*N)*((det_pu^2)/(det_p^2)*(1-(Fc*Fa)^N)-(det_pu/det_p)*((1+Fa)*((1-Fc)/(1-Fc*Fa)*(1-(Fc*Fa)^N))));
Variable q, Itotal;
Function Itotal=(R1+R2)/(4*pi*q^2);
Data;
q        Itotal
0.01128        0.05389
0.01292        0.04665
0.01544        0.03979
0.02042        0.02857
0.02628        0.02304
0.03208        0.01848
0.04002        0.016
0.05125        0.01713
0.05585        0.01943
0.05957        0.02038
0.06707        0.01695
0.07123        0.01465
0.07709        0.01142
0.08295        0.00894
0.09209        0.00629
0.09998        0.00476

2楼: Originally posted by dingd at 2013-11-08 14:38:30
Parameters d, beta1, fai, det_p, det_pu, N;
ComplexStr = i;
ConstStr Fc=exp(-i*fai*d*q)*exp(-0.5*(fai*beta1*d)^2*q^2),
         Fa=exp(-i*(1-fai)*d*q)*exp(-0.5*((1-fai)*beta1*d)^2*q^2),
        ...

3楼: Originally posted by 阿德 at 2013-11-08 15:59:02
非常感谢，但是参数不太对噢
我忘了设定限制范围了，d大于0，beta（0-1）, det_p(0-1),det_pu(0-1), N为大于零的自然数，谢谢...

5楼: Originally posted by dingd at 2013-11-08 20:24:04
均方差(RMSE):0.0025341163115579
残差平方和(SSE):0.000102747927688061
相关系数(R): 0.982569988533571
相关系数之平方(R^2): 0.965443782366861
确定系数(DC): 0.965413209513977
F统计(F-Statistic): 57. ...

7楼: Originally posted by dingd at 2013-11-10 11:23:31
仔细回看你自己前面的帖子，fai和beta的范围都是在0到1之间。