wangyuenter 发表于 2008-3-27 21:35

解微分方程组能直接用矩阵么

比如dx/dt=Ax+g(t)x(0)=0 中 A为300阶矩阵
用最简单的欧拉法 能直接编程为
x=zeros(300,1)
h=1/0.05 %从0到1 步长0.05
ht=0.05 % 步长
for t=0:0.05:1
    x=x+h*(A*x+subs(g,t))
end

还是得一个一个分量化?

sigma665 发表于 2008-3-27 22:20

回复 楼主 的帖子

这样编算出来的结果对不对

wangyuenter 发表于 2008-3-28 20:22

请教下 是这样出来的结果是错误的么

wr_9864 发表于 2008-3-28 20:24

matlab,提供左乘和右乘的乘法，你直接用“*”当然不行！

无水1324 发表于 2008-3-28 22:26

回复 楼主 的帖子

楼主为什么不用现有的ODE，算法呢，是不是有特殊的要求？

wangyuenter 发表于 2008-3-30 20:35

求教ode怎么导入矩阵呢?
页: [1]
查看完整版本: 解微分方程组能直接用矩阵么