for k=1:length(B)
BM=cell2mat(B(k));x=BM(1:length(BM),2);y=BM(1:length(BM),1);
plot(x,y,'r:.'), hold on
axis equal
abc=[x y ones(size(x))]\[-(x.^2+y.^2)];
a=abc(1);b=abc(2);c=abc(3);
xc=-.5*a;yc=-.5*b;R=sqrt((a^2+b^2)-4*c);
[xc yc R]
end
哪位大佬给我标注一下这段程序 |