MIMO-OFDM.zip_MIMO+OFDM_MIMO+OFDM+Matlab仿真_MIMO-OFDM_mimo_mimo o

Loop=1000;%循环次数 j=sqrt(-1);%虚数 N_Block=64;%子载波数 M_t=2;%发送天线数 R_t=2;%接收天线数 CP=8;%CP长度 SNR_dB=0:2:20;%信噪比 BER=zeros(1,length(SNR_dB)); Noise_Sigma=10.^(-SNR_dB./10); for ii=1:1:length(SNR_dB) for jj=1:1:Loop G=sqrt(0.5)*(randn(CP*R_t,M_t)+j*randn(CP*R_t,M_t)); H=zeros(N_Block*R_t,M_t); U=zeros(N_Block,R_t); V=zeros(R_t,N_Block); for m=1:1:N_Block h=zeros(2); for l=1:1:CP h=h+G((l-1)*R_t+1:l*R_t,:)*exp(-(j*2*pi*(m-1)*(l-1)/N_Block)); end H((m-1)*R_t+1:m*R_t,:)=h; Obj_fun1=h; [u,Eig_value1,Eig_vector1] =svd(Obj_fun1);%产生特征向量和特征值； [~, max_colum] = find( Eig_value1 == max(max(Eig_value1)) ); V(:,m)=Eig_vector1(:,max_colum); U(m,:)=u(:,max_colum)'; end Noise=sqrt(0.5*Noise_Sigma(ii))*(randn(R_t,(N_Block+2*CP-1))+1i*randn(R_t,(N_Block+CP*2-1))); B_stream=randi([0 1],1,N_Block); %%——————————BPSK编码——————————————% QAM_order =2; %QAM调制阶数:2 QAM_Es = sum(abs(qammod(0:QAM_order-1,QAM_order)).^2)/QAM_order; S_stream=qammod(B_stream,2)/sqrt(QAM_Es); %——————————————————————————————%%% X=V*diag(S_stream); X1=X(1,:); X2=X(2,:); X1_ifft=ifft(X1,N_Block); X2_ifft=ifft(X2,N_Block); X1temp=[X1_ifft(N_Block-CP+1:N_Block) X1_ifft]; X2temp=[X2_ifft(N_Block-CP+1:N_Block) X2_ifft]; G11=zeros(1,CP); G12=zeros(1,CP); G21=zeros(1,CP); G22=zeros(1,CP); for n=1:1:CP G11(n)=G((n-1)*2+1,1); G12(n)=G((n-1)*2+1,2); G21(n)=G((n-0)*2,1); G22(n)=G((n-0)*2,2); end Y1temp=conv(X1temp,G11)+conv(X2temp,G12)+Noise(1,:); Y1fft=Y1temp(CP+1:N_Block+CP); Y1=fft(Y1fft,N_Block); Y2temp=conv(X1temp,G21)+conv(X2temp,G22)+Noise(2,:); Y2fft=Y2temp(CP+1:N_Block+CP); Y2=fft(Y2fft,N_Block); Y=[Y1;Y2]; S_deqam=diag(U*Y)*sqrt(QAM_Es); B_bit=qamdemod(S_deqam,2); error=find((B_bit'-B_stream)==0); BER(ii)=BER(ii)+length(error); end BER(ii)=BER(ii)/(Loop*N_Block); end semilogy(SNR_dB,1-BER)