g_q_test.m

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

num_sample=10000;

for i=1:num_sample
PER=0.6;
Is_Succ=0;
num_re_tx=0;
max_num_re_tx=3;

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

        if(num_re_tx>=max_num_re_tx)
            break;
        end
    end
    num_re_tx_log(i)=num_re_tx;
end

p=PER;
Q=max_num_re_tx;
G_Q=(1-Q*p^Q+(p/(1-p))*(1-p^(Q-1)))+Q*p^Q;

disp(['simul: ',  num2str(mean(num_re_tx_log))]);
disp(['thy: ',  num2str(G_Q)])