//这是个画圆的程序
m=300;n=300;x=50;y=50;r=10 //m,n是图像大小,x,y是圆心坐标,r是半径
xy=[x y]; //xy为圆内点的坐标组成的矩阵,矩阵每行都是一个圆内的点,首先将圆心点加入
for i=1:m
for j=1:n
if (i-x)^2+(j-y)^2