uv2polar.m

一个研究声多普勒计程仪很好的工具箱

function [dir,spd] = uv2polar(u,v)

%  uv2polar:  direction and speed of (east, north) currents
%
%       u is eastward current component
%       v is northward current component
%       dir is compass direction of current
%             (90 for eastward flow, 270 for westward)
%       spd is magnitude of current
%
%       Fran Hotchkiss, February 20, 1998.

[dir,spd] = cart2pol(v,u);
i = find(dir < 0);
dir(i) =dir(i) + 2*pi;
dir = dir .* 180 ./ pi;