α β γ δ 2030 13 320 1 1850 12 320 1 1850 12 320 1 1970 13 320 1 2090 14 320 1 2210 15 320 1 3010 20 320 1 2540 17 320 1 2710 18 620 2p=[111110.909909910.090090090.5443037970.455696203];需要对图示的公式进行标定，式中a,b,c,d为需要进行拟合的参数。p，α_j……δ_j为已知的调查样本数据,如上面所示。需要将各参数进行标定计算，并评价标定结论。恳请各位高手帮忙，给出标定程序，最好能够有注释说明！谢了
%by dynamic%see also http://www.matlabsky.com%contact me matlabsky@gmail.com%2009.2.%程序编好了，但是初值真的很难选，我到现在还没有选这个一个合适的初值，你根据实际选择一个吧，否则程序是运行不起来的你的函数exp，初值稍微选的不好就会导致0或者∞function mmclcx=[2030 13 320 1 1850 12 320 1 1850 12 320 1 1970 13 320 1 2090 14 320 1 2210 15 320 1 3010 20 320 1 2540 17 320 1 2710 18 620 2];p=[1 1 1 1 1 0.90990991 0.09009009 0.544303797 0.455696203]; b0=[-.005 2 -0.05 8];b = nlinfit(x,p,@fun,b0)function yhat=fun(bb,x)alpha=x(:,1);beta=x(:,2);gamma=x(:,3);delta=x(:,4);a=bb(1);b=bb(2);c=bb(3);d=bb(4);num=exp(a*alpha+b*beta+c*gamma+d*delta);den=sum(num);yhat=num./den;
|
|
