本文共 447 字,大约阅读时间需要 1 分钟。
//使用XSLSF命名空间
数组xArray存放x的值,yArray存放y的值;ti为当前有效值的个数;tmax为ti对应的时间;tmin为起始时间。
xt(t:k:xArray,ti,tmax,tmin)函数定义:
k=(t-tmin)/(tmax-tmin)*ti-1
如果k小于0,则取k=0;如果k大于等于ti,则取k=ti-1;
返回xArray.getrai[k]的值。
类似地,yt(t:k:yArray,ti,tmax,tmin)函数定义:
k=(t-tmin)/(tmax-tmin)*ti-1
同上,返回yArray.getrai[k]的值。
微分方程组右端函数值计算函数f定义:
函数f(t,x,y,dx,dy):
dx = -(x³) - y
dy = x - y³
用连分式法对微分方程组进行积分,获得理论值。
函数f将被pbs1函数调用,用于微分方程组的积分步骤。
积分起点为t1,终点为t2;h和s为自动变量。
上述内容为对微分方程组进行连分式积分方法的实现基础。
转载地址:http://vfyfk.baihongyu.com/