récupérer les matrices

Version imprimable

Bonjour à tous!

J'ai un petit soucis avec mon code, je tente de trouver mon erreur mais impossible de mettre la main dessus. J'ai cette erreur qui survient :

Attempted to access v(2,:); index out of bounds because size(v)=[1,5].

Error in CreateAcq (line 47)
        xv = [xv v(i,:)]

voici mon 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
 % We want to create an acquisition
% We request for  BCD
% We know A(x,y)
%  a and b side
%  alpha is random number
%
%
 
k=input('Write the number of satellites you wish : ');
sat=[];
for j=1:k
x=input('Write the first coordinate (West-Est) : ');
y=input('Write the second coordinate (North-South) : ');
 
p(k,:)=[x y];
 
if j==1
sat=[p(k,:)];
else
    sat=[sat;p(k,:)];
end
 
end
xv=[];
yv=[];
 
for i=1:size(sat,1)
p1=p(k,:);
a = 15;
b = 5;
alpha=-45;
% Formula to find B and C
 pa = [x+b*cos(alpha) y+b*sin(alpha)];
 pb = [x-a*sin(alpha) y+a*cos(alpha)];
% D = C + (B-A);
pc = [pb(1)+(pa(1)-p1(1)) pb(2)+(pa(2)-p1(2))];
 
% Creating the area
v=[p1(1) pa(1) pc(1) pb(1) p1(1)];
w=[p1(2) pa(2) pc(2) pb(2) p1(2)];
 
%M(i,:)=[v w]
 
 if i==size(sat,1)
        % if the last satellite do this
        xv = [xv v(i,:)]
        yv = [yv w(i,:)]
 
 
    else
        % else add
        xv = [xv v(i,:) NaN]
        yv = [yv w(i,:) NaN]
 
 end
 
end
 
[X,Y,xq,yq] = Init_Region;
% showing the acquisition
in = inpolygon(X,Y,xv,yv);
figure
axis equal
hold on
plot(xv,yv,'LineWidth',2) % acquisition
plot(X(in),Y(in),'rp') % points inside
plot(X(~in),Y(~in),'b+') % points outside
hold off

13/03/2015, 17h40
Gooby

Bonjour,

C'est très simple, tu tentes de récupérer la 2eme ligne d'un vecteur à une ligne et 5 colonnes.
remplaces par v(:,2)
(ce qui dans ton cas reviens à faire v(2))

Ou alors, tu veux réellement récupérer toute la ligne et dans ce cas je pense que tu as oublié de faire ceci:

Code:

v=cat(1,v,[p1(1) pa(1) pc(1) pb(1) p1(1)]);
13/03/2015, 17h42
Gooby

Un petit conseil de programmation: Lorsque tu fais face à une erreur de ce type, je te conseille de mettre un point d'arrêt à la ligne problématique et de vérifier le contenu des variables pour être sûr qu'elles contiennent bien ce que tu veux qu'elles contiennent. Si elles ne contiennent pas ce que tu désires, tu peux vérifier ligne par ligne ton code afin de toujours vérifier le contenu de tes variables.
13/03/2015, 19h59
Korithov

Merci à vous 2 pour vos réponses. Ce problème est résolu en changeant v(i;:) par v(:;i) ;)mais un problème survient ensuite, lors de l'affichage des rectangle j'ai qu'une droite qui apparait, en remontant l'historique, je vois que p1=p(k; :) pause problème donc je change k par i, là j'ai juste une grille :?

Je voudrais que p1 prenne les valeurs de chaque pk = [x y]

Désolé de vous retenir mais je suis en gros newbie ^^ il y a encore 2jours je ne savais à peine déclarer une variable. J'apprends sur le tas
16/03/2015, 10h02
Gooby

Bonjour,

Il va falloir être un peu plus précis que cela:

Citation:

lors de l'affichage des rectangle j'ai qu'une droite qui apparait

De quels rectangles parles tu? Montre nous le code associé au tracé des rectangles.

Citation:

je vois que p1=p(k, :) pause problème donc je change k par i, là j'ai juste une grille

Qu'est-ce qui pose problème exaactement? quel est l'erreur? Qu'est-ce que tu entends par "j'ai juste une grille" ?

PS: v(i,:) par v(:,i)

2 pièce(s) jointe(s)

A la fin du code, je souhaite faire apparaitre les zones crées pour chaque p(k,:)

Le soucis c'est que j'obtiens la grille de base dès que je met 2 satellites ou plus :

Pièce jointe 171844

Avec 1 satellite, j'ai toujours la même figure :

Pièce jointe 171846

Théoriquement, j'ai un(des) rectangle(s) et à l'intérieur de ce(s) rectangle(s) les croix deviennent rouges.

Je remet le code que j'ai actuellement :

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
69
70
% We want to create an acquisition
% We request for  BCD
% We know A(x,y)
%  a and b side
%  alpha is random number
%
%
 
k=input('Write the number of satellites you wish : ');
sat=[];
for j=1:k
x=input('Write the first coordinate (West-Est) : ');
y=input('Write the second coordinate (North-South) : ');
 
p(k,:)=[x y];
 
if j==1
sat=[p(k,:)];
else
    sat=[sat;p(k,:)];
end
 
end
fg =size(sat,1);
xv=[];
yv=[];
p1=[];
 
for i=1:fg
p1= p(i,:);
a = 15;
b = 5;
alpha=45;
% Formula to find B and C
 pa = [p1(1)+b*cos(alpha) p(2)+b*sin(alpha)];
 pb = [p1(1)-a*sin(alpha) p(2)+a*cos(alpha)];
% D = C + (B-A);
pc = [pb(1)+(pa(1)-p1(1)) pb(2)+(pa(2)-p1(2))];
 
% Creating the area
v(:,i)=[p1(1);pa(1);pc(1);pb(1);p1(1)];
w(:,i)=[p1(2);pa(2);pc(2);pb(2);p1(2)];
 
%M(i,:)=[v w]
 
 if i==fg
        % if the last satellite do this
        xv = [xv v(:,i)]
        yv = [yv w(:,i)]
 
 
    else
        % else add
        xv = [xv v(:,i) NaN]
        yv = [yv w(:,i) NaN]
 
 end
 
end
 
[X,Y,xq,yq] = Init_Region %Je charge la grille
% showing the acquisition
in = inpolygon(X,Y,xv,yv); %C'est l'ensemble des rectangles
figure
axis equal
hold on
plot(xv,yv,'LineWidth',2) % acquisition
plot(X(in),Y(in),'rp') % points inside
plot(X(~in),Y(~in),'b+') % points outside
hold off

Dans ta boucle relatives à ton nombre de satellites, tu écrases chacune de tes variable à chaque itération. Ne devrais tu pas faire un truc de ce genre?
Code:

1 2 3 4 5 6 7 8 9 10 11 12 13 for j=1:k x(j)=input('Write the first coordinate (West-Est) : '); y(j)=input('Write the second coordinate (North-South) : '); p(j,:)=[x y]; if j==1 sat(j)=[p(j,:)]; else sat(j)=[sat;p(j,:)]; end end
Et du coup, adapter la suite du code en conséquence?

1 pièce(s) jointe(s)

J'ai enfin (plutôt à nouveau) mes rectangles :mrgreen: , j'ai un petit bug sur un des 2 qui n'est pas vraiment un rectangle ça doit venir de ma formule :)

Pièce jointe 171857

J'ai suivi ton conseil et j'ai simplifié mon code aussi :

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
69
70
 % We want to create an acquisition
% We request for  BCD
% We know A(x,y)
%  a and b side
%  alpha is random number
%
%
 
k=input('Write the number of satellites you wish : ');
sat=[];
xv=[];
yv=[];
for i=1:k
x(i)=input('Write the first coordinate (West-Est) : ');
y(i)=input('Write the second coordinate (North-South) : ');
 
p(i,:)=[x(i) y(i)];
 
%if j==1
%sat=[p(j,:)];
%else
%    sat=[sat;p(j,:)];
%end
 
%end
 
 
%for i=1:k
%p(i,:);
%p1= [
a = 15;
b = 5;
alpha=45;
% Formula to find B and C
 pa = [x(i)+b*cos(alpha) y(i)+b*sin(alpha)];
 pb = [x(i)-a*sin(alpha) y(i)+a*cos(alpha)];
% D = C + (B-A);
pc = [pb(1)+(pa(1)-p1(1)) pb(2)+(pa(2)-y(i))];
 
% Creating the area
v(:,i)=[x(i);pa(1);pc(1);pb(1);x(i)];
w(:,i)=[y(i);pa(2);pc(2);pb(2);y(i)];
 
%M(i,:)=[v w]
 
 if i==k
        % if the last satellite do this
        xv = [xv;v(:,i)]
        yv = [yv;w(:,i)]
 
 
    else
        % else add
        xv = [v(:,i);NaN]
        yv = [w(:,i);NaN]
 
 end
 
end
 
[X,Y,xq,yq] = Init_Region
% showing the acquisition
in = inpolygon(X,Y,xv,yv);
figure
axis equal
hold on
plot(xv,yv,'LineWidth',2) % acquisition
plot(X(in),Y(in),'rp') % points inside
plot(X(~in),Y(~in),'b+') % points outside
hold off