3f62d540 |
1 | %Dessin ondelettes |
2 | %script |
3 | figure(1) |
4 | courbe = X2(:,1)-mean(X2(:,1)); |
5 | [C,L] = wavedec(courbe, 4,'haar'); |
6 | |
7 | subplot(6,1,1) |
8 | plot(1:7,courbe(1:7),'b','LineWidth',2) |
9 | axis([0 48 -1.6 1.6]) |
10 | hold on |
11 | plot(8:48,courbe(8:48),'b','LineWidth',2) |
12 | %plot(waverec([C(1:12)',zeros(1,36)],L,'haar'),'b--','LineWidth',2) |
13 | %courbe2 = X2(:,1); |
14 | %[C2,L2] = wavedec(courbe2,4,'haar'); |
15 | %plot(waverec([zeros(1,3),C2(4:12)',zeros(1,36)],L2,'haar'),'b-.','LineWidth',2) |
16 | hold off |
17 | ylabel('z','FontSize', 30) |
18 | set(gca, 'FontSize', 20) |
19 | subplot(6,1,2) |
20 | A4=zeros(1,48); |
21 | Coeff = zeros(5,48); |
22 | for r=1:16 |
23 | A4(r)=C(1)/power(2,5/2); |
24 | A4(r+16)=C(2)/power(2,5/2); |
25 | A4(r+32) = C(3)/power(2,5/2); |
26 | end |
27 | Coeff(5,8)=C(1); |
28 | Coeff(5,24) = C(2); |
29 | Coeff(5,40) = C(3); |
30 | stairs(A4,'b','LineWidth',2) |
31 | axis([0 48 -0.5 0.5]) |
32 | ylabel('A_4','FontSize', 30) |
33 | set(gca, 'FontSize', 20) |
34 | subplot(6,1,3) |
35 | D4=zeros(1,48); |
36 | for r=1:8 |
37 | D4(r) = C(4)/power(2,4/2); |
38 | D4(r+8) = -C(4)/power(2,4/2); |
39 | D4(r+16) = C(5)/power(2,4/2); |
40 | D4(r+24) = -C(5)/power(2,4/2); |
41 | D4(r+32) = C(6)/power(2,4/2); |
42 | D4(r+40) = -C(6)/power(2,4/2); |
43 | end |
44 | Coeff(4,8) = C(4); |
45 | Coeff(4,24) = C(5); |
46 | Coeff(4,40) = C(6); |
47 | stairs(D4,'b','LineWidth',2) |
48 | axis([0 48 -0.9 0.9]) |
49 | ylabel('D_4','FontSize', 30) |
50 | set(gca, 'FontSize', 20) |
51 | subplot(6,1,4) |
52 | D3=zeros(1,48); |
53 | for k=1:12 |
54 | for r=1:4 |
55 | D3(r+4*(k-1)) = (-1)^(k+1) *C(7+floor((k-1)/2))/power(2,3/2); |
56 | end |
57 | end |
58 | stairs(D3,'b','LineWidth',2) |
59 | ylabel('D_3','FontSize', 30) |
60 | axis([0 48 -0.5 0.5]) |
61 | set(gca, 'FontSize', 20) |
62 | subplot(6,1,5) |
63 | D2=zeros(1,48); |
64 | for k=1:24 |
65 | for r=1:2 |
66 | D2(r+2*(k-1)) = (-1)^(k+1) *C(13+floor((k-1)/2))/power(2,2/2); |
67 | end |
68 | end |
69 | stairs(D2,'b','LineWidth',2) |
70 | ylabel('D_2','FontSize', 30) |
71 | axis([0 48 -0.8 0.8]) |
72 | set(gca, 'FontSize', 20) |
73 | subplot(6,1,6) |
74 | D1=zeros(1,48); |
75 | for k=1:48 |
76 | for r=1 |
77 | D1(r+1*(k-1)) = (-1)^(k+1) *C(25+floor((k-1)/2))/power(2,1/2); |
78 | end |
79 | end |
80 | plot(D1,'b','LineWidth',2); |
81 | axis([0 48 -0.9 0.9]) |
82 | ylabel('D_1','FontSize', 30) |
83 | set(gca, 'FontSize', 20) |