%给出五个角度的范围
Mix_Max1=deg2rad([-90 90]); %-90 to 90
Mix_Max2=deg2rad([-165 -15]); %-165 to -15
Mix_Max3=deg2rad([-60 -60]); %-60 to 60
Mix_Max4=deg2rad([30 150]); %30 to 150
Mix_Max5=deg2rad([-60 60]); %-60 to 60
for i=1:5000
%设置转角theta1 to theta5
theta1=Mix_Max1(1,1) + (Mix_Max1(1,2)-Mix_Max1(1,1))*rand(1);
theta2=Mix_Max2(1,1) + (Mix_Max2(1,2)-Mix_Max2(1,1))*rand(1);
theta3=Mix_Max3(1,1) + (Mix_Max3(1,2)-Mix_Max3(1,1))*rand(1);
theta4=Mix_Max4(1,1) + (Mix_Max4(1,2)-Mix_Max4(1,1))*rand(1);
theta5=Mix_Max5(1,1) + (Mix_Max5(1,2)-Mix_Max5(1,1))*rand(1);
%给出三个坐标值的表达式:
X(1,i)=423*cos(theta1)*cos(theta2) - 121*sin(theta1) - 168*cos(theta5)*sin(theta1) + 112*cos(theta4)*(cos(theta1)*cos(theta2)*sin(theta3) + cos(theta1)*cos(theta3