3楼: Originally posted by dingd at 2014-06-10 11:53:40
1stOpt可以很容易实现的。

3楼: Originally posted by dingd at 2014-06-10 11:53:40
1stOpt可以很容易实现的。

Constant a=5;
Parameters Cus,Ku,Kr,Kf,So;
Variables x,y;
Function y=Cus*Ku*(1+Kr)*sum(i=1:a+1)(x^i)/(Ku*(1+Kr)*sum(i=1:a+1)(x^i)+sum(i=1:a+1)(x^(i-a)))+So*Kf*x/(1+Kf*x);
Data;
0.56509 0
2.30222 0.444696164
4.03935 0.841178994
.....

6楼: Originally posted by dingd at 2014-07-30 14:44:02
求和的上限值a必须提前给定，不能当作未知参数来求，下面a取5：

Constant a=5;
Parameters Cus,Ku,Kr,Kf,So;
Variables x,y;
Function y=Cus*Ku*(1+Kr)*sum(i=1:a+1)(x^i)/(Ku*(1+Kr)*sum(i=1:a+1)(x^i)+sum( ...

IntParameter a=[1,50];
Parameters Cus,Ku,Kr,Kf,So;
Variables x,y;
StartProgram [Pascal];
Procedure MainModel;
var i, j, p: integer;
    temd, temd1, temd2, temd3: double;
Begin
    p := Round(a);
    for j := 0 to DataLength - 1 do begin
        temd1 := Cus*Ku*(1+Kr);
        temd2 := 0;
        for i := 1 to p + 1 do
            temd2 := temd2 + power(x[j],i);
        temd3 := 0;
        for i := 1 to p + 1 do
            temd3 := temd3 + power(x[j],i-p);
        y[j] := temd1*temd2/(Ku*(1+Kr)*temd2+temd3)+So*Kf*x[j]/(1+Kf*x[j]);
    end;
End;
EndProgram;
Data;
0.56509 0
2.30222 0.444696164
4.03935 0.841178994
5.96484 1.135857962
8.07869 1.532340792
10.98786 1.923468197
15.21557 2.368164361
20.21766 2.909308963
32.16827 4.195191377
39.68188 5.127456128
49.89535 6.365125209
60.10883 8.138568953
70.32231 10.50674057
80.13813 13.46427738
90.16325 16.22356915
94.60025 17.01117213

8楼: Originally posted by dingd at 2014-07-31 02:21:13
1stOpt中可以用编程模式解决：

IntParameter a=;
Parameters Cus,Ku,Kr,Kf,So;
Variables x,y;
StartProgram ;
Procedure MainModel;
var i, j, p: integer;
    temd, temd1, temd2, temd3: double;
Beg ...

8楼: Originally posted by dingd at 2014-07-31 02:21:13
1stOpt中可以用编程模式解决：

IntParameter a=;
Parameters Cus,Ku,Kr,Kf,So;
Variables x,y;
StartProgram ;
Procedure MainModel;
var i, j, p: integer;
    temd, temd1, temd2, temd3: double;
Beg ...