有关瑞利衰落信道的matlab编程，以及怎么改成莱斯分布呢？

Question

% 编写一个如瑞利多径的信道函数；
function y=Rayleigh_DuoJing(x,sigma,num_path,ray_factor,time_delay,p_path)

y=zeros(num_path,size(x,2));
for i=1:num_path
% 即使瑞利衰落因子是一样的，
% 对于每一行的衰落系数也要分别产生；
ray_number=raylrnd(ray_factor,size(x));
%ray_number=random('rician',0.5,ray_factor,size(x));
y(i,:)=sqrt(10^(p_path(i)/10))*(ray_number.*x);
temp=y(i,:);
y(i,:)=[zeros(1,time_delay(i)),temp(1:end-time_delay(i))];
randn('state',0);%设置初始状态
y(i,:)=y(i,:)+sigma*randn(size(x));
end;
y=sum(y,1);

    如上，p_path这个参数是什么呢？跟ray_factor有联系吗？
    sqrt(10^(p_path(i)/10))*(ray_number.*x);
    [zeros(1,time_delay(i)),temp(1:end-time_delay(i))];
    这两个行代码分别是什么意思呢？
    以及这段代码怎么改成莱斯分布呢？