24小时热门版块排行榜    

查看: 6164  |  回复: 13
【悬赏金币】回答本帖问题,作者wzxchars将赠送您 10 个金币
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

wzxchars

铁杆木虫 (著名写手)

[求助] 如何用matlab求解矩阵系数的二阶微分方程已有1人参与

以一个简单的有限元动力方程为例,M*X''+K*X=F(t),其中M和K为已知质量和刚度矩阵(16*16),F(t)为已知外部载荷输入(16*1),如何求解X呢?用ode45还是用simulink?希望高手们给出详细步骤,或相关例子

[ 来自科研家族 机器人 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

静坐常思自己过,闲谈莫论他人非
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jun_1129

铜虫 (初入文坛)

【答案】应助回帖

function z=zi(t,x)
global m
global k
z1=x(4);
z2=x(5);
z3=x(6);
z4=(k/m)*x(2)-3*(k/m)*x(1);
z5=(k/m)*x(1)-2*(k/m)*x(2)+(k/m)*x(3);
z6=(k/m)*x(2)-3*(k/m)*x(3);
z=[z1;z2;z3;z4;z5;z6];
end

……………………main.m………………………………
global m
global k
z1=1;
z2=1;
z3=-1;
z4=1;
z5=1;
z6=1;
m=2;
k=1;
t1=0;
t2=10;
[t,f]=ode45(@func,[t1 t2],[z1 z2 z3 z4 z5 z6])
plot(t,f)
8楼2013-10-30 17:59:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 14 个回答

csgt0

荣誉版主 (著名写手)

彩色挂图

【答案】应助回帖

★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
wzxchars: 金币+5, 有帮助 2012-10-24 19:30:25
xiegangmai: 金币+1, 谢谢参与 2012-10-24 19:57:25
比较麻烦,你搜搜帮助里的Several Differential Equations有个用dsolve的方法,但是估计不好解出来。
showmethemoney
2楼2012-10-24 17:07:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wzxchars

铁杆木虫 (著名写手)

引用回帖:
2楼: Originally posted by csgt0 at 2012-10-24 17:07:18
比较麻烦,你搜搜帮助里的Several Differential Equations有个用dsolve的方法,但是估计不好解出来。

ode45是能够解这种带有时间参数的常系数二阶微分方程的,但对于这种矩阵形式的解法我还没有用过,有一种最基本的方法是矩阵拆分后联立微分方程组求解,但是实际上我要解的是100*100的矩阵,光是函数表达就要有好几百个项,唉头大,希望能在这里遇到高手指点。谢谢版主的热心回复
静坐常思自己过,闲谈莫论他人非
3楼2012-10-24 19:30:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

eezyzyzzzee

木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
wzxchars: 金币+15, ★★★很有帮助, 谢谢了 2012-10-24 23:32:27
xiegangmai: 金币+1, 谢谢参与 2012-10-27 19:43:23
你搜搜——结构分析的有限元法与MATLAB程序设计_徐荣桥.pdf
看看最后一章
http://wenku.baidu.com/view/5cf31dc8050876323112124f.html

不懂你的领域,就是偶然看到过上面这个教程,有动力学响应这块的
4楼2012-10-24 21:25:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 关系导师与研究生的关系的一点看法 +12 平凡依旧 2025-01-01 38/1900 2025-01-04 16:07 by 平凡依旧
[基金申请] 报2停1的政策还是有道理的 +4 原因在哪里 2025-01-04 4/200 2025-01-04 13:56 by 六两废铜
[教师之家] 专利已受理,未授权,能做作为课题年度进展成果 +7 原因在哪里 2025-01-03 10/500 2025-01-04 13:41 by 六两废铜
[硕博家园] 江苏大学杨孙圣课题组 +9 姜一雪月 2024-12-31 10/500 2025-01-04 13:05 by qsx04152006
[论文投稿] 昨天online,今天发现单位顺序错了,还能改吗? +7 oooooo?o 2025-01-02 8/400 2025-01-04 12:44 by 一路禅修
[教师之家] 终于还是决定今年申报项目 +9 原因在哪里 2025-01-02 13/650 2025-01-04 11:28 by willbuilder
[基金申请] 国自然基金委主任窦贤康:下一步制定引导政策 +6 nandi2212 2024-12-31 8/400 2025-01-04 11:04 by xli1984
[教师之家] 我们组一个学生一直不明白研究生和本科生的区别 +9 苏东坡二世 2025-01-03 9/450 2025-01-04 10:58 by yuningning
[教师之家] 穷啊,准备考研辅导班当老师了 +4 guilin2018 2025-01-02 5/250 2025-01-04 09:46 by 六两废铜
[基金申请] 网传,基金委面上项目评审专家名单打听费7万,获批后再提10% +13 babu2015 2025-01-02 15/750 2025-01-04 09:28 by MonCrist
[论文投稿] Rare metals投稿求助 10+3 hngo 2025-01-02 11/550 2025-01-04 07:56 by steven_198377
[论文投稿] 期刊投稿咨询 +9 贝壳持儿 2025-01-01 9/450 2025-01-04 07:23 by steven_198377
[论文投稿] 工程科学学报期刊 5+3 巧巧玲 2024-12-30 5/250 2025-01-03 14:47 by 木木A123
[考博] 专家推荐信 +4 ℡额尔古纳 2025-01-02 5/250 2025-01-03 08:09 by ℡额尔古纳
[教师之家] 小教师身体第一,健康第一,保命最重要 +7 guilin2018 2024-12-30 10/500 2025-01-03 07:59 by yiran909
[论文投稿] 圣诞节国外放几天假啊,审稿状态都不动 +10 mollyzhang_2003 2024-12-30 14/700 2025-01-02 20:57 by frank1139
[考博] 考博综合面试,老师喜欢提什么样的问题? +8 爱喝风的龙卷茶 2025-01-02 9/450 2025-01-02 13:45 by 异形杀手
[论文投稿] Journal of Agricultural and Food Chemistry在哪里看文章详细状态 5+3 吃饭了嘎 2024-12-31 3/150 2025-01-02 09:17 by 北京莱茵润色
[教师之家] 现在的教改课题方案,平时真的会那么做吗? +4 akslis2024 2025-01-01 4/200 2025-01-02 09:10 by Quakerbird
[硕博家园] (常年有效)中外知名高校(QS排名前3)联培 - 招募大模型及强化学习博士生和博士后 +3 cowhorce 2024-12-29 3/150 2024-12-29 18:40 by 55bh@2u5
信息提示
请填处理意见