Частотная манипуляция

подскажите пожалуйста у меня есть модулированный по частоте сигнал



import matplotlib.pyplot as plt
import numpy as np
from math import pi
plt.close("all")

Fs=2000
t=np.arange(0,0.5,1/Fs)
fc=40
fm=7
b=3
fm=np.cos(2*pi*fc*t+b*np.sin(2*pi*fm*t))
plt.plot(t,fm)

и есть рандомные импульсы


import matplotlib.pyplot as plt
import numpy as np
plt.close('all')
sum_len=133
no_sum=32
rand_n=np.random.rand(no_sum)
rand_n[np.where(rand_n>=0.5)]=1
rand_n[np.where(rand_n<=0.5)]=0
plt.stem(rand_n)
sig=np.zeros(sum_len*no_sum)
id_n=np.where(rand_n==1)
for i in id_n[0]:
temp=int(i*sum_len)
sig[temp:temp+sum_len]=1
plt.figure()
plt.plot(sig) ;


как мне связать изменение частоты с двоичными импульсами
если проще написать другой код в котором частота изменялась согласно импульсам то лучше так (желательно чтоб он был с пояснениями)