平面六连杆机构的运动分析Matlab代码clc,clear%参数赋值l1=40;l2=55;l3=55;l4=22;M=-l;%装配模式omiga1=10;thetal=0:l:360;thetal=thetal*pi/l80;A=2*ll*l2*sin(thetal);B=2*l2*(ll*cos(thetal)-l4);C=11人2+12人2+14人2-13人2-2*l1*14*cos(theta1);E=2*ll*l3*sin(thetal);F=2*13*(11*cos(theta1)-14);G=12A2-11A2-13A2-14A2+2*11*14*cos(theta1);theta3=2*atan(E+M*sqrt(E.A2+F.A2-GA2)./(F-G);theta2=2*atan(A+M*sqrt(A.A2+B.A2-C.A2)./(B-C);omiga2=omiga1*1*sin(theta1-theta3)./(12*sin(theta3-theta2);omiga3=