求基本小波的变换矩阵的问题

mac_allen

我想使用基本小波的输入输出来得到它的变换矩阵，该小波变换函数应该是正交变换阵，满足T'T=I
我自己编了一个小波的变换函数：
function y=mbw(x,M)
y=[];
for i=1:M
    [a,b]=dwt(x,'db1');
    if rem(length(x),2)~=0
        b2=b(1:floor(length(x)/2));
        x=a;
    else
        x=a;
        b2=b;
    end
    y=[b2,y];
end
y=[a,y];

使得对于任意长度的输入均可以实现M-band的小波变换（自己的理解，不知道对不对）

通过9组1×９的行向量，得到的变换矩阵的T'T=
  Columns 1 through 4
    7.031250000040101e-002    1.546796083857247e-001   -1.308127261877218e-012   -4.063827085883133e-013
    1.546796083857247e-001    5.625000000028910e-001   -4.858158815950563e-012   -1.404424453952082e-012
   -1.308127261877218e-012   -4.858158815950563e-012    1.000000000003179e+000    1.162290659174126e-012
   -4.063827085883133e-013   -1.404424453952082e-012    1.162290659174126e-012    1.000000000000376e+000
   -8.789952890562265e-013   -3.308337867282228e-012    1.393302140334116e-012    6.018093082740653e-013
    1.097494387551132e-013    1.324561261763248e-013   -1.215662244491146e-013    1.202188392354075e-014
    1.345991309145983e-013    4.687104189290466e-013   -8.518183681299106e-013   -1.770805724277119e-013
   -3.339200091176912e-013   -1.419131500687587e-012    1.046477270669034e-012    3.851687939957369e-013
   -6.229395689807169e-013   -1.035803068438373e-012    2.905453655443967e-012    7.279206962356164e-013
  Columns 5 through 8
   -8.789952890562265e-013    1.097494387551132e-013    1.345991309145983e-013   -3.339200091176912e-013
   -3.308337867282228e-012    1.324561261763248e-013    4.687104189290466e-013   -1.419131500687587e-012
    1.393302140334116e-012   -1.215662244491146e-013   -8.518183681299106e-013    1.046477270669034e-012
    6.018093082740653e-013    1.202188392354075e-014   -1.770805724277119e-013    3.851687939957369e-013
    1.000000000000810e+000   -5.037830617207250e-013   -2.936568917801643e-013    1.318983109679414e-013
   -5.037830617207250e-013    1.000000000000255e+000   -4.087286437231033e-013    2.458193227644081e-013
   -2.936568917801643e-013   -4.087286437231033e-013    1.000000000000000e+000   -5.191835702181283e-013
    1.318983109679414e-013    2.458193227644081e-013   -5.191835702181283e-013    1.000000000000600e+000
    2.499056517279946e-012   -2.277136563486077e-013    2.376566191121050e-013    6.126714537493439e-013
  Column 9
   -6.229395689807169e-013
   -1.035803068438373e-012
    2.905453655443967e-012
    7.279206962356164e-013
    2.499056517279946e-012
   -2.277136563486077e-013
    2.376566191121050e-013
    6.126714537493439e-013
    1.000000000000178e+000
发现第一行第二列的元素不满足上面的条件，其他的都是可以的
不知道为什么
望各位解答