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;????????
???????