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

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

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

比如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怎么导入矩阵呢?

账号		自动登录	找回密码
密码			我要加入

[编程技巧] 解微分方程组能直接用矩阵么

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

回复楼主的帖子

回复楼主的帖子

[编程技巧] 解微分方程组能直接用矩阵么

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

回复 楼主 的帖子

回复 楼主 的帖子

回复楼主的帖子

回复楼主的帖子