minimax_gauss1.m

无线移动信道的matlab源码 用于对无线通信信道的编程实现

N=100;f_max=91;sigma_0_2=1;N_i=[5 10 15 20 25 30 35 40 50 60 80];PLOT=0;
beita=zeros(1,length(N_i));beitah=zeros(1,length(N_i));delta=zeros(1,length(N_i));
beita1=zeros(1,length(N_i));beitah1=zeros(1,length(N_i));delta1=zeros(1,length(N_i));
for i=1:length(N_i)
    Ni=N_i(i);
    f_c=sqrt(log(2))*f_max;
    kappa_c=f_max/f_c;
    F_list=[];
    save F_list F_list
    tau_max=Ni/(2*kappa_c*f_c);
    tau=linspace(0,0.1,N);
    r_mm=sigma_0_2*exp(-(pi*f_c/sqrt(log(2))*tau).^2);
    f_i_n=zeros(Ni,1);f_i_n1=zeros(Ni,1);
    c_i_n=sqrt(sigma_0_2)*sqrt(2/Ni)*ones(Ni,1);
    for ii=1:Ni
        f_i_n(ii)=f_c/sqrt(log(2))*erfinv((2*ii-1)/(2*Ni));
    end
    save data r_mm tau N_i c_i_n tau_max PLOT
    xo=f_i_n;
    options = optimset('MaxIter',10);
    [x,fval]=fminimax('fun_Gauss1',xo,[],[],[],[],[],[],[],options);
    load x
    f_i_n=sort(abs(x));

    for ii=1:Ni
        f_i_n1(ii)=f_c/sqrt(log(2))*erfinv(0.9999999); %MEA
    end
   beitah(i)=2*pi^2*sum((c_i_n.*f_i_n).^2);
   beita(i)=2*(pi*f_c*sqrt(sigma_0_2))^2/log(2);
   delta(i)=abs(beita(i)-beitah(i))/beita(i);
   beitah1(i)=2*pi^2*sum((c_i_n.*f_i_n1).^2);
   beita1(i)=2*(pi*f_c*sqrt(sigma_0_2))^2/log(2);
   delta1(i)=abs(beita1(i)-beitah1(i))/beita1(i);
end
plot(N_i,delta,'*',N_i,delta1,'o');

 xlabel('Ni');
%    ylabel('');
%    legend('minimax','MEA','MEDS');
legend('minimax','MEA');