figure
set(gcf,'defaultaxesnextplot','add'...
,'defaultaxesdataaspectratio',[1 1 1]...
,'defaultaxesDataAspectRatioMode','manual'...
,'defaultaxesPlotBoxAspectRatioMode','manual'...
,'defaultaxesXtickmode','manual'...
,'defaultaxesytickmode','manual'...
,'defaultaxesbox','on')
% set(gcf,'defaultaxesnextplot','add')
% set(gcf,'defaultaxesdataaspectratio',[1 1 1])
% set(gcf,'defaultaxesDataAspectRatioMode','manual')
% set(gcf,'defaultaxesPlotBoxAspectRatioMode','manual')
% set(gcf,'defaultaxesXtickmode','manual')
% set(gcf,'defaultaxesytickmode','manual')
% set(gcf,'defaultaxesbox','on')
get(gcf,'defaultaxesPlotBoxAspectRatio');
w=linspace(0,2*pi,100);
c=1;
a=5;
b=2;
N=200;
omega=rand(1,N);
wr=rand(1,N);
wr2=rand(1,N);
r=sqrt(rand(1,N));
r2=sqrt(rand(1,N));
%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%
subplot(231)
plot(c*cos(w),c*sin(w))
for k=1:4
xx=c*wr(k)*cos(2*pi*omega);
yy=c*wr(k)*sin(2*pi*omega);
plot(xx,yy,'r.')
end
%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%
subplot(232)
plot(c*cos(w),c*sin(w))
xx=c*wr.*cos(2*pi*omega);
yy=c*wr.*sin(2*pi*omega);
plot(xx,yy,'r.')
%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%
subplot(234)
plot(c*cos(w),c*sin(w))
xx=c*r.*cos(2*pi*omega);
yy=c*r2.*sin(2*pi*omega);
plot(xx,yy,'r.')
%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%
subplot(235)
plot(c*cos(w),c*sin(w))
xx=c*r.*cos(2*pi*omega);
yy=c*r.*sin(2*pi*omega);
plot(xx,yy,'r.')
%%%%%%%%%%%%%%%%%%%%%%%%%%%
%&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
%%%%%%%%%%%%%%%%%%%%%%%%
subplot(233)
plot(a*cos(w),b*sin(w))
xx=a*wr.*sin(2*pi*omega);
yy=b*wr2.*cos(2*pi*omega);
plot(xx,yy,'r.')
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%55
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
subplot(236)
plot(a*cos(w),b*sin(w))
xx=a*r.*cos(2*pi*omega);
yy=b*r.*sin(2*pi*omega);
plot(xx,yy,'r.')