%Filt=mwOpeninf('haar.txt',5);
%Filt=mwOpeninf('sym5.txt',5);
Filt=mwOpeninf('db3.txt',5);
x=mwOpeninf('Test_Sig.txt',8);
x1=Canal(x,1);
NS=100;
x1frag=Fragment(x1,1,NS);
LoD = Canal(Filt, 1);
HiD = Canal(Filt, 2);
LoR = Canal(Filt, 3);
HiR = Canal(Filt, 4);
%[cA2, cD1, cD2]=mwave2(x1frag,LoD,HiD);
[cA1, cD1] = mydwt(x1frag,LoD,HiD);
[cA2, cD2] = mydwt(cA1,LoD,HiD);
NS2=length(cA1);
A1=myidwt(cA2,cD2,LoR,HiR,NS2);
cA0=myidwt(A1,cD1,LoR,HiR,NS);

% A11=idwt(cA2,cD2,LoR,HiR,NS2);
% cA01=idwt(A11,cD1,LoR,HiR,NS);

cA1=myidwt(cA2,cD2,LoR,HiR,52);
cA0=myidwt(cA1,cD1,LoR,HiR,NS);
subplot (211); plot(x1frag); title('Исходный сигнал')
axis([1 NS min(x1frag) max(x1frag)]); 

subplot (212); plot(cA0); title('Аппроксимация уровня 2')
axis([1 NS min(cA0) max(cA0)]); 

% subplot (212); plot(cA01); title('Аппроксимация уровня 2')
% axis([1 NS min(cA01) max(cA01)]); 
