公平的席位分配(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)