d=1 %螺距
r=0.05 %螺线粗半径
th=linspace(0,10*pi,361)
ph=linspace(0,2*pi,37)
[t p]=meshgrid(th,ph)
x=(R+r*cos(p)).*cos(t)
y=(R+r*cos(p)).*sin(t)
z=d*t/(2*pi)+r*sin(p)
surf(x,y,z,'edgecolor','none','facecolor','b','facealpha',0.5)hold on
surf(x,-y,z,'edgecolor','none','facecolor','g','facealpha',0.5)hold off
axis equal
xlabel('X')
ylabel('Y')
zlabel('乱神Z')
这个程序可以绘制二维或三维的,你试试,俺枝做瞎琢磨的,好像还没什么问题。
clear all
n=10%圈数
v=0.1%时间变量t的间隔
t=0:v:n*2*pi
a=1%每一圈平面的间亏弯隔是a2pi左右
k=1%z最终上升多少个t,可计算销搭闷出每一圈z上升k2pi左右
m0=0.51*pi%初始相角矫正控制,控制结尾在哪个位置
x=a*(t+m0).*sin(t+m0)
y=a*(t+m0).*cos(t+m0)
z=k*t
figure(1)
plot3(x,y,z)grid on
xlabel('X'),ylabel('Y'),zlabel('Z')
axis equal
figure(2)
plot(x,y)grid on
xlabel('X'),ylabel('Y')
axis equal
假定简空C=1画r=1,2,3,4,5,6,7,8,9,10十条灶咐前线的程隐清序: clear clc close all C=1for r=1:10a=-pi/2:pi/100:pi/2k=tan(a)p=C*exp(r*k)plot(a,p) hold on end axis([-2 2 -10 20])欢迎分享,转载请注明来源:内存溢出
评论列表(0条)