map_bound.asv

短波信道抗多音干扰的性能分析及其仿真

function Pb=map_bound(BPH,number_of_states,rho_in_db)
%this function computes the MAP demodulation
%bit error bound

fanout=2^BPH;
rho=10^(rho_in_db/10);
beta=BPH*rho/2;
namda=(sqrt(1+6*beta+beta^2)-(1+beta))/2;
L=floor(log(number_of_states)/log(fanout));
Pb=(fanout/4)*(exp(-2*beta*(namda/(1+namda)))/(1-namda^2))^L;



%s2=BPH*n;
%sgma2=BPH/(2*u);
%Q=quad(@myfun,0,25,[],[],n,s2,sgma2);
%P=1-Q;
%Pb=(fanout^n/2)*P;
%g=inline('-exp(-(x./(1+x)).*rho)./(1-x.^2)','x');
%[namda,fval]=fminbnd(@myfun1,0,0.99);
%Pb=(M^n/4)*fval^n;
%g=inline('exp(namda.*x).*(1/gamma(n./2)).*x.^((n./2)-1).*exp(-x)','x','n','namda');
%h=inline('exp(-namda.*x).*(1/2).*((x./(n.*u)).^((n-2)./4)).*exp((-1/2).*(x+n.*u)).*besseli(n./2-1,sqrt(n.*u.*x))','x','n','u','namda');
%for i=0:1:99
    %Q(i+1)=quad(g,0,100,1e-3,[],n1,i/100);
    %J(i+1)=quad(h,0,100,1e-3,[],n1,u1,i/100);
    %end
%P=Q.*J;
%y=min(P);
%Pb=(fanout^n/2)*y;