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

佛自业障

程序介绍了数学建模中经典问题的两种解法，即席位分配问题！
3 I( @2 I) U2 t# A; x%适用于所有情况 BY Gu
' y0 F) B* Y7 U6 }" {0 K3 X+ Bclear all
) k6 Q$ F# q; m3 A: f' Q- _clc   %惯例Q值法分配席位，wy为席位数，ps为人数总和，R为分配方案
; B2 ~. P1 m" ]9 s. dwy=19;
P=[103 63 34]
%菜单选项
* {- d, g& J6 Y) M& a  jMENUN=menu('选择方法','惯例Q值法','dhondt法','都用');
%惯例Q值法
if (MENUN==1||MENUN==3)==1
n=size(P,2);
ps=sum(P);
RS=0;
% W+ N5 y7 K. G* v# Gfor i=1:n
9 Q8 i% F- O7 k9 c: S1 ~    R(i)=fix(P(i)/ps*wy);
    RS=RS+R(i);
/ E  ?& S  S0 ?* c3 B& O  t+ F0 l" jend
wy1=wy-RS;
for i=1:wy1
2 w- [: N& e3 |    for j=1:n
        Q(j,i)=P(j).^2/(R(j)*(R(j)+1));
9 k2 Y# A$ ?8 s2 R9 Z* w' Z    end
1 {0 g6 ~! e) L5 \' D6 |: C% q    t=Q(1,i);
5 G- m$ G1 J: C" M    a=1;
    for m=2:n
        if Q(m,i)>t
            t=Q(m,i);
6 B0 u; R4 P3 e* s            a=m;
        end
    end
    R(a)=R(a)+1;
* L7 i6 s' d3 T0 _1 [; ~2 `end
disp('惯例Q值法席位分配为：')
1 r4 r. s% N6 [R
end
%d'hondt法………………………………………………………………………………$$
if (MENUN==2||MENUN==3)==1
zs=P;
n=size(zs,2);
for i=1:n
0 P4 R. r5 F1 i. W) D$ j     for j=1:wy
. N" _( A( J; e+ E         r(i,j)=zs(i)/j;
7 `7 G1 B: h. g. T/ g" @) `     end
end
t=r(';
T=sort(t,'descend');
min=T(wy);
' u4 Z/ E7 k  w* }! ZR1=zeros(1,n);
for i=1:n
     for j=1:wy
: ^6 y) u* [7 w         if r(i,j)>min
             R1(i)=R1(i)+1;
$ I* e3 S+ k5 [  s         end      
     end
$ U2 r! p$ b$ e2 V& M. s( l# cend
rest=wy-sum(R1);%剩余席位
) c- V9 S. P5 ]$ t%对剩余席位的分配处理
# c5 W' V, P# B+ W3 e0 Dk=0;
for i=1:n
    if k>=rest%席位分配完毕，结束循环
, U5 A0 [" G5 @* @: S; b$ k        break;
    end
    for j=1:wy
        if r(i,j)==min
           R1(i)=R1(i)+1;
" ~+ }: E3 S7 I4 n           k=k+1;%对剩余席位分配计数
           break;
        end   
    end
end
disp ('dhondt法席位分配：')
R1
end
) i: ~# _6 S- ]4 ~' r+ {% \  @& _