【求助】matlab分析随机激励的动力方程

kb24lal · 发表于 2011-4-3 16:02:11

比如说很简单的一个单自由度弹簧振子系统:mx''+cx'+kx=f(t)
m、c、k均已知，在质量块上施加一个持续的随机激励f(t)（eg.wgn（1，10000，1）），这个f(t)由于是随机的，因此其显示表达式不能写出，让人困扰！
请问，版上各位高手这个问题怎么解决比较合适，前提是使用matlab。用什么方法可以得到系统在一段时间内（[0 10s]）的动力时程响应呢？

nwcwww · 发表于 2011-4-4 17:35:21

我的印象是对于非连续的输入函数，可以考虑先用spline去fit
不一定准确，五六年没用过ode45了

kb24lal · 发表于 2011-4-5 12:30:02

本帖最后由 kb24lal 于 2011-4-5 13:09 编辑

我的印象是对于非连续的输入函数，可以考虑先用spline去fit
不一定准确，五六年没用过ode45了
nwcwww 发表于 2011-4-4 17:35

呵呵，spline这个误差有些大啊。。还是三次厄密pchip好些～谢谢你了

账号		自动登录	找回密码
密码			注册

账号		自动登录	找回密码
密码			立即注册

【求助】matlab分析随机激励的动力方程

相关帖子