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

佛自业障

程序介绍了数学建模中经典问题的两种解法，即席位分配问题！
; ?* n. R/ A3 o% N%适用于所有情况 BY Gu
9 E1 \3 C# `3 Rclear all
clc   %惯例Q值法分配席位，wy为席位数，ps为人数总和，R为分配方案
wy=19;
& v0 j. w6 p7 lP=[103 63 34]
%菜单选项
MENUN=menu('选择方法','惯例Q值法','dhondt法','都用');
. k! G& U8 v1 @" S%惯例Q值法
if (MENUN==1||MENUN==3)==1
) {4 W- z/ d1 b( m2 Rn=size(P,2);
& g$ w2 l8 ]: ~4 r2 Wps=sum(P);
8 X! S6 V) r" O, }9 s) {- kRS=0;
) O6 W  a' \8 O4 ^2 dfor i=1:n
# M% B8 J6 b: u' `) V5 W7 J    R(i)=fix(P(i)/ps*wy);
    RS=RS+R(i);
! W6 I5 i' A3 c* \end
( L; [) U: N  cwy1=wy-RS;
0 f1 Z. x# p5 ]9 hfor i=1:wy1
# m: T& i3 d  z6 y! W    for j=1:n
+ S+ r& w7 R; |        Q(j,i)=P(j).^2/(R(j)*(R(j)+1));
    end
. Y9 X+ m, m9 ?* k    t=Q(1,i);
' x' s" ~2 I) E$ v' D- M    a=1;
    for m=2:n
! V; X$ h# k7 E+ ]/ m& D- e        if Q(m,i)>t
            t=Q(m,i);
            a=m;
& E; F) X! X% v$ D" E: [9 Q        end
    end
    R(a)=R(a)+1;
end
disp('惯例Q值法席位分配为：')
R
end
%d'hondt法………………………………………………………………………………$$
6 O5 Y. w  Z/ G1 K& Lif (MENUN==2||MENUN==3)==1
zs=P;
' g" {# @3 x8 An=size(zs,2);
0 u  D4 o; o' ^+ kfor i=1:n
     for j=1:wy
         r(i,j)=zs(i)/j;
8 m" S- D$ [! _' w- [     end
end
) ?( K# N! y" z" ~3 U. |& e* `' h: gt=r(';
+ n% l* m5 U  W4 {9 x$ mT=sort(t,'descend');
min=T(wy);
; ]# ?* J7 H7 s7 V1 ~R1=zeros(1,n);
for i=1:n
     for j=1:wy
         if r(i,j)>min
             R1(i)=R1(i)+1;
         end      
0 N$ ^2 M; S1 S8 u     end
end
! d( p' J9 z3 x4 a* c0 t7 Grest=wy-sum(R1);%剩余席位
%对剩余席位的分配处理
k=0;
for i=1:n
1 i2 Q* W. ^& r/ s5 P2 I3 i. [    if k>=rest%席位分配完毕，结束循环
8 {5 x' d2 U# a: X+ h& M        break;
    end
    for j=1:wy
, l4 t9 R( \5 e' I        if r(i,j)==min
           R1(i)=R1(i)+1;
( J! A. ]4 n' b1 e+ t+ t; k; K' T* @0 o           k=k+1;%对剩余席位分配计数
* E- I* }9 [6 h5 U+ c9 e% F0 ]           break;
$ R" A. A; [$ s% Z0 k+ S- ~        end   
    end
3 f6 b/ _3 [: t" ]  ^' @1 |6 [end
9 w0 D/ [, J# x; ^7 ldisp ('dhondt法席位分配：')
R1
end
9 t1 ~6 W2 g3 z4 \- D0 e
9 e% W$ i& q- v+ U% O/ X8 O
2 c3 [# R' z5 w( b1 q! ~! c