带有系数矩阵的二阶微分方程求解

chanelq

小弟现在需要求解一个二阶微分方程组y''=A*y'-2*y+B，其中A=[-2,1;998, -999],B=[2*sin(t);999*(cos(t)-sin(t))],为了求得数值解，采用ode45函数，编制程序如下：
function dy=f(t,Y)
A=[-2 1;998,-999];
B=[2*sin(t);999*(cos(t)-sin(t))];
Y=zeros(2,2);
dy1dt=Y(:,1);
dy2dt=A*Y(:,1)-2*Y(:,2)+B;
dy=[dy1dt;dy2dt];
在命令行输入ode45(@f,[0 10],[0,0;0,0])
得到的结果与我用dsolve测试的结果完全不同，我发现主要问题是在Y=zeros(2,2)这行上，可是我很想利用矩阵的形式来求解微分方程组，怎么办呢？恳请大虾指点，呵呵

eight

发帖前，先搜索。这类问题版面多的是，新手请先看置顶的会员守则

[ 本帖最后由 ChaChing 于 2010-6-15 14:52 编辑 ]