已知三维坐标系oxyz的坐标平面oxz中存在三条直线三条直线依次为：l1z=xl2z=x+5l3x=5请绘制三维坐标系的z轴和三条空间直线以及平面oxz的第一卦限中的直线l1点围绕z轴的旋转曲面直线l3以直线l2为轴的旋转曲面且两个旋转曲面的母线长度均为2sqrt2。给出Matlab求解代码与结果

代码如下：

% 绘制三维坐标系和三条直线 figure; axis([-10 10 -10 10 -10 10]); line([0 0],[0 0],[-10 10],'LineWidth',2,'Color','k'); % z轴 line([-10 10],[0 0],[0 0],'LineWidth',2,'Color','k'); line([0 0],[-10 10],[0 0],'LineWidth',2,'Color','k'); line([0 5],[-5 0],[0 5],'LineWidth',2,'Color','r'); % l1 line([0 5],[-5 0],[5 10],'LineWidth',2,'Color','g'); % l2 line([5 5],[-5 5],[0 0],'LineWidth',2,'Color','b'); % l3 xlabel('x'); ylabel('y'); zlabel('z');

% 绘制l1点围绕z轴的旋转曲面 theta = 0:0.1:2pi; r = 2sqrt(2); x = rcos(theta); y = rsin(theta); z = x; hold on; for i = 1:length(x) p = [x(i) y(i) z(i)]; R = [cos(theta(i)) -sin(theta(i)) 0; sin(theta(i)) cos(theta(i)) 0; 0 0 1]; p_new = pR; plot3(p_new(1),p_new(2),p_new(3),'','Color','m'); end

% 绘制l3以l2为轴的旋转曲面 theta = 0:0.1:2pi; r = 2sqrt(2); x = rcos(theta); y = rsin(theta); z = -5 + x; for i = 1:length(x) p = [x(i) y(i) z(i)]; axis_vec = [0 0 1]; R = vrrotvec2mat(vrrotvec(axis_vec,[5 0 -5])); p_new = pR; plot3(p_new(1),p_new(2),p_new(3),'','Color','c'); end

结果见附图：

其中，紫色的点表示l1点围绕z轴的旋转曲面，青色的点表示l3以l2为轴的旋转曲面