dyaddown

功能：对时间序列进行二元采样，每隔一个元素提取一个元素，得到一个降采样时间序列。

格式：

1.y = dyaddown(x, EVENODD)

当EVENODD=0时，从x中第二个元素开始采样（偶采样）；当EVENODD=1时，从x中第一个元素开始采样（奇采样）。

2.y = dyaddown(x)

EVENODD缺省，按EVENODD=0

dyadup

功能：对时间序列进行二元插值，每隔一个元素插入一个0元素，得到一个时间序列。

格式：

1.y = dyadup(x, EVENODD)

当EVENODD=0时，从x中第二个元素开始采样（偶采样）；当EVENODD=1时，从x中第一个元素开始采样（奇采样）。

2.y = dyadup(x)

EVENODD缺省，按EVENODD=0

interp

功能：对时间序列进行整数倍插值，使得时间序列曲线更光滑。

格式：

1.y = interp(x, r)

在x中插入一些数据，使得插值后的序列y的长度为x的r倍。

2.y = interp(x, r, l, alpha)

插值后得到的序列y的长度为x的r倍。

3.[y, b] = interp(x, r, l, alpha)

插值后同时得到一个低通插值滤波器的系数，长度为2rl+1.

说明：

x--时间序列

r--插入点的倍数

l--插值滤波器长度

alpha--滤波器的截止频率，0<alpha<=1，假设原序列的采样频率之半为1，缺省时l=4，alpha=0.5.

y--插值后得到的时间学列

b--低通插值滤波器的系数，长度为2rl+1

downsample

功能：对时间序列重采样，在原时间序列中等间隔地取出一些项，得到新序列。

格式：

1.y = downsample(x, n)

从第一项开始，等间隔n对x采样，得到的序列为y。

2.y = downsample(x, n, phase)

从第phase+1项开始，等间隔n对x采样，得到的序列为y，而0<=phase<n.

decimate

功能：对时间序列进行整数倍采样处理，使得时间序列的长度降低。

格式：

1.y = decimate(x, r)

将时间序列x的采样频率降低为原来的1/r，即length(y)=length(x)/r。在抽取之前，默认地采用了8阶chebyshevI型低通滤波器压缩频带。

2.y = decimate(x, r, n)

采用n阶chebyshevI型低通滤波器。

3.y = decimate(x, r, ‘fir’)

采用30阶的FIR型低通滤波器来压缩频带，对时间序列进行整数倍抽取。

4.y = decimate(x, r， n, ‘fir’)

指定当对时间序列进行整数倍抽取的时候，采用n点FIR型低通滤波器来压缩频带，对时间序列进行整数倍抽取。

说明：

x--时间序列

r--采样要降低的倍数

n--指定所采用的chebyshevI型低通滤波器的阶数

‘fir’--FIR滤波器

resample

功能：对时间序列进行重采样。

格式：

1.y = resample(x, p, q)

采用多相滤波器对时间序列进行重采样，得到的序列y的长度为原来的序列x的长度的p/q倍，p和q都为正整数。此时，默认地采用使用FIR方法设计的抗混叠的低通滤波器。

2.y = resample(x, p, q, n)

采用chebyshevIIR型低通滤波器对时间序列进行重采样，滤波器的长度与n成比例，n缺省值为10.

3.y = resample(x, p, q, n, beta)

beta为设置低通滤波器时使用Kaiser窗的参数，缺省值为5.

4.y = resample(x, p, q, b)

b为重采样过程中滤波器的系数向量。

5.[y, b] = resample(x, p, q)

输出参数b为所使用的滤波器的系数向量。

说明：

x--时间序列

p、q--正整数，指定重采样的长度的倍数。

n--指定所采用的chebyshevIIR型低通滤波器的阶数，滤波器的长度与n成比列。

beta--设计低通滤波器时使用Kaiser窗的参数，缺省值为5.

EVENODD=0

摘自新浪微博