通信算法之九：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)直接控制振荡器的频率使其按照调制信号的规律线性变化。调频信号的正交调制是间接法调频的一种实现方式.
三。4FSK调制解调算法设计流程
发射机设计：
接收机设计：
四。MATLAB仿真链路
%****************************************************************
% 内容概述：4FSK调制与4FSK解调 仿真链路
% 仿真链路：
% 创 建 人：
% 创建时间：2014年10月23日
% 修改时间：2014年10月28日 理想环境下调试通过
% 修改时间：2014年10月31日 高斯白噪声信道下调试通过
% 修改时间：2014年11月2日 高斯白噪声信道、有相差、有频偏直流分量环境下调试通过
% 参考文献：DMR通信系统各类工程实现文献
% 版权声明：不可随意复制、不可随意传播。
%****************************************************************
clc;
clear all;
close all;
%----------------发送端物理层基带算法---------------------------------------
%-1.输入基带信号-源比特序列--串并转换-映射-----------------------------------
%-2.成型滤波器-升余弦滤波器-------------------------------------------------
%-包含：八倍内插和平方根升余弦成型滤波器，减小基带信号带宽占用
Fsksignal = Squrootcosfilter(SourceData,Fd,Fs,Type,Rolloff,Delay);
%-----------------------以上是发送端物理层基带算法--------------------------
%--------------------发送端物理层中频算法数字上变频-------------------------
%-3.相位累加求和-----------------------------------------------------------
phase = Phasesum(Fsksignal);
%-符号速率2.4kbps，符号周期T ，调制系数H ，最大频偏D：D = 3h/2T
[II QQ] = phaseSinCos(phase,factor);
%-5.数字中频-数字上变频-内插滤波-CIC-HB-FIR---------------------------------
% 根据奈奎斯特抽样定理，采样频率FS应该大于2*FC+Bm.在实际应用中，一般要
% 取采样频率FS大于或等于四倍的模拟信号频率。FS >= 4*FC
% 【待开发数字中频算法：CIC-HB-FIR，符号速率：4*FC】
%-6.数字中频-数字上变频-载波体制-sin(Wc*t)--cos(Wc*t)------------------------
% 【待开发数字中频算法：正交，载波调制】
% s(t) = A*cos(2*pi*FC*t)*I(t) - A*sin(2*pi*FC*t)*Q(t)
%--------------------以上发送端物理层中频算法数字上变频----------------------
%-7.空间信道传输-高斯白噪声信道---------------------------------------------
title('过高斯白噪声信道之后的信号');
%-8.数字中频-数字下变频-正交解调算法----------------------------------------
%-9.数字中频-数字下变频-抽取滤波-------------------------------------------
% 【待开发，数字下变频算法抽取滤波】
%-10.空间信道传输-加相位差值--加频偏值----------------------------------------
PPHH = (........................);
NewSig = SigComplex.*PPHH;
II = real(NewSig);
QQ = imag(NewSig);
figure()
plot(II(500:800),'m-');
grid on;
title('过高斯白噪声信道/相位差/频率差之后的信号');
%-11.相位估计、相差，差分算法--消除相干解调中出现的本地载波与信号相差的问题
%-13.信号标准门限-动态范围确定------------------------
%-确定动态范围---计算标准门限---------------------------
[up down mid mean]=SignalDynamicRange(mx);
[up_level_std mid_level_std down_level_std]=SignalBiaozhunRange(up,down,mean);
%-16.门限判决，逆映射--------------------------------------
title('--4FSK调制解调仿真误码率曲线--');
disp('--4FSK调制，4FSK解调仿真结束--');
%---------------------------------------------
五。仿真图形
【详细资料，及相应MATLAB代码，C语言代码。咨询qq：1279682290】