calcul d'une formule

salut tous le monde

j'ai un formule que je peux le calculer sous Matlab ci joint

j'ai comme élément d'entré un pixel P( xi,yi,angle)

dans la formule il y a l'espérance et la variance, si possible de me signaler les fonctions qui les convient sous Matlab.

Merci d'avance

Salut!
Il te suffit de programmer ta formule plutôt que de chercher une fonction toute faite.
Jean-Marc Blanc

merci pour votre reponse sur le forum, mais vraiment je suis bloquee parceque j'ai comme element d'entré des points que j'ai extraire dans mon Image j'ai 5 point detecté P1(x1,y1)....P5(x5,y5)

et dans la formule il y a d'autre parametre que je sais pas comment je peux le trouver

si possible de me donner un coup de main

Citation:

---

Envoyé par hanane78

merci pour votre reponse sur le forum, mais vraiment je suis bloquee parceque j'ai comme element d'entré des points que j'ai extraire dans mon Image j'ai 5 point detecté P1(x1,y1)....P5(x5,y5)

et dans la formule il y a d'autre parametre que je sais pas comment je peux le trouver

si possible de me donner un coup de main

---

je pense que ca va etre difficile de vous aider car vous avez plus un probleme de comprehension de votre probleme que de programmation

A priori, les Ex et Ey correspondent aux moyenne par colonne et par ligne. Le sigma représente la variance.

Pour la moyenne

Code:

---

>> doc mean

---

pour la variance :

Code:

---

>> doc var

---

pour l'écart type :

Code:

---

1 2	>> doc std

---

Bonjour,
j'ai essayé de programmer la formule :

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12

I=[166 154 ; 240 100; 12 144 ; 12 252; 131 144 ; 131 252; 61 62 ; 88 88; 61 88 ; 88 62]; S1= std(I,0,1); S2= std(I,0,2); Ex = mean(I(:,1)); Ey = mean(I(:,2)); J = 1/ Ex*Ey; K = J( Ex^2*.S1 + Ey^2*.S2)

---

il me donne des erreurs

Citation:

---

??? Error using ==> plus
Matrix dimensions must agree.

Error in ==> GMC at 11
K = J( Ex^2.*S1 + Ey^2.*S2)

??? Error: File: GMC.m Line: 11 Column: 14
Unexpected MATLAB operator.

---

Merci d'avence

Bonjour,

Je pense qu'il manque un opérateur sur cette ligne:

Code:

---

1 2	K = J * ( Ex^2*.S1 + Ey^2*.S2)

---

Duf

toujours la meme probleme

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15

I=[166 154 ; 240 100; 12 144 ; 12 252; 131 144 ; 131 252; 61 62 ; 88 88; 61 88 ; 88 62];   Ex = mean(I,1); Ey = mean(I,2);   S1 = std(I,1); S2= std(I,2);   J = 1/ (Ex*Ey); K = J*( Ex^2*S1 + Ey^2*S1);

---

Citation:

---

??? Error using ==> var at 93
W must be a vector of nonnegative weights, or a scalar 0 or 1.

Error in ==> std at 32
y = sqrt(var(varargin{:}));

Error in ==> GMC at 12
S2= std(I,2);

---

Ce n'est plus la même erreur...

Que souhaites-tu calculer avec:

Code:

---

S2= std(I,2);

---

Car je suis pas sur que 2 soit une valeur valide pour le flag (voir aide)

Citation:

---

Ce n'est plus la même erreur...

Que souhaites-tu calculer avec:

Code :
S2= std(I,2);

---

je veux calculer l'ecart type selon l'axes y mais je pense que c'est faux

voila la formule

l1=sin(tetha)
l2=cos(tetha)
l3=xisin(tetha)-yicos(tetha)

avec x1 et y1 les position de pixels
std(I,1)= sigma(l1)
std(I,2)= sigma(l2)
std(I,3)= sigma(l3)

bonjour,

j'ai fait un essai pour programmer la formule en haut:

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18

clc; x = 240; y = 100; I = [x ,y ];   tetha = cart2pol( x,y);   J = [-sin(tetha); cos(tetha); x*sin(tetha) - y*cos(tetha)];   Ex = mean(I,1); Ey = mean(I,2); S1 = std(J(1)); S2 = std(J(2)); S3 = std(J(3));   K = 1 ./ (Ex .* Ey); R = J*( Ex.^2*S1 + Ey.^2*S2 + S3);

---

mais dans le resultat , il me donne des zero

Citation:

---

R =

0 0
0 0
0 0

---

avez vous une aidée ?

Citation:

---

Envoyé par hanane78

Code:

---

1 2 3 4

S1 = std(J(1)); S2 = std(J(2)); S3 = std(J(3));

---

---

Je ne suis pas un spécialiste en statistiques mais la déviation standard d'une valeur scalaire est toujours nulle.

La moyenne d'un scalaire est le scalaire lui-même, donc la variation de ce scalaire autour de sa moyenne est nulle

Citation:

---

la déviation standard d'une valeur scalaire est toujours nulle.

La moyenne d'un scalaire est le scalaire lui-même, donc la variation de ce scalaire autour de sa moyenne est nulle

---

oui ta raison parceque je travail seulement sur un seule point P(x,y)

mnt je travail sur un ensemble des point 10 point voila 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

clc; K =[]; I=[166 154;240 100; 12 144; 12 252; 131 144;131 252; 61 62; 88 88; 61 88;88 62 ]; [n l]= size(I); for i = 1:n for j = 1:l-1 tetha = cart2pol( I(i),I(j)); J = [-sin(tetha); cos(tetha); I(i)*sin(tetha) - I(j)*cos(tetha)]; K = [K,J]; end end   Ex = mean(I,1); Ey = mean(I,2);   S1 = std(K(1,:)); S2 = std(K(2,:)); S3 = std(K(3,:));   T = 1 ./ (Ex .* Ey);   R = T.*( Ex.^2*S1 + Ey.^2*S2 + S3);

---

voila l'erreur

Citation:

---

??? Error using ==> times
Matrix dimensions must agree.

Error in ==> GMC at 24
T = 1 ./ (Ex .* Ey);

---

je sais que j'ai question de taille de matrice Ex et Ey mais je dois pas changer la taille parceque se sont des moyenne selon x et Y de la matrice I

avez vous des remarques !!

Merci bcp