如题，业内人士都懂的吧？
figure('Name','overrun torque');
syms  m rl w r ra a f
m=0.14739;r=0.02475;w=(800*pi)/30;rl=0.266129;ra=linspace(0,4*pi,200);
h1=subs(-m.*r*(w^2).*(cos(ra)+rl.*cos(2.*(ra))).*r.*(sin(ra).*sqrt(1-(rl^2).*sin(ra).^2)+rl.*sin(ra).*cos(ra))./sqrt(1-(rl^2).*sin(ra).^2));
plot(ra*180/pi,h1,'r')
hold
SYMS  m r w rl ra
m=0.14739;r=0.02475;w=(800*pi)/30;rl=0.266129;ra=linspace(0,4*pi,200);
h2=subs((-m.*r*(w^2).*(cos(ra)+rl.*cos(2.*(ra))).*r.*(sin(ra).*(1-1/2.*(rl^2).*sin(ra).^2)+rl.*sin(ra).*cos(ra)))./(1-1/2.*(rl^2).*sin(ra).^2));
plot(ra*180/pi,h2)
xlabel('CA/°');
ylabel('N*m');
title('overrun torque');
grid on
hold off