msejake.m

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

Ni=[1 2 3 4 5 10 15 20 25 30 35 40 45 50];
Euu1=zeros(1,length(Ni));
Euu2=zeros(1,length(Ni));
for i=1:length(Ni)
    N_i=Ni(i);
    N=100;f_max=91;sigma_0_2=1;PLOT=0;
    tau=linspace(0,N_i/(2*f_max),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);
    f_i_n2=zeros(N_i,1);
    save  data Jo tau N_i c_i_n  PLOT  

    % Initial values:
    for ii=1:N_i
        f_i_n(ii)=f_max*sin((pi/(2*N_i))*(ii-1/2));
        f_i_n1(ii)=f_max*sin(pi*ii/(2*N_i)); %MEA
        
    end
    f_i_n2=f_i_n;
%     options = optimset('GradObj','on');
    xo=f_i_n;
    x=fminimax('fun_Jakes1',xo);
    load x
    f_i_n=x;
    r=acf_mue(f_i_n,c_i_n,tau);
    r1=acf_mue(f_i_n1,c_i_n,tau);
    r2=acf_mue(f_i_n2,c_i_n,tau);
    Euu1(i)=mse(abs(Jo-r).^2);
    Euu2(i)=mse(abs(Jo-r1).^2);
%     Euu3(i)=mse(abs(Jo-r2).^2);

end
%  plot(1:length(Euu1),Euu1,'*',1:length(Euu2),Euu2,'+',1:length(Euu3),Euu3,'o'); 
% plot(1:length(Euu1),Euu1,'*',1:length(Euu2),Euu2,'+');
plot(Ni,Euu1,'*',Ni,Euu2,'o');

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