求教非线性参数估计问题。

小弟需要估计下面函数的A，B两个参数吗，在论文看到1stOpt 1.5软件比较方便，就用了一下，但是运行下面程序提示“数据长度不够”，是什么意思？是数据个数不够吗？我应该修改软件参数的设置吗？Parameters A,B;

Variable FR,t;

Function FR=A*B*exp(-B*0.6/(0.8617*10^(-4)*(40+273)))*t^(B-1);

Data;

t=[14,14.5,15.5,16,16.5,17,17.5,18,18.5,19,19.5,20];

FR=[0.03965,0.04384,0.06766,0.06812,0.1219,0.1358,0.1964,0.2109,0.3078,0.3444,0.483,0.51];

如果各位有MATLAB 的参数估计程序也麻烦发一下，你之前用过的就行，不用针对我这个函数，我自己改一下试试。

Answer 1

[dojat]

加上过点约束即可：

Parameters A,B;
Variable FR,t;
Function FR=B*exp(-A*B*0.6/(0.8617*10^(-4)*(40+273)))*t^(B-1);
         14=B*exp(-A*B*0.6/(0.8617*10^(-4)*(40+273)))*0.038^(B-1);
Data;
14        0.03965
14.5        0.04384
15.5        0.06766
16        0.06812
16.5        0.1219
17        0.1358
17.5        0.1964
18        0.2109
18.5        0.3078
19        0.3444
19.5        0.483
20        0.51

均方差(RMSE):0.317369587414321

残差平方和(SSE):1.20868146018644

相关系数(R): 0.991036580671138

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

确定系数(DC): 0.970639316756605

卡方系数(Chi-Square): 0.0373432580787022

F统计(F-Statistic): 340.591529056107

约束条件: 14-(b*exp(0-a*b*0.6/0.02697121)*0.038^(b-1)) = 1.020081797E-009

参数 最佳估算

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

a -0.116966896721054

b 1.13775342075958