Variable appears to change size on every loop iteration - préallocation

**andryran** · 24/05/2013, 12h35

Bonjour à tous,je voudrais votre avis sur ce code.Nb j'utilise matlab R2012a.
Je veux transformer ce code matlab en langage C,pour cela je fais appel à l'outil coder de matlab.
Voici les deux codes:

Le fichier contenant la fonction

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
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
function [y,Y,E,xbar,Delta,Eder,Etrue,SNRspars,sigmanoise]= mygenerdatatest(M,N,nb_pic,NoiseLevel)
%#codegen
coder.inline('never');
% time samples
Tau = 0:N-1; Tau = Tau(:);
 
DeltaTau = 0.5;
Tau = min(max(0,Tau+DeltaTau*(rand(N,1)-1/2)),N-1);
 
%frequences sur la grille
Theta = (0:M-1)/M; Theta = Theta(:);
 
%position of peaks
if(nb_pic==4)
    pos = [100 150 200 400];
elseif(nb_pic==6)
    pos = [100 150 200 250 300 400];
else
    q = randperm(M);
	pos = q((1:nb_pic));
end
 
%perturbation
Delta = zeros(M,1);
if(nb_pic==4)
    Delta(pos) = [0.1 -0.4 -0.3 0.2]/M;
elseif(nb_pic==6)
    Delta(pos) = [0.1 -0.4 -0.3 0.2 -0.1 -0.2]/M;
else
    Delta(pos) = (0.1/M).*ones(M,1) ;      
end
 
%frequences
Nu = Theta + Delta;
o = Nu(pos);
 
%modulus of peaks
if(nb_pic==6)
    Al = [0.5 0.9 1 0.7 0.8 0.6]';
elseif(nb_pic==4)
    Al = [0.5 0.9 1 0.7]';
else
    Al = ones(nb_pic,1);
end
 
%phase of peaks
fi = rand(nb_pic,1).*(2*pi);
 
%exact sparse signal we want to estimate
xbar = zeros(M,1);
for ind = 1:length(o)
    ocur = o(ind);
    xbar(Nu == ocur) = Al(ind)*exp(1i*fi(ind));
end
 
%exact cisoid (before adding noise)
sig = 0;
%xx(k)=zeros(1,nb_pic);
for k =1:nb_pic;
    xx(k)= zeros(1,nb_pic);
end
for k = 1:nb_pic
    %xx(k)= zeros(1,nb_pic);
    xx(k) = Al(k)*exp(1i*fi(k));
    sig = sig + xx(k)*exp(1i*2*pi*Tau*o(k));
end
 
%IFFT linear operator
E = exp(1i*2*pi*Tau*Theta');
 
%derivative
Eder = ((1i*2*pi*Tau/M)*ones(1,M)) .* exp(1i*2*pi*Tau*Theta');
 
% exact frequency matrix
Etrue = exp(1i*2*pi*Tau*Nu');
 
%observation
Y = sig(:);
sigmanoise =  std(Y) * 10^(-NoiseLevel/20);
z = randn(2*N,1);
y = Y + (sigmanoise/sqrt(2)) * (z(1:N,1) + 1i * z(N+1:end,1));
 
SNRspars = 20*log10(min(Al)/sigmanoise);
 
end

Le fichier test ie fichier qui lance le programme en qlq sorte comme un main

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
%#codegen
 
M = 500;%Number of points in the fft
N = 50; %Number of observations
nb_pic = 4; %Number of peaks
NoiseLevel = 30;
% fonction
[y,Y,E,xbar,Delta,Eder,Etrue,SNRspars,sigmanoise]= mygenerdatatest(M,N,nb_pic,NoiseLevel);
 
% affichage
fprintf(1,'grid error = %g\n',norm(Y-E*xbar))
fprintf(1,'linearization error = %g\n',norm(Y-(E+M*(ones(N,1)*Delta').*Eder)*xbar))
fprintf(1,'modelling error = %g\n',norm(Y-Etrue*xbar))
fprintf(1,'SNR = %g dB\n',SNRspars)
fprintf(1,'sigmanoise = %g\n',sigmanoise)

Dans ce script il y a erreur que je n'arrive pas à résoudre pour lancer le coder de matlab.

Pour lancer le coder, matlab veut que tout les scripts sont transformer en fonction,j'ai essayé et c'est écrit là en haut,ensuite quand je fais appel à coder il y a encore une erreur qui dit que la variable xx(k) doit être définit mais je n'y arrive pas à le faire et la génération ne peut pas se faire.
Alors si quelqu'un a une idée s'il vous plaît aidez moi,merci beaucoup

**Merel** · 24/05/2013, 13h10

Salut,

il y a une erreur au niveau de la définition de xx.

Code :

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

1
2
3
for k=1:nb_pic;
    xx(k)= zeros(1,nb_pic);
end

Tu définies un scalaire comme étant égal à un vecteur.
Tu n'as pas besoin de boucle... Je suppose que tu souhaites définir un xx comme étant un vecteur nul:

Code :

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

xx(1:nb_pic)=zeros(1,nb_pic);

devrait convenir.

**Jerome Briot** · 24/05/2013, 13h30

Envoyé par Merel

Code :

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

xx(1:nb_pic)=zeros(1,nb_pic);

Plutôt

Code :

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

xx = zeros(1,nb_pic);

**andryran** · 24/05/2013, 14h51

Merel> il faut utiliser les indices car les fi aussi sont des vecteurs.
si on fait xx(1:nb_pic)=zeros(1,nb_pic); ça marche pas le compilateur de MATLAB n'accepte pas il souligne encore xx(k) en rouge.De plus si tu regardes bien en haut il y quelque code similaire à ce dernier en haut mais MATLAB accepte mais pour xx(k) je trouve pas de solution jusque là.

Invité · 24/05/2013, 15h02

Envoyé par andryran

il souligne encore xx(k) en rouge.

En rouge ou orange ? Et quel est le message ?
Si c'est l'avertissement

might be growing inside a loop

Le code de Merel / Dut règle justement ce "problème" et pour le comprendre je t'invite à lire ceci.

**andryran** · 24/05/2013, 15h06

Et oui il me répond the variable 'xx' appears to change size on every loop iteration. Consider preallocation for speed.Nb et il affiche en orange pour warning maintenant.

Invité · 24/05/2013, 15h10

Même réponse donc :

Envoyé par Winjerome

Le code de Merel / Dut règle justement ce "problème" et pour le comprendre je t'invite à lire ceci.

**andryran** · 24/05/2013, 15h17

J'ai déjà fait le préallocation mais pourquoi le compilateur marque en orange encore le xx(k)? c'est pas encore suffisant le bloc que j'ai alloué ou autres problèmes encore?

Invité · 24/05/2013, 15h19

As-tu fait cette préallocation avant ou après la ligne d'erreur ?

**andryran** · 24/05/2013, 15h32

avant le warning,mais j'ai lu sur le lien parlant de la pré-allocation que je cite " Utilisée avec des matrices de petites tailles, sa mise en place peut être plus longue que l'ensemble des affectations des éléments de la matrice."
Si tu as une idée le cette pré-allocation peut tu m'envoyer.

Invité · 24/05/2013, 15h41

Envoyé par andryran

Si tu as une idée le cette pré-allocation peut tu m'envoyer.

Ce n'est ni plus ni moins ce qui t'a été conseillé dès les début par Merel.

Envoyé par andryran

mais j'ai lu sur le lien parlant de la pré-allocation que je cite " Utilisée avec des matrices de petites tailles, sa mise en place peut être plus longue que l'ensemble des affectations des éléments de la matrice."

Tu ne l'as donc pas appliqué ? Alors c'est normal que tu obtiennes toujours ce message.
Si tu estimes que ce n'est pas nécessaire, garde à l'esprit que ce n'est qu'un avertissement.

**andryran** · 24/05/2013, 15h51

De ce fait ça n'affecte pas le résultat si je lance le coder avec cet avertissement?

Invité · 24/05/2013, 15h54

Non. Si nb_pic est grand cela peut juste engendrer un ralentissement, comme expliqué dans la FAQ.

**andryran** · 24/05/2013, 16h21

justement nb_pic = 4 c'est pas suffisant mais je vais augmenter cette valeur par une autre,merci pour ton aide.

**andryran** · 24/05/2013, 16h42

Code :

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

1
2
3
4
5
6
xx = complex(zeros(1,10));
for k = 1:nb_pic
    xx(1:nb_pic)=zeros(1,10);
    xx(k) = Al(k)*exp(1i*fi(k));
    sig = sig + xx(k)*exp(1i*2*pi*Tau*o(k));
end

Merci à vous c'est presque résolu avec la pré-allocation.

Variable appears to change size on every loop iteration - préallocation

MATLAB

Discussions similaires

Partager

Partager