5章MATLAB.rar_wireless green_中断概率_中断概率 认知_认知无线网络_认知网络

%Traditional case function pout=pro(es,dsd,dsr,dsu0,duu0,dur,dud,drd,dru0,rs,ru,p0,e0,N) es=10^(es/10); e0=10^(e0/10); ts=2^(rs)-1; tu=2^(2*ru)-1; t0=exp(-ts/(es*dsd))/(1-p0)-1; eu=es*dsd*max(t0,0)/(ts*dud); eu=min(eu,e0); t0=exp(-ts/(es*dsd))/(1-p0)-1; er=es*dsd*max(t0,0)/(ts*drd); er=min(er,e0); if eu==0 pout=1; else b=tu*(1+ts)/(1-tu*ts); t1=eu*dur*exp(-tu/(eu*dur))/(eu*dur+tu*es*dsr); t2=tu*es*dsr*exp(1/(es*dsr)-b/(eu*dur)-b/(tu*es*dsr))/(tu*es*dsr+eu*dur); t3=es*dsr*exp(-ts/(es*dsr)-tu/(eu*dur)-tu*ts/(es*dsr))/(es*dsr+ts*eu*dur); t4=es*dsr*exp(-ts/(es*dsr)-b/(eu*dur)-b*ts/(es*dsr))/(es*dsr+ts*eu*dur); if (tu*ts)<1 pd0=t1+t2+t3-t4; else pd0=t1+t3; end %Un-decode if (2-ts*tu)>0 a=ts*(2+tu)/(2-ts*tu); tem1=1-2*eu*duu0*exp(-tu/(2*eu*duu0))/(2*eu*duu0+tu*es*dsu0); tem2=ts*eu*duu0*exp(1/(eu*duu0)-a/(es*dsu0)-a/(ts*eu*duu0))... /(ts*eu*duu0+es*dsu0); tem3=2*eu*duu0*exp(-tu/(2*eu*duu0)-a/(es*dsu0)-a*tu/(2*eu*duu0))... /(2*eu*duu0+tu*es*dsu0); tem4=es*dsu0*exp(-ts/(es*dsu0)-tu/(2*eu*duu0)-ts*tu/(2*es*dsu0))... /(es*dsu0+ts*eu*duu0); pout0=tem1-tem2+tem3-tem4; else tem5=1-2*eu*duu0*exp(-tu/(2*eu*duu0))/(2*eu*duu0+tu*es*dsu0); tem6=es*dsu0*exp(-ts/(es*dsu0)-tu/(2*eu*duu0)-ts*tu/(2*es*dsu0))... /(es*dsu0+ts*eu*duu0); pout0=tem5+tem6; end pout_A=pout0*((1-pd0)^N); % pout=pout_A; %Decode for i=1:N y5=fun1(eu,er,es,ru,rs,duu0,dru0,dsu0,i); y6=fun2(eu,er,es,ru,rs,duu0,dru0,dsu0,i); y7=fun3(eu,er,es,ru,rs,duu0,dru0,dsu0,i); y8=fun4(eu,er,es,ru,rs,duu0,dru0,dsu0,i); pout_A1=nchoosek(N,i)*(y5+y6+y7+y8)*(pd0^i)*((1-pd0)^(N-i)); pout_A=pout_A+pout_A1; end pout_A; pout=pout_A; end end