利用Matlab编写循环码的调制解调程序；BPSK、QPSK

以下是一个基于 Matlab 的循环码调制解调程序示例，包括 BPSK 和 QPSK 两种调制方式：

% 循环码调制解调程序示例

% 参数设置
N = 7; % 循环码长度
K = 4; % 循环码信息位长度
n = 10000; % 发送比特数
SNR_dB = 5; % 信噪比（dB）

% 生成循环码
genPoly = cyclpoly(N,K);
trellis = poly2trellis(N,genPoly);
convEncoder = comm.ConvolutionalEncoder('TrellisStructure',trellis);

% 生成调制器和解调器
bpskMod = comm.BPSKModulator();
qpskMod = comm.QPSKModulator();
bpskDemod = comm.BPSKDemodulator();
qpskDemod = comm.QPSKDemodulator();

% 生成信道和误差率计算器
channel = comm.AWGNChannel('NoiseMethod','Signal to noise ratio (SNR)','SNR',SNR_dB);
errorRate = comm.ErrorRate();

% 生成发送数据和编码数据
txData = randi([0 1],n,1);
txCode = convEncoder(txData);

% BPSK 调制和解调
txSig_bpsk = bpskMod(txCode);
rxSig_bpsk = channel(txSig_bpsk);
rxCode_bpsk = bpskDemod(rxSig_bpsk);
rxData_bpsk = convDecoder(rxCode_bpsk);
err_bpsk = errorRate(txData,rxData_bpsk);

% QPSK 调制和解调
txSig_qpsk = qpskMod(txCode);
rxSig_qpsk = channel(txSig_qpsk);
rxCode_qpsk = qpskDemod(rxSig_qpsk);
rxData_qpsk = convDecoder(rxCode_qpsk);
err_qpsk = errorRate(txData,rxData_qpsk);

% 输出误差率
disp(['BPSK误码率：',num2str(err_bpsk(1)),'，QPSK误码率：',num2str(err_qpsk(1))]);

上述程序中，我们首先使用 `cyclpoly` 函数生成一个 $N=7$、$K=4$ 的循环码，并将其转换为卷积编码器的参数。然后，我们使用 `comm.BPSKModulator` 和 `comm.QPSKModulator` 分别生成 BPSK 和 QPSK 调制器，以及 `comm.BPSKDemodulator` 和 `comm.QPSKDemodulator` 分别生成 BPSK 和 QPSK 解调器。接着，我们生成一个 AWGN 信道，并使用 `comm.ErrorRate` 生成误差率计算器。

在程序的主体部分，我们首先生成一个随机的比特序列作为发送数据，并将其编码成循环码。然后，我们使用调制器将编码数据调制成 BPSK 和 QPSK 信号，并通过信道发送。接着，我们使用解调器将接收到的信号解调，并使用解码器将解调后的码字解码成比特序列。最后，我们使用误差率计算器计算误码率，并输出结果。

需要注意的是，在上述程序中，我们并没有包含循环码的解码过程。如果需要解码循环码，可以使用 `comm.ViterbiDecoder` 或 `convtraceback` 函数等方法。