导杆的角位移曲线和刨刀的位移曲线导杆的角位移曲线和刨刀的位移曲线导杆的角位移曲线和刨刀的位移曲线导杆的角位移曲线和刨刀的位移曲线
x1=linspace(0,2*pi,100); l1=100; l3=360; l4=200; l6=200; y=336; w1=2*pi; x3=atan((l6+l1*sin(x1))./(l1*cos(x1))); for i = 1 : 100; if x3(i)<0 x3(i)=x3(i)+pi; end end S3=l1*cos(x1)./cos(x3); V23=(-1)*w1*l1*sin(x1-x3); w3=w1*l1*cos(x1-x3)./S3; a23=w3.^2.*S3-w1.^2*l1*cos(x1-x3); A3=(w1.^2.*l1.*sin(x3-x1)-2.*w3.*V23)./S3; x4= atan((y-l3*sin(x3))./l4); Se=l3*cos(x3)+l4*cos(x4); w4=(-1)*w3.*l3.*cos(x3)./(l4.*cos(x4)); Ve=(-1)*w3.*l3.*sin(x3-x4)./cos(x4); A4=(x3.^2.*l3.*sin(w3)+x4.^2.*l4.*sin(x4)-A3.*l3.*cos(x3))./(l4.*cos(x4)); ae=(-1)*(A3.*l3.*sin(x3-x4)+w3.^2.*l3.*cos(x3-x4)-w4.^2.*l4)./cos(x4); plotyy(180*x1./pi,180*x3./pi,180*x1./pi,Se/1000); xlabel('x1/(°)'),ylabel('x3(°)'); title('导杆的角位移曲线和刨刀的位移曲线'); text(150,108,'x3'); text(200,60,'Se'); grid on; plotedit on; 导杆的角速度曲线和刨刀的速度曲线导杆的角速度曲线和刨刀的速度曲线导杆的角速度曲线和刨刀的速度曲线导杆的角速度曲线和刨刀的速度曲线
x1=linspace(0,2*pi,100); l1=100; l3=360; l4=200; l6=200; y=336; w1=2*pi; x3=atan((l6+l1*sin(x1))./(l1*cos(x1))); for i = 1 : 100; if x3(i)<0 x3(i)=x3(i)+pi; end end S3=l1*cos(x1)./cos(x3); V23=(-1)*w1*l1*sin(x1-x3); w3=w1*l1*cos(x1-x3)./S3; a23=w3.^2.*S3-w1.^2*l1*cos(x1-x3); A3=(w1.^2.*l1.*sin(x3-x1)-2.*w3.*V23)./S3; x4= atan((y-l3*sin(x3))./l4); Se=l3*cos(x3)+l4*cos(x4); w4=(-1)*w3.*l3.*cos(x3)./(l4.*cos(x4)); Ve=(-1)*w3.*l3.*sin(x3-x4)./cos(x4); A4=(x3.^2.*l3.*sin(w3)+x4.^2.*l4.*sin(x4)-A3.*l3.*cos(x3))./(l4.*cos(x4)); ae=(-1)*(A3.*l3.*sin(x3-x4)+w3.^2.*l3.*cos(x3-x4)-w4.^2.*l4)./cos(x4); plotyy(180*x1./pi,w3,180*x1./pi,Ve/1000); xlabel('x1/(°)'),ylabel('w3(rad/s)'); title('导杆的角速度曲线和刨刀的速度曲线'); grid on; text(100,-5,'Ve'); text(100,1.8,'w3'); axis([0 400 -7 4]); plotedit on; 导杆的角加速度曲线和刨刀的加速度曲线导杆的角加速度曲线和刨刀的加速度曲线导杆的角加速度曲线和刨刀的加速度曲线导杆的角加速度曲线和刨刀的加速度曲线
x1=linspace(0,2*pi,100); l1=100; l3=360; l4=200; l6=200; y=336; w1=2*pi; x3=atan((l6+l1*sin(x1))./(l1*cos(x1))); for i = 1 : 100; if x3(i)<0 x3(i)=x3(i)+pi; end end S3=l1*cos(x1)./cos(x3); V23=(-1)*w1*l1*sin(x1-x3); w3=w1*l1*cos(x1-x3)./S3; a23=w3.^2.*S3-w1.^2*l1*cos(x1-x3); A3=(w1.^2.*l1.*sin(x3-x1)-2.*w3.*V23)./S3; x4= atan((y-l3*sin(x3))./l4); Se=l3*cos(x3)+l4*cos(x4); w4=(-1)*w3.*l3.*cos(x3)./(l4.*cos(x4)); Ve=(-1)*w3.*l3.*sin(x3-x4)./cos(x4); A4=(x3.^2.*l3.*sin(w3)+x4.^2.*l4.*sin(x4)-A3.*l3.*cos(x3))./(l4.*cos(x4)); ae=(-1)*(A3.*l3.*sin(x3-x4)+w3.^2.*l3.*cos(x3-x4)-w4.^2.*l4)./cos(x4); plotyy(180*x1./pi,A3,180*x1./pi,ae/1000); xlabel('x1/(°)'),ylabel('A3(rad/s^2)'),grid; title