communications3.m

时隙alora和伪贝叶斯算法的形成仿真程序 matlab m文件

% lamida=input('input lamida=?  ');
% n=input('imput n=?  ');
% n_t=input('Guess n =?  ');

clear;
lamida=1/exp(1)-0.2;
n=50;
n_t=180;
N=800;
n1=n;
n1_t=n_t;
m=1000;
qa=1-exp(-lamida/m);


for j=1:N
    %j
       qr=min(1,1/n1_t);
        n=n1;
        n_t=n1_t;
%simulation
        u(1:n)=1; u(n+1:m)=0;  
        l(1:n)=0; l(n+1:m)=binornd(1,qa,m-n,1);
        y=u+l;

        d=length(find(y~=0));
        depart=binornd(1,qr,d,1);
        if sum(depart)==1
            n1=d-1;
        else sum(depart)==0
            n1=d;
        end
            
%           
%         for k=1:length(depart)
%         y(index)=y(index)-depart(k);
%         end
% 
%         n1=length(find(y~=0));
        N1(j)=n1;
        
 %theory
        if sum(depart)==1|sum(depart)==0
        n1_t=max(lamida,(lamida+n_t-1));
        else
        n1_t=n_t+lamida+(exp(1)-2)^(-1);
        end
        N1_t(j)=n1_t;

%---------
end

plot(1:N,N1,'b');
hold on;
plot(1:N,N1_t,'r');
hold on;