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

佛自业障

程序介绍了数学建模中经典问题的两种解法，即席位分配问题！
( C* q- t( K/ S9 T%适用于所有情况 BY Gu
clear all
# N) I$ l& R- t6 @2 ^clc   %惯例Q值法分配席位，wy为席位数，ps为人数总和，R为分配方案
wy=19;
+ u1 y( e2 c/ w9 v' fP=[103 63 34]
%菜单选项
MENUN=menu('选择方法','惯例Q值法','dhondt法','都用');
%惯例Q值法
- D9 \3 T* a  `" x: _/ aif (MENUN==1||MENUN==3)==1
n=size(P,2);
ps=sum(P);
RS=0;
for i=1:n
. L6 }3 p- W, m    R(i)=fix(P(i)/ps*wy);
    RS=RS+R(i);
- r, d3 e) F* M7 {end
wy1=wy-RS;
! B/ f  |+ R4 F5 z# jfor i=1:wy1
    for j=1:n
        Q(j,i)=P(j).^2/(R(j)*(R(j)+1));
9 K$ l) `- V9 ]  k    end
    t=Q(1,i);
    a=1;
    for m=2:n
  z$ \. _0 p- @, j/ d# r4 d7 H. K        if Q(m,i)>t
/ e& T; z3 Y) N3 Y, w/ d4 `            t=Q(m,i);
$ N1 A' r  K  ~- i+ `            a=m;
        end
# \* z" k$ ]! V& t$ t& E- z    end
    R(a)=R(a)+1;
end
. ~( t1 \, w, \# R  ]disp('惯例Q值法席位分配为：')
# q& M+ v6 _2 y) X, \R
end
: A8 p  ~1 e) z# t2 j% F& E- L%d'hondt法………………………………………………………………………………$$
6 i. G6 g. D0 {  c5 |) L- B- ]if (MENUN==2||MENUN==3)==1
zs=P;
n=size(zs,2);
$ ~, w' q, B" z( M  p' ]  Y. Lfor i=1:n
     for j=1:wy
         r(i,j)=zs(i)/j;
3 U) V2 ]0 j; a/ V     end
end
t=r(';
) |0 w" v! v# U5 i2 G- E4 ]( |T=sort(t,'descend');
min=T(wy);
; {  J; f" ~6 l5 A. z! t- L2 DR1=zeros(1,n);
for i=1:n
) w# {7 k+ P  Y0 @# `     for j=1:wy
; X# H* d2 x7 r" H         if r(i,j)>min
8 E  E6 o8 f/ N  h" o5 I             R1(i)=R1(i)+1;
/ _8 G) O" l! w         end      
1 \/ Y' O* u, {# @: h2 ~# v6 i9 V* h     end
end
) C. V2 h# ~- U7 n  \/ m9 brest=wy-sum(R1);%剩余席位
+ f0 u1 e1 G/ S, P+ S! k* s5 U%对剩余席位的分配处理
k=0;
* Y$ r8 y; O: m! x: ffor i=1:n
$ k3 S* j. C: N4 _6 L3 ?+ u) \    if k>=rest%席位分配完毕，结束循环
+ j* B9 [3 ~/ K        break;
& _* v+ I  h9 L/ s( a    end
9 z, \9 j# }+ I* b2 U; C) s    for j=1:wy
& \$ s; M5 @( [        if r(i,j)==min
- D$ c* f8 H. E  L3 I           R1(i)=R1(i)+1;
1 p2 `; K: @7 h4 p           k=k+1;%对剩余席位分配计数
: ?" E7 k5 P$ }1 K  q9 J$ \           break;
4 U( W7 D+ v0 T1 Q% S* F; @        end   
    end
end
* `5 A0 a+ S9 p1 p5 @. Edisp ('dhondt法席位分配：')
R1
end
- y3 W+ |- X2 H3 `  M8 K
, o9 @. P1 l% V' }! }
/ j) _3 J+ Y# w. e# F8 ^; \' l