席位分配问题——惯例Q值法和d'hondt法的MATLAB程序

佛自业障

程序介绍了数学建模中经典问题的两种解法，即席位分配问题！
, @. K' \( }% N) T" j0 I0 x5 ~%适用于所有情况 BY Gu
, a3 _: s8 _: C: D7 kclear all
/ b& n% p1 I6 @" w% w5 y" @clc   %惯例Q值法分配席位，wy为席位数，ps为人数总和，R为分配方案
wy=19;
P=[103 63 34]
! n$ Z6 q0 N0 L' ~2 p%菜单选项
MENUN=menu('选择方法','惯例Q值法','dhondt法','都用');
  i; M) V, w* T+ [( V( J6 T0 `2 i%惯例Q值法
if (MENUN==1||MENUN==3)==1
9 K! Y% h: m( |1 rn=size(P,2);
/ N0 c+ L' B' {7 ?" E1 F+ ~; R/ i$ tps=sum(P);
7 r6 n: w3 \! K7 U- f9 ?RS=0;
for i=1:n
    R(i)=fix(P(i)/ps*wy);
6 k6 }% Z, P+ s1 [: ~0 v    RS=RS+R(i);
, r/ _9 z! A' c: H- N% lend
  e1 y& g; @" d9 i0 l4 N( awy1=wy-RS;
1 W5 O9 V1 Z; Ifor i=1:wy1
    for j=1:n
        Q(j,i)=P(j).^2/(R(j)*(R(j)+1));
    end
    t=Q(1,i);
    a=1;
0 K3 R+ V  l  m% Z# ~    for m=2:n
8 E7 u9 o/ G2 \8 y        if Q(m,i)>t
. U5 h7 p+ K4 O8 y* C: p            t=Q(m,i);
            a=m;
        end
    end
    R(a)=R(a)+1;
end
disp('惯例Q值法席位分配为：')
R
: N+ y* q7 d+ P. E8 V) yend
%d'hondt法………………………………………………………………………………$$
* ]2 A$ _5 {% ~0 [0 W/ P& m4 v0 bif (MENUN==2||MENUN==3)==1
zs=P;
/ k6 e5 {1 W; L) C: Sn=size(zs,2);
" h8 e( u7 n# p; ]& x" U5 D  ofor i=1:n
     for j=1:wy
         r(i,j)=zs(i)/j;
     end
0 v$ j% m9 Y5 X; K( @8 r) q# ]; Qend
t=r(';
T=sort(t,'descend');
min=T(wy);
R1=zeros(1,n);
for i=1:n
     for j=1:wy
: O7 A0 M" P- m         if r(i,j)>min
             R1(i)=R1(i)+1;
         end      
     end
end
2 g  R* `7 J7 Y- Rrest=wy-sum(R1);%剩余席位
  }0 e1 |6 X1 P! u5 G%对剩余席位的分配处理
k=0;
for i=1:n
    if k>=rest%席位分配完毕，结束循环
        break;
    end
' z# }- Q5 j' L1 A0 t    for j=1:wy
        if r(i,j)==min
) l( ^7 c! X. O" |1 z           R1(i)=R1(i)+1;
           k=k+1;%对剩余席位分配计数
           break;
        end   
( q$ `4 R% m5 z% t$ \' [$ [    end
end
disp ('dhondt法席位分配：')
% ?" u  Q+ c8 s5 Y- CR1
end

/ w; ^2 [1 s: x) h