c.m

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

clear all;
N_1=25;
N_3=15;
sigma_0_2=0.7697;
kappa_0=0.4045;
theta_0=164*pi/180;
f_max=200;
sigma_3=0.0062;
m_3=-0.3861;
rho=1.567;
theta_rho=127*pi/180;
f_c=sqrt(log(2))*f_max;
T_s=1/144e+3;
T_sim=5;
PLOT=0;
eta_t=Suzuki_Type_II(N_1,N_3,sigma_0_2,kappa_0,theta_0,...
                            f_max,sigma_3,m_3,rho,theta_rho,f_c,...
                            T_s,T_sim,PLOT);
mean(eta_t)
r=0.8;
a=zeros(1,length(eta_t)/2);
b=zeros(1,length(eta_t)/2);
c=1;
d=1;
for i=1:length(eta_t)
    if eta_t(i)<=r
        a(c)=i;
        c=c+1;
         
    else b(d)=i;
        d=d+1;
    end
end
a1=diff(a);
b1=diff(b);
a2=find(a1>1);
b3=a1(a2)-1;             %error-free recorder
b2=find(b1>1);
a3=b1(b2)-1;             %error recorder
pe=[];
te=min(a3):max(a3);
for k=1:length(te)
    pe(k)=length(find(a3==k));
     p_re(k)=sum(pe)/length(a3);
end    
 p_red=pe/length(a3);
 al=find(p_red>0.99);
 pe(al)
%  plot(te,p_re)
 pf=[];
 tf=min(b3):max(b3);
for k=1:length(tf)
    pf(k)=length(find(b3==k));
     p_rf(k)=sum(pf)/length(eta_t);
end    
 p_rfd=pf/length(eta_t);
  bl=find(p_rfd>0.99);
 pf(b1)


%  subplot(211),plot(te,p_re);
%  subplot(212),plot(tf,p_rf);
% N_1=10;
% N_2=11;
% p=[0.02 0.2 0.4 0.8 1];
% for i=1:1
% sigma_1_2=0.069417;
% sigma_2_2=0.069417;
% % sigma_1_2=1;
% % sigma_2_2=p(i)*sigma_1_2;
% % sigma_1_2=0.10391;
% % sigma_2_2=0.030488;
% 
% pp=sigma_1_2/sigma_2_2;
% kappa_0=sqrt(pp*(2-pp));
% theta_0=pi/2;
% f_max1=23.19;
% f_max2=f_max1;
% f_max=f_max1;
% BEITA1=1251.96;
% % BEITA1=1103.4298;
% % BEITA2=1091.5206;
% % f_max2=sqrt((sigma_1_2*beita2)/(sigma_2_2*BEITA1))*f_max1;
% % f_max2=42.58;
% rho=0;
% theta_rho=0;
% 
% T_s=0.0005;
% T_sim=0.5;
% PLOT=1;
% % [xi_t,thita]=gen_Nakagami_proc(N_1,N_2,sigma_1_2,sigma_2_2,kappa_0,...
% %                             theta_0,rho,theta_rho,f_max1,f_max2,T_s,...
% %                             T_sim,0);    
% xi_t=gen_Rice_proc(N_1,N_2,sigma_1_2,sigma_2_2,kappa_0,...
%                             theta_0,0,0,f_max,T_s,...
%                             T_sim,0);
% 
% % z=power(xi_t,xi_t);
% % plot(xi_t)
% % hist(xi_t,z)
% % p_z=pdf(sigma_1_2,sigma_2_2,z,PLOT); 
% % p_thita=pdf_phase(sigma_1_2,sigma_2_2,thita,PLOT);
% r=-30:5:5;
% % N_r=lcr_sim(20*log10(xi_t),r,f_max,PLOT);
% N_r=lcr_simch(sigma_1_2,BEITA1,r,f_max,PLOT);
% hold on;
% ee=[0.3168 0.5507  0.9115 1.2987 1.2046 0.4052 0.0140 3.5546e-007 4.1924e-021];
% plot(-30:5:10,ee)
% 
% end