怎么来对已经求出来的速度曲线进行求导？

猴王

我现在通过状态空间法，求解出来一个结构的速度和位移曲线，现在想通过对速度求导的方式来得到加速度，请各位好手指点，具体的表达式是怎么样的，这是个动力系统，我希望能通过求导来画出加速度的曲线

[ 本帖最后由 eight 于 2007-6-22 09:52 编辑 ]

eight

help diff

[ 本帖最后由 ChaChing 于 2010-1-12 17:26 编辑 ]

VibrationMaster

中心差分法:
0.5*(diff(y(1:end-1)+diff(2:end))/dt
这样比原始数据的头尾少两个点,可以自己补一下

TurtleRock

可以试试根据数据用polyfit求出多项式拟合，然后再diff得到加速度表达式，再作图

cheers1986

可以试试这个，用polyfit拟合出多项式，然后再用polyder求导

beyondhxf

如果有速度的表达式，可以用diff求导，如果只是数值解，可以使用VibrationMaster 提出的差分表示微分的方法，对于TurtleRock和cheers1986提出的先拟合再求导的方法，本人深表反对，虽然加了分，但是这种方法显然是太轻浮，首先拟合会造成误差，对拟合的结果求导会带来非常高的误差，在某些点上甚至和真实解相差甚远！

上当了，这怎么是2007年的帖子？

[ 本帖最后由 ChaChing 于 2010-1-12 17:06 编辑 ]

ChaChing

原帖由 beyondhxf 于 2010-1-12 12:53 发表

                               
登录/注册后可看大图

...对于TurtleRock和cheers1986提出的先拟合再求导的方法，本人深表反对，虽然加了分，但是这种方法显然是太轻浮，首先拟合会造成误差 ...

LS说的好, 但个人以为可不一定. 直接diff或先拟合再求导, 其实各有利弊, 应该视实际工程数据而定! 而且个人以为说diff也是一种拟合, 可能也不为过!

还有老帖没差吧! 前人走过的经验就是宝贝!
个人是较喜欢与回覆老帖的朋友讨论, 至少表示有自行搜索并看过资料!

beyondhxf

一般来说，如果速度曲线比较平滑，那么不管你用什么办法求加速度，都是很准确的，可是你看看楼主说的，动力系统，使用的是状态向量法，就是说基本变量不是一个量，而是一个向量，包含系统的多个状态，如位移速度等。一般来说动力系统的速度曲线不平滑，有大幅度的改变，就是说加速度变化比较快。对于这样的情况，使用polyfit拟合，凶多吉少，对本来的数据还可以，若是求导，肯定算飞了！