矩阵维度必须一致。 出错 Untitled6 (第 47 行) mse(i) = mean(abs(symbols - rxSymbols).^2, 'all');

非常抱歉，我之前的回答中存在错误。在计算均方误差（MSE）时，由于接收到的符号和发送的符号大小不匹配，会导致矩阵维度不一致的错误。

要解决这个问题，我们可以修改代码以确保接收到的符号大小与发送的符号大小相匹配。可以通过以下方式修改代码：

mse(i) = mean(abs(symbols - rxSymbols(:)).^2, 'all');

修改后的完整代码如下所示：

% 设置参数
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);
    
    % 计算误码率（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变化的曲线。如果还有其他问题，请随时提问！