MATLAB实现OFDM MIMO技术的仿真研究

clear all;
close all;                              %初始化
N=2048;                                 %FFT点数
SINR=10;                                %信噪比
CP=N/4;                                 %CP长度
GI=1/32*(CP+N);                         %保护间隔长度
carrier_count=300;                      %LTE在5M带宽配置下，15k间隔的载波数
SymbolsPerCarier=12;                    %每个子载波包含的OFDM符号数
bitPerSymbol=4;                         %每个OFDM符号携带的二进制信号bit数，等于log(16),16为采用16QAM
disp('************************仿真开始*****************')
carriers =(1:carrier_count) +(floor(N/4)-floor(carrier_count/2));           %共轭对称子载波映射  复数数据对应的IFFT点坐标,有用的FFT点数（信道数）
conjugate_carriers =N - carriers + 2;                                       %共轭对称子载波映射  共轭复数对应的IFFT点坐标
source=randint(1,carrier_count*SymbolsPerCarier*bitPerSymbol);              %以下为待传送二进制码流生成过程

modSource=qam16(source);                                                    %16QAM调制（符号映射）
figure(1)
subplot(221)
plot(modSource,'*k')
grid on;
axis([-4 4 -4 4])
title('传送前的直角坐标系星座图')
figure(1)
subplot(222)
polar(angle(modSource),abs(modSource),'kd');
title('传送前的极坐标系星座图')
modSource=reshape(modSource',carrier_count,SymbolsPerCarier)';               %变换矩阵，行代表子载波数所含符号数，列代表子载波所包含的的符号数

preIFFT=zeros(SymbolsPerCarier,N);
preIFFT(:,carriers ) = modSource ;                                           %代替导频信号映射点
preIFFT(:,conjugate_carriers ) = conj(modSource);                            %共轭复数映射

IFFTmodSource=ifft(preIFFT,N,2);                                              %时域波形矩阵，行为每载波所含符号数，列ITTF点数，N个子载波映射在其内，每一行即为一个OFDM符号

%***************添加循环前缀CP******************%

preCPadd=zeros(SymbolsPerCarier,N+CP+GI);
for k=1:SymbolsPerCarier
    for i=1:N
       preCPadd(k,i+CP)=IFFTmodSource(k,i);
    end
    for i=1:CP
        preCPadd(k,i)=IFFTmodSource(k,i+N-CP);                              %添加循环前缀
    end
     for j=1:GI
        preCPadd(k,N+CP+j)=IFFTmodSource(k,j);%添加循环后缀
    end
end

CPadd=preCPadd;                                                             %加入CP后的信号

%***********************加窗(带通滤波)并串转换***************************%

tr=zeros(1,N+CP+GI);
for i = 1:SymbolsPerCarier 
tr(i,:) = real(CPadd(i,:)).*rcoswindow(1/32,N+CP)';                         %添加一个滚降系数为1/32的
end  

Tx=zeros(1,SymbolsPerCarier*(N+CP)+GI);
Tx(1:N+CP+GI)=tr(1,:);
for i = 1:SymbolsPerCarier-1 ;
    Tx((N+CP)*i+1:(N+CP)*(i+1)+GI)=tr(i+1,:);