1stopt拟合复杂曲线求救

uglydu

实在不知道怎么表达了，还请各位高手指点下，其中变量Y只有86，132，164，184四种.而b和Y之间的关系在表达式下面。
谢谢了 鞠躬。。。
Variable x,y,z ;
Parameter a1[0,1],a2[0,],a3[0,],a4[0,],a5[0,],a6[0,],a7[0,],a8[0,],a9[0,],a10[0,],
Constant y=86,b=0.02271,y=132,b=0.03461,y=164,b=0.01487,y=184,b=0.04593;
Function z=(1-a1*(y/0.0001)^a2*x^a3*(x/b)^a10)*(a4*x+a5*x^2+a6*a7*y*(1-exp(-x/(y*a7)))+a8*a9*y*(1-exp(-x/(y*a9))));
data；
0，86，0
2E-05，86，0
4E-05，86，0.25023
0.00926，86，2.41996
0.00934，86，2.43096
0.00942，86，2.4422
0.02051，86，4.43271
0.0206，86，4.45378
0.02069，86，4.47387
0.02077，86，4.4916
0.02621，86，3.28568
0.02627，86，3.23328
0.02632，86，3.10092
.....
0，132，0
3E-05，132，0.32901
6E-05，132，0.37437
0.0001，132，0.43106
0.00014，132，0.48718
......
0,164,0
2E-05,164,0.2341
4E-051640.2598
0.00302,164,1.7294
0.00331,164,1.80163
0.00361,164,1.87243

........
0,184,0
2E-05,184,0.77194
4E-05,184,0.84217
7E-05,184,0.93125
0.0001,184,1.01095
0.00792,184,3.02387
0.00837,184,3.06403
0.00881,184,3.10349
0.00926,184,3.14383
0.00969,184,3.18393

......


[ 本帖最后由 uglydu 于 2008-12-17 15:09 编辑 ]

shamohu

可以象下面在编程模式下实现：

Variable x,y,z ;
Parameter a1[0,1],a2[0,],a3[0,],a4[0,],a5[0,],a6[0,],a7[0,],a8[0,],a9[0,],a10[0,];
StartProgram [Pascal];
Procedure MainModel;
var j: integer;
    b: double;
Begin
    for j := 0 to DataLength - 1 do begin
        if y[j] =86 then b := 0.02271
        else if y[j] = 132 then b := 0.03461
        else if y[j] = 164 then b := 0.01487
        else b := 0.04593;
        z[j] := (1-a1*(y[j]/0.0001)^a2*x[j]^a3*(x[j]/b)^a10)*(a4*x[j]+a5*x[j]^2+
                a6*a7*y[j]*(1-exp(-x[j]/(y[j]*a7)))+a8*a9*y[j]*(1-exp(-x[j]/(y[j]*a9))));
    end;
End;
EndProgram;
Data;
0        86        0
2E-05        86        0
4E-05        86        0.25023
0.00926        86        2.41996
0.00934        86        2.43096
0.00942        86        2.4422
0.02051        86        4.43271
0.0206        86        4.45378
0.02069        86        4.47387
0.02077        86        4.4916
0.02621        86        3.28568
0.02627        86        3.23328
0.02051        164        5.54367
0.02066        164        5.59151
0.02081        164        5.63396
......

uglydu

好的 谢谢 我立刻去试试

uglydu

好的 谢谢 我立刻去试试


可是按下F9运行后没反应呀。。。

shamohu

如果用的是网上的PJ版，是无法运行上面代码的。

uglydu

那怎么办呢？这些放进matlab又出错，恩，我知道可能两者代码不通用吧，哎。

wanglei5201118

迭代数: 1653
计算用时(时:分:秒:微秒): 00:00:09:891
优化算法: 最大继承法
计算结束原因: 达到收敛判定标准
均方差(RMSE): 0.0688542932138693
残差平方和(SSE): 0.0711137054097224
相关系数(R): 0.999463892097194
相关系数之平方(R^2): 0.998928071606071
决定系数(DC): 0.998693996385168
F统计(F-Statistic): 12114.675760456

参数                  最佳估算
--------------------        -------------
a1                 0.999999999317417
a2                 9.56176694986926
a3                 36.141540391217
a4                 5.1975731720829E-8
a5                 2891.5050234386
a6                 346.710381325455
a7                 0.000127847356618798
a8                 3.07679055320162E-9
a9                 1.28167522926204
a10                1.41117341450273E-10