应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 帮忙用MATLABode45求解2阶微分方程-----新问题！
51/1返回列表
查看: 1656  |  回复: 4
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

老花

木虫 (著名写手)

老花

[求助] 帮忙用MATLABode45求解2阶微分方程-----新问题！

帮忙用ode45求解2阶微分方程。方程如下！


求解q的特征值和特征性向量！


求解 q值！

ps：m，k 请用矩阵表示，m为对角阵！因为实际上m，k都是56行，56列的矩阵！
赋值的时候，能否用矩阵表示？
只写出过程来就可以！
那个帖子已经无法追加金币，故再次开贴！见谅！


[ Last edited by 老花 on 2011-5-23 at 22:13 ]
回复此楼

» 收录本帖的淘帖专辑推荐

〖经典科研帖子〗

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
平凡中不平凡。我们都希望长寿，但是只有有意义的人生和生命中那些闪光的瞬间才是最重要的。用我们的精神而不是钟表来度量时间的流逝吧。....
1楼 2011-05-23 21:37:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

老花

木虫 (著名写手)

老花
引用回帖:
Originally posted by xiegangmai at 2011-05-24 00:53:16:
楼主要做什么东东的自由振动分析？
建议看看《Vibration Simulation Using  MATLAB and ANSYS》，里面有MATLAB的求解过程和示例代码。


质量矩阵很好办:
m = [ m1, ……, m56 ]; %先将质量表示成一行或一 ...

谢谢！  你回应助贴 给你点币币了！
赞一下 回复此楼
平凡中不平凡。我们都希望长寿，但是只有有意义的人生和生命中那些闪光的瞬间才是最重要的。用我们的精神而不是钟表来度量时间的流逝吧。....
3楼2011-05-25 11:44:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 5 个回答

xiegangmai

版主 (职业作家)

我没头衔

优秀版主优秀版主优秀版主
★ ★
sunyang1988(金币+2): 谢谢交流 2011-05-24 15:24:00
引用回帖:
Originally posted by 老花 at 2011-05-23 21:37:11:
帮忙用ode45求解2阶微分方程。方程如下！


求解q的特征值和特征性向量！


求解 q值！

ps：m，k 请用矩阵表示，m为对角阵！因为实际上m，k都是56行，56列的矩阵！
赋值的时候，能 ...

楼主要做什么东东的自由振动分析？
建议看看《Vibration Simulation Using  MATLAB and ANSYS》，里面有MATLAB的求解过程和示例代码。


质量矩阵很好办:
m = [ m1, ……, m56 ]; %先将质量表示成一行或一列
m = diag( m ); % 构造对角矩阵

刚度矩阵可以先把对角元素按上面的方式构造，非对角的元素单独赋值。
赞一下(1人) 回复此楼
明德厚学、求是创新
2楼2011-05-24 00:53:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

catqq

铜虫 (小有名气)

【答案】应助回帖

★ ★
老花(金币+150): 谢谢应助，那50给上面那个斑竹吧，他私下也给我解决了点问题！ 2011-05-25 17:40:50
sunyang1988(金币+2): 谢谢交流 2011-05-27 14:43:47
只要改这个文件
CODE:
function dy = func(t,y)
% 假设y是行向量
m=[]; %填m
k=[];  %填k

y1=y(size(m,2)+1:size(y,2));
y2=-inv(m)*k*y(1:size(m,1));
dy=[y1,y2];
赞一下(1人) 回复此楼
4楼2011-05-25 17:17:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

老花

木虫 (著名写手)

老花
引用回帖:
Originally posted by catqq at 2011-05-25 17:17:50:
只要改这个文件
CODE:
function dy = func(t,y)
% 假设y是行向量
m=[]; %填m
k=[];  %填k

y1=y(size(m,2)+1:size(y,2));
y2=-inv(m)*k*y(1:size(m,1));
dy=[y1,y2];


画图的时候怎么画啊？
请在http://muchong.com/bbs/viewthread.php?tid=3246430 回答！
现在是画不出来图了！ 解的图！谢谢！！
赞一下 回复此楼
平凡中不平凡。我们都希望长寿，但是只有有意义的人生和生命中那些闪光的瞬间才是最重要的。用我们的精神而不是钟表来度量时间的流逝吧。....
5楼2011-05-30 21:55:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 5 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿985，本科211，0817化学工程与技术319求调剂 +5 Liwangman 2026-03-15 5/250 2026-03-16 17:10 by 我的船我的海
[考研] 一志愿211 0703方向310分求调剂 +3 努力奋斗112 2026-03-15 3/150 2026-03-16 16:44 by houyaoxu
[考研] 0703化学调剂 +6 妮妮ninicgb 2026-03-15 9/450 2026-03-16 16:40 by houyaoxu
[考研] 312求调剂 +3 陌宸希 2026-03-16 4/200 2026-03-16 15:06 by peike
[考研] 材料与化工专硕调剂 +3 heming3743 2026-03-16 3/150 2026-03-16 15:05 by peike
[考研] 277材料科学与工程080500求调剂 +3 自由煎饼果子 2026-03-16 3/150 2026-03-16 14:10 by 运气yunqi
[考研] 0856求调剂 +3 刘梦微 2026-03-15 3/150 2026-03-16 10:00 by houyaoxu
[考研] 材料工程327求调剂 +3 xiaohe12w 2026-03-11 3/150 2026-03-14 20:20 by ms629
[考研] 考研材料与化工，求调剂 +8 戏精丹丹丹 2026-03-09 8/400 2026-03-14 01:14 by JourneyLucky
[考研] 271求调剂 +10 生如夏花… 2026-03-11 10/500 2026-03-14 00:35 by 卖报员小雨
[考研] 一志愿华中农业大学071010，总分三百二,求调剂 +3 困困困困坤坤 2026-03-10 3/150 2026-03-14 00:35 by JourneyLucky
[考研] 一志愿湖师大化学289求调剂 +6 XMCMM3.14159 2026-03-10 6/300 2026-03-14 00:28 by JourneyLucky
[考研] 材料与化工（0856）304求B区调剂 +6 邱gl 2026-03-12 7/350 2026-03-13 23:24 by 邱gl
[考研] 材料专硕288分求调剂 一志愿211 +4 在家想你 2026-03-11 4/200 2026-03-13 22:49 by JourneyLucky
[考研] 0856材料与化工301求调剂 +5 奕束光 2026-03-13 5/250 2026-03-13 22:00 by 星空星月
[考研] 332求调剂 +3 Zz版 2026-03-13 3/150 2026-03-13 20:36 by 18595523086
[考研] 工科278分求调剂 +5 周慢热啊 2026-03-12 7/350 2026-03-13 15:49 by JourneyLucky
[论文投稿] 投稿问题 5+4 星光灿烂xt 2026-03-12 6/300 2026-03-13 14:17 by god_tian
[考研] 08食品或轻工求调剂，本科发表3篇sci一区top论文，一志愿南师大食品科学与工程 +3 我是一个兵， 2026-03-10 3/150 2026-03-13 10:21 by Yuyi.
[考研] 420求调剂 +4 莫向外求11 2026-03-10 6/300 2026-03-12 14:41 by ruiyingmiao
信息提示
关闭
请填处理意见
关闭