| 查看: 2658 | 回复: 21 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
shengjunjie木虫 (正式写手)
|
[求助]
如何求解描述振动的二阶微分方程
|
||
那位大侠给看看这种描述振动的二阶微分方程，困扰我一阵子了 其中C1=0.08；C2=2.3；e=0.3*(sin(335*t))^2 如何根据上面的公式求解λ和时间的关系，用matlab画出类似下面的图形  |
回复此楼» 猜你喜欢
真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗
已经有8人回复
-
寻求一种能扛住强氧化性腐蚀性的容器密封件
已经有5人回复
-
论文投稿，期刊推荐
已经有6人回复
-
请问哪里可以有青B申请的本子可以借鉴一下。
已经有4人回复
-
孩子确诊有中度注意力缺陷
已经有14人回复
-
请问下大家为什么这个铃木偶联几乎不反应呢
已经有5人回复
-
请问有评职称，把科研教学业绩算分排序的高校吗
已经有5人回复
-
2025冷门绝学什么时候出结果
已经有3人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有4人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
shengjunjie
木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 3613
- 散金: 208
- 红花: 2
- 帖子: 320
- 在线: 128.7小时
- 虫号: 988929
- 注册: 2010-04-04
- 性别: GG
- 专业: 橡胶及弹性体
|
我在公式里面加了阻尼那个选项，和别人的公式一样，但是感觉结果和他的差距很大啊（位移曲线一样，但是振幅曲线差距很大），不知道什么原因 这是原图  然后我根据你给我程序的结果是 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*410*t)).^2)./(1+0.08*x(1).^3)-(5*10^3)*x(1).^3*x(2)./(25+2*x(1).^3)]; [t,x]=ode45(odefun,[0 0.06],[1 0]);% 自变量的范围和初始值 figure(1); plot(t,x(:,1)) Y=x(:,1) pn=floor(log2(length(Y))); N=2^pn; S=Y(1:N); Y1 = fft(S);%进行N个点的fft变换 figure; f1 = 1/(t(2)-t(1))*(0:N/2)/N; Yabs=abs(Y1(1:N/2+1)); Ym=find(diff(sign(diff(Yabs)))<0)+1;%求幅值极大值的索引 Ymax=max(Yabs(Ym)); Yabs(Yabs plot(f1,Yabs(1:N/2+1)) hold on; plot(f1(Ym1),Yabs(Ym1),'r*'); axis([0 f1(Ym1(end)) 0 1.1*Ymax]); xlabel('频率 (Hz)') ylabel('振幅|A| (m)') 结果   |
13楼2012-04-17 13:53:19
dbb627
荣誉版主 (著名写手)
-
专家经验: +4 - 仿真EPI: 6
- 应助: 289 (大学生)
- 贵宾: 0.589
- 金币: 24640.4
- 散金: 551
- 红花: 61
- 沙发: 1
- 帖子: 1246
- 在线: 1794.8小时
- 虫号: 149791
- 注册: 2005-12-29
- 性别: GG
- 专业: 污染控制化学
- 管辖: 计算模拟
2楼2012-04-16 15:02:28
shengjunjie
木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 3613
- 散金: 208
- 红花: 2
- 帖子: 320
- 在线: 128.7小时
- 虫号: 988929
- 注册: 2010-04-04
- 性别: GG
- 专业: 橡胶及弹性体
3楼2012-04-16 20:34:29
shengjunjie
木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 3613
- 散金: 208
- 红花: 2
- 帖子: 320
- 在线: 128.7小时
- 虫号: 988929
- 注册: 2010-04-04
- 性别: GG
- 专业: 橡胶及弹性体
4楼2012-04-16 20:43:50