matlab_蜂窝系统中小区用户通信调度算法的程序，包含了三种经典调度算法，比如比例调度算法

clc; clear; num_user=10; %用户数 ntimeslot=100; %时隙数 nt=2; %发送天线数 nr=2; %接收天线数 snr_db=10; Inr=eye(nr,nr); snr=10^(snr_db/10); for k=1:num_user %瑞利信道 for m=1:ntimeslot h(:,:,k,m)=(1/sqrt(2))*(randn(nr,nt)+1i*randn(nr,nt)); end end for m=1:num_user %capa为虚数，公平性的计算中需要用绝对值来计算 for k=1:ntimeslot capa(m,k)=log2(det((Inr+snr*h(:,:,m,k)*h(:,:,m,k)'/nt))); abscapa(m,k)=abs(capa(m,k)); end end user_capa1=zeros(num_user,ntimeslot); for k=0:(ntimeslot/num_user-1) for m=1:num_user n=m+k*num_user; capa1(n)=capa(m,n); %轮询调度算法，每个时隙依次取用户 user_capa1(m,n)=capa(m,n); %计算公平性需要每个用户的速率 end end s11=0; %求所有时隙吞吐量的和，虚数 for k=1:length(capa1) s11=s11+abs(capa1(k)); end sum_user1=sum(abs(user_capa1),2); %每个用户获得的吞吐量之和，虚数 s12=0; for k=1:num_user s12=s12+sum_user1(k).^2; end J1=s11.^2/(num_user*s12);%Jain's准则 user_capa2=zeros(num_user,ntimeslot);%设置这样一个矩阵，记录被调度的用户和瞬时容量，方便求每个用户获得的容量和 %最大载干比调度算法，怎么求最大载干比? for k=1:ntimeslot [max_capa,m]=max(abscapa(:,k)); %取用户容量最大的情况,这里比较是取成最大值 capa2(k)=capa(m,k); user_capa2(m,k)=capa(m,k); end s21=0; for k=1:length(capa2) s21=s21+abs(capa2(k)); end sum_user2=sum(abs(user_capa2),2); s22=0; for k=1:num_user s22=sum_user2(k).^2+s22; end J2=s21.^2/(num_user*s22); for k=1:num_user R(k,1)=1; %设置初始值为1,用户瞬时速率和平均速率都取成吞吐量的绝对值 for t=1:ntimeslot p(k,t)=abscapa(k,t)/R(k,t); R(k,t+1)=(1-1/ntimeslot)*R(k,t)+abscapa(k,t)/ntimeslot;%平均速率更新算法 end end user_capa3=zeros(num_user,ntimeslot); for k=1:ntimeslot [pmax,m]=max(p(:,k)); capa3(k)=capa(m,k); user_capa3(m,k)=capa(m,k);%记录每个用户被调度的历史 end s31=0; for k=1:ntimeslot s31=s31+abs(capa3(k)); %公平性的计算中用户速率取成系统吞吐量的绝对值 end s32=0; sum_user3=sum(abs(user_capa3),2); for k=1:num_user s32=s32+sum_user3(k).^2; end J3=s31.^2/(num_user*s32); k=1:ntimeslot plot(k,abs(capa1),'kd-'); hold on; plot(k,abs(capa2),'bd--'); hold on; plot(k,abs(capa3),'rd-.'),xlabel('时隙'),ylabel('吞吐量'),title('轮询，最大载干比及比例调度算法');