这段程序不是画圆的,是求点到圆心的距离;
p
中保存的是点的组数(多个点),然后批量求这些点和圆心的距离。
d=sqrt((p(:,1)-xc).^2+(p(:,2)-yc).^2)-r;
上面这句就是求距离d
d=[d,d];
这句没有看出有什么意义,可能和函数在具体应用有关吧!
下面是一个运行的例子:
clc
clear all
close all
x=3;y=4;R=5;
m=x-R:0.1:x+R;
n=[]
p=length(m);
for i=1:p
n(i)=sqrt(R^2-(m(i)-x)^2)+y;
n2(i)=-sqrt(R^2-(m(i)-x)^2)+y;
end
plot(m,n,'r',m,n2,'r')