席位分配问题的MATLAB程序
说明：
1. 本程序用三种方法，分别是惯例法、d’honht分配法和Q值法。 2. 可以模拟出任意一种分配情况，即可以推广到N种情形。 3. 三种分配方案供你选择，相互比较。 4. 请务必阅读注意事项。
注意：
1. 以下包含两个程序，下载完后把程序拷贝到matlab的M文件中， 2. 第一个程序可以任意命名，只要符合规范就可以（本人以”xiweifenpei”命名，
这样便于查看），第二个程序一定要命名为“xiwei”，因为程序中要用到函数。 3. 下载完后先把程序拷贝到txt文件中，再从txt拷贝到M文件中，这样可以避免乱
码。
程序一： clear all clc
disp('席位分配：') P=1000
p=[235 333 432] N=10
[x,y]=size(p); zu=x*y;
disp('惯例分配方法:') for i = 1:zu
n(i) =p(i)*N/P; end n;
m=n-fix(n); for i=1:zu
if n(i)==max(m)+fix(n(i)) n(i)=fix(n(i))+1; else
n(i)=fix(n(i)); end end
disp('惯例分配人数：') n
disp('d’honht方法:') pp=[]; for i=1:N pi=p/i;
pp=[pp; pi]; end pp
m=zeros(1,zu); for i=1:N
[x,y]=find(pp==max(pp(:))); pp(x,y)=0; m(y)=m(y)+1; end pp
disp('d’honht分配人数:') m
disp('Q值法分配方法:')
q=ones(1,zu); Q=[]; p;
for i=1:zu
Q(i)=p(i)*p(i)/(q(i)*(q(i)+1)); end Q;
xiwei(p,q,Q,N,zu) 程序二：
再次提醒，以“xiwei”为文件名保存， function xiwei(p,q,Q,N,zu) if sum(q)==N
disp('Q值法分配人数：') q
return; else
for i=1:zu
if Q(i)==max(Q) q(i)=q(i)+1;
Q(i)=p(i)*p(i)/(q(i)*(q(i)+1)); break; end end end
xiwei(p,q,Q,N,zu)