clc
clear all
x0=[15,210];
[x,fval]=fmincon('objfun22',x0,[],[],[],[],[],[],'confun22')
%%%优化结束后,计算优化值所得的结果
L=4430;
M=2150;
M1=2360; %主销轴线与地面交点的距离,实际的“轮距”
gam=x(1)*pi/180;
m=x(2);
for i=1:35
afa(i)=i*pi/180;%外转角循环
N=sqrt(M^2+m^2-2*M*m*sin(gam+afa(i)));
deta1=asin(m*cos(gam+afa(i))/N);
b=M-2*m*sin(gam);
beta=gam+deta1-asin((N^2+m^2-b^2)/(2*N*m)); %实际nei转角
beta0=atan(1/(M1/L+1/tan(afa(i)))); %理论nei转角
c(i)=abs(beta-beta0)*180/pi;
bt0(i,:)=beta0.*180/pi;
bt(i,:)=beta.*180/pi;
ai(i,:)=afa(i).*180/pi;
end
plot(ai,bt0,ai,bt,'r*');grid on;
figure(2);plot(ai,c);grid on;
评论2