带有系数矩阵的二阶微分方程求解
小弟现在需要求解一个二阶微分方程组y''=A*y'-2*y+B，其中A=[-2,1;998, -999],B=,为了求得数值解，采用ode45函数，编制程序如下：function dy=f(t,Y)
A=[-2 1;998,-999];
B=;
Y=zeros(2,2);
dy1dt=Y(:,1);
dy2dt=A*Y(:,1)-2*Y(:,2)+B;
dy=;
在命令行输入ode45(@f,,)
得到的结果与我用dsolve测试的结果完全不同，我发现主要问题是在Y=zeros(2,2)这行上，可是我很想利用矩阵的形式来求解微分方程组，怎么办呢？恳请大虾指点，呵呵 发帖前，先搜索。这类问题版面多的是，新手请先看置顶的会员守则
[ 本帖最后由 ChaChing 于 2010-6-15 14:52 编辑 ]
页:
[1]