matlab 机械原理,机械原理matlab编程

导杆的角位移曲线和刨刀的位移曲线导杆的角位移曲线和刨刀的位移曲线导杆的角位移曲线和刨刀的位移曲线导杆的角位移曲线和刨刀的位移曲线

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

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值