将二阶方程
化为一阶方程组。取h=0.1,用四阶龙格-库塔法求y(0.2)的近似值,保留5位有效数字。
用罚函数法求解问题
(1)写出ck=0,1,10时相应的增广目标函数,并画出它们对应的图形;
(2)取ck=k-1(k=1,2,...)求出近似最优解的迭代点列;
(3)利用(2)求问题的最优解。
式中: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]上的数值解。
采用位移法解图中所示连续梁时,基本未知量Δ1和Δ2,分别为B和C结点的角位移,基本方程为
试用式(13-38)和式(13-39)求kij和Fip并与用第7章的方法得出的结果加以比较。设各杆EI=常数。