小波变换中的信号扩展(延拓)问题

题目：小波变换中的信号扩展(延拓)问题

        本篇主要是做为上一篇《压缩感知稀疏基之离散小波变换》的一个补充说明。

        首先给出一个卷积的例子：

       

        从上面的卷积计算过程中可以看出，在计算y(0)和y(4)时，x(n)和h(n)只有一个点对应，在计算其它值时有两个点对应，以上计算还可以用矩阵表示如下：

        我们可以对x(n)按某种规律进行扩展，以使求卷积过程中x(n)和h(n)一直有两个点（h的长度）对应，这就是信号的扩展。

        观察可知，只要对x(n)左右各扩展一个点即可，在Matlab中有函数wextend对信号进行扩展，扩展的模式如下：

        注意最后一种’per’模式，因为我们在生成N×N的小波变换矩阵时用的是此模式，所以重点阐述这种模式。

        在Matlab分别运行两条命令，可得输出结果如下：

>> x=[1 2 34];wextend(1,’per’,x,1)

ans =

    4     1     2    3     4     1

>> x=[1 2 3];wextend(1,’per’,x,1)

ans =

    3     1     2    3     3     1

        命令wextend(1,’per’,x,1)中的输入参数中第一个“1”表示是一维扩展，最后一个“1”表示左右各扩展一个点。从上面的两个结果基本可知得到’per’模式的扩展规律：

        1）对于长度为偶数的序列，直接按周期进行扩展；

        2）对于长度为奇数的序列，先把序列延长一个点变为偶数后再按周期进行扩展。

        在上面的例子中，若对信号按’per’模式进行扩展，则计算过程如下：

        

        观察可以发现由h构成的矩阵各行实际上是一个循环移位的结果。

        我们可以在’per’扩展模式下写出一个较为复杂、但代表性更强的表达式，设x(n)长度为8，h(n)的长度为4，若对x(n)按’per’扩展模式扩展，则卷积表达式的矩阵形式为：

在Matlab中编程得到这个变换矩阵时，可以分成几个步骤：

1）参数初始化：N=length(x)，L=length(h)，在这里N=8，L=4

2）得到h=[h(0),h(1),h(2),h(3)]的反转序列，即h_1=[h(3),h(2),h(1),h(0)]，这一操作可使用命令“h_1 = fliplr(h);”实现

3）由h_1得到一个长为N的基本行向量：p1_0 = [h_1zeros(1,N-L)];

这里得到的p1_0实际上即是上面变换矩阵的第4行

4）通过对p1_0循环移位得到变换矩阵中的各行：

Matlab中有循环移位函数circshift，但这个函数是对行之间循环移位，而这里的行向量p1_0是1行N列，为了实现循环移位，可以对其先进行转置变为列向量再进行循环移位，移位结束后再转置变回行行量;

若要得到变换矩阵中的第一行：circshift(x’,-(L-1))，如：

>> x=[4 3 2 10 0 0 0]

x =

     4    3     2     1    0     0     0    0

>>circshift(x’,-(4-1))’

ans =

     1    0     0     0    0     4     3    2

同理可以得到其它各行。

                转载自彬彬有礼的专栏