Discussion :

[compaq610]

Salut tout le monde
j'essaie de faire un programme en MATLAB qui fait la décomposition d'un signal par ondelettes, mais ça ne marche pas.
Voilà mon code source
j'attends votre réponse, merci

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
clear all
fid1 = fopen('1.txt', 'r');
x=fscanf(fid1,'%f', [1 Inf]);
N=256;
figure(1)
plot(x)
title('signale originale')
    a=2;
  for b=0:N-1
    for t=1:N  
w(t)=(1/sqrt(a)).*exp(-((t-b)./a).^2).*cos(5.*(t-b)/a);
    end
    figure(6)
    plot(w)
    g=w.*x;
         G(1,b+1)=0;
         for i=1:N
         G(1,b+1)=G(1,b+1)+g(i);
         end  
end
figure(2)
plot(G(1,:))
 a=4;
  for b=0:N-1
    for t=1:N  
w(t)=(1/sqrt(a)).*exp(-((t-b)./a).^2).*cos(5.*(t-b)/a);
    end
    figure(6)
    plot(w)
    g=w.*G(1);
         G(2,b+1)=0;
         for i=1:N
         G(2,b+1)=G(2,b+1)+g(i);
         end  
end
figure(3)
plot(G(2,:))
a=8;
  for b=0:N-1
    for t=1:N  
w(t)=(1/sqrt(a)).*exp(-((t-b)./a).^2).*cos(5.*(t-b)/a);
    end
    figure(6)
    plot(w)
    g=w.*G(2);
         G(3,b+1)=0;
         for i=1:N
         G(3,b+1)=G(3,b+1)+g(i);
         end  
end
figure(4)
plot(G(3,:))

[littledrummer]

beaucoup de chose qui ne vont pas dans ton programme.

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
 
clear all
close all;
x=rand(1,256)
N=256;
figure(1)
plot(x)
title('signale originale')
    a=2;
for b=0:N-1
    for t=1:N  
        w(t)=(1/sqrt(a)).*exp(-((t-b)./a).^2).*cos(5.*(t-b)/a);
    end
    g=w.*x;
    figure(6)
    subplot(311)
    plot(w)
    subplot(312)
    plot(x)
    subplot(313)
    plot(g) 
    G(1,b+1)=sum(g);
 end
figure(2)
plot(G(1,:))
 a=4;
  for b=0:N-1
    for t=1:N  
        w(t)=(1/sqrt(a)).*exp(-((t-b)./a).^2).*cos(5.*(t-b)/a);
    end
    figure(6)
    subplot(311)
    plot(w)
    subplot(312)
    plot(x)
    subplot(313)
    plot(g) 
    g=w.*x;
    G(2,b+1)=sum(g);
 
end
figure(3)
plot(G(2,:))
a=8;
  for b=0:N-1
    for t=1:N  
        w(t)=(1/sqrt(a)).*exp(-((t-b)./a).^2).*cos(5.*(t-b)/a);
    end
    figure(6)
    subplot(311)
    plot(w)
    subplot(312)
    plot(x)
    subplot(313)
    plot(g) 
    g=w.*x;
    G(3,b+1)=sum(g);
  end
 
figure(4)
subplot(311)
plot(G(1,:))
subplot(312)
plot(G(2,:))
subplot(313)
plot(G(3,:))

tu ne peux pas décomposer les détails de la décomposition précédente. c'est ton signal que tu décompose. Regarde de plus près les équations de la CWT.

y a encore moyen d'optimiser le code. par exemple tu calcules a chaque fois ton ondelette alors qu'elle est juste translatée. tu la calcules au début une seul fois et dans les boucles tu la translates. Encore mieux sert toi du calcule matriciel. La ton signal est court, mais si sa taille augmente tu va avoir des temps de calcul vraiment longs.

[sigma10]

ila faut faire l'etape de seuillage apres la decomposition et enfin le codage.