1.算法描述

在光通信领域，多址技术主要有WDM(波分复用)、TDM(时分复用)及OCDMA(光码分多址)三种方式。OCDMA技术从70年代中期开始出现，现在引起了人们的广泛重视。将CDMA(码分多址)技术用于光通信中，不仅保持了CDMA技术在无线电中的抗干扰、保密、软容量及网络协议简单等特性，而且可以利用光载体的巨大带宽。

因此，CDMA技术和光通信的结合存在着优势互补。本文围绕OCDMA系统及其关键技术展开研究工作。在本文中，我们对OCDMA的关键技术进行了理论分析，提出了几种新的地址码，提出了一些新方案来改善系统性能。同时，开展了时域非相干OCDMA的实验工作。

编写无线光通信-CDMA中闭环链路的功率控制误差（Power control error）程序，其定义参考文献An Analytical Approach for Closed-Loop Power Control Error Estimations in CDMA Cellular Systems中详细说明了。仿真出功率控制误差（Power control error）在OOK 2PPM 4PPM 8PPM 16PPM的不同调制方式下；以及功率控制误差（Power control error）在1/4code rate, 1/2 code rate, 3/4 code rate不同码率下与BER的关系。通过给每个用户指定互相关为0的两个地址码分别发送1和0信号以及在接收端采用差分检测技术，避免了OOK(开关键控)光通信系统中0不发送光脉冲和OOK常规检测方式下最佳判决门限必须根椐用户数等外界环境而改变，从而提高了系统性能和降低了系统的复杂度。通过引入OPPM(重叠脉冲脉位调制)方式，进一步提高了系统的传输效率。

这里我们首先建立一个cdma小区通信环境，然后根据光通信代码在CDMA的基础上，使用光通信进行通信，然后在此基础上加入功率工作，不同的调制方式以及不同的码率等进行仿真。

加入功率控制PCE模块，其基本结构如下所示：

对OOK，2PPM，4PPM，8PPM，16PPM五种方法进行仿真。

光CDMA方式是基于OTDM和WDM之上的一种多路存取方式，它可以使相同带宽和比特率的光信号携带更多的信息。在光纤CDMA（OCDMA)中，每一位数据都被一个序列编码，每个用户都有一个单独的序列，在发送端，对要传输的数据的序列码进行光正交编码，然后实现多个用户共享同一光纤信道；在接收端，用与发送端相同的地址码进行光正交解码，恢复原用户数据。
CDMA技术的原理是基于扩频技术，即将需传送的具有一定信号带宽信息数据，用一个带宽远大于信号带宽的高速伪随机码进行调制，使原数据信号的带宽被扩展，再经载波调制并发送出去。CDMA使用带扩频技术的模数转换（ADC).输入音频首先数字化为二进制元。传输信号频率按指定类型编码，因此只有频率响应编码一致的接收机才能拦截信号。

2.仿真效果预览

matlab2022a仿真结果如下：

3.MATLAB部分代码预览

PACKET_LENGTH  = 1000;  %Packet size 
PEEK_POWER     = 1;
SYMBOL_LENGTH  = 36;   %here symbol length= packet length 
SLOT_LENGTH    = 4;  
TRANS_LENGTH   = 9*SLOT_LENGTH*PACKET_LENGTH;
pulse_drtn     = 4;
CHANNELS       = 3;
ELECTRON       = 1.6e-19;
INTERFERENCE   = 1e-3;
RECEIVER_RES   = 0.53;
COEF_RED       = 0.1550;
COEF_GREEN     = 1.7309;
COEF_BLUE      = 1.1142;
RES_RED        = 0.5200;
RES_GREEN      = 0.4800;
RES_BLUE       = 0.4000;
%Derived parameters 
bit_rate       =  16e8;                          % bit rate [100Mbps]  
rms_delay      =  2e-9;                         % delay spread (1 ns)                                                      
multipath      =  3;                            % multi-path dispersion
sample_rate    =  bit_rate / 3.0 * SLOT_LENGTH;
%struct built
for i=1:CHANNELS
    sentbit(i,:)  = zeros(1,3*PACKET_LENGTH);
    signal(i,:)   = zeros(1,12*PACKET_LENGTH);
    recivbit(i,:) = zeros(1,3*PACKET_LENGTH);
end
 
dwTransmitPower           = 7;        % 发送功率
dShadowFadingStd          = 4;        % 阴影衰落标准差BS2UE dB
dUserVelocity             = 0;        % 用户移动速度0km/h，这里假设是静止的
 
MAX_Dis                   = 7;
STEP                      = 0.2;
MIN_Dis                   = 1;
 
DIS                       = [MIN_Dis:STEP:MAX_Dis];%定义用户和基站之间的距离
MTKL                      = 20;%多次仿真，计算平均值
Avgtime                   = 200;
Delay                     = 5;
 
startpow                  = -70;    
endpow                    = -40;  
powstep                   =  5;     
%**************************************************************************
 
%%
%主题函数，通过for循环，仿真不同距离下的误码率曲线
%为了对比，我们设置三组不同的SNR下的AWGN进行仿真
%由于原先的程序是基于数据组进行仿真，而基于PCE估计的算法，需要基于时间流模型，所以必须将
%仿真的方法修改下，即对每一个数据位单独的循环仿真，从而模拟出实时的效果
 
E1 = zeros(PACKET_LENGTH,1);
E2 = zeros(PACKET_LENGTH,1);
E3 = zeros(PACKET_LENGTH,1);
cnt = 0;
for jj = 0.1:0.05:0.5   
    cnt = cnt + 1;
    for mm = 1:MTKL
        jj
        mm
        %根据路径距离计算路径损耗
        PASS_LOSS = -25;
        %计算光三组元色对应的SNR值
        snr_red   = 10*log10(power((power(10,PASS_LOSS/10)*1e-3)*(COEF_RED  /(COEF_RED+COEF_GREEN+COEF_BLUE)),2)/((ELECTRON*INTERFERENCE/RES_RED)  *(bit_rate*(SLOT_LENGTH/pulse_drtn)/3.0)));
        snr_green = 10*log10(power((power(10,PASS_LOSS/10)*1e-3)*(COEF_GREEN/(COEF_RED+COEF_GREEN+COEF_BLUE)),2)/((ELECTRON*INTERFERENCE/RES_GREEN)*(bit_rate*(SLOT_LENGTH/pulse_drtn)/3.0)));
        snr_blue  = 10*log10(power((power(10,PASS_LOSS/10)*1e-3)*(COEF_BLUE /(COEF_RED+COEF_GREEN+COEF_BLUE)),2)/((ELECTRON*INTERFERENCE/RES_BLUE) *(bit_rate*(SLOT_LENGTH/pulse_drtn)/3.0)));
        VarR      = calc_snr(snr_red);
        VarG      = calc_snr(snr_green);
        VarB      = calc_snr(snr_blue);
 
 
        %产生随机数作为发送数据
        for i=1:CHANNELS
            sentbit(i,:) = round(rand(1,3*PACKET_LENGTH)); 
        end
 
 
        for times = 1:3*PACKET_LENGTH%每一时刻，发送一个数据包
            S1 = 1*(times-1) + 1;
            F1 = 1*times;
            S2 = 4*(times-1) + 1;
            F2 = 4*(times);
            
            %**********************************************************************
            %根据上一次反馈得到的E进行功率的调整
            %**********************************************************************
            %对E取倒数
            %对应公式Pt(t) = Pt(t-D)/E(t-D) = 1/r(t-D)
            Coff(1,times) = 1/jj;% 之后，根据前一次的E进行调整发送功率
            Coff(2,times) = 1/jj;% 之后，根据前一次的E进行调整发送功率
            Coff(3,times) = 1/jj;% 之后，根据前一次的E进行调整发送功率  
            
            %调制,OOK 
            for i=1:CHANNELS
                signal(i,S2:F2) = func_1in4(sentbit(i,S1:F1),Coff(i,times));
            end
 
            %加多径
            if multipath~=1 
               for i=1:CHANNELS
                   signa2(i,S2:F2) = add_multipath_dispersion1(signal(i,S2:F2),4,sample_rate,rms_delay);
               end
            end
 
            
            %加噪声
            signal(1,S2:F2)=add_gauss_noise(signal(1,S2:F2),length(S2:F2), VarR);      
            signal(2,S2:F2)=add_gauss_noise(signal(2,S2:F2),length(S2:F2), VarG);
            signal(3,S2:F2)=add_gauss_noise(signal(3,S2:F2),length(S2:F2), VarB); 
        
 
            %解调
            for i=1:CHANNELS
                recivbit(i,S1:F1) = func_de1in4(signal(i,S2:F2),Coff(i,times));
            end
        end
 
        for i=1:CHANNELS
            error(i,:) = size(find(recivbit(i,:)~=sentbit(i,:)),2);
        end
        ERR(mm) = (error(1,:)+error(2,:)+error(3,:))/(CHANNELS*3*PACKET_LENGTH);
    end
    ber(cnt) = mean(ERR);
end
01_057_m