公平的席位分配(MATLAB程序)

席位分配问题的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)。