[示例]如何求解微分、积分方程的几个例子

2016-10-13 09:42| 发布者: aspen| 查看: 3685| 评论: 0|原作者: xjzuo|来自: 声振论坛

摘要: 看到有不少人问过二阶动力微分方程的求解问题,现举一个简单的例子, 其余的情形希望读者能举一反三, 自己多思考.

一、关于二阶动力微分方程的求解问题

看到有不少人问过二阶动力微分方程的求解问题,
现举一个简单的例子, 其余的情形希望读者能举一反三, 自己多思考.
%%%=========================================%%%
clear all
n=3;
F=[25;24;20];
m1=31.2;
m2=31.2;
m3=31.2;
k1=67.51;
k2=67.51;
k3=67.51;
c1=0.01;
c2=0.01;
c3=0.01;
M=[m1,0,0;0,m2,0;0,0,m3];
B=[c1+c2,-c2,0;-c2,c2+c3,-c3;0,-c3,c3];
K=[k1+k2,-k2,0;-k2,k2+k3,-k3;0,-k3,k3];
DL=inline('[x(n+1:end,1); inv(M)*(F-B*x(1:n,1)-K*x(1:n,1))]',...
      't','x','flag','n','M','K','F','B');
options = odeset('RelTol',1e-4,'AbsTol',[1e-4 1e-4 1e-5]);
[t,x]=ode45(DL,[0,3],rand(n,1),options,n,M,K,F,B);
plot(t,x(:,1:n))
%%%==========================================%%%

讨论、回复专贴： “二阶动力微分方程的求解问题”回复专贴

[ 本帖最后由 eight 于 2007-6-1 00:13 编辑 ]

12 3 / 3 页下一页

收藏分享邀请

上一篇：Matlab图像处理小结下一篇：Matlab中龙格-库塔(Runge-Kutta)方法原理及实现

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

[示例]如何求解微分、积分方程的几个例子

最新评论

相关分类