DCSK的实用仿真和理论误码率分析

clear all;
beta=5;
n_symbol=5000;
Data1=randi（[01]n_symbol1）;                                             
%%%%%%%%%%%%%%%%%%  产生混沌序列
chaos=chaos_generate（n_symbol*beta0.1）;
Eb_av=sum（chaos.^2）/n_symbol;
%发送接收开始
dB_max=15;
BER_simulation=zeros（1dB_max+1）;
BER_GA=zeros（1dB_max+1）;
%BER_numerical=zeros（1dB_max+1）;
Eb=zeros（1n_symbol）;

for snr=0:dB_max
    %信息比特开始发送
    error_bit=0;
    for j=1:n_symbol
        
        ref=chaos（（j-1）*beta+1:（j-1）*beta+beta）; 
        Eb（j）=sum（ref.^2）;
        sigma=sqrt（2*Eb_av/（10^（snr/10））/2）; 
        temp=ref+randn（beta1）*sigma;                                          %加入AWGN
        
        
        t=（2*Data1（j）-1）*ref;