erreur d'utilisation de la quantification VQ

Version imprimable

09/11/2009, 17h10
kadjuv

erreur d'utilisation de la quantification VQ
Bonjour,

Je veux quantifier un signal audio de type (wav) et pour cela j'ai utilisé la quantification vectorielle (VQ)
Code:

1 2 3 4 5 6 7 8 9 10 11 12 13 fid= fopen('D:\1_MES COURS\oooooooooooooo\[MON MéMO\simulation\sondict.wav') fid = 3 >> s = fread(fid,'uint16'); >> length(s) ans = 185243
Le problème que MATLAB m'affiche une erreur lorsque exécute (VQ)
Code:

1 2 [codebook,distortion]=VQ(s,4,16,0.004)
Code:

1 2 3 4 5 ??? Error using ==> unknown Matrix dimensions must agree. Error in ==> VQ at 29 distort=sum((training_seq_block-updated_codebook(j,:)).^2);
Je suis un nouveau utilisateur de MATLAB :( .....

S'il y a quelqu'un qui pourrait m'aider?

Merci :)
09/11/2009, 21h08
Jean Dumoncel

Salut,

d'ou vient la fonction VQ?

voici le code

Code:

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
67
68
function [codebook,distortion]=VQ(training_seq,dimension,codebook_size,tolerance)
%VQ.m 	Vector quatizer design using K-means algorithm
%	[codebook,distortion]=VQ(training_seq,dimension,codebook_size,tolerance)
%	training_seq = training sequence
%	dimension = dimension of the quantizer
%	codebook_size = size of the codebook (rate=log2(codebook_size)/dimension)
%	tolerance = desired relative distortion (default=0.001)
%	Length of training_seq must be a multiple of the dimension of quantizer
 
if (nargin==3)
    tolerance=0.001;
end
m=round(length(training_seq)/dimension);
if (m*dimension-length(training_seq)<0)
    error('length of training_seq is not a multiple of dimension')
end
%	initialize the codebook
initial=training_seq(1:dimension*codebook_size);
initialcodebook=(reshape(initial,dimension,codebook_size))';
updated_codebook=initialcodebook;
%	first update
newdistortion=0;
distortion=0;
for i=1:m;
    training_seq_block=training_seq((i-1)*dimension+1:i*dimension);
    training_matrix(i,:)=training_seq_block;
    distortion_block=[ ];
    for j=1:codebook_size;
        distort=sum((training_seq_block-updated_codebook(j,:)).^2);
        distortion_block=[distortion_block distort];
    end
    [distortion_min,ind]=min(distortion_block);
    newdistortion=newdistortion+distortion_min;
    index(i)=ind;
end
for l=1:codebook_size;
    partition=(index==l);
    if sum(partition)>0
        updated_codebook(l,:)=partition*training_matrix/sum(partition);
    end
end
newdistortion=newdistortion/m;
%	furthur updated until the desired tolerance is met
while(abs(distortion-newdistortion)/newdistortion>tolerance)
    distortion=newdistortion;
    newdistortion=0;
    for i=1:m;
        training_seq_block=training_seq((i-1)*dimension+1:i*dimension);
        training_matrix(i,:)=training_seq_block;
        distortion_block=[ ];
        for j=1:codebook_size;
            distort=sum((training_seq_block-updated_codebook(j,:)).^2);
            distortion_block=[distortion_block distort];
        end
        [distortion_min,ind]=min(distortion_block);
        newdistortion=newdistortion+distortion_min;
        index(i)=ind;
    end
    for l=1:codebook_size;
        partition=(index==l);
        if sum(partition)>0
            updated_codebook(l,:)=partition*training_matrix/sum(partition);
        end
    end
    newdistortion=newdistortion/m;
end
codebook=updated_codebook;
distortion=newdistortion/dimension;

10/11/2009, 16h26
Jean Dumoncel

Ca vient peut-être de là :

Citation:

Length of training_seq must be a multiple of the dimension of quantizer

or tu as un training_seq de longueur 185243, et 185243/4 n'est pas un nombre entier, donc 185243 n'est pas un multiple de 4...

essaie en prenant : s = s(1:185240);
11/11/2009, 14h14
kadjuv

j'ai essayé avec un multiple de '4' mais sa marche pas, toujours le même problème
11/11/2009, 16h15
Jean Dumoncel

Essaie ceci :

Code:

[codebook,distortion]=VQ(s',4,16,0.004)

J'ai transposé s car apparemment, s doit être une matrice ligne et non une matrice colonne.

Pour trouver toi-même comment résoudre ce genre d'erreur, tu peux mettre un point d'arrêt dans ton mfile à la ligne ou se situe l'erreur (ici ligne 29 de ta fonction VQ). Quant tu lances la commande, le programme s'arrêtera à cette ligne, et comme le message d'erreur indique un problème de dimensions, tu pourras plus facilement regarder les dimensions des variables et ainsi résoudre le problème.
12/11/2009, 11h54
kadjuv

merci le problem est réglé j'ai essayé avec s' est ça marche
16/11/2009, 12h44
Jerome Briot

Citation:

Envoyé par kadjuv

merci le problem est réglé j'ai essayé avec s' est ça marche

Au passage, comme les valeurs ne sont pas complexes => :faq: : Quelle est la différence entre ' et .' ?