2022B第(1)问

污浊的双黑

已于 2023-02-13 10:14:07 修改

阅读量90

点赞数

于 2022-11-25 23:24:23 首次发布

本文链接：https://blog.csdn.net/qq_53011270/article/details/128046233

版权

matlab

% N=2
alpha = 2*pi/9
for n = 2:10    
    if alpha1 <= alpha2
%        sin (alpha1)/R == sin (pi - theta - alpha1)/p;
%        sin (alpha2)/R == sin (pi - (theta - (n-1)*alpha + alpha2))/p;
       syms R p alpha1 alpha3 theta;
       S = solve(sin (alpha1)/R == sin (pi - theta - alpha1)/p,sin (alpha2)/R == sin (pi - (theta - (n-1)*alpha + alpha2))/p,p,theta);
       pp=S.p;theta=S.theta;
    else
%        sin (alpha1)/R == sin (pi - theta - alpha1)/p;
%        sin (alpha2)/R == sin (pi - (2*pi - theta - (n-1)*alpha + alpha2))/p;
       syms R p alpha1 alpha3 theta;
       S = solve(sin (alpha1)/R == sin (pi - theta - alpha1)/p,sin (alpha2)/R == sin (pi - (2*pi - theta - (n-1)*alpha + alpha2))/p,p,theta);
       pp=S.p;theta=S.theta;
    end
end


alpha = 2*pi/9; N=9
for n = 2:10
    if alpha1 <= alpha2
%        sin(alpha1)/R == sin(pi - theta - alpha1)/p;
%        sin(alpha2)/R == sin(pi - (theta - (10-n)*alpha + alpha2))/p;
         syms R p alpha1 alpha3 theta;
         S = solve(sin(alpha1)/R == sin(pi - theta - alpha1)/p,sin(alpha2)/R == sin(pi - (theta - (10-n)*alpha + alpha2))/p,p,theta);
         pp=S.p;theta=S.theta;
      else
%        sin(alpha1)/R == sin(pi - theta - alpha1)/p;
%        sin(alpha2)/R == sin(pi - (2*pi - theta - (10-n)*alpha + alpha2))/p;
         syms R p alpha1 alpha3 theta;
         S = solve(sin(alpha1)/R == sin(pi - theta - alpha1)/p,sin(alpha2)/R == sin(pi - (2*pi - theta - (10-n)*alpha + alpha2))/p,p,theta);
         pp=S.p;theta=S.theta;
    end
end


alpha = 2*pi/9;
for N = 3:5;n = 2:10
  if alpha3 <pi/2     
          alpha1 + alpha3 == alpha2
%          sin(alpha1)/R = sin(pi - theta - alpha1)/p;
%          sin(alpha2)/R = sin(pi - (theta - (n-1)*alpha + alpha2))/p;
           syms R p alpha1 alpha3 theta;
           S = solve(sin(alpha1)/R == sin(pi - theta - alpha1)/p,sin(alpha2)/R == sin(pi - (theta - (n-1)*alpha + alpha2))/p,p,theta);
           pp=S.p;theta=S.theta;       
       
       elseif alpha1 + alpha2 == alpha3
%          sin(alpha1)/R == sin(pi - theta - alpha1)/p;
%          sin(alpha2)/R == sin(pi - ( - theta + (n-1)*alpha + alpha2))/p;
           syms R p alpha1 alpha3 theta;
           S = solve(sin(alpha1)/R == sin(pi - theta - alpha1)/p,sin(alpha2)/R == sin(pi - ( - theta + (n-1)*alpha + alpha2))/p,p,theta);
           pp=S.p;theta=S.theta;
       else
%          sin(alpha1)/R == sin(pi - theta - alpha1)/p;
%          sin(alpha2)/R == sin(pi - (2*pi - theta - (n-1)*alpha + alpha2))/p;
           syms R p alpha1 alpha3 theta;
           S = solve(sin(alpha1)/R == sin(pi - theta - alpha1)/p,sin(alpha2)/R == sin(pi - (2*pi - theta - (n-1)*alpha + alpha2))/p,p,theta);
           pp=S.p;theta=S.theta;
  end
end


alpha = 2*pi/9;
for N = 6:8
   if alpha3 <pi/2  
         if alpha1 +  alpha3 == alpha2
%        sin(alpha1)/R == sin(pi - theta - alpha1)/p;
%        sin(alpha2)/R == sin(pi - (theta - (10-n)*alpha + alpha2))/p;
         syms R p alpha1 alpha3 theta;
         S = solve(sin(alpha1)/R == sin(pi - theta - alpha1)/p,sin(alpha2)/R == sin(pi - (theta - (10-n)*alpha + alpha2))/p,p,theta);
         pp=S.p;theta=S.theta;
         end
       
      elseif alpha1 +  alpha2 == alpha3
%        sin(alpha1)/R == sin(pi - theta - alpha1)/p;
%        sin(alpha2)/R == sin(pi - ( - theta + (10-n)*alpha + alpha2))/p;
         syms R p alpha1 alpha3 theta;
         S = solve(sin(alpha1)/R == sin(pi - theta - alpha1)/p,sin(alpha2)/R == sin(pi - ( - theta - (10-n)*alpha + alpha2))/p,p,theta);
         pp=S.p;theta=S.theta;
      else
%        sin(alpha1)/R == sin(pi - theta - alpha1)/p;
%        sin(alpha2)/R == sin(pi - (2*pi - theta - (10-n)*alpha + alpha2))/p;
         syms R p alpha1 alpha3 theta;
         S = solve(sin(alpha1)/R == sin(pi - theta - alpha1)/p,sin(alpha2)/R == sin(pi - (2*pi - theta - (10-n)*alpha + alpha2))/p,p,theta);
         pp=S.p;theta=S.theta;
   end
end