lcrgauss.m

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

% Ni=[1 2 3 4 5 10 15 20 25 30 35 40 45 50];
% Euu=zeros(1,length(Ni));
% Euu1=zeros(1,length(Ni));
% Euu2=zeros(1,length(Ni));
% for i=1:length(Ni)
    N=100;f_max=91;sigma_0_2=1;N1=5;N2=6;PLOT=0;
    f_c=sqrt(log(2))*f_max;kappa_c=f_max/f_c;
    F_list=[];
    save F_list F_list
    tau_max1=N1/(2*kappa_c*f_c);
    tau1=linspace(0,tau_max1,N);
    r_mm1=sigma_0_2*exp(-(pi*f_c/sqrt(log(2))*tau1).^2);
    f_i_n1=zeros(N1,1);f_i_n2=zeros(N2,1);
    c_i_n1=sqrt(sigma_0_2)*sqrt(2/N1)*ones(N1,1);
    c_i_n2=sqrt(sigma_0_2)*sqrt(2/N2)*ones(N2,1);
    for ii=1:N1
        f_i_n1(ii)=f_c/sqrt(log(2))*erfinv((2*ii-1)/(2*N1));
    end 
    N_i=N1;r_mm=r_mm1;tau=tau1;c_i_n=c_i_n1;
    save data r_mm tau N_i c_i_n  PLOT
    xo=f_i_n1;
    options = optimset('MaxIter',10);
    [x,fval]=fminimax('fun_Gauss1',xo,[],[],[],[],[],[],[],options);
    load x
    f_i_n1=sort(abs(x));
    F_list=[];
    save F_list F_list
    tau_max2=N2/(2*kappa_c*f_c);
    tau2=linspace(0,tau_max2,N);
    r_mm2=sigma_0_2*exp(-(pi*f_c/sqrt(log(2))*tau2).^2);
     for ii=1:N2
        f_i_n2(ii)=f_c/sqrt(log(2))*erfinv((2*ii-1)/(2*N2));
     end 
    N_i=N2;r_mm=r_mm2;tau=tau2;c_i_n=c_i_n2;
    save data r_mm tau N_i c_i_n  PLOT
    xo=f_i_n2;
    options = optimset('MaxIter',10);
    [x,fval]=fminimax('fun_Gauss1',xo,[],[],[],[],[],[],[],options);
    load x
    f_i_n2=sort(abs(x));
    th1=rand(N1,1)*2*pi;th2=rand(N2,1)*2*pi;
    T_s=1/30*max(f_i_n1);K=1e6;T_sim=K*T_s;
    xi_t=Rice_proc(f_i_n1,c_i_n1,th1,f_i_n2,c_i_n2,th2,0,0,0,...
                        T_s,T_sim,0);
    r=0:9;rou=0;
    [p_z,r]=pdf(sigma_0_2,sigma_0_2,r,rou,f_max,0);
    beita=2*(pi*f_c*sqrt(sigma_0_2))^2/log(2);
    N_realr=sqrt(beita/(2*pi))*p_z;
    N_realr=N_realr/max(N_realr);
    plot(r,N_realr);
    hold on;
    beitah=2*pi^2*sum((c_i_n1.*f_i_n1).^2);
    N_r=sqrt(beitah/(2*pi))*p_z/f_max;
    N_r=N_r/max(N_r);
    plot(r,N_r,'*');
    xlabel('r')
    ylabel('N(r)/fmax')  
    legend('reference model','simulation model');