标签：倾斜 text sqrt MATLAB rho plot3 theta circle 绘制

绘制倾斜圆

绘制倾斜圆，设半径为\(\rho\)；在x，y，z轴上的投影分别为x，y，z，则\(\rho\)² = x²+y²+z²。

圆平面绕着y轴倾斜\(\theta\)=45\(^{\circ}\)，则tan\(\theta\)=\(\frac{z}{y}\)，x²+y²+{ytan\(\theta\)}²=\(\rho\)²，求解得到\(y = \sqrt{ \frac{ \rho^2 - x^2}{ 1 + tan^2\theta} }\)

假设\(\rho\)=5，\(\theta\)=45\(^{\circ}\)，则令x=linspace(-\(\rho\),0.1,\(\rho\))，可求得y，进而求得z。

附录

附录1:demo_for_circle_theta

clc,clear all
% 三维图窗
figure
r=5;
scatter3(0,0,0,'r','filled');
xlabel('x');
ylabel('y');
zlabel('z');
% 绘制xy平面圆和绕x轴倾斜45度圆
hold on
circle_theta(r,pi/4);
hold on
circle_theta(r,0);
hold on
% 绘制重要的点
scatter3(3,0,0,'r','filled');
scatter3(3,4,0,'r','filled');
scatter3(3,2*sqrt(2),2*sqrt(2),'r','filled');
scatter3(3,2*sqrt(2),0,'r','filled');
% 绘制线段
plot3([-r,r],[0,0],[0,0]);
plot3([3,3],[0,4],[0,0]);
plot3([3,3],[0,2*sqrt(2)],[0,2*sqrt(2)]);
plot3([0,3],[0,2*sqrt(2)],[0,2*sqrt(2)]);
plot3([3,3],[2*sqrt(2),2*sqrt(2)],[0,2*sqrt(2)]);
% 标注
text(3,1,0.5,'\theta=45^\circ');
text(2,2,2,'\rho');
text(1.5,0,0,'x');
text(3,sqrt(2),0,'y');
text(3,2*sqrt(2),sqrt(2),'z');

附录2：circle_theta

function [] = circle_theta(radius,theta)
% 绘制绕x轴旋转theta角度得到的圆
% 不能等于PI/2
x=linspace(-radius,radius,50);
y=sqrt((radius^2-x.^2)/(1+tan(theta)^2));
z=y.*tan(theta);
plot3(x,y,z,'k');
hold on
plot3(-x,-y,-z,'k');
hold off
end

标签：倾斜,text,sqrt,MATLAB,rho,plot3,theta,circle,绘制
来源： https://www.cnblogs.com/xingdu/p/15584036.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。