A.interp1(x,y,4.4,’spline’)
B.interp1(x,y,’spline’)
C.interp1(x,y,4.4,’cubic’)
D.interp1(x,y,4.4)
A.y=4 1 2 3
B.y=3 2 2 3
C.y=3 2 3 2
D.y=4 2 1 1
其中x0是给定的x(t)的初始值,xp0是任意给定的x(1)的初始值,fixed_:x0和fixed_xp0是与xp0同维数的列向量,其分量为1表示需要保留的初值,为0表示需要求解的初始值。若fixed_x0和fixed_xp0等于空矩阵[],表示允许所有的初值分量可以发生变化。分别用显式和隐式解法求下列微分方程的数值解
式中:T≥0为时滞常数。在Matlab中提供了命令dde23来直接求解时滞微分方程。其调用格式为801=dde23(ddefun,lags,history,tspan,options),
其中,ddfun为描述时滞微分方程的函数;lags为时滞常数向量;history为描述t≤to时的状态变量值的函数;tspan为求解的时间区间;options为求解器的参数设置。该函数的返回值sol是结构体数据,其中sol.x成员变量为时间向量l,sol.y成员变量为各个时刻的状态向量构成的矩阵,其每一个行对应着一个状态变量的取值。求解如下时滞微分方程组:
已知,在i≤0时,x(t)=5,x2(t)=0,x(1)=1,试求该方程组在[0,40]上的数值解。
利用现成的数学通用软件(如MATLAB、Mathematica、Maple等), 对于N = 32, 64,128 :
(5) 改变δ0的值,重复(4),分析不同的δ0对逆变换所得到的数据的影响。
用解析法和MATLAB法相结合求解差分方程。设系统差分方程为
输入为r(m)=cos(πn/3)u(n)边界条件为:y(-1)=-2.y(-2)=-3.x(-1)=1.x(-2)=1