%---------------------------------------------------------------------------------------
%  exa060402.m,   例 6.4.2 和 圖 6.4.2  
%  利用半帶濾波器及對應的高通濾波器對一個含有噪聲的信號進行分解。
%  半帶濾波器h0來自程序exa070501
%----------------------------------------------------------------------------------------
clear;

load h0;
h1=qmf(h0,1)

NN=1000;w1=pi/4;
u=0.8*rand(1,NN);
for i=1:NN
   x(i)=sin(w1*i)+u(i);
end
subplot(311)
plot(1:150,x(1:150))
title(' test signal x(n)')

x10=filter(h0,1,x);
x11=filter(h1,1,x);
subplot(312)
plot(1:150,x10(101:250))
title('lowpass signal v0(n)')
subplot(313)
plot(1:150,x11(101:250))
title('highpass signal v1(n)')