pmean.m v1.1 (Sep 2009)

MEAN and STD of angles, polar, cylindrical, spherical or cartesian coordinates!

This function calculates the mean and standard deviation of vectors in polar coordinates (and others), that is, from angles and magnitudes.

The usage is very simple:

For polar or cylindrical: Direction, Radius, Height
[mD,eD] = pmean(D);
[mD,mR,eD,eR] = pmean(D,R);
[mD,mR,mH,eD,eR,eH] = pmean(D,R,H);

For spherical coordinates: Azimuth, Elevation, Radius
[mA,eA] = pmean(A);
[mA,mE,eA,eE] = pmean(A,E,'s');
[mA,mE,mR,eA,eE,eR] = pmean(A,E,R,'s');

For cartesian coordinates:
[mX,eX] = pmean(X,'c');
[mX,mY,eX,eY] = pmean(X,Y,'c');
[mX,mY,mZ,eX,eY,eZ] = pmean(X,Y,Z,'c');

where "m" stands for MEAN's, and "e" for STD's.

Input angles should be in radians.

It ignores NaN's.

Enjoy it!
Any comments or bug reports will be very appreciated!