用MATLAB求三圆交点（条件如下）新手没分，跪求好心人

2025-03-28 07:37:20

推荐回答（3个）

回答1：

clear;
%------------圆相交----------%
%-----------k打酱油k---------%
x1=input('圆心1坐标x1=?');
y1=input('圆心1坐标y1=?');
r1=input('圆1半径r1=?');
x2=input('圆心2坐标x2=?');
y2=input('圆心2坐标y2=?');
r2=input('圆2半径r2=?');
x3=input('圆心3坐标x3=?');
y3=input('圆心3坐标y3=?');
r3=input('圆3半径r3=?');
syms X1 Y1 R1 X2 Y2 R2 X3 Y3 R3;
eq1=sym('(x-X1)^2+(y-Y1)^2=R1^2');
eq2=sym('(x-X2)^2+(y-Y2)^2=R2^2');
eq3=sym('(x-X3)^2+(y-Y3)^2=R3^2');
eq1=subs(eq1,{X1 Y1 R1},{x1 y1 r1});
eq2=subs(eq2,{X2 Y2 R2},{x2 y2 r2});
eq3=subs(eq3,{X3 Y3 R3},{x3 y3 r3});
[A1 B1]=solve(eq1,eq2);
[A2 B2]=solve(eq1,eq3);
[A3 B3]=solve(eq2,eq3);
X1=eval(A1);
Y1=eval(B1);
X2=eval(A2);
Y2=eval(B2);
X3=eval(A3);
Y3=eval(B3);
n1=length(A1)*(imag(X1(1))==0);
switch n1
case 0
fprintf('\n')
disp('圆1和圆2没有交点')
otherwise
intersection=['圆1和圆2的交点有', num2str(n1),'个'];
fprintf('\n')
disp(intersection)
for i=1:n1
showX=['X' num2str(i) '=' num2str(X1(i))];
showY=['Y' num2str(i) '=' num2str(Y1(i))];
disp([showX ' ' showY])
end
end
n2=length(A2)*(imag(X2(1))==0);
switch n2
case 0
fprintf('\n')
disp('圆1和圆3没有交点')
otherwise
intersection=['圆1和圆3的交点有', num2str(n2),'个'];
fprintf('\n')
disp(intersection)
for i=1:n2
showX=['X' num2str(i) '=' num2str(X2(i))];
showY=['Y' num2str(i) '=' num2str(Y2(i))];
disp([showX ' ' showY])
end
end
n3=length(A3)*(imag(X3(1))==0);
switch n3
case 0
fprintf('\n')
disp('圆2和圆3没有交点')
otherwise %有实根的情况
intersection=['圆2和圆3的交点有', num2str(n3),'个'];
fprintf('\n')
disp(intersection)
for i=1:n3
showX=['X' num2str(i) '=' num2str(X3(i))];
showY=['Y' num2str(i) '=' num2str(Y3(i))];
disp([showX ' ' showY])
end
end

%本来还想设计个GUI界面给你，参数可以直接在对话框输入，并且还可以作图画出3个圆，这样会比较直观一点。不过实在没时间了，明天还得上班。你有兴趣的话自己去做吧。

回答2：

我ca 楼上牛人啊

回答3：

5分还是太少了点

最新问答

三星960EVO M.2固态装win7系统开机启动就蓝屏怎么回事？求解

嘉合小镇周边环境怎么样？生活便利吗？

猫和老鼠一共有多少集

美式箱变里的ZGS11变压器，和普通的S11油变是一样的吗？不一样，有什么差别？

拼多多几百块钱一部6g运存加128g内存是真的吗

华为手机屏幕突然变黄,连拍照镜头看到的也是黄色,怎么办?

火影忍者oljjc六道匹配哪些段位

赣州市大余县池江镇到村下古上下有多远

有谷的猫粮对猫咪肠胃不友好？是这样吗？

能帮忙查一下从日本发的邮件CD098995848JP到哪里了吗