通信算法之二十：DMR系統 4FSK調制解調仿真鍊路

一.  4FSK原理簡介

         DMR系統，數字叢集通信系統。DMR協定采用的調制方式4FSK，是一種恒包絡調制，調制時每秒發送2400個符号，每個符号攜帶兩比特的資訊。最大頻偏D定義如下：D = 3h/2T，h代表每個調制的頻偏系數，0.6。T标示符号周期，為1/2400。D=2160。

       4FSK調制器由一個平方根升餘弦濾波器級聯一個頻率調制器組成。第一部分是成型濾波器子產品，産生四電平的基帶帶限信号作為調制信号。第二部分是跳頻部分FM。

      調頻是角度調制的一種，角度調制一般表示為：Sm(t) = Acos(Wc(t)+ )。頻率調制是利用基帶信号載波頻率得到的，即瞬時頻偏随調制信号m(t)線性變化。滿足： =Kf.m(t) 。其中，Kf為頻偏常數。

二. 4FSK調制信号産生

       産生調制信号主要有兩種，直接法和間接法。間接法實作成本和複雜度高。直接法調頻是調制信号m(t)直接控制振蕩器的頻率使其按照調制信号的規律線性變化。調頻信号的正交調制是間接法調頻的一種實作方式.

三. MATLAB仿真鍊路

%****************************************************************

% 内容概述：4FSK調制與4FSK解調 仿真鍊路

% 仿真鍊路：

% 創 建 人：lee gang

% 建立時間：2016年10月23日

% 參考文獻：DMR通信系統

% 版權聲明：不可随意複制、不可随意傳播。

%****************************************************************

clc;

clear all;

close all;

%系統參數設計---------------------------------

Fs = 30000;%采樣頻率.30k，256k

Ts = 0.1;%符号周期.輸入符号速率10Hz,2.4k

SampleNum = Fs * Ts;%采樣點                     

A = 1;   %ampilitude is 1

w1 = 300;%2160,%初始化載波頻率

w2 = 600;%720

w3=900;  %-720

w4=1200; %-2160 

%輸入源比特序列基帶信号-------------------------------------

Sign_Set= randint(1,200);

LenthSourceBit = length(Sign_Set);

%信号分離成IQ信号-------------------------------------------

LenthQp = LenthSourceBit/2; 

t = 0......................;

for I=1:LenthSourceBit/2                            

    Sign_Set1(I)= Sign_Set(2*I-1);

    Sign_Set2(I)=Sign_Set(2*I);

end

%産生基帶信号---------------------------------------------------

[sign_orign,st] = Fsk4Modulation(Sign_Set1,Sign_Set2,LenthQp,SampleNum,t,w1,w2,w3,w4);

%調制之後的4FSK-畫信号圖----------------------------------

figure;su......;

axis([0,.......]);

title('輸入源比特基帶信号');grid;

subplot....);

axis([0,......]);

title('4FSK調制信号');grid;

%高斯白噪聲信道-------------------------------------------------------

st = awgn(st,-2);

figure;plot(t,st);

title('接收到的4FSK信号');grid;

%帶通低通濾波器-------------------------------------------------------

[....1,..2,......t3,...t4] = Fsk4Demodulation(s.....,F...s,w...,w...,w...,w....);

%畫圖---部分代碼省略（）

figure;subplot(2,2,1)....................

axis([0,...........................................]);

title('F1分量低通濾波後的波形');grid;

subplot(....,....,2);plot(t,........);

axis([0,Ts..............]);

title('F2分量低通濾波後的波形');grid;

subplot(2.............................

axis([0,T.................................)]);

title('F3分量低通濾波後的波形');

....................................................

title('F4分量低通濾波後的波形');grid;

%抽樣判決---------------------------------------------------

[...,..] = Fsk4DeSampsetence(d...,dt...,.....,...,....,...);

%輸出源比特序列----------------------------------------------

.............................

......................................

[errNum,errRate] = biterr(Sign_Set,FskDataOut);

disp(['誤碼率 = ',num2str(errRate),', 錯誤比特數 = ',num2str(errNum)]);

四. 仿真圖形

[詳細資料，及相應MATLAB代碼，C語言代碼。咨詢qq：1279682290 ]