提升小波及其算法-基本程序-有说明

1

层的近似系数

ca1=lwtcoef('ca',xDec,lsnew,2,1);

%

提取第

2

层的近似系数

ca2=lwtcoef('ca',xDec,lsnew,2,2);

%

提取第

1

层的细节系数

cd1=lwtcoef('cd',xDec,lsnew,2,1);

%

提取第

2

层的细节系数

cd2=lwtcoef('cd',xDec,lsnew,2,2);

subplot(311);

plot(x);

title('

原始信号

');

subplot(323);

plot(ca1);

title('

第一层近似信号

');

subplot(324);

plot(ca2);

title('

第二层近似信号

');

subplot(325);

plot(cd1);

title('

第一层细节信号

');

subplot(326);

plot(cd2);

title('

第二层细节信号

');

%

获得

d b2

小波的提升方案

lsdb2=liftwave('db2');

%

显示提升方案

displs(lsdb2);

%

获得

sym2

小波的提升方案

lssym2=liftwave('sym2');

%

显示提升方案

displs(lssym2);

%

获得整数变换提升方案

lsdb2=liftwave('db2','Int2Int');

x=[1:10];

lwtx=lwt(x,lsdb2)

%

获得

H aar

小波的提升方案

lshaar=liftwave('haar');

% < /p>

将提升步骤

ELS

加入到提升方案中

els={'p',[-0.125 0.125],0};

lsnew=addlift(lshaar,els);

%

进行单层提升小波分解

load noisdopp;

x=noisdopp;

subplot(211);

plot(x);

ylabel('x');

%

实施提升小波变换

[cA,cD]=lwt(x,lsnew);

xRec=ilwt(cA,cD,lsnew);

err=max(max(abs(x-xRec)))

subplot(212);

plot(xRec);

ylabel('xRec');

els={'p',[-0.125 0.125],0};

lsnew=addlift(lshaar,els);

% 2

层提升小波分解

load noisdopp;

x=noisdopp;

xDec=lwt(x,lsnew,2);

%

重构近似信号和细节信号

a1=lwtcoef('a',xDec,lsnew,2,1);

a2=lwtcoef('a',xDec,lsnew,2,2);

d1=lwtcoef('d',xDec,lsnew,2,1);

d2=lwtcoef('d',xDec,lsnew,2,2);

%

检查重构误差

err=max(abs(x-a2-d2-d1))

subplot(311);

plot(x);

title('

原始信号

');

subplot(323);

plot(a1);

title('

重构第一层近似信号

');

subplot(324);

plot(a2);

title('

重构第二层近似信号

');

subplot(325);

plot(d1);

title('

重构第一层细节信号

');

subplot(326);

plot(d2);

title('

重构第二层细节信号

');

%

得到

Haar

小波的提升方案

lshaar=liftwave('haar');

% < /p>

将提升步骤

ELS

添加到提升方案中

els={'p',[-0.125 0.125],0};

lsnew=addlift(lshaar,els);

% 2

层提升小波分解

load noisdopp;

x=noisdopp;

subplot(211);

plot(x);

ylabel('x');

%

对信号实施整数提升小波变换

lshaarInt=liftwave('haar','int2int');

lsnewInt=addlift(lshaarInt,els);

[cAint,cDint]=lwt(x,lsnewInt);

%

实施提升小波变换

xRecInt=ilwt(cAint,cDint,lsnewInt);

errInt=max(max(abs(x-xRecInt)))

subplot(212);

plot(xRecInt);

ylabel('xRecInt');

%

使用

Haar

小波，得到相应的提升方案

lshaar=liftwave('haar');

% < /p>

添加

ELS

到提升方案

< br>

els={'p',[-0.125 0.125],0}

lsnew=addlift(lshaar,els);

%

对于简单信号，尺度为

1

进行