麻烦高手给我看下协作通信的程序，急，谢谢了

wptx032

%this code is for fixed DF performance evaluation, compare the diffirent SR distance
%该仿真用来评估不同的SR距离下的固定DF转发性能

%for S-R-D three nodes
%the SD and RD channel is Rayleigh fading, and the average SNR is varied.

clear all;

%the times of Monte Carlo
%蒙特卡罗仿真的次数
Monte_MAX=10^(7);

%the length of data frame
M2=10^(0);

MIN_SNR=0;
MAX_SNR=30;
INTERVAL=2;

%功率划分因子，在协同情况下，为了保证总的功率一定，每个节点使用1/2的功率发送
POW_DIV=1/2;

%open a new record file
%打开一个新的文件，如果文件已经存在，那么清除以前的内容
fid=fopen('record.txt','w');
fprintf(fid,'%%%s\n',datestr(now));
fclose(fid);

%set the SD distance and RD distance
%设置SD距离和RD的距离，归一化为1
sd_distance=1;
rd_distance=1;
%变化不同的SR距离
sr_distances=[0.1,0.5,0.7,0.9];

%使用循环的方式仿真不同的SR距离性能
for sr_INDEX=1:2
   
    %通过sr_INDEX变量获得真正的SR距离
    sr_distance=sr_distances(sr_INDEX);
     
    %snrcount这个变量用来对下面的SNR循环进行计数
    snrcount=0;
   
    %变化SNR循环
    for SNR=MIN_SNR:INTERVAL:MAX_SNR
        sig=10^(SNR/10);
        
        %just display something on the screen.
        %仅仅是在计算机上显示一些信息，因为仿真往往很漫长，只是想知道仿真到哪里了^.^
        datestr(now)
        M2
        SNR
        
        snrcount=snrcount+1;
        
        %用来对错误的bit计数
        err_num_dsd=0;
        err_num_coop=0;
        
        for tries=0:Monte_MAX
            
            %flag indicate if the cooperate is used, if 1, use cooperation,else direct tran is perferred.
            %一个标志，表明是否转发
            tx_coop=1;
            
            %generate the source BPSK signal, 0,1
            %产生一个长度为M2的0，1随机数
            X1=rand(1,M2)>0.5;
            %generate the source BPSK signal,+1,-1
            %产生一个长度为M2的BPSK随机信号
            Xs=X1*2-1;

            %generate the source-relay channel, the channel keep constant in a frame
            %产生SR信道，假设信道在一帧内保持不变
            CH_sr=xy_RayleighCH(1)/(sr_distance)^2;
            
            %the received signal of source-relay transmission
            %中继接收到的源的信号，y=hx+n
            R_dsr=CH_sr*sqrt(POW_DIV*sig).*Xs+xy_noise(M2);
            %the decode signal of source-relay transmission
            %中继对信道进行估计，匹配后硬判决得到的源的信号，根据MRC（最大比合并）的基本原理
            Xr=(conj(CH_sr)*sqrt(POW_DIV*sig).*R_dsr>0)*2-1;
            
            %如果中继解码错误，则不转发
            if (sum(Xs~=Xr)>0) tx_coop=0; end
            %如果使用固定DF，我们强制设置tx_coop为1
            tx_coop=1;
            
            %generate the source-destination channel, the channel keep constant in a frame
            %产生SD信道，假设信道在一帧内保持不变
            CH_sd=xy_RayleighCH(1)/(sd_distance)^2;  

            %the received signal of direct source-destination transmission
            %假设不使用协同，直接传输情况下目的节点接收到的源的信号，y=hx+n，这里，源节点使用全功率发送
            R_dsd=CH_sd*sqrt(sig).*Xs+xy_noise(M2);
            %the decode signal of direct source-destination transmission
            %目的节点对信道进行估计，匹配后硬判决得到的源的信号，根据MRC（最大比合并）的基本原理
            Y_dsd=(conj(CH_sd)*sqrt(sig).*R_dsd>0)*2-1;


            %the received signal of cooperative source-destination transmission
            %在使用协同情况下，目的节点接收到的源的信号，y=hx+n，这里，源节点使用半功率发送
            R_csd=CH_sd*sqrt(POW_DIV*sig).*Xs+xy_noise(M2);

            %generate the relay-destination channel, the channel keep constant in a frame
            %产生RD信道，假设信道在一帧内保持不变
            CH_rd=xy_RayleighCH(1)/(rd_distance)^2;
            
           
            if (tx_coop==1)
            %the received signal of cooperative relay-destination transmission
            %在使用协同情况下，目的节点接收到的来自中继的信号
            R_crd=CH_rd*sqrt(POW_DIV*sig).*Xr+xy_noise(M2);
            %MRC的合并信号
            R_combine=conj(CH_sd)*sqrt(POW_DIV*sig).*R_csd+conj(CH_rd)*sqrt(POW_DIV*sig).*R_crd;
            else
            %if direct transmission is chosen, the source just retransmit the signal again.
            %如果不使用协同，源节点重复发送自己的数据，假设两次重复发送期间信道保持不变
            R_crd=CH_sd*sqrt(POW_DIV*sig).*Xs+xy_noise(M2);
            %MRC的合并信号
            R_combine=conj(CH_sd)*sqrt(POW_DIV*sig).*R_csd+conj(CH_sd)*sqrt(POW_DIV*sig).*R_crd;
            end
                       
            %the decode signal of of cooperative MRC combining
            %MRC后的判决信号
            Y_combine=(R_combine>0)*2-1;

            %the number of error bit
            %统计每一帧里面错误的比特数目
            err_num_dsd=sum(Xs~=Y_dsd)+err_num_dsd;
            err_num_coop=sum(Xs~=Y_combine)+err_num_coop;

          end %try=0:Monte_MAX

            %the real ber
            %计算每个SNR下的平均错误BER
            ber_dsd(snrcount)=err_num_dsd/(M2*Monte_MAX);  
            ber_coop(snrcount)=err_num_coop/(M2*Monte_MAX);
           

    end % SNR=MIN_SNR:MAX_SNR  
     
     %以追加的方式打开文件，将数据写入记录文件
     fid=fopen('record.txt','a');
     fprintf(fid,'\n%%ber_dsd is\nber_dsd%d=[ ',sr_INDEX);
     fprintf(fid,'%d ',ber_dsd);
     fprintf(fid,'];\n');
     fprintf(fid,'\n%%fixed ber_coop is\nber_coop%d=[ ',sr_INDEX);
     fprintf(fid,'%d ',ber_coop);
     fprintf(fid,'];\n');
     fclose(fid);
     
end %Sr index



SNR1=MIN_SNR:INTERVAL:MAX_SNR;

%瑞利信道下的直接发送的理论BER表达式
sig1=10.^(SNR1/10);
u1=sqrt(sig1./(1+sig1));
Pe10=1/2*(1-u1);

%瑞利信道下的二阶分集理论BER表达式，使用半功率
sig2=1/2*10.^(SNR1/10);
u2=sqrt(sig2./(1+sig2));
Pe20=1/4*(1-u2).^(2).*(2+u2);

%使用半对数画出SNR和BER的曲线
semilogy(SNR1,ber_dsd,'-o',SNR1,ber_coop,'-+',SNR1,Pe10,'-+',SNR1,Pe20,'-*');
legend('direct','cooperative','L=1, therotic','L=2 therotic');
grid on;
ylabel('The AVERAGE BER');
xlabel('SNR(dB)');
title('the cooperative BER of BPSK in Rayleigh fading channel');
axis([MIN_SNR,MAX_SNR,10^(-6),1]);
我想问的就是我现在要把多径衰落的程序加进去，好像是matlab自带的berfading文件吧，但我不知道应该怎么用那个berfading文件，里面的几个参数也不知道该取什么值合适，还有这个berfading程序应该放在我上面程序中的哪个位置啊，那位高手帮我解答下啊，本人要写论文，很急。谢谢了