g_q_test.asv

来自「MATSNL is a package of MATLAB M-files fo」· ASV 代码 · 共 64 行

PER=0.7;
Is_Succ=0;
num_re_tx=0;

while 1
    Is_Succ=bern_gen(1- PER);
    if Is_Succ==1
        brea
    
    if(num_re_tx>max_num_re_tx)
        break;
    end
    num_re_tx=num_re_tx+1;

end
     


    p=PER;
    Q=max_num_re_tx;
%     G_Q=(1/(1-p))*(1+p^(Q-1))+Q*p^(2*Q-1);
G_Q=(1-Q*p^Q+(p/(1-p))*(1-p^(Q-1)))+Q*p^Q;
    N=floor(G_Q*(n/TDMA_cycle)*data_size*(pkt_hdr_size+unit_pkt_size)/(unit_pkt_size*tx_rate));

    if N>0
        Z=G_Q*data_size*(pkt_hdr_size+unit_pkt_size)/(unit_pkt_size*tx_rate*N);
        L=TDMA_cycle/(N)+(n-1)*TDMA_cycle/n;
        wait_time_thy(h)=N*L;
        tx_time_thy(h)=N*Z;
    else
        N=1;L=TDMA_cycle;
        Z=G_Q*data_size*(pkt_hdr_size+unit_pkt_size)/(unit_pkt_size*tx_rate);
        wait_time_thy(h)=N*L;
        tx_time_thy(h)=N*Z;
    end

    disp(['wait time(thy): ', num2str(N*L)]);
    disp(['tx time(thy): ', num2str(N*Z)]);


    wait_time_sim(h)=mean(wait_time_log);
    wait_time_std(h)=std(wait_time_log);

    tx_time_sim(h)=mean(tx_time_log(k));
    tx_time_std(h)=std(tx_time_log(k));
    N_log(h)=N;
end

plot( TDMA_cycle_log,wait_time_sim,'b',TDMA_cycle_log,wait_time_thy,'.-r');
% hold on;
% plot( TDMA_cycle_log,tx_time_sim,'b',TDMA_cycle_log,tx_time_thy,'.-r');
set(gca,'FontSize',24);
xlabel('TDMA cycle length');
ylabel('Time duraiton, (sec)');


% h1=errorbar((P),mean_num_sim, L,U,'+');
% hold on;
% h2=plot( (P),mean_num_sim);
% set(gca,'FontSize',18);
% xlabel('The maximum number of backoff, P');
% ylabel('Average listening delay, L(sec)');