function aRefSig = MBSFNRefSig(Ns,NDLRB,Nidmbsf,cMBSFN) % 將生成的MBSFN參考信號，映射至虛擬資源塊aRefSig；
% Ns 時(shí)隙號；
% NDLRB 下行鏈路占用的資源塊數(shù)；
% Nidmbsf  MBSFN（多播廣播單頻網(wǎng)絡(luò)）的號，用于PMCH；
aRefSig = zeros( 6,12*NDLRB);
for i = 1:6 % MBSFN參考信號只在天線端口4上傳輸，與PMCH共用一個(gè)天線端口。MBSFN參考信號只用于擴(kuò)展型循環(huán)前綴。
    Cinit(i) = 2 ^ 9 * (7 * (Ns + 1) + (i - 1) + 1) * (2 * Nidmbsf + 1) + Nidmbsf ; % 偽隨機(jī)序列要在每個(gè)OFDM符號的開頭初始化，Cinit(i) = 2 ^ 9 * (7 * (Ns + 1) + l + 1) * (2 * Nidmbsf + 1) + Nidmbsf;
end
cMBSFN = MBSFNpseudoranseq(Nidmbsf,Ns);% 調(diào)用偽隨機(jī)序列
NmaxDLRB = 110 ; % 協(xié)議規(guī)定下行鏈路中所占用的最大資源塊數(shù)是110；
c1 = cMBSFN(1, : ); % 將偽隨機(jī)序列分開，每個(gè)對應(yīng)著一個(gè)OFDM符號
c2 = cMBSFN(2, : );
c3 = cMBSFN(3, : );
c4 = cMBSFN(4, : );
c5 = cMBSFN(5, : );
c6 = cMBSFN(6, : );
for m = 1:(6 * NmaxDLRB) % MBSFN參考信號序列rlNs(m) = (1 / sqrt(2)) * (1 - 2 * c(2 * m)) + j * (1 / sqrt(2)) * (1 - 2 * c(2 * m + 1))
    r1Ns(m) = (1 / sqrt(2)) * (1 - 2 * c1(2 * m)) + j * (1 / sqrt(2)) * (1 - 2 * c1(2 * m + 1)); 
    r2Ns(m) = (1 / sqrt(2)) * (1 - 2 * c2(2 * m)) + j * (1 / sqrt(2)) * (1 - 2 * c2(2 * m + 1));
    r3Ns(m) = (1 / sqrt(2)) * (1 - 2 * c3(2 * m)) + j * (1 / sqrt(2)) * (1 - 2 * c3(2 * m + 1));
    r4Ns(m) = (1 / sqrt(2)) * (1 - 2 * c4(2 * m)) + j * (1 / sqrt(2)) * (1 - 2 * c4(2 * m + 1));
    r5Ns(m) = (1 / sqrt(2)) * (1 - 2 * c5(2 * m)) + j * (1 / sqrt(2)) * (1 - 2 * c5(2 * m + 1));
    r6Ns(m) = (1 / sqrt(2)) * (1 - 2 * c6(2 * m)) + j * (1 / sqrt(2)) * (1 - 2 * c6(2 * m + 1));
end
if mod(Ns,2) == 0;  %如果是偶時(shí)隙，參考信號占用的是第三個(gè)OFDM符號；
    a3 = zeros(1,12*NDLRB); % 先預(yù)定第三個(gè)OFDM符號的大小；
    for m = 1:6*NDLRB  % 參考信號在載波上從低層向高層開始映射，當(dāng)l=3時(shí)，映射到奇數(shù)載波上；
        k = 2 * m - 1 ;
        mm = m + 3 * (NmaxDLRB - NDLRB) ;  %mm是MBSFN參考信號的計(jì)數(shù)
        a3(k) = r3Ns(mm) ; % 將MBSFN參考信號映射至資源塊上；
    end
    aRefSig(3,:) = a3(1,:);
end
if mod(Ns,2) == 1 ; % 如果是奇時(shí)隙，參考信號占用的是第一個(gè)和第五個(gè)OFDM符號；
    a5 = zeros(1,12*NDLRB);
    a1 = zeros(1,12*NDLRB);
    for m = 1:6*NDLRB % 對于第五個(gè)OFDM符號，參考信號映射到奇數(shù)載波上；
        k = 2 * m - 1 ;
        mm = m + 3 * (NmaxDLRB - NDLRB) ; 
        a5(k) = r5Ns(mm) ;  % 將MBSFN參考信號映射至資源塊上；
    end
    for m = 1:6*NDLRB  % 對于第一個(gè)OFDM符號，參考信號映射到偶數(shù)載波上；
        k = 2 * m ;
        mm = m + 3 * (NmaxDLRB - NDLRB) ;
        a1(k) = r1Ns(mm) ;  % 將MBSFN參考信號映射至資源塊上；
    end
    aRefSig(5,:) = a5(1,:);
    aRefSig(1,:) = a1(1,:);
end
% 輸出的是MBSFN參考信號在資源塊上的映射規(guī)則。