f_s_k_p.m

《移动衰落信道》Mobiel_Fading_Channels一书后面的相关仿真程序。

%-------------------------------------------------------------------- 
% F_S_K_p.m --------------------------------------------------------- 
% 
% Program for the generation of the matrices used in F_S_K.m. 
% 
% Used m-file: pCOST207.m 
%-------------------------------------------------------------------- 
% [C1,F1,TH1,C2,F2,TH2,F01,F02,RHO,F_RHO,q_l,T]= 
%                      F_S_K_p(N_1,AREA,f_max) 
%-------------------------------------------------------------------- 
% Explanation of the input parameters: 
% 
% N_1: minimum number of discrete Doppler frequencies 
% AREA: according to COST 207, 4 types of channels are specified: 
%            1) Rural Area:    'ra' 
%            2) Typical Urban: 'tu' 
%            3) Bad Urban:     'bu' 
%            4) Hilly Terrain: 'ht' 
% f_max: maximum Doppler frequency 
 
function [C1,F1,TH1,C2,F2,TH2,F01,F02,RHO,F_RHO,q_l,T]=... 
         F_S_K_p(N_1,AREA,f_max) 
 
% The greatest common divisor of the discrete propagation delays 
% defines the sampling interval T_s: 
T_s=0.2E-6; 
 
if     all(lower(AREA)=='ra'), 
       a_l=[1,0.63,0.1,0.01]; 
       tau_l=[0,0.2,0.4,0.6]*1E-6; 
       DOPP_KAT=['RI';'JA';'JA';'JA']; 
elseif all(lower(AREA)=='tu'), 
       a_l=[0.5,1,0.63,0.25,0.16,0.1]; 
       tau_l=[0,0.2,0.6,1.6,2.4,5]*1E-6; 
       DOPP_KAT=['JA';'JA';'G1';'G1';'G2';'G2']; 
elseif all(lower(AREA)=='bu'), 
       a_l=[0.5,1,0.5,0.32,0.63,0.4]; 
       tau_l=[0,0.4,1.0,1.6,5.0,6.6]*1E-6; 
       DOPP_KAT=['JA';'JA';'G1';'G1';'G2';'G2']; 
elseif all(lower(AREA)=='ht'), 
       a_l=[1,0.63,0.4,0.2,0.25,0.06]; 
       tau_l=[0,0.2,0.4,0.6,15,17.2]*1E-6; 
       DOPP_KAT=['JA';'JA';'JA';'JA';'G2';'G2']; 
end 
 
% Generate the parameters and assign them to the matrices: 
num_of_taps=length(DOPP_KAT); 
F1=zeros(num_of_taps,N_1+2*num_of_taps-1); 
F2=F1;C1=F1;C2=F1;TH1=F1;TH2=F1; 
F01=zeros(1,num_of_taps);F02=F01; 
RHO=zeros(1,num_of_taps);F_RHO=RHO; 
NN1=N_1+2*(num_of_taps-1):-2:N_1; 
for k=1:num_of_taps, 
    [f1,f2,c1,c2,th1,th2,rho,f_rho,f01,f02]=... 
    pCOST207(DOPP_KAT(k,:),NN1(k)); 
    F1(k,1:NN1(k))=f1; 
    C1(k,1:NN1(k))=c1*sqrt(a_l(k)); 
    TH1(k,1:NN1(k))=th1; 
    F2(k,1:NN1(k)+1)=f2; 
    C2(k,1:NN1(k)+1)=c2*sqrt(a_l(k)); 
    TH2(k,1:NN1(k)+1)=th2; 
    F01(k)=f01;F02(k)=f02; 
    RHO(k)=rho;F_RHO(k)=f_rho; 
end 
 
% Determine indices of the delay elements of the FIR filter: 
q_l=tau_l/T_s+1; 
 
% Initialization of the delay elements of the FIR filter: 
T=zeros(1,max(q_l));