求基本小波的变换矩阵的问题
我想使用基本小波的输入输出来得到它的变换矩阵,该小波变换函数应该是正交变换阵,满足T'T=I我自己编了一个小波的变换函数:
function y=mbw(x,M)
y=[];
for i=1:M
=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=;
end
y=;
使得对于任意长度的输入均可以实现M-band的小波变换(自己的理解,不知道对不对)
通过9组1×9的行向量,得到的变换矩阵的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
发现第一行第二列的元素不满足上面的条件,其他的都是可以的
不知道为什么
望各位解答
页:
[1]