微分方程式：

x'=dx/dt=a*0.0321*(b-x)-d*x-dy/dt,

y'= dy/dt=0.25*p1*exp(-p1*t)*x ,

四个待求参数：a、b、d、p1

t、x、y数据见下面：

//0 0 0  //这是初值

0,0,0,

0.1,0.486966799,0.048018378,

0.167,1.6657,0.05823,

0.2,0.860306078,0.060834243,

0.3,1.156255213,0.064254733,

0.4,1.390856542,0.065167644,

0.5,1.67518,0.06638,

0.6,1.724247244,0.065476325,

0.7,1.841108525,0.065493681,

0.8,1.93374543,0.065498314,

0.9,2.007179471,0.06549955,

1,1.92438,0.05641,

1.1,2.111536196,0.065499968,

1.2,2.148115682,0.065499991,

1.3,2.177112544,0.065499998,

1.4,2.200098596,0.065499999,

1.5,2.218319829,0.0655,

1.6,2.232763952,0.0655,

1.7,2.244213928,0.0655,

1.8,2.253290418,0.0655,

1.9,2.260485427,0.0655,

2,2.16156,0.07359,

2.1,2.270710216,0.0655,

2.2,2.274294246,0.0655,

2.3,2.277135335,0.0655,

2.4,2.27938749,0.0655,

2.5,2.281172792,0.0655,

2.6,2.282588016,0.0655,

2.7,2.283709875,0.0655,

2.8,2.284599183,0.0655,

2.9,2.285304144,0.0655,

3,2.44976,0.15449,

3.1,2.286305961,0.0655,

3.2,2.286657121,0.0655,

3.3,2.286935489,0.0655,

3.4,2.287156153,0.0655,

3.5,2.287331076,0.0655,

3.6,2.287469738,0.0655,

3.7,2.287579657,0.0655,

3.8,2.287666791,0.0655,

3.9,2.287735862,0.0655,

4,2.287790616,0.0655,

4.1,2.287834019,0.0655,

4.2,2.287868426,0.0655,

4.3,2.2878957,0.0655,

4.4,2.287917321,0.0655,

4.5,2.287934459,0.0655,

4.6,2.287948045,0.0655,

4.7,2.287958815,0.0655,

4.8,2.287967352,0.0655

用MATLAB拟合不了，在本版搜索到一个帖子，用forcal写，但是不太会

我的Forcal代码如下：

i: OutVector(p:k,i)= k=FCDLen(p),printff{"\r\n"},i=0,(i

{"\r\n"};    //输出一维数组

!using["XSLSF"]; //使用命名空间XSLSF

f(t,x,y,dx,dy::p1,p2,p3,p4)={ //函数定义，连分式法对微分方程组积分一步函数pbs1中要用到

dx=p1*0.0321*(p2-x)-p3*x-dy,

dy=0.25*p4*exp(-p4*t)*x

};

//用连分式法对微分方程组进行积分，获得理论值。

//t1,t2为积分的起点和终点。

//h,s为自动变量。

//模块变量：hf为函数f的句柄，要预先获得该句柄；Array为工作数组；step为积分步长；eps为积分精度。

t_i_2(hf,Array,step,eps,t1,t2,x_1,x_2:x1,x2,h,i)=    //用于计算目标函数

{

h=(t2-t1)/step,

{   pbs1[hf,t1,Array,h,eps],  //连分式法对微分方程组积分一步函数pbs1，hf为函数f的句柄

t1=t1+h

}.until[abs(t1-t2)

Array.getra(0,&x1,&x2),

(x1-x_1)^2+(x2-x_2)^2

};

//目标函数定义，自变量_p1,_p2,_p3,_p4为需要优化的参数，需要将这些参数传递给对应的模块变量p1,p2,p3,p4。

//模块变量：Array为工作数组；DataArray存放实验数据t,x,y；max为数组DataArray的列长度；a,b,d,p1为优化变量。

//用于优化

Opt0(_p1,_p2,_p3,_p4:t1,s,i:hf,Array,step,eps,max,p1,p2,p3,p4,DataArray)={

p1=_p1,p2=_p2,p3=_p3,p4=_p4,    //传递优化变量，函数f中要用到p1,p2,p3,p4

t1=0,Array.setra(0,0,0,0),//设置积分初值，通过模块变量传递，Array是一个数组

s=0,i=0,

(i

s=s+t_i_2[hf,Array,step,eps:&t1,DataArray.getra(i,0):DataArray.getra(i,1),DataArray.getra(i,2)],

i++

},

s

};

//验证函数定义

验证(_p1,_p2,_p3,_p4:t1,s1,s2,i:hf,Array,step,eps,max,p1,p2,p3,p4,DataArray)={

p1=_p1,p2=_p2,p3=_p3,p4=_p4,

t1=0,Array.setra(0,0,0),

i=0,printff{"\r\n  Time        实验数据x              计算x             实验数据y              计算y \r\n\r\n"},

i

t_i_2[hf,Array,step,eps:&t1,DataArray.getra(i,0):DataArray.getra(i,1),DataArray.getra(i,2)],

Array.getra(0,&s1,&s2),

printff{"{1,r,6.3}{2,r,22.16}{3,r,22.16}{4,r,22.16}{5,r,22.16}\r\n",DataArray.getra(i,0),DataArray.getra(i,1),s1,DataArray.getra(i,2),s2},

i++

}

};

main(:d1,u,v,x,_eps,k,xx,g,i,s1,s2:hf,Array,step,eps,max,DataArray)=

{

hf=HFor("f"),           //模块变量hf保存函数f的句柄，预先用函数HFor获得该句柄

step=0.01,eps=1e-7,    //积分步长step要合适，积分精度eps越小越精确，用于对微分方程组积分一步函数pbs1

max=50,

DataArray=new[rtoi(real_s),rtoi(max),rtoi(3),rtoi(EndType),

0,0,0,

0.1,0.486966799,0.048018378,

0.167,1.6657,0.05823,

0.2,0.860306078,0.060834243,

0.3,1.156255213,0.064254733,

0.4,1.390856542,0.065167644,

0.5,1.67518,0.06638,

0.6,1.724247244,0.065476325,

0.7,1.841108525,0.065493681,

0.8,1.93374543,0.065498314,

0.9,2.007179471,0.06549955,

1,1.92438,0.05641,

1.1,2.111536196,0.065499968,

1.2,2.148115682,0.065499991,

1.3,2.177112544,0.065499998,

1.4,2.200098596,0.065499999,

1.5,2.218319829,0.0655,

1.6,2.232763952,0.0655,

1.7,2.244213928,0.0655,

1.8,2.253290418,0.0655,

1.9,2.260485427,0.0655,

2,2.16156,0.07359,

2.1,2.270710216,0.0655,

2.2,2.274294246,0.0655,

2.3,2.277135335,0.0655,

2.4,2.27938749,0.0655,

2.5,2.281172792,0.0655,

2.6,2.282588016,0.0655,

2.7,2.283709875,0.0655,

2.8,2.284599183,0.0655,

2.9,2.285304144,0.0655,

3,2.44976,0.15449,

3.1,2.286305961,0.0655,

3.2,2.286657121,0.0655,

3.3,2.286935489,0.0655,

3.4,2.287156153,0.0655,

3.5,2.287331076,0.0655,

3.6,2.287469738,0.0655,

3.7,2.287579657,0.0655,

3.8,2.287666791,0.0655,

3.9,2.287735862,0.0655,

4,2.287790616,0.0655,

4.1,2.287834019,0.0655,

4.2,2.287868426,0.0655,

4.3,2.2878957,0.0655,

4.4,2.287917321,0.0655,

4.5,2.287934459,0.0655,

4.6,2.287948045,0.0655,

4.7,2.287958815,0.0655,

4.8,2.287967352,0.0655

],

Array=new[rtoi(real_s),rtoi(2*max)], //申请工作数组

x=new[rtoi(real_s),rtoi(3)],                 //申请工作数组

xx=new[rtoi(real_s),rtoi(2),rtoi(3)],        //申请工作数组

g=new[rtoi(real_s),rtoi(3)],                 //申请工作数组

_eps=1e-100, d1=1,u=1.6,v=0.4,k=800,           //变换d1、u、v进一步求解，k为允许的最大迭代次数

i=jsim[HFor("Opt0"),d1,u,v,x,_eps,k,xx,g],       //求n维极值的单形调优法

printff{"\r\n实际迭代次数={1,r}\r\n",i},     //输出实际迭代次数

OutVector[x],                                //输出最优参数值及目标函数终值

x.getra(0,&s1,&s2),

验证[s1,s2],

delete[x],delete[xx],delete[g],delete[Array],delete[DataArray] //销毁申请的对象

}

但是总是提示，验证[s1,s2],调用实数表达式时参数不匹配。

，请问该怎么改呢？新学中，谢谢帮助！