How can I plot a circle use given point?

4 visualizzazioni (ultimi 30 giorni)
Seong-Won Ham
Seong-Won Ham il 25 Gen 2021
Commentato: weikang zhao il 25 Gen 2021
I want to plot a circle, use given data.
but, I just know the some of the point in circle, but do not midpoint and radius.
Then, How Can I plot the circle using a given point?

Risposte (2)

KSSV
KSSV il 25 Gen 2021
If you have points.....you can find the radius of circle. It depends on how many points you know.
Frame the distance formula and solve the equations.

weikang zhao
weikang zhao il 25 Gen 2021
If you have coordinates of three points, the radius and center of the circle can be obtained by solving a positive definite equation set. If you have more points, there will be a overdetermined set, the midpoint and radius have a least squares estimation.
  2 Commenti
Seong-Won Ham
Seong-Won Ham il 25 Gen 2021
Thanks for your ccomments.
Now, I have a lot of point data and will try least squares estimation to find the equation.
I wonder that the best code of calculate radius.
If you have an idea of code, Please comment for me.
I'm sorry that my poor english.
thanks.
weikang zhao
weikang zhao il 25 Gen 2021
give you a demo:
clear
c=[5;5];%real coordinate of midpoint
r=3;%real radius
pointNum=20;%num of all the points
thetatrain=unifrnd(0,2*pi,1,pointNum);
rtrain=r+0.01*randn(1,pointNum);
xytrain=c+[rtrain.*cos(thetatrain);rtrain.*sin(thetatrain)];%Generate a set of points perturbed around the circle
plot(xytrain(1,:),xytrain(2,:),'.');
A=[xytrain.',ones(pointNum,1)];
b=sum(xytrain.^2).';
est=(A.'*A)\A.'*b;
rhat=sqrt(est(1)^2/4+est(2)^2/4+est(3));%LS estimate of radius
xhat=est(1)/2;%LS estimate of the x coordinate of midpoint
yhat=est(2)/2;%LS estimate of the y coordinate of midpoint
hold on;
plot(xhat+rhat*cos(0:pi/180:2*pi),yhat+rhat*sin(0:pi/180:2*pi));
If you are interested in mathematical principles, email me. my email address

Accedi per commentare.

Categorie

Scopri di più su 2-D and 3-D Plots in Help Center e File Exchange

Tag

Prodotti

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by