【求助】关于离散信号求导的问题

diying813

我是课题是做拖拉机主动减振的，然后建立的数学模型涉及到路面输入的一阶导数和二阶导数。用simulink建立的路面模型仿真图如下，然后也把仿真的数据导入TXT文件了，就是相当于是一连串的离散点。然后需要对这个路面输入求导，只要低频的，10HZ以下就行，我应该怎么做呢，步长取得1e-3，之前直接用相邻点的差除以步长，得到的结果一个是会出现许多很大的峰值，再一个肯定非常不准确不能用。我应该怎样得到精确的导数呢？

江树空

diff函数不可以么？

diying813

江树空 发表于 2014-3-27 18:14
diff函数不可以么？

问题是只是一系列离散的点，都没有函数呢

跳动的杯子

diying813 发表于 2014-3-27 18:29
问题是只是一系列离散的点，都没有函数呢

diff也可用于离散点，叫做差分，求导的话可以用diff（x）./dy

马社

http://forum.vibunion.com/thread-130768-1-2.html
看牛主任的解释，挺对你的

simon21

这种情况要先对数据进行平滑处理，然后在对数据进行求导

江树空

diying813 发表于 2014-3-27 18:29
问题是只是一系列离散的点，都没有函数呢

请在matlab的help中查看完diff函数的用法后再提问

diying813

跳动的杯子 发表于 2014-3-27 22:39
diff也可用于离散点，叫做差分，求导的话可以用diff（x）./dy

我想问下，diff对于离散点，是不是就是相当于相邻两项相减再除以dt，可是这样不行的呀，会出现许多很大的尖峰，结果也不准确

diying813

simon21 发表于 2014-3-28 05:44
这种情况要先对数据进行平滑处理，然后在对数据进行求导

我想只留下低频的信号，想求教下是不是在这之前再滤波下呢？

diying813

江树空 发表于 2014-3-28 09:12
请在matlab的help中查看完diff函数的用法后再提问

之前没看到它用于离散点的情况，但是我想求教下，diff对于离散点，是不是就是相当于相邻两项相减再除以dt，可是这样不行的呀，会出现许多很大的尖峰，结果也不准确。

diying813

马社 发表于 2014-3-27 22:50
http://forum.vibunion.com/thread-130768-1-2.html
看牛主任的解释，挺对你的

可是直接差分得到的结果不能用啊，我之前虽然没有用diff这个函数，但是我写的语句应该差不多就是那个意思，我把相邻两项相减再处于dt，但是会出现许多大峰值，结果肯定也不对

江树空

diying813 发表于 2014-3-28 11:17
之前没看到它用于离散点的情况，但是我想求教下，diff对于离散点，是不是就是相当于相邻两项相减再除以dt ...

嗯，你这样求导从算法上来讲是没错的，结果不好只能是因为数据本身就不好了。你那个图是时域信号对吧？再就是你说只要10Hz以下的的这里我不是太明白，如果说你是要以求导结果为输入，那我觉得你应该在求完导之后，还要再进行低通滤波。

马社

diying813 发表于 2014-3-28 11:23
可是直接差分得到的结果不能用啊，我之前虽然没有用diff这个函数，但是我写的语句应该差不多就是那个意思 ...

还是要认真读一下牛主任的意思，要领会其中的意思

diying813

马社 发表于 2014-3-28 18:28
还是要认真读一下牛主任的意思，要领会其中的意思

求点拨

diying813

江树空 发表于 2014-3-28 17:13
嗯，你这样求导从算法上来讲是没错的，结果不好只能是因为数据本身就不好了。你那个图是时域信号对吧？再 ...

这个数据是根据标准路面谱的公式生成的，嗯是时域信号。然后我仿真的时候是这个数据，然后它的导数，二阶导数都需要用到，作为系统的输入。然后说是10HZ以下，是因为系统是分析路面输入对拖拉机振动的影响，只考虑低频振动。现在就是不知道怎么得到那个导数输入