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

24小时热门版块排行榜    

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

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的回帖
查看全部 5 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 304求调剂 +4 ahbd 2026-03-14 4/200 2026-03-16 16:48 by 我的船我的海
[考研] 材料专硕306英一数二 +4 z1z2z3879 2026-03-16 4/200 2026-03-16 13:53 by laoshidan
[基金申请] 国自科面上基金字体 +5 iwuli 2026-03-12 6/300 2026-03-16 13:13 by Kamiu_MK
[考研] 0703 物理化学调剂 +3 我可以上岸的对� 2026-03-13 5/250 2026-03-16 10:50 by 我可以上岸的对�
[考研] 290求调剂 +5 孔志浩 2026-03-12 10/500 2026-03-16 09:01 by 余晖&
[考研] 070305求调剂 +3 mlpqaz03 2026-03-14 4/200 2026-03-15 11:04 by peike
[考研] 中科大材料与化工319求调剂 +3 孟鑫材料 2026-03-14 3/150 2026-03-14 20:10 by ms629
[考研] 复试调剂 +3 呼呼？~+123456 2026-03-14 3/150 2026-03-14 16:53 by WTUChen
[考研] 308 085701 四六级已过求调剂 +7 温乔乔乔乔 2026-03-12 14/700 2026-03-14 10:49 by JourneyLucky
[考研] 295复试调剂 +5 简木ChuFront 2026-03-09 5/250 2026-03-14 01:29 by JourneyLucky
[考研] 考研材料与化工，求调剂 +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 卖报员小雨
[考研] 311求调剂 +5 牛乳糖的卡卡 2026-03-10 5/250 2026-03-14 00:05 by JourneyLucky
[考研] 341求调剂 +3 番茄头--- 2026-03-10 3/150 2026-03-13 23:07 by JourneyLucky
[考研] 材料工程调剂 +4 咪咪空空 2026-03-11 4/200 2026-03-13 19:57 by JourneyLucky
[考研] 289求调剂 +3 李政莹 2026-03-12 3/150 2026-03-13 11:02 by 求调剂zz
[考研] 321求调剂（食品/专硕） +3 xc321 2026-03-12 6/300 2026-03-13 08:45 by xc321
[考博] 26读博 +4 Rui135246 2026-03-12 10/500 2026-03-13 07:15 by gaobiao
[考博] 2026年博士申请 +3 QwQwQW10 2026-03-11 3/150 2026-03-12 17:58 by gxch43
[考研] 大连大学化学专业研究生调剂 +3 琪久. 2026-03-10 8/400 2026-03-11 10:02 by 琪久.
信息提示
关闭
请填处理意见
关闭