%%?跳頻信號的產(chǎn)生（調(diào)制方式為BPSK）
%%%%信號參數(shù)%%%%
%fs:采樣率
%Rs:跳速
%time_delay:時(shí)間延遲
%f1-f2:載波頻段
%data:跳頻序列（未碼元擴(kuò)展的）
%Rd:碼速率
%dataout:跳頻序列引入
%period:跳頻碼周期，小于L
%M:模式選擇
function?v=signal_bpsk（fsRstime_delayf1f2RdNFdataoutperiodMInput_Duty_ratio）
[Um?L]=size（dataout）;????????????????????????%其中Um是用戶數(shù)，L為跳頻碼的個(gè)數(shù)
for?n?=?1:Um
s=randsrc（round（period*Rd/Rs）1[0:M-1]）;????%隨機(jī)產(chǎn)生長度為period*Rd/Rs的0?1
s_mod=pskmod（sM）;
%%?計(jì)算出每個(gè)碼元在一個(gè)矩形脈沖內(nèi)的點(diǎn)數(shù)
m?=?0;
y?=?[];
N?=?zeros（1length（s_mod））;
for?k?=?1?:?length（s_mod）
????N（k）?=?ceil（k*fs/Rd）-m;
????m?=?m+N（k）;
????y_k?=?rectpulse（s_mod（k）N（k））;
????y_k?=?y_k‘;
????y?=?[y?y_k];
end
%%?低通濾波
%?a?=?0.35;
L_zeros?=?round（（1-Rd/fs）*length（y）/2）;
L_ones?=?length（y）-L_zeros*2;
Window?=?[zeros（1?L_zeros）?ones（1?L_ones）?zeros（1?L_zeros）];
s_base?=?ifft（fftshift（fftshift（fft（y））.*Window））;
signal_base（n:）=s_base;?????????%基帶信號
%?figure（1）
%?plot（real（signal_base））;
%?此處這個(gè)signal_base就是原來程序里的signal_base
end
%%?BPSK調(diào)制
Nh=round（fs/Rs）;??????????????????????%跳頻碼的一個(gè)周期采樣點(diǎn)數(shù)
deta=（f2-f1）/（NF-1）;
for?k=1:Um
%?????x（k:）?=?signal_base（k:）;
????s=sign（rand（1round（Rd/Rs*period））-0.5）;
????x（k:）=rectpulse（sround（fs/Rd））;
????data（k:）=rectpulse（dataout（k:）Nh）;
????if?k~=1
%?????????data（k:）=[zeros（1time_delay*（k-1））?data（k1:end-time_delay*（k-1））];%延遲
????????x（k:）=[zeros（1time_delay*（k-1））?x（k1:end-time_delay*（k-1））];
????end
end
if?length（data）?<=?length（x）
???x=x（:1:length（data））;
else
???data=data（:1:length（x））;?
end
F=data*deta+f1;???????????????????????%跳頻頻率表
for?k=1:Um
????if?k~=1
????????F（k:）=[zeros（1time_delay*（k-1））?F（k1:end-time_delay*（k-1））];%延遲
????end
end
L=length（F）;
t=（1/fs）*（0:L-1）;
for?k=1:Um
????v（k:）=x（k:）.*exp（1i*2*pi*（F（k:）+rand）.*t）.*exp（1i*2*pi*rand）;%跳頻信號生成
end
nh=round（Nh*（1-Input_Duty_ratio））;
for?k=1:Um
????for?i=1:period
????????????v（kNh*i-nh+time_delay*（k-1）:Nh*i+time_delay*（k-1））=0;%跳頻間隔
????end
end
v?=?v（:1:round（period*fs/Rs））;
%?for?k=1:Um
%?????if?k==1
%????????v（kNh-4*nh:Nh）=0;
%?????else
%????????v（kNh-nh*4+time_delay*（k-1）:Nh+time_delay*（k-1））=0;
%?????end
%?end?
clearvars?-except?v
end

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----

?????文件???????2323??2018-05-21?17:13??跳頻信號-BPSK調(diào)制\signal_bpsk.m

?????目錄??????????0??2018-05-21?17:13??跳頻信號-BPSK調(diào)制

-----------?---------??----------?-----??----

?????????????????2323????????????????????2