cpolspc.m

这是一个关于hht变换很有用的工具箱

function q=cpolspc(a,x,u)

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

close
polar2(a);
hold on
ss=max(max(u));
u=flipud(rot90(u))-2*ss;
%q=surf(x,x,u);
q=pcolor(x,x,u);
shading interp;
colormap(hsv);
hold off