1.算法概述
那么FM调制端的基本原理图如下所示：
FM解调
FM解调的基本原理框图为：
在不考虑频偏的情况下，FM解调运算就相当于FM调制的逆运算，任务的第一步比较简单，所以这里关于FM解调的原理就不在叙述。
对该系统进行系统仿真，在仿真之前，要将系统做进一步的模块化从而有利于仿真分析。
2.仿真效果预览
matlab2013b
上图是系统的发送信号，第二章图是通过FM调制以后的信号。系统发送出去的就是这个信号。
这是系统接收到的信号，我们可以看到，接收到的信号在初始有一个抖动的过程，这个就是PLL的锁定过程，通过这个锁定过程，系统就可以得到正确的解调信号，上图第三章图是PLL的相位输出，通过一开始的抖动，相位误差就锁定在0附近，这说明相位已经锁定了。此后的解调信号都是正确的。
以上是一个正弦信号做为发送信号，下面我们用一个方波来作为发送信号检测系统。
当发送发波的时候也可以锁定，但是在波形的边沿会产生误差信号。
3.核心MATLAB代码预览
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
