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

佛自业障

程序介绍了数学建模中经典问题的两种解法，即席位分配问题！
9 d" q) U' q3 H. c( z%适用于所有情况 BY Gu
clear all
/ b$ I/ S* j! T5 Uclc   %惯例Q值法分配席位，wy为席位数，ps为人数总和，R为分配方案
wy=19;
/ Q0 ^, \; Z. f' k! c5 N; v% vP=[103 63 34]
%菜单选项
: a  ]) [0 ]! H8 \MENUN=menu('选择方法','惯例Q值法','dhondt法','都用');
%惯例Q值法
if (MENUN==1||MENUN==3)==1
n=size(P,2);
ps=sum(P);
RS=0;
. D: o& s' b1 jfor i=1:n
  ^7 n: S1 @% i: `- u" j9 v    R(i)=fix(P(i)/ps*wy);
2 T1 c$ P4 {0 [5 l+ P8 j5 A    RS=RS+R(i);
# G* A- P( b6 Dend
wy1=wy-RS;
for i=1:wy1
    for j=1:n
        Q(j,i)=P(j).^2/(R(j)*(R(j)+1));
. P0 V2 T$ h8 u    end
    t=Q(1,i);
    a=1;
    for m=2:n
/ j! N+ I+ |+ E3 |* Z        if Q(m,i)>t
" r; G; n9 R9 \! y, e            t=Q(m,i);
            a=m;
        end
    end
    R(a)=R(a)+1;
2 d4 r6 ?3 x0 B8 A4 aend
% e0 x$ G( k: hdisp('惯例Q值法席位分配为：')
* C- \4 R' [$ T# x5 DR
) y* R& w& j7 W  F. J" Wend
%d'hondt法………………………………………………………………………………$$
if (MENUN==2||MENUN==3)==1
/ q  `  A+ p4 u5 d% e/ s" }zs=P;
n=size(zs,2);
for i=1:n
     for j=1:wy
         r(i,j)=zs(i)/j;
: X* |9 \, h% Z- `0 O8 g" O     end
# D& {" j" }$ E+ l, jend
. C2 w8 v' z2 w8 L) W& Ct=r(';
! u! K( c% z4 t2 E# J# XT=sort(t,'descend');
min=T(wy);
5 Q+ J* f# v$ R7 X1 C1 B" rR1=zeros(1,n);
, T: G9 Y/ M. h# R7 k' O* W$ Qfor i=1:n
     for j=1:wy
4 _7 D. F& \  f$ k* H         if r(i,j)>min
             R1(i)=R1(i)+1;
         end      
     end
end
rest=wy-sum(R1);%剩余席位
%对剩余席位的分配处理
k=0;
for i=1:n
; F4 d, {4 G4 g" ?( y: r    if k>=rest%席位分配完毕，结束循环
        break;
5 M6 D1 l3 j1 r+ s6 w/ B    end
- I% `: w% c7 Z  d8 b1 D( A    for j=1:wy
. y: H5 [- \' y! H/ v6 S- F        if r(i,j)==min
           R1(i)=R1(i)+1;
           k=k+1;%对剩余席位分配计数
# a: o' V8 F4 r! q$ S           break;
        end   
    end
end
disp ('dhondt法席位分配：')
4 B8 |4 M5 n/ T  @9 j2 O+ \R1
end