错误使用 biterr (第 162 行) Input column vector must contain as many elements as there are rows in the input matrix. 出错 Untitled5 (第 40 行) [~, ber(i)] = biterr(bits, rxBits);

抱歉，我在示例代码中犯了一个错误。biterr函数要求输入的比特向量大小必须与解调后的比特向量大小相同。在多天线接收信号处理部分，我们使用了qamdemod函数对接收信号进行解调，这将返回一个与发送比特序列大小不同的比特序列。

要解决这个问题，我们需要对接收到的比特序列进行重新调整，使其大小与发送比特序列相同。这可以通过添加以下代码来实现：

rxBits = rxBits(:); % 将接收到的比特序列转换为列向量
rxBits = rxBits(1:numBits); % 裁剪接收到的比特序列为与发送比特序列相同的大小

修改后的代码如下所示：

% 设置参数
M = 4; % 调制阶数
K = 2; % 发送天线数
N = 2; % 接收天线数
numBits = 1e5; % 总比特数
SNRdB = -10:2:20; % 信噪比范围
SNR = 10.^(SNRdB/10); % 信噪比转换为线性单位

% 初始化误码率（BER）和均方误差（MSE）
ber = zeros(length(SNR), 1);
mse = zeros(length(SNR), 1);

% 开始循环计算每个信噪比下的BER和MSE
for i = 1:length(SNR)
    % 生成随机比特序列
    bits = randi([0, 1], numBits, 1);
    
    % 调制
    symbols = qammod(bits, M);
    
    % 平均功率归一化
    symbols = symbols / sqrt(mean(abs(symbols).^2));
    
    % 多天线发送信号
    txSignal = reshape(symbols, [], K);
    
    % 添加高斯白噪声
    noise = sqrt(1/(2*SNR(i))) * (randn(size(txSignal)) + 1j*randn(size(txSignal)));
    
    % 接收信号
    rxSignal = txSignal + noise;
    
    % 多天线接收信号处理
    rxSymbols = reshape(rxSignal, [], N);
    
    % 解调
    rxBits = qamdemod(rxSymbols, M);
    
    % 调整接收比特序列大小
    rxBits = rxBits(:); % 将接收到的比特序列转换为列向量
    rxBits = rxBits(1:numBits); % 裁剪接收到的比特序列为与发送比特序列相同的大小
    
    % 计算误码率（BER）
    [~, ber(i)] = biterr(bits, rxBits);
    
    % 计算均方误差（MSE）
    mse(i) = mean(abs(symbols - rxSymbols).^2, 'all');
end

% 绘制BER和MSE随SNR变化的曲线
figure;
semilogy(SNRdB, ber, 'bo-', 'LineWidth', 2);
hold on;
semilogy(SNRdB, mse, 'rx-', 'LineWidth', 2);
grid on;
xlabel('SNR (dB)');
ylabel('BER / MSE');
legend('BER', 'MSE');
title('BER and MSE vs. SNR');

现在，代码应该能正确运行并绘制BER和MSE随SNR变化的曲线。如果还有其他问题，请随时提问！