1.算法概述

那么FM调制端的基本原理图如下所示：

FM解调

FM解调的基本原理框图为：

在不考虑频偏的情况下，FM解调运算就相当于FM调制的逆运算，任务的第一步比较简单，所以这里关于FM解调的原理就不在叙述。

对该系统进行系统仿真，在仿真之前，要将系统做进一步的模块化从而有利于仿真分析。

2.仿真效果预览

matlab2013b

上图是系统的发送信号，第二章图是通过FM调制以后的信号。系统发送出去的就是这个信号。

这是系统接收到的信号，我们可以看到，接收到的信号在初始有一个抖动的过程，这个就是PLL的锁定过程，通过这个锁定过程，系统就可以得到正确的解调信号，上图第三章图是PLL的相位输出，通过一开始的抖动，相位误差就锁定在0附近，这说明相位已经锁定了。此后的解调信号都是正确的。

以上是一个正弦信号做为发送信号，下面我们用一个方波来作为发送信号检测系统。

当发送发波的时候也可以锁定，但是在波形的边沿会产生误差信号。

3.核心MATLAB代码预览

...............................
N  = 4000;                         % 采样点数
Fs = 40000;                        % 采样频率
Ts = 1/Fs;                         % 采样周期


%注意，这里不同的频域，不同改变的太夸张，否则系统也无法正常工作。46000~54000

%5000~15000,25000~35000,45000~55000,65000~75000,85000~95000,105000~115000..
%.......................
k=input('请输入你要的电台序号：0,1，2，3，4，5，6，7，8，9\n\n\n');
 

fc1= k*20000+8000+4000*rand(1)                     %发送载波频率

fc2 =k*20000+10000                                 %接收载波频率


t = [0:Ts:(N*Ts)- Ts];             
start_point = 1;
end_point = 1000; %查看仿真的持续时间周期

%锁相环的初始化                                              



TestFreq = 500;                                                %FM调制频率
msg = 0.5*square(2*pi*TestFreq*t);                             %要发送的信号
[fmmsg, phase] = fm_modulation(msg, fc1, Fs, 1.0, 0.05, 0);     %调用FM调制函数进行调制

%以上就是FM的调制过程
%=============================================================
%=============================================================
%一下是PLL-FM解调过程
[pd_out,lfsum,lf_out,vco_phase,vco_out] = pll_function(fmmsg,Fs,2.0,1.0,1.0,fc2,50000,1.0);



figure(1);
subplot(3,1,1);
plot(fmmsg(start_point:end_point));
title('10K的载波对500HZ的信号进行调制');
grid;
subplot(3,1,2);
plot(lf_out(start_point:end_point));
title('PLL的环路滤波输出信号');
grid;

subplot(3,1,3);
plot(pd_out(start_point:end_point));%   fc1>fc2,值<0，
title('PLL相位输出');                %   fc1<fc2,值>0，
%                                       Rfc1=fc2,值-5~5
A01-11