bat3.m

编程实现基于小波变换尺度相关性去噪算法

%非正交的二次样条二进小波
h = [0 0 0.125 0.375 0.375 0.125 0 0].*sqrt(2);
g = [0 0 0 -1 1 0 0 0].*sqrt(2);
g1 = [0 0.015625 0.109375 0.34375 -0.34375 -0.109375 -0.015625 0].*sqrt(2);

%正交的二次样条二进小波
% h = [0 0 0 0 0.375 0.125 0.373].*sqrt(2);
% g = [0 0 0 0 0.5798 0.0869 0.0061].*sqrt(2);
%g1 = [0 0.015625 0.109375 0.34375 -0.34375 -0.109375 -0.015625 0].*sqrt(2);

%正交的三次样条二进小波
% h = [0 0.0625 0.2500 0.3750 0.2500 0.0625 0].*sqrt(2);
% g = [-0.00008 -0.01643 -0.10872 -0.59261 0 0.59261 0.10872 0.01643 0.00008].*sqrt(2);
%g1 = [0 0.015625 0.109375 0.34375 -0.34375 -0.109375 -0.015625
%0].*sqrt(2);

l = 512;
x = [1:l];
n = 7;
sa = 0;
sn = 5;
dsa = 2.5;
snum = 12;
noise = rand(1,l)-0.5;

% so(1:l/4) = 2;
% so(l/4+1:l*1/2) = 5;
% so(l*1/2+1:l*3/4) = 1;
% so(l*3/4+1:l) = 4;

signalName = 'doppler';
signalName = 'quadchirp';
signalName = 'blocks';
signalName = 'heavy sine';
signalName = 'bumps';
[so,s] = wnoise(signalName,9,2);
noise = s-so;

SSNFtest(so,noise,sa,sn,dsa,snum,h,g,g1,n,signalName);
% SSNFtest(so,noise,sa,sn,dsa,snum,h,g,g1,n,'pic_recw');

% so = sin(2*pi*x/l*4).*sin(2*pi*x/l*8).*sin(2*pi*x/l)+2.5;
% % 
% SSNFtest(so,noise,sa,sn,dsa,snum,h,g,g1,n,'pic_sinw');

% sa = 2.5
% so(1:l/4) = [1:l/4]./(l/4) + 2;
% so(l/4+1:l*1/2) = so(l/4:-1:1);
% so(l/2+1:l) = so(1:l/2);
% SSNFtest(so,noise,sa,sn,dsa,snum,h,g,g1,n,'pic_triw');