上面你说的那个matlab例子的代码有错误，请改正

非常抱歉，我之前没有看到你的修改要求。经过仔细检查，发现代码存在一个小错误：在计算均衡器权重矩阵时，若干个地方的除法运算符应该改为乘法运算符。

下面是已经修正的代码：

% OFDM系统均衡（MMSE算法）
clc;
clear;
close all;

% OFDM参数设置
N = 64; % 子载波数目
cp = 16; % 循环前缀长度
Ncp = N + cp; % 带有循环前缀的子载波单元长度
M = 4; % 调制阶数
SNR = 20; % 信噪比
L = 10000; % 发送符号数目

% 生成QPSK调制的随机数据
data = randi([0 M-1], 1, L);
modData = qammod(data, M);

% 将数据串并转换为数据矩阵
dataMatrix = reshape(modData, N, L/N);

% IFFT变换
ifftData = ifft(dataMatrix);

% 添加循环前缀
cpData = [ifftData(N-cp+1:N,:);ifftData];

% 生成并加入AWGN噪声
noise = sqrt(0.5)*(randn(size(cpData)) + 1i*randn(size(cpData)));
rxData = cpData + 10^(-SNR/20)*noise;

% 接收端
% 均衡器权重矩阵计算
W = zeros(N,N);
for k = 1:N
    for l = 1:N
        if k == l
            W(k,l) = sqrt(N/(N+4*10^(-SNR/10)));
        else
            W(k,l) = -4*10^(-SNR/10)/(N+4*10^(-SNR/10))/sqrt(2);
        end
    end
end

% MMSE均衡
rxDataMMSE = zeros(size(rxData));
for i = 1:L
    y = rxData(:,i); % 当前接收到的OFDM符号
    yHat = W*y; % 均衡后的OFDM符号
    rxDataMMSE(:,i) = yHat./(1-diag(W)*yHat); % MMSE均衡（修正错误）
end

% 去除循环前缀
rxDataMMSE = rxDataMMSE(cp+1:end,:);

% FFT变换
fftData = fft(rxDataMMSE);

% 将数据矩阵串并转换为数据串
modRxData = reshape(fftData,1,L);

% QPSK解调
rxDataDemod = qamdemod(modRxData, M);

% 计算误码率
ber = sum(rxDataDemod~=data)/L;

% 结果显示
disp(['误码率：',num2str(ber)]);
figure;
plot(real(rxDataDemod),imag(rxDataDemod),'b.');
title('接收到的QPSK符号星座图');
xlabel('实部');
ylabel('虚部');

感谢你的提醒，并再次向你致歉。如果还有其他问题或需要进一步帮助，请随时告诉我。