跳频扩频系统的MATLAB仿真

此系统是对跳频扩频系统所作的MATLAB仿真，取得了良好的效果.可作为毕业设计和教学使用

clc;
clear all;
N=8;
begin=0;
 picture=[0 1 0 1 0 0 1 1 0 1 1 1 0 0 0 1 0 1 1 0 1 0 0 1 0];         
if begin==0                                                       
    temp=25*rand;
    begin=floor（temp）+1;
    if （begin<19）
        n=picture（begin:begin+7）;
    else
      n=[picture（begin:25）picture（1:begin-18）];
    end
else                                                          
    begin=begin+1;
    if begin==26
        begin=1;
    end
    if （begin<19）
        n=picture（begin:begin+7）;
    else
      n=[picture（begin:25）picture（1:begin-18）];
    end
end
PN=n;
%PN=[0 1 2 3 0 1 2 3];
b=PN;
z=zeros（N101）;
%信息码元
figuresubplot（211）
for i=1:8
    if（b（i）==0）
        m=zeros（1101）;
        t=i-1:0.01:i;
        plot（tm）
        hold on;
    else
        m=b（i）*ones（1101）;
        t=i-1:0.01:i;
        plot（tm）
        hold on;
    end;   
end;
xlabel（‘信息码元‘）;
axis（[0 8 -2 2]）;
grid on;
hold off;
%----------------------------------------------------
  %FSK调制
 b=PN;
subplot（212）
for i=1:8
    if（b（i）==0）
        m=zeros（1101）;
        t=i-1:0.01:i; 
        y=sin（t.*（2）*pi）;
        plot（ty）;
        hold on;
    else
        m=b（i）*ones（1101）;
        t=i-1:0.01:i;
         y=sin（t.*（2+4*m）*pi）;
        plot（ty）;
        hold on;
    end;
   %y=sin（t.*（2+2*m）*pi）; 
   z（i:）=y;
  % plot（ty）;   
end;
xlabel（‘FSK调制‘）;
axis（[0 8 -2 2]）;
grid on;
hold off;
y1=z;
yfsk=[];

for i=1:8
    yfsk=[yfsky1（i1:100）];
end
figure  
yfsk=[yfsk0];
lnx=length（yfsk）;
nfft=lnx+1;
f = lnx*（0:nfft/2）/nfft;
%freq=-pi:2*pi/（lnx-1）:pi;% the frequency vector频率向量，长度为num
X=fft（yfsknfft）;
%plot（freq/pi abs（X））;ylabel（‘|X|‘）;axis（[0 pi/pi min（abs（X）） max（abs（X））]）;
plot（fabs（X（1:nfft/2+1）））;
title（‘FSK信号频谱图‘）;
N=8;
 picture=[0 1 2 3 0 1 2 3 0 1 2 3 0 1 2 3 0 1 2 3 0 1 2 3 0];          %跳频图案
if begin==0                                           %判断是否是跳频开始时刻
    temp=25*rand;
    begin=floor（temp）+1;
    if （begin<19）
        n=picture（begin:begin+7）;
    else
      n=[picture（begin:25）picture（1:begin-18）];
    end
else                                                     %按跳频图案顺序跳变
    begin=begin+1;
    if begin==26
        begin=1;
    end
    if （begin<19）
        n=picture（begin:begin+7）;
    else
      n=[picture（begin:25）picture（1:begin-18）];
    end  
end
PN=n;
%PN=[0 1 2 3 0 1 2 3];
b=PN;
z=zeros（N101）;
%----------------------------------------------------
  %画PN码的数字直观图
%axes（e1）;
figuresubplot（211）
for i=1:8
    if（b（i）==0）
        m=zeros（1101）;
        t=i-1:0.01:i;
        plot（tm）
        hold on;
    else
        m=b（i）*ones（1101）;
        t=i-1:0.01:i;
        plot（tm）
        hold on;
    end;
end;
xlabel（‘跳频码‘）;
axis（[0 8 -0.2 6.2]）;
grid on;
hold off;
%---------------------------------------------------------
   %画频率合成后的图
subplot（212）
for i=1:N  
a=b（i）;
%为频率合成部分
if（a==0）
    m

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       9639  2008-01-03 14:29  fhss.m

----------- ---------  ---------- -----  ----

                 9639                    1