跳频通信系统MATLAB仿真

主要是针对跳频通信整个流程的原理性仿真，包含主程序和两个子函数，附带一份自制的PPT，PPT中的结果基于该代码，供初学者参考学习使用。

clc;clear;
g=40;fs=100000; 
r=-10;delay=0;  
sig1=round（rand（1g））;     %产生随机信号源 
signal1=[];   
for k=1:g                  %离散点化     
if sig1（1k）==0          
sig=-ones（11000）;         % bit 0设置1000个样点     
else          
sig=ones（11000）;          % bit 1设置1000个样点     
end      
signal1=[signal1 sig];
end  
figure（1）  
plot（signal1‘b‘‘linewidth‘1）;           %信源波形
grid on;
axis（[-100 1000*g -1.5 1.5]）; 
title（‘信号源‘）

T0=200; f0=1/T0; 
T1=400; f1=1/T1;  
u0=gensig（‘sin‘T01000*g-11）;
u0=rot90（u0）; u1=gensig（‘sin‘T11000*g-11）;
u1=rot90（u1）; y0=u0.*sign（-signal1+1）; 
y1=u1.*sign（signal1+1）;  
SignalFSK=y0+y1;          % 生成的FSK信号 
figure（2）;
% subplot（211）;  
plot（SignalFSK）           % FSK信号的时域波形 
axis（[-100 1000*g -3 3]）; 
title（‘SignalFSK‘） 

% % FSK信号频谱 
% subplot（212）  
% Plot_f（SignalFSK fs）; 
% title（‘FSK调制后的频谱‘）; 

t1=（0:100*pi/999:100*pi）;            
t2=（0:110*pi/999:110*pi）;                 
t3=（0:120*pi/999:120*pi）;                            
t4=（0:130*pi/999:130*pi）;                   
t5=（0:140*pi/999:140*pi）;                
t6=（0:150*pi/999:150*pi）;        
t7=（0:160*pi/999:160*pi）;   
t8=（0:170*pi/999:170*pi）;   
c1=cos（t1）;                      
c2=cos（t2）;  
c3=cos（t3）;                      
c4=cos（t4）;
c5=cos（t5）;
c6=cos（t6）; 
c7=cos（t7）; 
c8=cos（t8）; 
adr1=Mcreate（1001203）;  
adr1=[adr1adr1（1）adr1（2）];      %用户地址为初始m序列
fh_seq1= []; 
for k=1:g   
seq_1=adr1（3*k-2）*2^2+adr1（3*k-1）*2+adr1（3*k）;   
fh_seq1=[fh_seq1 seq_1];              %生成用户载波序列 
end

spread_signal1=[];           %用户一载波
fhp=[]; 
for k=1:g      
c=fh_seq1（k）;     
switch（c）         
case（0）              
spread_signal1=[spread_signal1 c8];         
case（1）              
spread_signal1=[spread_signal1 c1];                %形成随机载频序列         
case（2）              
spread_signal1=[spread_signal1 c2];         
case（3）              
spread_signal1=[spread_signal1 c3];         
case（4）              
spread_signal1=[spread_signal1 c4];         
case（5）                      
spread_signal1=[spread_signal1 c5]; 
case（6）              
spread_signal1=[spread_signal1 c6];         
case（7）              
spread_signal1=[spread_signal1 c7];                
end      
fhp=[fhp （500*c+5000）]; 
end

figure（3） %跳频图案
plot（fhp‘s‘‘markerfacecolor‘‘b‘‘markersize‘12）; 
grid on;

freq_hopped_sig1=SignalFSK.*spread_signal1;   %跳频扩频调制（类似于幅度调制）
figure（4）;  
% subplot（211）;  
plot（（1:1000*g）freq_hopped_sig1）;    %跳频扩频后的时域信号
axis（[-100 1000*g -2 2]）;  
title（‘跳频扩频后的时域信号‘）; 
% % 扩频调制后的频谱 
% subplot（212）;  
% Plot_f（freq_hopped_sig1fs）; 
% title（‘扩频调制后的频谱‘）;

% figure（5）%时域跳频信号
% plot（spread_signal1）;
% axis（[0 inf -5 5]）;
% 
% figure（6）%时域跳频信号频谱
% Plot_f（spread_signal1fs）;

% 加多径 
s1=freq_hopped_sig1;  
s=[zeros（1delay） s1（1:（1000*g-delay））]; 
freq_hopped_sig1=freq_hopped_sig1+s;    
% 加高斯白噪声  
awgn_signal=awgn（freq_hopped_sig1r1/2）;%信噪比为r；

figure（7）; 

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

     文件       5241  2017-12-11 11:37  跳频通信原理\FH_practice2.m

     文件        449  2017-12-11 11:37  跳频通信原理\Mcreate.m

     文件        242  2017-11-28 16:54  跳频通信原理\Plot_f.m

     文件    1756924  2017-12-11 11:43  跳频通信原理\跳频通信原理.pdf

     目录          0  2017-12-11 11:46  跳频通信原理

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

              1762856                    5