新手求助：如何用matlab对数组求导？

晗菲

在下刚刚用matlab不久，需要用希尔伯特黄变换对信号进行处理，已经把信号分解为一系列的固有模态函数（IMF），接下来要进行Hibert分析，信号的瞬时振幅和瞬时相位都可以由matlab自带的hilbert()求出，但是求瞬时频率时要对瞬时相位进行求导。我要处理的信号是在workspace加载进来的1000*1的数组，求出的瞬时相位也是1000*1的数组，没有明确的y关于x的表达式，请问各位高手如何求出这1000*1数组的导数啊？

[ 本帖最后由 晗菲 于 2006-10-29 11:35 编辑 ]

晗菲

:'(:'(:'(第一次到这个论坛发帖求助都末有人理，是不是偶的问题太简单了所以大家都不屑回答啊？同学说用插值试试，但是这么多插值法偶又不知道用哪种，烦呐，谁有好方法就告诉俺吧

eight

原帖由 晗菲 于 2006-10-29 11:22 发表
:'(:'(:'(第一次到这个论坛发帖求助都末有人理，是不是偶的问题太简单了所以大家都不屑回答啊？同学说用插值试试，但是这么多插值法偶又不知道用哪种，烦呐，谁有好方法就告诉俺吧

diff函数就可以完成求导，事实上这个是向前差分（还是向后，忘了：P），如果要中心差分，那么就自己写一个吧。另外，求hilbert变换后的瞬时频率可以直接利用flandrin的instfreq函数，具体上flandrin网页下载相关程序（参考我在 信号处理 板块写的文章：送给搞EMD或HHT但.....）

晗菲

谢谢eight，偶去找找看

visstic

谢谢！

hyl2323

提醒，直接用diff两点差分公式不好，我试过了，误差大，建议用5点差分公式，得到的结果明显好多了。