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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 帮忙用MATLABode45求解2阶微分方程-----新问题！
51/1返回列表
查看: 1658  |  回复: 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的回帖

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的回帖
查看全部 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的回帖

老花

木虫 (著名写手)

老花
引用回帖:
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的回帖

老花

木虫 (著名写手)

老花
引用回帖:
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 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0854控制工程 359求调剂 可跨专业 +3 626776879 2026-03-14 9/450 2026-03-16 17:42 by 626776879
[考研] 梁成伟老师课题组欢迎你的加入 +8 一鸭鸭哟 2026-03-14 9/450 2026-03-16 17:35 by 沐霖12138
[考研] 0703化学调剂 ，六级已过，有科研经历 +7 曦熙兮 2026-03-15 7/350 2026-03-16 16:34 by houyaoxu
[考研] 085600材料与化工 求调剂 +13 enenenhui 2026-03-13 14/700 2026-03-16 15:19 by 了了了了。。
[考研] 274求调剂 +4 时间点 2026-03-13 4/200 2026-03-15 15:29 by Rambo13
[考研] 211本，11408一志愿中科院277分，曾在中科院自动化所实习 +3 Losir 2026-03-12 3/150 2026-03-14 12:11 by 热情沙漠
[考研] 331求调剂（0703有机化学 +5 ZY-05 2026-03-13 6/300 2026-03-14 10:51 by Jy?
[考研] 求调剂 +6 yfihxh 2026-03-09 6/300 2026-03-14 01:18 by JourneyLucky
[考研] 307求调剂 +7 超级伊昂大王 2026-03-10 7/350 2026-03-14 00:49 by JourneyLucky
[考研] 271求调剂 +10 生如夏花… 2026-03-11 10/500 2026-03-14 00:35 by 卖报员小雨
[考研] 一志愿湖师大化学289求调剂 +6 XMCMM3.14159 2026-03-10 6/300 2026-03-14 00:28 by JourneyLucky
[考研] 0805，333求调剂 +3 112253525 2026-03-10 3/150 2026-03-13 23:42 by JourneyLucky
[考研] 材料与化工（0856）304求B区调剂 +6 邱gl 2026-03-12 7/350 2026-03-13 23:24 by 邱gl
[考研] 一志愿中科院，化学方向，295求调剂 +4 一氧二氮 2026-03-11 4/200 2026-03-13 22:35 by JourneyLucky
[考研] 工科，求调剂 +3 我887 2026-03-11 3/150 2026-03-13 21:39 by JourneyLucky
[考研] 301求调剂 +6 Liyouyumairs 2026-03-11 6/300 2026-03-13 20:11 by JourneyLucky
[考研] 求调剂 +7 18880831720 2026-03-11 7/350 2026-03-13 16:10 by JourneyLucky
[考研] 277求调剂 +4 anchor17 2026-03-12 4/200 2026-03-13 11:15 by 白夜悠长
[考研] 一志愿河海大学085900土木水利专硕279求调剂不挑专业 +4 SunWwWwWw 2026-03-10 8/400 2026-03-13 02:23 by SunWwWwWw
[考博] 2026年博士申请 +3 QwQwQW10 2026-03-11 3/150 2026-03-12 17:58 by gxch43
信息提示
关闭
请填处理意见
关闭