ieeechannel.txt

run in matlab 6.0 or above UWB channel simulink .

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),'-*')