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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 帮忙用MATLABode45求解2阶微分方程-----新问题！
51/1返回列表
查看: 1804  |  回复: 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返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 u1xkze8ybe 2026-07-30 7/350 2026-07-31 12:10 by syhhffkrym
[论文投稿] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 u1xkze8ybe 2026-07-30 3/150 2026-07-31 12:08 by syhhffkrym
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 u1xkze8ybe 2026-07-30 6/300 2026-07-31 12:05 by syhhffkrym
[教师之家] 基础研究怎么拉横向，学校到款任务越来越多，难以完成 拉横向，都有哪些途径啊 +6 锦衣卫寒战 2026-07-28 6/300 2026-07-31 11:27 by syl200707
[高分子] HXDI做水性聚氨酯乳液，是不是特别容易出渣 15+3 yuyusuv 2026-07-29 3/150 2026-07-31 09:21 by huizingga
[基金申请] 小木虫看见有人已经知道结果了 +15 1234567wang 2026-07-25 19/950 2026-07-31 08:30 by ZJTJZ
[基金申请] 2026年国自然面上资助率 +8 布布和一二 2026-07-30 8/400 2026-07-31 01:51 by alongwaytogo
[基金申请] 系统今天又提示维护了，估计离放榜不远了 +11 winnerche 2026-07-29 15/750 2026-07-30 23:08 by jnhyjjm
[基金申请] 微信指数没变化，科研之友没阅读 +14 wangze12014 2026-07-28 17/850 2026-07-30 20:10 by kissu88
[考博] 申请2027年材料/化学类博士 10+3 考研老狗? 2026-07-25 5/250 2026-07-30 19:33 by 白色木鱼
[基金申请] 你们的时间戳变了吗 +3 archvillain 2026-07-30 4/200 2026-07-30 18:53 by levinzhwen
[基金申请] 今年的WR进展到哪一步了？ +5 wsgjhwz 2026-07-25 9/450 2026-07-30 12:59 by 可淡不可忘
[有机交流] 产物和副产物价值比较 30+4 小汤02 2026-07-27 5/250 2026-07-30 12:44 by czyzsu
[基金申请] 时间戳他又来了 +14 晓晓爱翠翠 2026-07-26 17/850 2026-07-29 16:30 by 呵呵二哥
[有机交流] 这个自发加氧反应的机理是什么？ 50+4 YUAN2273 2026-07-25 7/350 2026-07-29 12:27 by xiaobenpu
[基金申请] 没消息就是被刷了呗 +13 lambert2014 2026-07-24 14/700 2026-07-29 10:38 by msh145800
[基金申请] 这种情况还有戏吗 +5 drbart 2026-07-27 11/550 2026-07-29 05:50 by drbart
[基金申请] 同事接到电话了，我却没有 +4 1234567wang 2026-07-27 4/200 2026-07-28 10:11 by GOODLUCKER
[基金申请] 准备明年的基金了 +5 Tide man 2026-07-27 5/250 2026-07-28 09:25 by newfuzzy1
[分析] 水蒸气蒸馏法乳化现象 +3 月光娇娇 2026-07-26 3/150 2026-07-26 23:25 by caiyun
信息提示
关闭
请填处理意见
关闭