untitled4.m

这些代码包括ofdm的实现和信道估计,可以运行,有指导程序的说明文档.

%%%%%%%%%%%%% mmse initialization %%%%%%%%%%%

clc;
clear all;
fftlen=64;
gilen=16;
fftlen1=52;
ebn0=3;
L=gilen;
bita=1;% the number of qpsk 
SNR=exp(ebn0*log(10)/10);
tao_rms=L/4;
for m=1:fftlen1
   for n=1:fftlen1
  %      R_HH(m,n)=(1-exp(-L*((1/tao_rms)+2*pi*i*(m-n)/fftlen)))/(tao_rms*(1-exp(-(L/tao_rms)))*(1/tao_rms+2*pi*i*(m-n)/fftlen));
        if(m==n) 
            R_HH(m,n)=1;
            R_HH1(m,n)=1;
        else 
        R_HH(m,n)=(1-exp(-2*pi*i*L*(m-n)/fftlen))/(2*pi*i*L*(m-n)/fftlen);
        R_HH1(m,n)=(1-exp(-2*pi*i*(m-n)/fftlen))/(2*pi*i*L*(m-n)/fftlen);
    end
end
end
B_inter=zeros(52,52);
for i=1:52
    B_inter(i,i)=bita/SNR;
end
B_inter1=R_HH+B_inter;
B_inter11=R_HH1+B_inter;
B_inter2=inv(B_inter1);
B_inter21=inv(B_inter11);
B_inter3=R_HH*B_inter2;

B_inter31=R_HH1*B_inter21;
hls=ones(52,1);
H=B_inter3*hls;
H1=B_inter31*hls;