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