08.10.2016 Views

Foundations of Data Science

2dLYwbK

2dLYwbK

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

Exercise 11.10 Compute an approximation to the scaling function that comes from the<br />

dilation equation<br />

φ(x) = 1 + √ 3<br />

φ(2x) + 3 + √ 3<br />

φ(2x − 1) + 3 − √ 3<br />

φ(2x − 2) + 1√ 3<br />

φ(2x − 3).<br />

4<br />

4<br />

4<br />

4<br />

Exercise 11.11 Consider f(x) to consist <strong>of</strong> the semi circle (x − 1 2 )2 + y 2 = 1 and y ≥ 0<br />

4<br />

for 0 ≤ x ≤ 1 and 0 otherwise.<br />

1. Using precision j = 4 find the coefficients for the scale functions and the wavelets<br />

for D 4 defined by the dilation equation<br />

φ(x) = 1 + √ 3<br />

4<br />

φ(2x) + 3 + √ 3<br />

4<br />

φ(2x − 1) + 3 − √ 3<br />

4<br />

2. Graph the approximation to the semi circle for precision j = 4.<br />

Exercise 11.12 What is the set <strong>of</strong> all solutions to the dilation equation<br />

φ(2x − 2) + 1√ 3<br />

φ(2x − 3)<br />

4<br />

φ(x) = 1 + √ 3<br />

φ(2x) + 3 + √ 3<br />

φ(2x − 1) + 3 − √ 3<br />

φ(2x − 2) + 1√ 3<br />

φ(2x − 3)<br />

4<br />

4<br />

4<br />

4<br />

Exercise 11.13 Prove that if scale functions defined by a dilation equation are orthogonal,<br />

then the sum <strong>of</strong> the even coefficients must equal the sum <strong>of</strong> the odd coefficients in the<br />

dilation equation. That is, ∑ c 2k = ∑ c 2k+1 .<br />

k<br />

k<br />

function<br />

= wavelets<br />

acc=32; %accuracy <strong>of</strong> computation<br />

phit=[1:acc zeros(1,3*acc)];<br />

c1=(1+3^0.5)/4; c2=(3+3^0.5)/4; c3=(3-3^0.5)/4; c4=(1-3^0.5)/4;<br />

for i=1:10<br />

temp=(phit(1:2:4*acc)+phit(2:2:4*acc))/2;<br />

phi2t=[temp zeros(1,3*acc)];<br />

phi2tshift1=[ zeros(1,acc) temp zeros(1,2*acc)];<br />

phi2tshift2=[ zeros(1,2*acc) temp zeros(1,acc)];<br />

phi2tshift3=[ zeros(1,3*acc) temp ];<br />

phit=c1*phi2t+c2*phi2tshift1+c3*phi2tshift2+c4*phi2tshift3;<br />

plot(phit)<br />

figure(gcf)<br />

pause<br />

end plot(phit) figure(gcf) end<br />

374

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!