|
回复 楼主 的帖子
无聊建个函数如下:
function chilunplot(r,N)
a=2/5*2*pi*r/N;b=3/5*2*pi*r/N;
R=((a/2)^2+(r+b)^2)^(1/2);
A=atan(a/2/(r+b));
for n=1:N
p=2*pi/N*(n-1);q=2*pi/N*n;
theta=2*pi/N*(n-1)+7/10*2*pi/N;
tao=2*pi/N*(n-1)+2/5*2*pi/N;
afa=theta-A;beta=theta+A;
ma1=linspace(p,tao,20);x1=r*cos(ma1);y1=r*sin(ma1);plot(x1,y1);hold on;
plot([r*cos(tao),R*cos(afa)],[r*sin(tao),R*sin(afa)]);
plot([r*cos(q),R*cos(beta)],[r*sin(q),R*sin(beta)]);
plot([R*cos(afa),R*cos(beta)],[R*sin(afa),R*sin(beta)]);
end
axis equal;
axis([-2*r 2*r -2*r 2*r]);
调用格式:chilunplot(r,N),r代表齿轮半径,N代表一共有几个齿;
例如:chilunplot(10,15) |
|