程序如下
gama = 2*pi/9;
h = 300;
Rr = 400;
Dr = 850;
R1 = h/(2*sin(gama/2));
theta1 = linspace(0,gama/8,10);
theta2 = linspace(gama/8,gama/4,10);
theta3 = linspace(gama/4,3*gama/8,10);
theta4 = linspace(3*gama/8,gama/2,10);
theta5 = linspace(gama/2,5*gama/8,10);
theta6 = linspace(5*gama/8,3*gama/4,10);
theta7 = linspace(3*gama/4,7*gama/8,10);
theta8 = linspace(7*gama/8,gama,10);
theta = [theta1 theta2 theta3 theta4 theta5 theta6 theta7 theta8];
Y1 = -80.9186*theta1.^3+279.4697*theta1.^2-153.9023*theta1+474.3163;
Y2 = -80.9186*(theta2-gama/8).^3+258.2853*(theta2-gama/8).^2-106.7218*(theta2-gama/8)+462.9603;
Y3 = -407.0496*(theta3-gama/4).^3+237.1009*(theta3-gama/4).^2-63.7437*(theta3-gama/4)+455.5382;
Y4 &#
