simulation_data.m

PPM调制的超宽带系统RAKE接收机的误码率仿真程序

clc;
clear;


% 仿真参数：
% Pow = -30;       % average transmitted power (dBm)
% fc = 50e9;      % sampling frequency
% numbits =1000;     % number of bits generated by the source
% Ts = 50e-9;      % frame time, i.e., average pulse repetition period [s]
% Ns = 1;          % number of pulses per bit
% Tc = 1e-9;       % chip time [s]
% Nh = 5;          % cardinality of the TH code
% Np = 2000;          % periodicity of the TH code
% 
% Tm = 0.5e-9;       % pulse duration [s]
% tau = 0.2e-9;      % shaping factor for the pulse [s]   
% dPPM = 0.5e-9;   % time shift introduced by the PPM [s]
% 仿真结果：
%  第一次仿真结果：
   ABER(1,:)=[2.1100000e-001  1.1100000e-001  4.0000000e-002  5.0000000e-003];
   BBER(1,:)=[2.6400000e-001  1.7400000e-001  8.1000000e-002  3.1000000e-002];
   CBER(1,:)=[3.0800000e-001  2.4500000e-001  1.6000000e-001  9.6000000e-002];
   DBER(1,:)=[3.3000000e-001  2.6600000e-001  1.7000000e-001  1.1500000e-001];
   EBER(1,:)=[3.9800000e-001  3.7600000e-001  3.0000000e-001  2.4900000e-001];
% 第二次仿真结果：  
  ABER(2,:)=[1.9000000e-001  9.3000000e-002  3.0000000e-002  5.0000000e-003];
  BBER(2,:)=[1.9000000e-001  1.1600000e-001  4.9000000e-002  9.0000000e-003];
  CBER(2,:)=[2.1700000e-001  1.3800000e-001  6.4000000e-002  1.0000000e-002];
  DBER(2,:)=[2.3700000e-001  1.5300000e-001  7.6000000e-002  1.9000000e-002];
  EBER(2,:)=[2.4800000e-001  1.7000000e-001  9.3000000e-002  2.6000000e-002];
%  第三次仿真结果：  
  ABER(3,:)=[2.0600000e-001  1.2500000e-001  4.2000000e-002  8.0000000e-003];
  BBER(3,:)=[2.5000000e-001  1.7000000e-001  9.5000000e-002  2.8000000e-002];
  CBER(3,:)=[3.5000000e-001  2.4400000e-001  1.9700000e-001  1.0400000e-001];
  DBER(3,:)=[3.0900000e-001  2.6300000e-001  1.8900000e-001  9.7000000e-002];
  EBER(3,:)=[4.6700000e-001  4.4600000e-001  3.8600000e-001  3.3800000e-001];  
%  第四次仿真结果：  
  ABER(4,:)=[1.2300000e-001  5.2000000e-002  1.2000000e-002  1.0000000e-003];
  BBER(4,:)=[1.5500000e-001  8.4000000e-002  2.3000000e-002  1.0000000e-003];
  CBER(4,:)=[2.2000000e-001  1.4200000e-001  4.1000000e-002  1.7000000e-002];
  DBER(4,:)=[2.6000000e-001  1.8300000e-001  8.8000000e-002  3.5000000e-002];
  EBER(4,:)=[3.4100000e-001  2.5200000e-001  1.6800000e-001  1.2100000e-001]; 
%  第五次仿真结果：
   ABER(5,:)=[2.2900000e-001  1.1700000e-001  5.5000000e-002  1.3000000e-002];
   BBER(5,:)=[2.3200000e-001  1.3600000e-001  7.3000000e-002  2.1000000e-002];
   CBER(5,:)=[2.5500000e-001  1.5700000e-001  8.5000000e-002  2.9000000e-002];
   DBER(5,:)=[2.4700000e-001  1.5900000e-001  9.2000000e-002  3.0000000e-002];
   EBER(5,:)=[4.4700000e-001  4.2400000e-001  4.2000000e-001  3.2900000e-001];
[W,L] =size(ABER);
aber=zeros(1,L);bber=zeros(1,L);cber=zeros(1,L);dber=zeros(1,L);eber=zeros(1,L);
exno=[0 3 6 9];
for j=1:L
    for i=1:W
        aber(j)=aber(j)+ ABER(i,j);
        bber(j)=bber(j)+ BBER(i,j);
        cber(j)=cber(j)+ CBER(i,j);
        dber(j)=dber(j)+ DBER(i,j);
        eber(j)=eber(j)+ EBER(i,j);
    end
    aber(j)=aber(j)/W;
end
aber 
figure(11);
semilogy(exno,aber,'b-o',exno,bber,'g-*',exno,cber,'r-d',exno,dber,'c-+',exno,eber,'y-s');
legend('理想RAKE','选择性RAKE(5个支路)','选择性RAKE(2个支路)','部分RAKE(5个支路)','部分RAKE(2个支路)');
xlabel('Ex/N0/dB');
ylabel('Prb');
title('MRC合并RAKE性能仿真');
grid;

% Pow = -30;       average transmitted power (dBm)
% fc = 50e9;      sampling frequency
% numbits =1000;     number of bits generated by the source
% Ts = 10e-9;      frame time, i.e., average pulse repetition period [s]
% Ns = 1;          number of pulses per bit
% Tc = 1e-9;       chip time [s]
% Nh = 5;          cardinality of the TH code
% Np = 2000;          periodicity of the TH code
% 
% Tm = 0.5e-9;       pulse duration [s]
% tau = 0.2e-9;      shaping factor for the pulse [s]   
% dPPM = 0.5e-9;   time shift introduced by the PPM [s]


%  第一次仿真结果：
   ABER(1,:)=[3.2200000e-001  2.6000000e-001  1.8400000e-001  1.5100000e-001];
   BBER(1,:)=[3.4300000e-001  2.9100000e-001  2.1700000e-001  1.9300000e-001];
   CBER(1,:)=[3.9500000e-001  3.5100000e-001  2.9400000e-001  2.7700000e-001];
   DBER(1,:)=[2.8900000e-001  2.3200000e-001  1.5100000e-001  7.6000000e-002];
   EBER(1,:)=[3.3000000e-001  2.9300000e-001  2.3100000e-001  1.7900000e-001];

% 第二次仿真结果：  
  ABER(2,:)=[3.5800000e-001  2.9850000e-001  2.6000000e-001  2.0000000e-001];
  BBER(2,:)=[3.4200000e-001  3.0050000e-001  2.5000000e-001  1.8350000e-001];
  CBER(2,:)=[4.2050000e-001  4.0300000e-001  3.9000000e-001  3.4700000e-001];
  DBER(2,:)=[2.8050000e-001  2.2350000e-001  1.4300000e-001  8.0500000e-002];
  EBER(2,:)=[3.4800000e-001  2.8750000e-001  2.3000000e-001  1.8200000e-001];
%  第三次仿真结果：  
  ABER(3,:)=[3.1250000e-001  2.5450000e-001  1.7800000e-001  8.6500000e-002];
  BBER(3,:)=[3.1750000e-001  2.3700000e-001  1.5900000e-001  8.4000000e-002];
  CBER(3,:)=[3.0900000e-001  2.4250000e-001  1.5200000e-001  8.2500000e-002];
  DBER(3,:)=[3.0200000e-001  2.3750000e-001  1.4900000e-001  7.8000000e-002];
  EBER(3,:)=[3.0900000e-001  2.4250000e-001  1.5200000e-001  8.2500000e-002];  
%  第四次仿真结果：  
  ABER(4,:)=[2.2700000e-001  1.6760000e-001  9.8200000e-002  4.7800000e-002];
  BBER(4,:)=[2.1580000e-001  1.5040000e-001  7.8200000e-002  3.1200000e-002];
  CBER(4,:)=[2.4200000e-001  1.6980000e-001  9.3400000e-002  3.2200000e-002];
  DBER(4,:)=[2.2660000e-001  1.4700000e-001  7.5200000e-002  2.6800000e-002];
  EBER(4,:)=[3.5820000e-001  3.2080000e-001  2.6320000e-001  1.9500000e-001]; 
%  第五次仿真结果：
   ABER(5,:)=[2.8620000e-001  2.3560000e-001  1.8720000e-001  1.3100000e-001];
   BBER(5,:)=[2.6500000e-001  1.9680000e-001  1.3660000e-001  7.9800000e-002];
   CBER(5,:)=[3.0300000e-001  2.4260000e-001  1.8320000e-001  1.1220000e-001];
   DBER(5,:)=[2.6580000e-001  1.8440000e-001  1.1740000e-001  5.2800000e-002];
   EBER(5,:)=[3.4340000e-001  2.8740000e-001  2.3260000e-001  1.7080000e-001];
%  第6次仿真结果：
   ABER(5,:)=[2.6600000e-001  1.9960000e-001  1.3320000e-001  9.7200000e-002];
   BBER(5,:)=[2.4800000e-001  1.8760000e-001  1.1380000e-001  7.1600000e-002];
   CBER(5,:)=[2.8340000e-001  2.0120000e-001  1.3420000e-001  8.4400000e-002];
   DBER(5,:)=[2.6780000e-001  1.9160000e-001  1.0040000e-001  4.9200000e-002];
   EBER(5,:)=[3.1400000e-001  2.5280000e-001  1.8340000e-001  1.2900000e-001];
  
%  第7次仿真结果：
   ABER(5,:)=[2.8860000e-001  2.0400000e-001  1.5220000e-001  9.7200000e-002];
   BBER(5,:)=[2.7160000e-001  1.8900000e-001  1.2360000e-001  6.5600000e-002];
   CBER(5,:)=[3.1460000e-001  2.5500000e-001  1.8020000e-001  1.1960000e-001];
   DBER(5,:)=[2.7140000e-001  1.9300000e-001  1.2140000e-001  6.3400000e-002];
   EBER(5,:)=[4.3300000e-001  4.1420000e-001  3.7880000e-001  3.4580000e-001];
  
%  第8次仿真结果：
   ABER(5,:)=[4.3600000e-001  4.1400000e-001  3.9380000e-001  4.0300000e-001];
   BBER(5,:)=[4.4260000e-001  4.3060000e-001  4.1440000e-001  4.0680000e-001];
   CBER(5,:)=[4.5000000e-001  4.5520000e-001  4.3660000e-001  4.3320000e-001];
   DBER(5,:)=[3.8240000e-001  3.3700000e-001  2.6340000e-001  1.9920000e-001];
   EBER(5,:)=[4.6400000e-001  4.3620000e-001  4.2560000e-001  4.0040000e-001];
  
  

[W,L] =size(ABER);
aber=zeros(1,L);bber=zeros(1,L);cber=zeros(1,L);dber=zeros(1,L);eber=zeros(1,L);
exno=[0 3 6 9];
for j=1:L
    for i=1:W
        aber(j)=aber(j)+ ABER(i,j);
        bber(j)=bber(j)+ BBER(i,j);
        cber(j)=cber(j)+ CBER(i,j);
        dber(j)=dber(j)+ DBER(i,j);
        eber(j)=eber(j)+ EBER(i,j);
    end
    aber(j)=aber(j)/W;
end
aber 
figure(12);
semilogy(exno,aber,'b-o',exno,bber,'g-*',exno,cber,'r-d',exno,dber,'c-+',exno,eber,'y-s');
legend('理想RAKE','选择性RAKE(5个支路)','选择性RAKE(2个支路)','部分RAKE(5个支路)','部分RAKE(2个支路)');
xlabel('Ex/N0/dB');
ylabel('Prb');
title('MRC合并RAKE性能仿真');
grid;