f_s_k_p.m

继续上传我自己遍的信道！谢谢！ 这个我没有完全的调试成功

%--------------------------------------------------------------------
% 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       离散Doppler频率的最小数
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']; %PDF形式 
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));