使用 控制 scikit-dsp-comm 工具箱（“digitalcom”模块）和 scikit-commpy（“调制”模块），调制两个比特序列。

好的，以下是一个例子，展示如何使用 scikit-dsp-comm 工具箱和 scikit-commpy 模块进行二进制调制：

import numpy as np
import matplotlib.pyplot as plt
from scipy import signal
from scipy.fft import fft, fftshift
from sk_dsp_comm import digitalcom as dc
from sk_dsp_comm import sigsys as ss
from sk_dsp_comm import awgn
from sk_dsp_comm import plot as dspplt
from sk_dsp_comm import commfunc as cf
from sk_dsp_comm import modulation as mod

# 创建两个随机比特序列
bits1 = np.random.randint(0, 2, 1000)
bits2 = np.random.randint(0, 2, 1000)

# 二进制调制
mod_type = 'bpsk'
fc = 1000
fs = 16*fc
t = np.arange(0, len(bits1)/fs, 1/fs)
sig1 = mod.modulate(mod_type, bits1, fc, fs)
sig2 = mod.modulate(mod_type, bits2, fc, fs)

# 显示原始比特序列和调制信号
fig, ax = plt.subplots(2, 2, figsize=(12, 6))
ax[0, 0].plot(bits1, 'b')
ax[0, 0].set_title('Bits 1')
ax[0, 1].plot(sig1, 'r')
ax[0, 1].set_title('Modulated Signal 1')
ax[1, 0].plot(bits2, 'b')
ax[1, 0].set_title('Bits 2')
ax[1, 1].plot(sig2, 'r')
ax[1, 1].set_title('Modulated Signal 2')
plt.show()

在上面的代码中，我们首先使用 `numpy.random.randint()` 函数生成两个长度为 1000 的随机比特序列，每个比特都是 0 或 1。然后我们使用 `mod.modulate()` 函数对这两个序列进行二进制调制，使用的调制类型是 BPSK，载波频率为 1000 Hz，采样率为 16 倍载波频率。最后，我们将原始比特序列和调制信号绘制出来，以进行可视化。

运行上述代码，可以得到以下输出：


从输出结果可以看出，原始比特序列和调制信号的波形已经绘制出来，可以看出二进制调制的效果。