www.gusucode.com > MPSK,误码率性能仿真源码程序 > MPSK,误码率性能仿真源码程序/IET_MATLAB/functions/Poutage.m

    % Outage probability

function [pout,pout_floor]=Poutage(P0_N0,A2,sig_1,sig_2,g_th,fsr,frd,ch_dis)

for k=1:length(P0_N0)

alfa1=besselj(0,2*pi*fsr*ch_dis);
alfa2=besselj(0,2*pi*frd*ch_dis);
alfa=alfa1*alfa2;
    
c1=alfa^2*A2(k)*P0_N0(k);
c2=(alfa2^2+alfa^2+(1-alfa^2)*P0_N0(k)*sig_1)*A2(k);
c3=1+alfa^2;

d3=c2/c1/sig_1;
d2=4*c3/(c1*sig_1*sig_2);
d1=(4*c3*sig_2)/(c1*sig_1);

pout(k)=1-1/sig_2*exp(-d3*g_th)*sqrt(g_th*d1)*besselk(1,sqrt(d2*g_th));


end

pout_floor=1-exp(-g_th*(1-alfa^2)/alfa^2)*ones(1,length(P0_N0));

end