【无标题】MATLAB 关于圆的绘制和平面的相切

最新推荐文章于 2023-03-01 22:32:41 发布

302号民工

最新推荐文章于 2023-03-01 22:32:41 发布

阅读量498

点赞数 1

文章标签：算法

本文链接：https://blog.csdn.net/weixin_45388992/article/details/123749477

版权

本文详细介绍了如何使用MATLAB进行圆的绘制，并探讨了如何模拟平面与圆相切的情况，提供了相关的算法实现。

摘要由CSDN通过智能技术生成

% n=[0.00660129 -0.00955624 -0.999933]; %法向量n
% r=194.159; %圆的半径为1
% c=[812.699 -386.513 1.68927]; %圆心的坐标
% theta=(0:2*pi/100:2*pi)'; %theta角从0到2*pi
% a=cross(n,[1 0 0]); %n与i叉乘，求取a向量
% if ~any(a) %如果a为零向量，将n与j叉乘
%     a=cross(n,[0 1 0]);
% end
% b=cross(n,a); %求取b向量
% a=a/norm(a); %单位化a向量
% b=b/norm(b); %单位化b向量
% 
% c1=c(1)*ones(size(theta,1),1);
% c2=c(2)*ones(size(theta,1),1);
% c3=c(3)*ones(size(theta,1),1);
% 
% x=c1+r*a(1)*cos(theta)+r*b(1)*sin(theta);%圆上各点的x坐标
% y=c2+r*a(2)*cos(theta)+r*b(2)*sin(theta);%圆上各点的y坐标
% z=c3+r*a(3)*cos(theta)+r*b(3)*sin(theta);%圆上各点的z坐标

plot3(x,y,z)
hold on;
xlabel('x轴')
ylabel('y轴')
zlabel('z轴')


x4=600:0.1:1200;
y4=-600:0.1:0;

z4=((y4+157245)/-0.999952)*(-0.00258615)-405.677;
plot3(x4,y4,z4,'r');
hold on;


x5=600:0.1:1200;
y5=-600:0.1:0;
z5=((y5+126172)/-0.999959)*(-0.002