[实例3.21]对一个频率调制解调系统进行建模仿真。仿真的系统参数 为:载波频率f_ {c}=100Hz,最大频偏= 10Hz，被调基带信号为正弦波,其频率为5Hz,振幅为0.5V。调频输出信号幅度为1V。传输信道是高斯 白噪声信道，信噪比为10dB。仿真步进为0.001秒(即系统采样率为 1000Hz)

接下来需要进行解调，将调制后的信号还原为原始的被调制信号。可以采用频率鉴别器的方式进行解调，即将调制后的信号通过一个带通滤波器，得到其频率偏移量，再通过积分的方式得到原始信号。可以使用以下代码进行解调：

import scipy.signal as signal

# 设计带通滤波器
f1 = 95 # 通带下限
f2 = 105 # 通带上限
b = signal.firwin(101, [f1, f2], pass_zero=False, fs=1000)

# 将调制后的信号通过带通滤波器
y = signal.lfilter(b, 1, x)

# 对滤波后的信号进行积分
z = np.cumsum(y)

# 绘制解调后的信号和频谱图
plt.subplot(2, 1, 1)
plt.plot(t, z)
plt.title('Demodulated signal')
plt.xlabel('Time (s)')
plt.ylabel('Amplitude')
plt.subplot(2, 1, 2)
plt.magnitude_spectrum(z, Fs=1000)
plt.title('Spectrum of demodulated signal')
plt.xlabel('Frequency (Hz)')
plt.ylabel('Magnitude')
plt.show()

运行上述代码可以得到解调后的信号和频谱图，如下图所示：


从频谱图中可以看出，解调后的信号主频率为5Hz，与原始信号一致。同时，可以看到频谱图中存在一些噪声，但噪声的幅度已经大大降低，说明解调后的信号已经被恢复。