minimax_jakes1.m

来自「无线移动信道的matlab源码 用于对无线通信信道的编程实现」· M 代码 · 共 38 行

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)
   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(i))*ones(N_i(i),1);
   f_i_n=zeros(N_i(i),1);f_i_n1=zeros(N_i(i),1);
   Ni=N_i(i);
   save  data Jo tau Ni c_i_n  PLOT
   % Initial values:
   for ii=1:N_i(i)
       f_i_n(ii)=f_max*sin((pi/(2*N_i(i)))*(ii-1/2));
   end
   xo=f_i_n;
   options = optimset('MaxIter',6);
   [x,fval]=fminimax('fun_Jakes1',xo,[],[],[],[],[],[],[],options);
   load x
   f_i_n=x;
    for ii=1:Ni
        f_i_n1(ii)=f_max*sin(pi*ii/(2*Ni)); %MEA
    end
   beitah(i)=2*pi^2*sum((c_i_n.*f_i_n).^2);
   beita(i)=2*(pi*f_max*sqrt(sigma_0_2))^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_max*sqrt(sigma_0_2))^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');