麻烦高手帮忙编下这段加窗滑动滤波的程序

ambelmm

计算对数谱差
先对数据提取谱信息(1261920个点,变为239*5280的矩阵,再用周期图法加矩形窗,窗宽120点,数据间距150点,每次滑动10个点)
我的理解是 每次用这样一个单元移动,120个点的窗+150点数据+120个点的窗 也就是单元包括两个窗和中间空出的点
该单元每次移动10个点 一直移动到这一列的最后 对每次移动中的两个窗做对数谱差,最后每列得到(5280-390)/10=489个点

大侠们帮忙看看思路正确不啊
还有我实在不知道该怎么遍这段程序,觉得有以下几个难点
1.把点变为矩阵后 怎么样分别对每列进行一系列处理啊
2.periodogram函数可以加窗 但加两个窗,并把两个窗数据相减

谢谢各位

[ 本帖最后由 eight 于 2007-7-23 12:46 编辑 ]

花如月

matlab里专门对矩阵进行邻域和块处理的函数nlfilter（滑动邻域操作）colfilt（列处理）roifilt2（区域滤波）