03-连杆机构的运动精度综合-2.m

MATLAb 机械源码文件好东西啊

disp '            ****** 曲柄滑块机构的等影响法精度综合 *******'
N=input('            输入机构运动精度影响尺度数目       N = ');
H=input('            输入滑块行程的均值(mm)             H = ');
P=input('            输入曲柄轴心至滑销最远距离(mm)     P = ');
DH=input('            输入滑块位置允许误差(mm)          DH = ');
R=H/2;
fprintf(1,'              曲柄长度的均值(mm)               R = %3.3f \n',R) 
L=P-R;
fprintf(1,'              连杆长度的均值(mm)               L = %3.3f \n',L) 
theta=0:10:360;
hd=theta.*pi/180;
% 计算曲柄长度和滑块长度的影响系数(偏导数的最大绝对值)
CR=1-cos(hd);
CL=0.5.*sin(hd).^2;
CRM=max(abs(1-cos(hd)));
CLM=max(abs(0.5.*sin(hd).^2));
fprintf(1,'              曲柄长度影响系数的最大绝对值   CRM = %3.6f \n',CRM) 
fprintf(1,'              连杆长度影响系数的最大绝对值   CLM = %3.6f \n',CLM) 
% 计算曲柄长度和滑块长度的最大允许偏差
DRM=DH/sqrt(N)/CRM;
DLM=DH/sqrt(N)/CLM;
fprintf(1,'              曲柄长度允许的最大偏差(mm)     DRM = %3.6f \n',DRM) 
fprintf(1,'              连杆长度允许的最大偏差(mm)     DLM = %3.6f \n',DLM) 
plot(theta,CR,'r')
gtext('曲柄长度影响系数曲线')
title('\bf 机构尺度影响系数线图')
xlabel('\bf 曲柄转角\theta(度)')
ylabel('\bf 尺度影响系数')
hold;
plot(theta,CL,'k')
gtext('连杆长度影响系数曲线')