一. 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 ]