机械原理大作业牛头刨床-matlab编程

1、子程序：function F=myfun(x,theta1,l1,l3,l4,h,h1,h2)F=l4*cosd(x(2)+x(1)*cosd(x(3)-l1*cosd(theta1)-h2; l4*sind(x(2)+x(1)*sind(x(3)-l1*sind(theta1)-h1; l4*cosd(x(2)+l3*cosd(x(3)-x(4); l4*sind(x(2)+l3*sind(x(3)-h;主程序：theta1=0:360;s3=0:360;theta3=0:360;theta4=0:360;sc=0:360;v3=0:360;omega3=0:360;omega4=0:360

2、;vc=0:360;a3=0:360;alpha3=0:360;alpha4=0:360;ac=0:360;m=0.5;0;70;0.47;l1=160;l3=960;l4=input('Please input Lde ');h=900;h1=460;h2=120; for i=1:1:361 x = fsolve(x) myfun(x,theta1(i),l1,l3,l4,h,h1,h2),m(1);m(2);m(3);m(4); m=x; s3(i)=x(1); theta4(i)=x(2); theta3(i)=x(3); sc(i)=x(4);endfor i=1:1

3、:361V1=cosd(theta3(i) -l4*sind(theta4(i) -s3(i)*sind(theta3(i) 0; sind(theta3(i) l4*cosd(theta4(i) s3(i)*cosd(theta3(i) 0; 0 -l4*sind(theta4(i) -l3*sind(theta3(i) -1; 0 l4*cosd(theta4(i) l3*cosd(theta3(i) 0;V2=-l1*sind(theta1(i);l1*cosd(theta1(i);0;0;V3=V1V2;v3(i)=V3(1);omega4(i)=V3(2);omega3(i)=V3(

4、3);vc(i)=V3(4); A111= omega3(i)*sind(theta3(i) l4*omega4(i)*cosd(theta4(i) v3(i)*sind(theta3(i)+s3(i)*omega3(i)*cosd(theta3(i) 0; -omega3(i)*cosd(theta3(i) l4*omega4(i)*sind(theta4(i) -v3(i)*cosd(theta3(i)+s3(i)*omega3(i)*sind(theta3(i) 0; 0 l4*omega4(i)*cosd(theta4(i) l3*omega3(i)*cosd(theta3(i) 0;

5、 0 l4*omega4(i)*sind(theta4(i) l3*omega3(i)*sind(theta3(i) 0;A112=v3(i);omega4(i);omega3(i);vc(i);A11=A111*A112;A12=-l1*cosd(theta1(i); -l1*sind(theta1(i); 0; 0;A1=A11+A12;A2=cosd(theta3(i) -l4*sind(theta4(i) -s3(i)*sind(theta3(i) 0; sind(theta3(i) l4*cosd(theta4(i) s3(i)*cosd(theta3(i) 0; 0 -l4*sin

6、d(theta4(i) -l3*sind(theta3(i) -1; 0 l4*cosd(theta4(i) l3*cosd(theta3(i) 0;A3=A2A1;a3(i)=A3(1);alpha4(i)=A3(2);alpha3(i)=A3(3);ac(i)=A3(4);endfigure(1);plot(theta1,theta3,'k',theta1,theta4,'k:');grid on;xlabel('theta_1( o )');ylabel('theta_3 theta_4 ( o )');legend(

7、9;theta_3','theta_4');figure(2)plot(theta1,omega3,'k',theta1,omega4,'k:');grid on;xlabel('theta_1( o )');ylabel('omega_3 omega_4 (r a d / s)');legend('omega_3','omega_4'); figure(3)plot(theta1,alpha3,'k',theta1,alpha4,'k:');

8、grid on;xlabel('theta_1( o )');ylabel('alpha_3 alpha_4 (r a d / s2)');legend('alpha_3','alpha_4');figure(4) plot(theta1,sc,'k');grid on;xlabel('theta_1( o )');ylabel('s_c (mm)');figure(5) plot(theta1,vc,'k');grid on;xlabel('theta_1( o )');ylabel(