Commit 8959174b authored by Riccardo La Placa's avatar Riccardo La Placa
Browse files

Geometric dist added

parent e39b6875

power_dist_test.m

+48 −8
Original line number Diff line number Diff line
@@ -126,14 +126,34 @@ for i = 1:length(mario) 
            meanpow/vsuetot

            meanpow/vsuedat

            extrtot4 = ptotk(p4_1,extrpoiss);



            ncr = 1;

            epscr = 1 - (length(xtemp(xtemp == 1))/N)/(mi*ptot0);

            pcr = 1- (1-epscr)^(1/ncr);

            expec1 = 1+ epscr + (3*(ncr-1)/(2*ncr))*epscr^2;

            var1 = epscr + (epscr^2)*(7-9/ncr)/2;

            enf = 1 + var1/(expec1^2);

            expectot = expec1*mi

            vartot = var1*mi + mi*expec1^2

            qcr = (1-epscr)^(1/ncr);

            p1_1 = zeros(1,maxcount);

            p1_1(1:maxcount) = (pcr.^(0:maxcount-1)).*(1-pcr);

            vard = var(xtemp)

            vsuetot = vartot/expectot

            vsuedat = vard/mean(xtemp)

            meanpow = mean(Y0(100*Tseg+1:end))/2

            meanpow/vsuetot

            meanpow/vsuedat

            extrtot1 = ptotk(p1_1,extrpoiss);



            tiledlayout(2,1)

            nexttile

            semilogy([1:maxcount],extrtot4,[1:maxcount],extrtot8,[1:maxcount],datadist,'+')

            legend('extrtot4','extrtot8','datadist')

            semilogy([1:maxcount],extrtot1,[1:maxcount],extrtot4,[1:maxcount],extrtot8,[1:maxcount],datadist,'+')

            legend('extrtot1','extrtot4','extrtot8','datadist')

            nexttile

            semilogy([1:maxcount],abs(extrtot4-datadist)./datadist,[1:maxcount],abs(extrtot8-datadist)./datadist)

            semilogy([1:maxcount],abs(extrtot1-datadist)./datadist,[1:maxcount],abs(extrtot4-datadist)./datadist,[1:maxcount],abs(extrtot8-datadist)./datadist)

            ylim([0.0001 100])

            legend('abs(extrtot4-datadist)/datadist','abs(extrtot8-datadist)/datadist')

            legend('abs(extrtot1-datadist)/datadist','abs(extrtot4-datadist)/datadist','abs(extrtot8-datadist)/datadist')

            yline([0.01 0.1 0.5 1.0])

        else

            disp('Too many counts!')
@@ -244,14 +264,34 @@ for m = 1:M 
        meanpow/vsuetot

        meanpow/vsuedat

        extrtot4 = ptotk(p4_1,extrpoiss);



        ncr = 1;

        epscr = 1 - (length(xtemp(xtemp == 1))/N)/(mi*ptot0);

        pcr = 1- (1-epscr)^(1/ncr);

        expec1 = 1+ epscr + (3*(ncr-1)/(2*ncr))*epscr^2;

        var1 = epscr + (epscr^2)*(7-9/ncr)/2;

        enf = 1 + var1/(expec1^2);

        expectot = expec1*mi

        vartot = var1*mi + mi*expec1^2

        qcr = (1-epscr)^(1/ncr);

        p1_1 = zeros(1,maxcount);

        p1_1(1:maxcount) = (pcr.^(0:maxcount-1)).*(1-pcr);

        vard = var(xtemp)

        vsuetot = vartot/expectot

        vsuedat = vard/mean(xtemp)

        meanpow = mean(Y0(100*Tseg+1:end))/2

        meanpow/vsuetot

        meanpow/vsuedat

        extrtot1 = ptotk(p1_1,extrpoiss);



        tiledlayout(2,1)

        nexttile

        semilogy([1:maxcount],extrtot4,[1:maxcount],extrtot8,[1:maxcount],datadist,'+')

        legend('extrtot4','extrtot8','datadist')

        semilogy([1:maxcount],extrtot1,[1:maxcount],extrtot4,[1:maxcount],extrtot8,[1:maxcount],datadist,'+')

        legend('extrtot1','extrtot4','extrtot8','datadist')

        nexttile

        semilogy([1:maxcount],abs(extrtot4-datadist)./datadist,[1:maxcount],abs(extrtot8-datadist)./datadist)

        semilogy([1:maxcount],abs(extrtot1-datadist)./datadist,[1:maxcount],abs(extrtot4-datadist)./datadist,[1:maxcount],abs(extrtot8-datadist)./datadist)

        ylim([0.0001 100])

        legend('abs(extrtot4-datadist)/datadist','abs(extrtot8-datadist)/datadist')

        legend('abs(extrtot1-datadist)/datadist','abs(extrtot4-datadist)/datadist','abs(extrtot8-datadist)/datadist')

        yline([0.01 0.1 0.5 1.0])

    else

        disp('Too many counts!')