matlab齿轮转动代码,matlab编程齿轮 - 程序语言 - MATLAB/Mathematica - 小木虫论坛-学术科研互动平台...

function chilun(m,z,ha,c)

r=m*z/2;

rb=r*cos(pi/9);

ra=r+ha*m;

rf=r-(ha+c)*m;

sita=2*pi/z;

almax=acos(rb/ra);

alg=acos(rb/r);

ct1=tan(alg)-alg;

ct2=tan(almax)-almax;

bita=sita/4+ct1-ct2;

for i=0:1:z

fai=-bita:0.001:bita;

xa=ra*cos(fai+i*sita);

ya=ra*sin(fai+i*sita);

hold on;

plot(xa,ya,'b');

fail=bita+ct2:0.01:sita-bita-ct2;

xb=rb*cos(fail+i*sita);

yb=rb*sin(fail+i*sita);

hold on;

plot(xb,yb,'r');

for al=almax:-0.001:0;

ct=tan(al)-al;

ri=rb/cos(al);

alfa=-ct+(i+0.25)*sita+ct1;

x2=ri*cos(alfa);

y2=ri*sin(alfa);

plot(x2,y2,'b');

end

for al2=0:0.001:almax

ct=tan(al2)-al2;

ri=rb/cos(al2);

beta=ct+(i-0.25)*sita-ct1;

x3=ri*cos(beta);

y3=ri*sin(beta);

plot(x3,y3,'g');

end

end