minimax_jakes.m

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

N=100;f_max=91;sigma_0_2=1;Ni=25;PLOT=1;
% Ni=[5 10 15 20 25 30 35 40 50];
Euu1=zeros(1,length(Ni));
Euu2=zeros(1,length(Ni));
for ii=1:length(Ni)
    N_i=Ni(ii);
    tau=linspace(0,0.2,N);
    Jo=sigma_0_2*besselj(0,2*pi*f_max*tau);
    c_i_n=sqrt(sigma_0_2)*sqrt(2/N_i)*ones(N_i,1);
    f_i_n=zeros(N_i,1);
    f_i_n1=zeros(N_i,1);
    save  data Jo tau N_i c_i_n  PLOT
    tic
%     Initial values:
   for i=1:N_i
       f_i_n(i)=f_max*sin((pi/(2*N_i))*(i-1/2));
   end
%    [f_i_n0,dummy1,dummy2]=parameter_Jakes('es_j',N_i,...
%                       sqrt(sigma_0_2),f_max,'rand',0);
   xo=f_i_n;
   options = optimset('MaxIter',5);
   [x,fval]=fminimax('fun_Jakes1',xo,[],[],[],[],[],[],[],options);
   load x
   f_i_n=x;

 
%     o=foptions;
%     o(1)=1;
%     o(1)=0;
%     o(2)=1e-9;
%     o(14)=N_i/10*200;
%     o(9)=0;
% 
%     xo1=f_i_n0;
%     yo1=c_i_n;
%     w1o=0.25*ones(1,length(tau));
%     w2o=0.75*ones(1,length(tau));
% 
%     xx=fmincon('fun_Jakes',xo1,yo1,w1o,w2o,o,'grad_Jakes');
% 
%     load xx
% 
%     f_i_n1=xx;
%     
%     r1=acf_mue(f_i_n,c_i_n,tau);
%     r2=acf_mue(f_i_n1,c_i_n,tau);
%     Euu1(ii)=mse(abs(Jo-r1).^2);
%     Euu2(ii)=mse(abs(Jo-r2).^2);
    
    toc
end
% plot(Ni,Euu1,'*',Ni,Euu2,'o');
% 
%  xlabel('Ni');
%    ylabel('Erui');
% %    legend('minimax','MEA','MEDS');
% legend('minimax','Lp-norm');