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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 帮忙用MATLABode45求解2阶微分方程-----新问题！
北京石油化工学院2026年研究生招生接收调剂公告
51/1返回列表
查看: 1682  |  回复: 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的回帖

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的回帖

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的回帖
相关版块跳转 我要订阅楼主 老花 的主题更新
51/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 265求调剂 +8 小木虫085600 2026-04-06 8/400 2026-04-06 17:20 by swb0218
[考研] 338求调剂 +4 wxygxsaaaaa 2026-04-06 4/200 2026-04-06 17:14 by 尚水阁主
[考研] 0857大类环境工程B区求调剂 +3 龚禹铭 2026-04-05 3/150 2026-04-06 10:22 by 蓝云思雨
[考研] 0817化学工程与技术求调剂，一志愿中海洋319 +14 lv945 2026-04-04 14/700 2026-04-06 10:20 by 蓝云思雨
[考研] 319求调剂 +3 handrui 2026-04-05 3/150 2026-04-06 09:33 by jp9609
[考研] 一志愿武汉理工大学-085601材料工程（专硕）-总分353求调剂 +3 2626262626li 2026-04-02 3/150 2026-04-06 09:08 by 无际的草原
[考研] 308求调剂 +3 终不似从前 2026-04-05 3/150 2026-04-05 22:23 by hemengdong
[考研] 求调剂 +11 熊二想上岸 2026-04-04 11/550 2026-04-05 22:21 by 醉翁wl
[考研] 377求调剂 +6 by.ovo 2026-04-05 6/300 2026-04-05 22:18 by dongzh2009
[考研] 调剂 +5 好好读书。 2026-04-01 5/250 2026-04-05 17:54 by liucky
[考研] 0854求调剂 +4 assdll 2026-04-04 4/200 2026-04-05 09:44 by zhq0425
[考研] 353求调剂 +10 MayUxw1 2026-04-03 10/500 2026-04-05 09:23 by 无际的草原
[考研] 男生，一志愿沪9生物学071000，初试308求调剂 +3 刘墨墨 2026-04-04 3/150 2026-04-05 08:26 by barlinike
[考研] 266分，一志愿电气工程，本科材料，求材料专业调剂 +11 哇呼哼呼哼 2026-04-01 12/600 2026-04-04 23:17 by 永字号
[考研] 材料专硕306英一数二 +8 z1z2z3879 2026-03-31 8/400 2026-04-04 22:08 by hemengdong
[考研] 341求调剂 +3 洛多罗 2026-04-02 4/200 2026-04-04 21:36 by 智能智慧
[考研] 一志愿085404，总分291，四级已过，求调剂 +5 阿俊阿俊阿俊 2026-04-04 7/350 2026-04-04 13:23 by 莲菜就是藕吧
[考研] 求调剂，一志愿南京航空航天大学 ，080500材料科学与工程学硕 +10 @taotao 2026-04-03 10/500 2026-04-04 09:01 by T可可西里T
[考研] 279求调剂 +6 qazplm0852 2026-04-02 6/300 2026-04-03 10:03 by 蓝云思雨
[考研] 296求调剂 +4 汪！？！ 2026-03-31 7/350 2026-04-01 22:04 by 客尔美德
信息提示
关闭
请填处理意见
关闭