求曲柄滑块机构的matlab程序，已知杆长以及滑块行程，要优化一下各个参数，好的话会另外加分

%Matlab仿真程序肆纯如下：

r1=36 %单位mm

r2=140

omiga1=60%单位d/sec

x11=1:720

for i=1:720

x1(i)=i*pi/180

%sin(x2(i))=-r1/r2*sin(x1(i))

x2(i)=asin(-r1/r2*sin(x1(i)))

x22(i)=x2(i)*180/pi

r3(i)=r1*cos(x1(i))+r2*cos(x2(i))

B=[-r1*omiga1*sin(x1(i))

r1*omiga1*cos(x1(i))]

A=[r2*sin(x2(i)) 1-r2*cos(x2(i)) 0] X=inv(A)*B

omiga2(i)=X(1,1)

v3(i)=X(2,1)

end

plot(x11/60,0.5*r1*sin(x1))

xlabel('时间 t/sec')

ylabel('谈中连杆质裂侍咐心在Y轴上位置/mm')

figure(2)

plot(x11/60,r3)

xlabel('时间 t/sec')

ylabel('滑块位移r3/mm')

figure(3)

plot(x11/60,omiga2)

xlabel('时间 t/sec')

ylabel('连杆角速度omiga2/rad/sec')

figure(4)

plot(x11/60,v3*pi/180)

xlabel('时间 t/sec')

ylabel('滑块速度v3/mm/sec')

用matlab进行机构分析的步骤，桥姿可以根据机械原理的方法来进行敏竖绝。

1、分析机构的综合，其机构由曲柄导杆机构和偏心曲柄滑块机构组成。

2、建立其平面运动方程，求其滑块C在n1转速下的位移、速度、加速度。

3、建立其平面力平衡方程，求其滑块C的摩擦阻力。

如需要，可能要花费比较长时间去做。纤大

---