matlab 分式拟合,如何用matlab拟合微分方程

//这里是代码窗口，请将Forcal代码写在下面

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

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

//函数定义，用于计算微分方程组中各方程右端函数值，连分式法对微分方程组积分一步函数pbs1将调用该函数f。

//t为自变量，x为函数值，dx为右端函数值(即微分方程的值)。

//该函数被调用时将用到参数A,k，这2个参数正好是拟合参数，用模块变量传递这2个参数。

f(t,x,dx::A,k)={

dx=A*x/(k+x)

};

//计算相邻两组实验数据的理论值与实验值的误差：第一组数据作为理论值的起点，用以计算出终点处的理论值，与第二组数据代表的实验值做比较。

//hf为函数f的句柄，预先获得了该句柄；a为工作数值；step为积分步数；eps为积分精度；t1,t2为积分的起点和终点；

//x1是对应t1的实验数据，即积分起始点处的函数值；x_1是对应t2的实验数据，与积分终点处的函数值做比较。

//h,i为自动变量。

t_i_j(hf,a,step,eps,t1,t2,x1,x_1:h,i)=

{

a.setra(0,x1),                  //工作数组赋初值

h=(t2-t1)/step,