Diminuer le temps d'exécution d'un programme

**air_doux** · 31/05/2008, 00h31

salut,j'ai fait ce petit programme sous matlab.

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
Mat=imread('imag.bmp');
[H,W]=size(Mat)
 for i=1:1:H
      for j=1:1:W
      s=(i-1)*W+j;
      x(s,1)=i;
      x(s,2)=j;
      u(s)=1;
      if Mat(i,j)<210
      u(s)=-1;
      end
      end
 end
for  i=1:1:H*W
     for j=1:1:W*W
      G(i,j)=u(i)*u(j)*(x(i,:)*x(j,:)'+1)^4;
     A(i,j)=0;
     if i==j
         A(i,j)=-1;
     end
     end
  Aeq(i)=u(i); 
  e(i)=1;
  B1(i)=0;
  x1(i)=0;
  Ib(i)=0;
  Ub(i)=1000;
 end
 B=B1';
Ib=Ib';
 Ub=Ub';
Beq=0;
x0=x1';
[a,fval] = fmincon(@(a)1/2*a'*G*a-e*a,x0,A,B,Aeq,Beq,Ib,Ub)
n=0;
for i=1:1:100
    n=n+a(i)*u(i)*(x(i,:)*x(1,:)'+1)^4;
end
b=u(1)-n;
for  h=1:1:10     
A=0;
B=0;
C=0;
D=0;
E=0;
for i=1:1:100
     A=A+a(i)*u(i)*x(i,2)^4;
     B=B+a(i)*u(i)*3*(h*x(i,1)+1)*x(i,2)^3;
     C=C+a(i)*u(i)*6*(h*x(i,1)+1)^2*x(i,2)^2;
     D=D+a(i)*u(i)*((h*x(i,1)+1)^2*x(i,2)+3*(h*x(i,1)+1)^3*x(i,2)+(h*x(i,1)+1)*x(i,2)^3);
     E=E+a(i)*u(i)*(h*x(i,1)+1)^4;
end
    E=E+b;
    P=[A B C D E];
    r=roots(P);
    for i=1:1:4
        if imag(r(i))==0 
            if r(i)>=0
    figure(1),plot(r(i),h,'og'),hold on;
            end
        end
    end
end
figure(2),imshow(Mat);

mon problème suppose sur la taille de l'image pour des valeurs de H*W>100
SVP quelqu'un peux m'orienter et m'aider ???

**Jerome Briot** · 31/05/2008, 09h11

Quel est le problème ?
Que fait ce code ?

**air_doux** · 31/05/2008, 16h05

salut,ce code fait la segmentation des images.
je l'ai executé pour H*W=100 et ça marche trés bien mais pour H*W superieur le temps d'execution devenu trés tres grand

et j'ai pas pu le minimisé.

**Caro-Line** · 31/05/2008, 17h05

Pour commencer essayes d'utiliser à fond les possibilités de l'editeur de MATLAB : par exemple l'outil M-Lint : Menu Tools > M-Lint.
Il te donnera plein d'informations, comme par exemple le fait qu'il faut initialiser tes variables en dehors des boucles (pré-allocation).
Pour cela tu peux t'aider des fonctions comme ZEROS ou ONES (ou EYE).
Par exemple tu as à un moment :

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
for  i=1:1:H*W
    for j=1:1:W*W
        G(i,j)=u(i)*u(j)*(x(i,:)*x(j,:)'+1)^4;
        A(i,j)=0;
        if i==j
            A(i,j)=-1;
        end
    end
    Aeq(i)=u(i);
    e(i)=1;
    B1(i)=0;
    x1(i)=0;
    Ib(i)=0;
    Ub(i)=1000;
end

Tu peux le modifier en :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
%pre-allocations
A=eye(H*W,W*W);
e=ones(1,H*W);
B1=zeros(1,H*W);
x1=B1;
Ib=B1;
Ub=1000*e;
Aeq=u;
for  i=1:1:H*W
    for j=1:1:W*W
        G(i,j)=u(i)*u(j)*(x(i,:)*x(j,:)'+1)^4;
    end
end

Il est même possible que le calcul de G soit "vectorisable" mais là il faut réfléchir un peu plus.