polvec.m

一种新的时频分析方法的matlab源程序。

function q=polvec(a,x,ux,uy)

% The function POLVEC draws a vector plot of two 2-D data arrays in polar coordinates
% with 0 degrees located at the 12 o'clock position and 
% degrees increasing counter-clockwise.
% To add labels to the contours do: 
%	hold on; eclabel(q,'fontsize',6); .
% Non MATLAB Library routine called is: POLARH.
%
% Calling sequence-
% q=polvec(a,x,ux,uy)
%
% Input-
%	a	- the radius
%	x	- the x-y axis column vector
%	ux	- the matrix (issued from ro2xy.cpp)
%	uy	- the matrix (issued from ro2xy.cpp)
% Output-
%	q	- the handle of the figure
%
% See also-     
% POLARS, RO2XY.cpp.

close
r=a*ones(200,1);
d=linspace(0,2*pi,200)';
polarh(d,r);
%polar(d,r);
hold on
u=rot90(ux);
v=rot90(uy);
%quiver(x,x,flipud(ux),-flipud(uy));
quiver(x,x,flipud(v),flipud(u));
hold off