polar1.m

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

function h1=polar1(a)
%
% polar1(a):
%
% Function to create polar coordinates in a new figure with
% 0 degrees located at the 3 o'clock position and 
% degrees increasing counter-clockwise.
%
% Input-
%	a	- the radius
% Output-
%	h1	- the handle of the figure
%
% Notes-
%	(1) do first:  h=econtour(x,x,matrix);
% 	you can then find a propriate radius.
%	(2) do:        polar1(r);
%	(3) do again:  h=econtour(x,x,matrix);
%	(4) do:        eclabel(h,'fontsize',6); hold off;     
%	This program is related to ro2xy.cpp, econtour.m, eclabel.m.  

d=linspace(0,2*pi,200)';
r=a*ones(200,1);
polar(d,r,'w');    
h=get(gca,'children');
set(h(1),'linewidth',1.);
n=length(h);
for i=14:19,
   set(h(i),'linestyle',':','linewidth',.25);
end
for i=21:2:n,
   set(h(i),'linestyle',':','linewidth',.25);
end
hold on