Discussion :

[Globoxx]

Bonjour,

Je cherche à calculer un vecteur à partir d'une norme. Alors oui les possibilités de vecteurs sont infinies mais la boucle s'arrête dès qu'il a trouvé un vecteur correspondant.

Code :

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

1
2
3
4
5
6
7
8
9
10
11
 
% 9000 est la plus grand valeur de Distance_from_Origine donc la plus grande norme
 
Vector = zeros(N,3);
 
for i = 1:N
    Vector(i,1:3) = 9000*rand(1);
    while sqrt(Vector(i,1)^2+Vector(i,2)^2+Vector(i,3)^2) ~= Distance_from_Origine(i)
        Vector(i,1:3) = 9000*rand(1);
    end
end

Evidemment ce script est beaucoup trop lent... Y'aurait-il un autre moyen ?
Merci d'avance !

[shaiHulud]

Bonjour,

Je cherche à calculer un vecteur à partir d'une norme

Il faut préciser un peu car le vecteur (N,0,0,0,...) a pour norme N et réponds à ta question

[Globoxx]

Je vais faire un exemple.
Prenons l'univers, les normes seraient la distance entre la Terre et d'autres corps célestes.
Et bien je voudrais déterminer les coordonnées de ces corps célestes, mais pour se faire je dois d'abord les disséminer aléatoirement, avec seule information la distance qui doit les séparer de la Terre.
J'espère avoir été plus clair

[shaiHulud]

J'espère avoir été plus clair

Pas vraiment... cherches-tu à générer des vecteurs aléatoires ayant une norme donnée, ou à retrouver dans un ensemble de vecteurs ceux qui ont une norme donnée ?

[Globoxx]

Je cherche à générer des vecteur ayant des normes données oui voilà.

[shaiHulud]

Je cherche à générer des vecteur ayant des normes données

il suffit de générer des vecteurs quelconques, de les diviser par leurs normes, et de remultiplier par la norme souhaitée. Par exemple si tu veux 1000 vecteurs de coordonnées gaussiennes dans R^3

Code :

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

1
2
3
4
5
 
x= randn(1000,3);
norms= sqrt(sum(x.^2,2));
x= bsxfun(@times,x,1./norms); % vecteurs normés
x= x* norme_voulue;

Mais bien sur rien ne t'oblige à prendre des coordonnées gaussiennes, ni à respecter la symétrie sphérique

[Globoxx]

Merci beaucoup !