mmsejake.m

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

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=200;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);
    p=2;
    save data Jo tau N_i c_i_n p PLOT

    % Initial values: 
    [f_i_n,dummy1,dummy2]=parameter_Jakes('lp_j',N_i,...
                      sqrt(sigma_0_2),f_max,'rand',0);
    o=foptions;
    o(1)=1;
    o(1)=0;
    o(2)=1e-9;
    o(14)=N_i/10*200;
    o(9)=0;

    xo=f_i_n;
    yo=c_i_n;
    w1o=0.25*ones(1,length(tau));
    w2o=0.75*ones(1,length(tau));

    x=fminu('fun_Jakes',xo,yo,w1o,w2o,o,'grad_Jakes');

    load x

    f_i_n=x;
end
 plot(1:length(Euu1),Euu1,'*',1:length(Euu2),Euu2,'+');