ZF/MMSE检测的MIMO误码率对比matlab仿真

clc;
clear;
N=1000;
frame=100;
EbN0=（0:5:25）;
N0=zeros（1length（EbN0））;
H=zeros（22N）;
H1=zeros（22N）;
H2=zeros（22N）;
W_ZF=zeros（22N）;
W_MMSE=zeros（22N）;
xzf_hat=zeros（21N）;
xmmse_hat=zeros（21N）;
z1=zeros（N1）;
z2=zeros（N1）;
y1=zeros（N1）;
y2=zeros（N1）;
y=zeros（21N）;
x1zf_hat=zeros（N1）;
x2zf_hat=zeros（N1）;
x1zfhat=zeros（N1）;
x2zfhat=zeros（N1）;
x1mmse_hat=zeros（N1）;
x2mmse_hat=zeros（N1）;
x1mmsehat=zeros（N1）;
x2mmsehat=zeros（N1）;
num1ZF（6）=0;
num2ZF（6）=0;
BER_ZF（6）=0;
num1MMSE（6）=0;
num2MMSE（6）=0;
BER_MMSE（6）=0;
h11=（randn（1N）+1j*randn（1N））/sqrt（2）;
h12=（randn（1N）+1j*randn（1N））/sqrt（2）;
h21=（randn（1N）+1j*randn（1N））/sqrt（2）;
h22=（randn（1N）+1j*randn（1N））/sqrt（2）;
x1=round（rand（N1））;
x2=round（rand（N1））;
u1=2*x1-1;
u2=2*x2-1;
for t=1:N
    H（::t）=[h11（t）h12（t）;h21（t）h22（t）];
end

for i=1:length（EbN0）
    berZF=0;
    berMMSE=0;
    berML=0;
    N0（i）=1/（10^（（EbN0（i））/10））;
    for m=1:frame
        z1=sqrt（N0（i）/（2））*randn（size（x1））;
        z2=sqrt（N0（i）/（2））*randn（size（x2））;
        y1=u1.*（h11）.‘+u2.*（h12）.‘+z1;
        y2=u1.*（h21）.‘+u2.*（h22）.‘+z2;
        for q=1:N
            y（::q）=[y1（q）;y2（q）];
        end
        %ZF================================================================
        for n=1:N
            W_ZF（::n）= inv（conj（H（::n）.‘）*H（::n））*（conj（H（::n）.‘））;
        end
        for p=1:N
            xzf_hat（::p）=W_ZF（::p）*y（::p）;
        end
        for o=1:N
            x1zf_hat（o）=xzf_hat（11o）;
            x2zf_ha