Delta1=3.4;%簇幅度的分布的方差
Delta2=3.4;%脈沖幅度的分布的方差
IndexOfCluDecay=1;
IndexOfRayDecay=1;

%IndexOfDecayOfClu=7.1;
%IndexOfDecayOfRay=4.3;
%POfClu=0.0233;
%POfImpulse=0.0233*2.5;

%POfClu=0.4;
%POfImpulse=0.5*0.4;
%IndexOfDecayOfClu=5.5;
%IndexOfDecayOfRay=6.7;

%POfClu=0.0667;
%POfImpulse=2.1*0.0667;
%IndexOfDecayOfClu=14;
%IndexOfDecayOfRay=7.9;

POfClu=0.0667;
POfImpulse=2.1*0.0667;
IndexOfDecayOfClu=24;
IndexOfDecayOfRay=12;

Impulse=[1:1:300];%存儲每個位置的脈沖幅度
for i=1:300
  Impulse(i)=0;
end;
for j=1:1000%1000次隨機
  DecayOfClu=1;
  StartOfClu=1;
  for i=1:300 %考查范圍300ns
    Pclu=unifrnd(0,1);
    if (Pclu<0.0233) %產生了一個新簇
      IndexOfCluDecay=10^(normrnd(0,Delta1)/10);
      DecayOfClu=IndexOfCluDecay*exp(-(i-1)/IndexOfDecayOfClu);%更新簇衰落指數
      StartOfClu=i;%紀錄簇的起始位置
    end;

    Pimpulse=unifrnd(0,1);
    if(Pimpulse<(0.0233*2.5))%簇內產生脈沖
      IndexOfRayDecay=10^(normrnd(0,Delta2)/10);
      DecayOfImp=IndexOfRayDecay*DecayOfClu*exp(-(i-StartOfClu)/IndexOfDecayOfRay);%簇內衰落
      Impulse(i)=Impulse(i)+DecayOfImp;
    end;  
   end;%考查范圍300ns
end;%1000次隨機 



%%%%%%%%求碼間干擾
Imp=Impulse1;
E=0;

for i=1:300
  E=E+Imp(i)^2;
end;
Imp=Imp/sqrt(E);%能量歸一化


u=91.5180/400;% 高斯脈沖的平均干擾
Tf=[1:1:10];
Tf=20*Tf;%脈沖重復周期 40，80，120，160，200
a=[1:1:10];%記錄對應不同周期的干擾
aNoth=[1:1:10];%記錄對應不同周期的干擾
for i=1:10%不同的Tf
    sum0=0;
    sum0Noth=0;
    for j=1:10000
        x=unifrnd(0,1);
        y=unifrnd(0,1);
        z=Tf(i)*(x-y);%采用跳時方案，前后脈沖之間隨機距離
        %%%前后脈沖之間的距離取整
        for jj=-Tf(i):Tf(i)
          if(((jj-1)<z)&&(z<=jj))            
            z=jj;
            break;
          end;
        end;
        %%%
         
        zNoth=1;%不采用跳時方案
        
        %跳時干擾
        for k=1:50%前面50個周期可能有干擾
          if((z+Tf(i)*k)>300)
            break;
          end; 
          sum0=sum0+Imp(z+Tf(i)*k);  
        end;
        %非跳時干擾
        for k=1:50%前面50個周期可能有干擾 
          if((zNoth+Tf(i)*k)>300)
            break;
          end; 
          sum0Noth=sum0Noth+Imp(zNoth+Tf(i)*k);  
        end;        
    end;
    a(i)=sum0/10000;%平均干擾
    aNoth(i)=sum0Noth/10000;%平均干擾
end;
ICI1=a*u;
ICI1Noth=aNoth*u;

plot(Tf,10*log10(ICI1),'-o',Tf,10*log10(ICI1Noth),'-*')