张家谔 发表于 2009-8-14 15:40

求解微分方程时遇到这种情况怎么办


M,K为常数对角矩阵,T为常数
A,B为时间函数矩阵
用matlab能否求出其数值解?
{X(t-T)}应该怎么处理

无水1324 发表于 2009-8-14 17:20

help dde23



DDE23Solve delay differential equations (DDEs) with constant delays.
    SOL = DDE23(DDEFUN,LAGS,HISTORY,TSPAN) integrates a system of DDEs
    y'(t) = f(t,y(t),y(t - tau_1),...,y(t - tau_k)). The constant, positive
    delays tau_1,...,tau_k are input as the vector LAGS. DDEFUN is a function
    handle. DDEFUN(T,Y,Z) must return a column vector corresponding to
    f(t,y(t),y(t - tau_1),...,y(t - tau_k)). In the call to DDEFUN, a scalar T
    is the current t, a column vector Y approximates y(t), and a column Z(:,j)
    approximates y(t - tau_j) for delay tau_j = LAGS(J).The DDEs are
    integrated from T0=TSPAN(1) to TF=TSPAN(end) where T0 < TF. The solution
    at t <= T0 is specified by HISTORY in one of three ways: HISTORY can be
    a function handle, where for a scalar T, HISTORY(T) returns a column
    vector y(t). If y(t) is constant, HISTORY can be this column vector.
    If this call to DDE23 continues a previous integration to T0, HISTORY
    can be the solution SOL from that call.

张家谔 发表于 2009-8-14 17:49

回复 沙发 无水1324 的帖子

非常感谢无水1324
正在学习中
页: [1]
查看完整版本: 求解微分方程时遇到这种情况怎么办