untitled.m

基于matlab上进行扩频技术的仿真研究（M文件实现）

%==============平移对应法==============
%==============基于Logistic映射==============
clear;
q=8;%频隙数
L=60;
N=L*q;
d=10;%跳频间隔
x0(1)=0.35;
for i=2:N
  x0(i)=4*x0(i-1)*(1-x0(i-1));
end
%==============将混沌序列变换为二进制序列==============
for i=1:N
    if  (x0(i)>=0.5)
        x0(i)=1;
    else
        x0(i)=0;
    end
end
x0=reshape(x0,N/6,6);
%==============将产生的二值符号序列转换为q元跳频序列==============
%==============进行T位截断，得到频隙为q=2^k=2^6=64的混沌序列==============
for i=1:N/6
        if x0(i,1)==1
            x0(i,1)=32;
        else
            x0(i,1)=0;
        end
    end
for i=1:N/6
    if x0(i,2)==1
            x0(i,2)=16;
        else
            x0(i,2)=0;
        end
    end
for i=1:N/6
    if x0(i,3)==1
            x0(i,3)=8;
        else
            x0(i,3)=0;
        end
    end
for i=1:N/6
    if x0(i,4)==1
            x0(i,4)=4;
        else
            x0(i,4)=0;
        end
    end
for i=1:N/6
    if x0(i,5)==1
            x0(i,5)=2;
        else
            x0(i,5)=0;
        end
    end
x0=(x0)';
x0=sum(x0);
x1=x0;
%==============进行宽间隔改进==============
for i=1:(N/6)-1
    if  (abs(x1(i+1)-x1(i))>=(d+1))
        x1(i+1)=x1(i+1);
    else if (abs(x1(i+1)-x0(i))<(d+1)&(x1(i+1)>=x1(i)))
        x1(i+1)=mod((x1(i+1)+d+1),q);
    else (abs(x1(i+1)-x0(i))<(d+1)&x1(i+1)<x1(i))
           x1(i+1)=mod((x1(i+1)-d-1),q);
       end
   end
end
%汉明自相关
y1=[x1 x1];
for j=1:N/6
y2(j,:)=y1(j:(N/6)-1+j);
end
for i=1:N/6
    for j=1:N/6
    x1(i);
    y2(j,:);
if (x1(i)==y2(j,i))
    h(j,i)=1;
else
    h(j,i)=0;
end
end
end
H=sum(h,2);
for i=2:N/6
    R(i)=H(i);
end
R1=max(R);%R2为自相关主峰，R1为自相关得最大旁瓣
R2=N/6;     %R1/R2可以决定某一跳频序列的抗干扰能力
A=R1/R2