%% 系统参数设置% OFDM参数nSubcarrier = 64; % 子载波数目nBitPerSymbol = 4; % 每个符号的比特数nSymbol = 1000; % 发送符号数目nCP = 16; % 循环前缀长度% MIMO参数nTxAntenna = 2; % 发射天线数目nRxAntenna = 2; % 接收天线数目% 信道参数nSNRdB = 10; % 信噪比%% 数据生成txData = randi([0, 1], nTxAntenna, nSubcarrier * nBitPerSymbol * nSymbol);%% 调制txSymbol = qammod(txData, 2^ nBitPerSymbol);%% IFFT变换txSignal = ifft(txSymbol, nSubcarrier, 2);%% 添加循环前缀txSignalPrefix = [txSignal(:, end-nCP+1:end), txSignal];%% 信道传输H = randn(nRxAntenna, nTxAntenna) + 1i * randn(nRxAntenna, nTxAntenna); % 随机信道rxSignalPrefix = H * txSignalPrefix;%% 去除循环前缀rxSignal = rxSignalPrefix(:, nCP+1:end);%% FFT变换rxSymbol = fft(rxSignal, nSubcarrier, 2);%% 解调rxData = qamdemod(rxSymbol, 2^ nBitPerSymbol);%% 比较发送和接收数据ber = sum(sum(xor(txData, rxData))) / (nTxAntenna * nSubcarrier * nBitPerSymbol * nSymbol);%% 显示结果fprintf('误码率：%f\n', ber);在上面的代码中增加预编码，预编码方式为LT编码，接收端的信号检测为MMSE检测，并绘制误比特率图像

好的，以下是增加预编码和接收端的信号检测，绘制误比特率图像的代码：

%% 系统参数设置
% OFDM参数
nSubcarrier = 64; % 子载波数目
nBitPerSymbol = 4; % 每个符号的比特数
nSymbol = 1000; % 发送符号数目
nCP = 16; % 循环前缀长度

% MIMO参数
nTxAntenna = 2; % 发射天线数目
nRxAntenna = 2; % 接收天线数目

% 信道参数
nSNRdB = 10; % 信噪比

%% 数据生成
txData = randi([0, 1], nTxAntenna, nSubcarrier * nBitPerSymbol * nSymbol);

%% 调制
txSymbol = qammod(txData, 2^ nBitPerSymbol);

%% IFFT变换
txSignal = ifft(txSymbol, nSubcarrier, 2);

%% 添加循环前缀
txSignalPrefix = [txSignal(:, end-nCP+1:end), txSignal];

%% 预编码
precodeMatrix = dftmtx(nTxAntenna); % LT编码
txSignalPrecode = precodeMatrix * txSignalPrefix;

%% 信道传输
H = randn(nRxAntenna, nTxAntenna) + 1i * randn(nRxAntenna, nTxAntenna); % 随机信道
rxSignalPrecode = H * txSignalPrecode;

%% 接收端信号检测
Hest = rxSignalPrecode(:, 1:nSymbol) / txSignalPrecode(:, 1:nSymbol); % MMSE估计信道
rxData = zeros(nTxAntenna, nSubcarrier * nBitPerSymbol * nSymbol);
for i = 1:nSymbol
    y = rxSignalPrecode(:, i);
    xHat = (Hest' * Hest + 10^(nSNRdB / 10) * eye(nTxAntenna)) \ Hest' * y;
    rxData(:, (i-1)*nSubcarrier*nBitPerSymbol+1:i*nSubcarrier*nBitPerSymbol) = qamdemod(xHat, 2^ nBitPerSymbol);
end

%% 比较发送和接收数据
ber = sum(sum(xor(txData, rxData))) / (nTxAntenna * nSubcarrier * nBitPerSymbol * nSymbol);

%% 绘制误比特率图像
snrRange = 0:2:20;
berVec = zeros(1, length(snrRange));
for i = 1:length(snrRange)
    nSNRdB = snrRange(i);
    rxData = zeros(nTxAntenna, nSubcarrier * nBitPerSymbol * nSymbol);
    for j = 1:nSymbol
        y = rxSignalPrecode(:, j);
        xHat = (Hest' * Hest + 10^(nSNRdB / 10) * eye(nTxAntenna)) \ Hest' * y;
        rxData(:, (j-1)*nSubcarrier*nBitPerSymbol+1:j*nSubcarrier*nBitPerSymbol) = qamdemod(xHat, 2^ nBitPerSymbol);
    end
    berVec(i) = sum(sum(xor(txData, rxData))) / (nTxAntenna * nSubcarrier * nBitPerSymbol * nSymbol);
end
semilogy(snrRange, berVec, 'o-', 'LineWidth', 2);
grid on;
xlabel('信噪比(dB)');
ylabel('误比特率');
title('2x2 MIMO OFDM系统误比特率曲线');

以上代码中增加了预编码和接收端的信号检测，并绘制了误比特率图像，可以根据需要进行调整和改进。