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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » MATLAB » 如何求解描述振动的二阶微分方程
北京石油化工学院2026年研究生招生接收调剂公告
221/3返回列表123下一页
查看: 2902  |  回复: 21
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

shengjunjie

木虫 (正式写手)

[求助] 如何求解描述振动的二阶微分方程

那位大侠给看看这种描述振动的二阶微分方程，困扰我一阵子了

其中C1=0.08；C2=2.3；e=0.3*(sin(335*t))^2
如何根据上面的公式求解λ和时间的关系，用matlab画出类似下面的图形
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
好好学习，多运动，多发Paper...
1楼 2012-04-16 14:30:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

dbb627

荣誉版主 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
shengjunjie: 金币+10, ★★★很有帮助 2012-04-17 09:35:06
臭水沟: 金币+1, 谢谢交流~~ 2012-04-17 10:11:24
引用回帖:
5楼: Originally posted by shengjunjie at 2012-04-16 20:54:38:
画出来了，如下图，你好，能否赐教如何画振幅，还有，我该如何给你赠送金币呢
ec/75/988929_1334580862_122.jpg

做下傅立叶变化画出幅频曲线求相应频率的振幅
赞一下(1人) 回复此楼
The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.
6楼2012-04-16 20:56:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

dbb627

荣誉版主 (著名写手)

【答案】应助回帖

★ ★
感谢参与，应助指数 +1
臭水沟: 金币+2, 谢谢交流~~ 2012-04-16 15:55:10
需要给初值条件
假设λ（0）=0.1 一阶导初值为0.01，
odefun=@(t,x)[x(2);3*x(2)./(2*x(1)*(1+0.08*x(1).^3))-2.3*(x(1).^4-x(1)+0.3*(sin(335*t)).^2)./(1+0.08*x(1).^3)];
[t,x]=ode45(odefun,[0 1],[0.1 0.01]);
plot(t,x(:,1))
赞一下 回复此楼
The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.
2楼2012-04-16 15:02:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shengjunjie

木虫 (正式写手)
引用回帖:
2楼: Originally posted by dbb627 at 2012-04-16 15:02:28:
需要给初值条件
假设λ（0）=0.1 一阶导初值为0.01，
odefun=@(t,x);
=ode45(odefun,,);
plot(t,x(:,1))

非常感谢。你好，我今天下午用了一个类似的方法计算出来了。忘记写初始条件了，λ（0）=1，λ（0）的导数为0.
你知道怎么根据这个方程求振幅不：其中C1=0.08；C2=1.342*10^7；e=9.23*(sin(2*pi*f*t))^2.
赞一下 回复此楼
好好学习，多运动，多发Paper...
3楼2012-04-16 20:34:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shengjunjie

木虫 (正式写手)
引用回帖:
2楼: Originally posted by dbb627 at 2012-04-16 15:02:28:
需要给初值条件
假设λ（0）=0.1 一阶导初值为0.01，
odefun=@(t,x);
=ode45(odefun,,);
plot(t,x(:,1))

非常感谢你的答复，太好啦。

下面是初始条件和各个参数

怎么根据他们求振幅
赞一下 回复此楼
好好学习，多运动，多发Paper...
4楼2012-04-16 20:43:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shengjunjie

木虫 (正式写手)
引用回帖:
2楼: Originally posted by dbb627 at 2012-04-16 15:02:28:
需要给初值条件
假设λ（0）=0.1 一阶导初值为0.01，
odefun=@(t,x);
=ode45(odefun,,);
plot(t,x(:,1))

画出来了，如下图，你好，能否赐教如何画振幅，还有，我该如何给你赠送金币呢
赞一下 回复此楼
好好学习，多运动，多发Paper...
5楼2012-04-16 20:54:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shengjunjie

木虫 (正式写手)
引用回帖:
6楼: Originally posted by dbb627 at 2012-04-16 20:56:55:
做下傅立叶变化画出幅频曲线求相应频率的振幅

你好，该如何做傅里叶变换，直接针对求出的lambda，还是怎么样，能否给个基本的程序，感谢啦
赞一下 回复此楼
好好学习，多运动，多发Paper...
7楼2012-04-17 09:36:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dbb627

荣誉版主 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
shengjunjie: 金币+10, ★★★很有帮助 2012-04-17 10:23:04
xiegangmai: 金币+1, 谢谢参与！ 2012-04-17 12:48:56
引用回帖:
5楼: Originally posted by shengjunjie at 2012-04-16 20:54:38:
画出来了，如下图，你好，能否赐教如何画振幅，还有，我该如何给你赠送金币呢
ec/75/988929_1334580862_122.jpg

把你画这个图的代码，贴一下，由于你的条件没有给f，我画不出这样的图
得到数据后，我可以给你做下fourier分析
赞一下 回复此楼
The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.
8楼2012-04-17 10:16:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shengjunjie

木虫 (正式写手)
引用回帖:
6楼: Originally posted by dbb627 at 2012-04-16 20:56:55:
做下傅立叶变化画出幅频曲线求相应频率的振幅

你好，这样编程对不对
%其中lambda（0）=1，D（lambda（0））=0
odefun=@(t,x)[x(2);3*x(2).^2/(2*x(1)*(1+0.08*x(1).^3))-1.342*10^7*(x(1).^4-x(1)+(61.98/6.71)*(0.185*sin(2*pi*335*t)).^2)./(1+0.08*x(1).^3)];
[t,x]=ode45(odefun,[0 0.06],[1 0]);% 自变量的范围和初始值
figure(1);
plot(t,x(:,1))
figure(2);
y=fft(x(:,1));
plot(abs(y)); % 画幅频曲线
figure(3);
plot(angle(y)); % 画相频曲线

在你的基础上直接加了这句话，
figure(2);
y=fft(x(:,1));
plot(abs(y)); % 画幅频曲线
figure(3);
plot(angle(y)); % 画相频曲线

这样的结果好像不对
但是函数中的2*pi*f=335是固定的
赞一下 回复此楼
好好学习，多运动，多发Paper...
9楼2012-04-17 10:21:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shengjunjie

木虫 (正式写手)
引用回帖:
8楼: Originally posted by dbb627 at 2012-04-17 10:16:54:
把你画这个图的代码，贴一下，由于你的条件没有给f，我画不出这样的图
得到数据后，我可以给你做下fourier分析

频率f的范围是0-1000hz
赞一下 回复此楼
好好学习，多运动，多发Paper...
10楼2012-04-17 10:22:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 shengjunjie 的主题更新
221/3返回列表123下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085600，321分求调剂 +9 大馋小子 2026-04-04 10/500 2026-04-04 23:42 by lqwchd
[考研] 材料调剂 +4 dxy调剂 2026-04-04 4/200 2026-04-04 22:04 by hemengdong
[考研] 材料专硕 调剂 +13 CXN123456 2026-04-03 13/650 2026-04-04 21:59 by hemengdong
[考研] 数一英一274机械调剂 +4 星陨流霞 2026-04-04 5/250 2026-04-04 09:35 by liucky
[考研] 283分材料与化工求调剂 +29 罗KAKA 2026-04-02 29/1450 2026-04-03 23:56 by userper
[考研] 一志愿中国石油大学化学工程323分求调剂 +4 化工专硕323分 2026-04-03 6/300 2026-04-03 22:12 by dongzh2009
[考研] 兽医调剂 +3 wh119216 2026-04-02 3/150 2026-04-03 19:34 by zrongyan
[考研] 求调剂机会 +5 意染ivy 2026-04-03 5/250 2026-04-03 15:13 by qoooooo614
[考研] 建环，能源，土木老师路过看一看！！！ +5 嘿嘿uu 2026-04-01 5/250 2026-04-03 11:47 by znian
[考研] 土木水利328分求调剂 +6 疾风知劲草666 2026-04-02 6/300 2026-04-03 11:38 by znian
[考研] 325分化学调剂 +5 15771691647 2026-04-02 5/250 2026-04-03 09:58 by ChemPharm
[考研] 一志愿山东大学，085600，344 +7 魏子per 2026-04-02 8/400 2026-04-02 21:12 by 百灵童888
[考研] 一志愿上海海洋大学083200食品学硕，求调剂，接受其他专业 +6 what张 2026-04-01 7/350 2026-04-02 16:48 by zzsw+
[考研] 318求调剂，计算材料方向 +10 吸喵有害笙命 2026-04-01 11/550 2026-04-02 16:29 by oooqiao
[考研] 354求调剂 +4 lxb598 2026-03-31 5/250 2026-04-02 09:55 by Jaylen.
[考研] 303分 0807学硕求调剂 +3 TYC3632 2026-04-01 3/150 2026-04-01 19:24 by lwk2004
[考研] 生物与医药考研调剂 +5 铁憨憨123425 2026-03-31 5/250 2026-04-01 18:01 by syh9288
[考研] 调剂 +3 好好读书。 2026-04-01 3/150 2026-04-01 17:06 by zhouyuwinner
[考研] 福建理工大学材料学院先进合金团队招收考研调剂学生 +3 大华金商都 2026-03-30 4/200 2026-03-31 01:04 by 方英俊602
[考研] 293求调剂 +3 末未mm 2026-03-30 5/250 2026-03-30 17:23 by 王保杰33
信息提示
关闭
请填处理意见
关闭