polspc.m

来自「一种新的时频分析方法的matlab源程序。」· M 代码 · 共 36 行

function q=polspc(a,x,u,cn)

% The function POLSPC draws a contour plot of a 2-D data 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: POLARS
%
% Calling sequence-
% q=polspc(a,x,u,cn)
%
% Input-
%	a	- the radius
%	x	- the x-y axis column vector
%	u	- the matrix (issued from ro2xy.cpp)
%	cn	- number of contours
% Output-
%	q	- the handle of the figure
%
% See also-     
% POLARS, RO2XY.cpp.

d=linspace(0,2*pi,200)';
r=a*ones(200,1);
polars(d,r,'w');

h=get(gca,'children');
set(h(1),'linewidth',1.);
n=length(h);

hold on
u=flipud(rot90(u));
q=contour(x,x,u,cn);
hold off