| 查看: 2609 | 回复: 5 | ||||
[求助]
matlab二阶微分方程求解程序
|
||||
|
%已知二阶微分方程组，求解Mx’’+Kx=Qe 其中M，K均为6*6矩阵，x为6*1向量，每一个x代表一个位移， 希望求出每一个x关于t的曲线 %主程序 global k1 m k1=[5.31645e9 0 0 -5.31645e9 0 0; 0 319.4473e9 252.3634e9 0 -319.4473e9 252.3634e9; 0 252.3634e9 265.8228e9 0 -252.3634e9 132.9114e9; -5.31645e9 0 0 5.31645e9 0 0; 0 -319.4473e9 -252.3634e9 0 319.4473e9 319.4473e9; 0 252.3634e9 132.9114e9 0 -252.3634e9 168.2423e9;]; m=[ 1.9333 0 0 0.9667 0 0; 0 2.1543 0.48 0 0.7457 0.2836; 0 0.48 0.1379 0 0.2836 -0.1034; 0.9667 0 0 1.9333 0 0; 0 0.7457 0.2836 0 2.1543 -0.48; 0 -0.2836 -0.1034 0 -0.48 0.1379;]; y0=zeros(1,12); [t,x]=ode45(@order2,[0 5],y0); plot(t, x);%绘制t,q图形 %子程序 function sol=order2(t,x) global k1 m a=0.06*t^3; B=[sin(a) cos(a) 0 0 0 0; -cos(a) sin(a) 0 0 0 0; 0 0 1 0 0 0; 0 0 0 sin(a) cos(a) 0; 0 0 0 -cos(a) sin(a) 0; 0 0 0 0 0 1;]; F=[-2.842*cos(a);-4.263*sin(a)+99.114;-1.648*sin(a)-55.614; -2.842*sin(a);-2.842*cos(a);1.123*sin(a)+15.66;]; M=B'*m*B; K=B'*k1*B; Qe=B'*F; sol=zeros(12,1); sol(7：12)=x(1：6); sol(1:6)=M\(Qe-K*x(7:12)); end 能够运行，但是结果是发散的，急，急，急 |
回复此楼» 收录本帖的淘帖专辑推荐
 机械科研经验专辑 | 科研资料 |
» 猜你喜欢
26材料博士申请自荐
已经有3人回复
-
高温合金金属球（直径10mm）加工
已经有2人回复
-
工程热物理与能源利用论文润色/翻译怎么收费?
已经有167人回复
-
定向工程博士如何申请
已经有2人回复
-
紧急招收2026年秋季入学博士生1名（湘潭大学 固体废弃物低碳利用湖南省工程研究中心）
已经有20人回复
-
测面比电阻，600-850℃
已经有0人回复
-
【急招】“双一流”高校-新能源材料课题组招收2026年秋季入学博士生1名(湘潭大学)
已经有22人回复
-
爱尔兰都柏林大学工学院全额奖学金招聘博士研究生
已经有1人回复
-
请问下国内的中文EI学报有哪些投稿周期比较快的
已经有7人回复
-
祈祷面上基金能中
已经有36人回复
-
含铝高熵合金 真空感应熔炼
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 高手进来用matlab求解常微分方程，十分感谢！
已经有17人回复
- 如何求解描述振动的二阶微分方程
已经有21人回复
- matlab解微分方程
已经有10人回复
- Matlab求解偏微分方程组
已经有22人回复
- matlab数值求解边界条件微分方程组
已经有7人回复
- 用MATLAB ode45求解2阶微分方程
已经有4人回复
- 【求助】matlab 求解微分方程中的未知参数
已经有20人回复
- 【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题
已经有21人回复
- 【求助】MATLAB 有限差分法（FDM）求解偏微分方程
已经有22人回复
- 【求助】积分微分方程matlab求解
已经有6人回复
- 【求助】用MATLAB求解偏微分方程的数值解 悬赏20
已经有8人回复
huaijuliu
铁杆木虫 (正式写手)
- 应助: 7 (幼儿园)
- 金币: 5653.8
- 红花: 2
- 帖子: 598
- 在线: 296.3小时
- 虫号: 583924
- 注册: 2008-08-01
- 性别: GG
- 专业: 传动机械学
2楼2012-05-20 06:23:52
ssttrroonngg
金虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 2570.6
- 帖子: 334
- 在线: 89.3小时
- 虫号: 1336640
- 注册: 2011-07-03
- 性别: GG
- 专业: 机构学与机器人
3楼2012-05-22 11:14:05
huaijuliu
铁杆木虫 (正式写手)
- 应助: 7 (幼儿园)
- 金币: 5653.8
- 红花: 2
- 帖子: 598
- 在线: 296.3小时
- 虫号: 583924
- 注册: 2008-08-01
- 性别: GG
- 专业: 传动机械学
4楼2012-05-22 21:33:47
5楼2012-05-22 21:58:28
6楼2013-03-24 09:47:55
100