凸轮Matlab程序

r0=50

r1=10

e=0

q=2.617994

w=3.1415926

t=5.235988

syms d

s0=sqrt((r0).^2-e^2)

h=30*((d/2.617994)-sin(2*pi*d/2.617994)/(2*pi))

j=30

k=30*(1-(d-pi)/(2*pi/3)+sin(2*pi*(d-pi)/(2*pi/3))/(2*pi))

l=0

s=h

x=(s0+s)*sin(d)+e*cos(d)

y=(s0+s)*cos(d)-e*sin(d)

dx=diff(x,'d')

dy=diff(y,'d')

a=dx/sqrt((dx).^2+(dy).^2)

c=-dy/sqrt((dx).^2+(dy).^2)

m1=x+r1*c

n1=y+r1*a

m2=x-r1*c

n2=y-r1*a

ezplot(n1,m1,[0,q])

hold on

ezplot(y,x,[0,q])

ezplot(n2,m2,[0,q])

s=j

x=(s0+s)*sin(d)+e*cos(d)

y=(s0+s)*cos(d)-e*sin(d)

dx=diff(x,'d')

dy=diff(y,'d')

a=dx/sqrt((dx).^2+(dy).^2)

c=-dy/sqrt((dx).^2+(dy).^2)

m1=x+r1*c

n1=y+r1*a

m2=x-r1*c

n2=y-r1*a

ezplot(n1,m1,[q,w])

ezplot(y,x,[q,w])

ezplot(n2,m2,[q,w])

s=k

x=(s0+s)*sin(d)+e*cos(d)

y=(s0+s)*cos(d)-e*sin(d)

dx=diff(x,'d')

dy=diff(y,'d')

a=dx/sqrt((dx).^2+(dy).^2)

c=-dy/sqrt((dx).^2+(dy).^2)

m1=x+r1*c

n1=y+r1*a

m2=x-r1*c

n2=y-r1*a

ezplot(n1,m1,[w,t])

ezplot(y,x,[w,t])

ezplot(n2,m2,[w,t])

s=l

x=(s0+s)*sin(d)+e*cos(d)

y=(s0+s)*cos(d)-e*sin(d)

dx=diff(x,'d')

dy=diff(y,'d')

a=dx/sqrt((dx).^2+(dy).^2)

c=-dy/sqrt((dx).^2+(dy).^2)

m1=x+r1*c

n1=y+r1*a

m2=x-r1*c

n2=y-r1*a

ezplot(n1,m1,[t,2*pi])

ezplot(y,x,[t,2*pi])

ezplot(n2,m2,[t,2*pi])

x=r0*cos(d)

y=r0*sin(d)

ezplot(y,x,[0,2*pi])

grid on

gtext('外包络线')

gtext('理论廓线')

gtext('实际廓线')

gtext('基圆²')