数字集群MDR中4FSK的matlab仿真

1 . 基本概念 物理层的定义h=0.27；T 是采样周期，1/4800 秒。而 D=3h/2T=1.944kHz。D 的物理含义是Maximum deviation，最大频偏。 4FSK 是恒包络调制，所以在终端功率上，应该比较省电。 4 进制的随机信号产生，并把信号强度调整到规范规定的3,1,-1,-3： M=4; Ts=1/4800; % 4.8K 的波特率； x=randint(1000,1,M); x=x.*2-3; %调整到正确的symbol； plot(x); 部分信号绘制如下： 2 1 02 2 02 3 02 4 02 5 02 6 02 7 02 8 02 9 0-3-2-10123 由于是随机信号，所以是宽带信号。频谱如下。只绘制到信号的采样频率，绘制频谱： lx=length(x); Freqs=1/Ts; freq=[-Freqs/2:Freqs/lx:Freqs/2-Freqs/lx]; S=fftshift(abs(fft(x))); Sdb=10*log10(S); plot(freq,Sdb);： -2 5 0 0-2 0 0 0-1 5 0 0-1 0 0 0-5 0 005 0 01 0 0 01 5 0 02 0 0 02 5 0 0-5051 01 52 02 5 进行基带滤波，滤波后的信号和频谱为下图。注意信号带宽被限制为根升余弦滤波的3dB带宽以内。这个带宽|cos(f*pi/1920)|=1/2，f>1920，得到 f=640+1920=2560Hz。这是一半带宽： % 滤波器的基本参数：sqrt，0.2。Fs 至少要是 2Fd，delay 用默认的3。 y=rcosflt(x,1,2,'fir/sqrt',0.2,3,0); %绘制频谱： ly=length(y); Freqs=2/Ts; freq=[-Freqs/2:Freqs/ly:Freqs/2-Freqs/ly]; Sy=fftshift(abs(fft(y))); Sdby=10*log10(Sy); figure; plot(freq,Sdby); -5 0 0 0-4 0 0 0-3 0 0 0-2 0 0 0-1 0 0 001 0 0 02 0 0 03 0 0 04 0 0 05 0 0 0-2 0-1 5-1 0-5051 01 52 02 5 6 2 06 4 06 6 06 8 07 0 07 2 07 4 07 6 0-3-2-10123 接下去，对 4FSK 调制进行仿真。 2. 4FSK 调制仿真 4FSK 有 3 种方式，方式 1，用前面余弦滤波后的基带数据进行频率振荡器的控制。这种方式实现容易，但是频率的稳定性差。方式2，键控法。键控分频器的分频比，以控制 2 种不同分频脉冲的输出。在经过整形，就可以得到 FSK 信号。这种在 DSP 的实现中常用到。方法 3，用 fskmod 函数。 不管用何种方式，DMR 规范中，对 4FSK 的调制的结果给出了要求，要求在给定序列下，输出的调制信号，应该和 1.2KHz 的sine 信号，通过最大偏差2.794KHz 的FM 调制器输出的信号一致。 方法 1： fmmod。 Fmmod 是模拟调制函数。它以 Fs 采样...