上面是我的结果,我把程序和部分数据拿出来。。
Parameters a,b,c,d,e,f,g,h,w;
Variable s1,x1,x2,p,s2,x3,x4 ;
Function p=a*ln(s1/s2) +b*(x1-x3)+c*(x2-x4)+ln(1+exp(-w*(a*ln(s1)+d*x1+e*x2-f*s1-g*x1-h*x2))/(1+exp(-w*(a*ln(s2)+d*x3+e*x4-f*s2-g*x3-h*x4)))) ;
data;
1169.71 951852 3.6 -0.89942 6693.38 956620 0
1540.08 504066 4.0 -1.89169 6693.38 956620 0
507.44 331209 6.5 -2.66831 6693.38 956620 0
1142.49 618975 7.8 -2.40685 6693.38 956620 0
832.16 434457 5.8 -4.52618 6693.38 956620 0
68.76 135101 10.3 -4.98759 6693.38 956620 0
643.36 295955 12.1 -3.46358 6693.38 956620 0
910.57 2253856 9.1 -3.27470 6693.38 956620 0
193.39 505752 15.1 -4.43021 6693.38 956620 0
810.51 562840 16.0 -6.87665 6693.38 956620 0
599.77 351826 14.3 -5.68216 6693.38 956620 0
1019.88 703859 20.0 -5.70528 6693.38 956620 0
1464.57 590200 22.6 -6.06363 6693.38 956620 0
2157.56 1599877 20.5 -4.98759 6693.38 956620 0
1475.31 741276 23.5 -5.55113 6693.38 956620 0
996.82 287472 15.2 -5.13486 6693.38 956620 0
884.73 671137 13.8 -4.38080 6693.38 956620 0
757.39 467801 14.3 -4.30881 6693.38 956620 0
1336.28 764069 23.0 -4.99474 6693.38 956620 0
415.16 641804 21.0 -5.88061 6693.38 956620 0
527.16 640279 32.7 -7.77889 6693.38 956620 0
419.46 391710 15.5 -6.51621 6693.38 956620 0
293.71 135434 26.6 -6.15560 6693.38 956620 0
6693.38 956620 3.6 -1.72870 1169.71 951852 0
1540.08 504066 6.0 -3.01030 1169.71 951852 0
507.44 331209 7.0 -3.66981 1169.71 951852 0
1142.49 618975 9.7 -3.61823 1169.71 951852 0
832.16 434457 8.8 -5.40551 1169.71 951852 0
68.76 135101 8.0 -5.33563 1169.71 951852 0
643.36 295955 13.2 -5.40551 1169.71 951852 0
910.57 2253856 7.8 -4.08143 1169.71 951852 0