洛伦兹曲线参数估计求​助,知道x,y求参数​b1、b2、b3。

x=[0.1:0.1:0.9,0.95];

y=[0.025, 0.0673,0.1221,0.1882,0.2663,0.3569,0.4631,0.5901,0.7485,0.8493];

plot(x,y,'*'),xlabel('p'),ylabel('L'),title('A地区，年份之一的洛伦兹曲线')

hold on;

b0=[0.2,0.2,0.5]; %b初始值

fun=inline('(1-(1-x)^b(1))/(1-b(2)*log(x))^b(3)','b','x');

%将b看成参考量，b(i)为分量，本意是求：(1-（1-p）^b1)/(1-(1/b2)*log(x))^b3

[b,r,j]=nlinfit(x,y,fun,b0);

b %最佳参数

R=sun(r.^2) %误差平方和

Answer 1

[yoxer]

均方差(RMSE):0.00104754738219215

残差平方和(SSE):1.09735551793764E-5

相关系数(R): 0.999993102934655

相关系数之平方(R^2): 0.999986205916879

确定系数(DC): 0.999985292738451

卡方系数(Chi-Square): 5.65156168961001E-5

F统计(F-Statistic): 237977.681190137

参数 最佳估算

---------- -------------

b1 0.688826282349568

b2 0.36058678259219

b3 1.60970899512608