unction plot_polygon(n,a)
%n为多边形的边数
%a为多边形的边长,默认值为1
if nargin<2
a=100;n=4;
end
t=(0:n)*2*pi/n;
theta=(n-2)*pi/2/n;
r=a/2/cos(theta);
x=r*cos(t);
y=r*sin(t);
%plot(x,y,'r-');
axis equal;
axis([0 100 0 100]);
title('半径为10的圆所截节点图')
hold on
for m=-9:9
for k=-5:5
t=0:pi/200:2*pi;
plot((10*cos(t)+50)+10*sqrt(3)*k-5*sqrt(3)*m,(10*sin(t)+50)+15*m,'r')
end
end
|