瑞利衰落信道

[那年下雪]

我利用多普勒滤波器Sf1来产生瑞利衰落（网上找的），具体步骤如下。
首先产生多普勒滤波器的频响：
f=1:2*fm-1; %通频带长度
y=0.5./((1-((f-fm)/fm).^2).^(1/2))/pi; %多普勒功率谱(基带)
Sf=zeros(1,len_data);
Sf1=y;%多普勒滤波器的频响（通频带长度限制？？？？）
产生复高斯信号,并与多普勒滤波器的开方作用：
x1=randn(1,len_data);
x2=randn(1,len_data);
nc1=conv(x1+i*x2,Sft); %同相分量
nc=nc1(1:len_data);

x3=randn(1,len_data);
x4=randn(1,len_data);
ns1=conv(x3+i*x4,Sft); %正交分量
ns=ns1(1:len_data);
得到瑞利信号：
r0=(real(nc)+j*real(ns)); %瑞利信号
r=abs(r0); %瑞利信号幅值
但是我想加到我的基带仿真中，只知道和数据点乘即可，但怎样和数据匹配，才能得到真正的瑞利衰落信号。
如：BPSK，码元速率sr，一次仿真数据长度是100；过采样点数8；
成型滤波器长度21*8；如此要送入瑞利信道的数据长度是800+168-1；
而
谢谢啦！！！！！！！！