MIMO-OFDM仿真.rar_MIMO-OFDM_MQP2_mimo ofdm_mimo仿真_ofdm

NgType=1; % 对于CP或ZP，NgType=1或2 if NgType==1 nt='CP'; elseif NgType==2 nt='ZP'; end Ch=0; % 对于AWGN/多径瑞利信道，channelCh=0/1 if Ch==0 chType='AWGN'; Target_neb=100; else chType='CH'; Target_neb=500; end PowerdB=[0 -8 -17 -21 -25]; % 信道抽头功率特性（dB） Delay=[0 3 5 6 8]; % 信道时延（采样点） Power=10.^(PowerdB/10); % 信道抽头功率特性（线性尺度） Ntap=length(PowerdB); % 信道抽头数 Lch=Delay(end)+1; % 信道长度 Nbps=4; M=2^Nbps; % 调制阶数=2/4/6：QPSK/16-QAM/64-QAM Nfft=64; % FFT大小 Ng=3; % GI(保护间隔)长度，没保护间隔时，Ng=0 Nsym=Nfft+Ng; % 符号周期 Nvc=Nfft/4; % Nvc=0 ：没有VC(虚拟载波) Nused=Nfft-Nvc; EbN0=(0:3:30); % Eb/N0 N_iter=1e4; % 对于每一个EbN0的迭代次数 Nframe=3; % 每一帧的符号数 sigPow=0; % 初始信号功率 file_name=['OFDM_BER_' chType '_' nt '_' 'GL' num2str(Ng) '.dat']; fid=fopen(file_name, 'w+'); norms=[1 sqrt(2) 0 sqrt(10) 0 sqrt(42)]; % BPSK 4-QAM 16-QAM for i=0:length(EbN0) randn('state',0);rand('state',0); %Ber2=ber(); 初始化BER Neb=0; % 初始化错误比特数 Ntb=0; % 初始化总比特数 for m=1:N_iter % Tx_______________________________________________________________________________ X= randi([0 M-1],1,Nused*Nframe); % bit: 整数向量(问题：randint怎么换成randi) Xmod= qammod(X,M,'gray')/norms(Nbps); if NgType~=2 x_GI=zeros(1,Nframe*Nsym); elseif NgType==2 x_GI= zeros(1,Nframe*Nsym+Ng); % 用Ng个零扩展OFDM符号 end kk1=(1:Nused/2); kk2=(Nused/2+1:Nused); kk3=1:Nfft; kk4=1:Nsym; for k=1:Nframe if Nvc~=0 X_shift= [0 Xmod(kk2) zeros(1,Nvc-1) Xmod(kk1)]; else X_shift= [Xmod(kk2) Xmod(kk1)]; end x= ifft(X_shift); x_GI(kk4)= guard_interval(Ng,Nfft,NgType,x); kk1=kk1+Nused; kk2= kk2+Nused; kk3=kk3+Nfft; kk4=kk4+Nsym; end if Ch==0 y= x_GI; % 没有信道 else % 多径衰落信道 channel=(randn(1,Ntap)+1i*randn(1,Ntap)).*sqrt(Power/2); h=zeros(1,Lch); %信道脉冲响应 h(Delay+1)=channel; y = conv(x_GI,h); end if i==0 % 只测量增加的AWGN噪声的信号功率 y1=y(1:Nframe*Nsym); sigPow = sigPow + y1*y1'; continue; end % 加AWGN噪声________________________________________________ snr = EbN0(i)+10*log10(Nbps*(Nused/Nfft)); % SNR vs. Eb/N0 by Eq.(4.28) noise_mag = sqrt((10.^(-snr/10))*sigPow/2); y_GI = y + noise_mag*(randn(size(y))+1i*randn(size(y))); % Rx_____________________________________________________________ kk1=(NgType==2)*Ng+(1:Nsym); kk2=1:Nfft; kk3=1:Nused; kk4=Nused/2+Nvc+1:Nfft; kk5=(Nvc~=0)+(1:Nused/2); if Ch==1 H= fft([h zeros(1,Nfft-Lch)]); % 信道频率响应 H_shift(kk3)= [H(kk4) H(kk5)]; end for k=1:Nframe Y(kk2)= fft(remove_GI(Ng,Nsym,NgType,y_GI(kk1))); Y_shift=[Y(kk4) Y(kk5)]; if Ch==0 Xmod_r(kk3) = Y_shift; else Xmod_r(kk3)=Y_shift./H_shift; % 均衡器 end kk1=kk1+Nsym; kk2=kk2+Nfft; kk3=kk3+Nused; kk4=kk4+Nfft; kk5=kk5+Nfft; end X_r=qamdemod(Xmod_r*norms(Nbps),M,'gray'); Neb=Neb+sum(sum(de2bi(X_r,Nbps)~=de2bi(X,Nbps))); % 判决 Ntb=Ntb+Nused*Nframe*Nbps; %[Ber,Neb,Ntb]=ber(bit_Rx,bit,Nbps); if Neb>Target_neb,break;end end if i==0 sigPow= sigPow/Nsym/Nframe/N_iter; else Ber = Neb/Ntb; fprintf('EbN0=%3d[dB], BER=%4d/%8d =%11.3e\n', EbN0(i), Neb,Ntb,Ber) fprintf(fid, '%d\t%11.3e\n', EbN0(i), Ber); if Ber<1e-6, break; end end end if (fid~=0) fclose(fid); end disp('Simulation is finished') ber_AWGN = ber_QAM(EbN0,M,'AWGN'); ber_Rayleigh = ber_QAM(EbN0,M,'Rayleigh'); semilogy(EbN0,ber_AWGN,'r:'),hold on semilogy(EbN0,ber_Rayleigh,'m-') a= load(file_name); semilogy(a(:,1),a(:,2),'b--d'), grid on xlabel('EbN0[dB]'), ylabel('BER') axis([a(1,1) 30 1e-4 1]) Ch=1; % 对于AWGN/多径瑞利信道，channelCh=0/1 file_name=['OFDM_BER_' chType '_' nt '_' 'GL' num2str(Ng) '.dat']; fid=fopen(file_name, 'w+'); for i=0:length(EbN0) randn('state',0);rand('state',0); %Ber2=ber(); 初始化BER Neb=0; % 初始化错误比特数 Ntb=0; % 初始化总比特数 for m=1:N_iter % Tx______________________________________________________________ X= randi([0 M-1],1,Nused*Nframe); % bit: 整数向量(问题：randint怎么换成randi) Xmod= qammod(X,M,'gray')/norms(Nbps); if NgType~=2 x_GI=zeros(1,Nframe*Nsym); elseif NgType==2 x_GI= zeros(1,Nframe*Nsym+Ng); % 用Ng个零扩展OFDM符号 end kk1=(1:Nused/2); kk2=(Nused/2+1:Nused); kk3=1:Nfft; kk4=1:Nsym; for k=1:Nframe if Nvc~=0 X_shift= [0 Xmod(kk2) zeros(1,Nvc-1) Xmod(kk1)]; else X_shift= [Xmod(kk2) Xmod(kk1)]; end x= ifft(X_shift); x_GI(kk4)= guard_interval(Ng,Nfft,NgType,x); kk1=kk1+Nused; kk2= kk2+Nused; kk3=kk3+Nfft; kk4=kk4+Nsym; end if Ch==0 y= x_GI; % 没有信道 else % 多径衰落信道 channel=(randn(1,Ntap)+1i*randn(1,Ntap)).*sqrt(Power/2); h=zeros(1,Lch); %信道脉冲响应 h(Delay+1)=channel; y = conv(x_GI,h); end if i==0 % 只测量增加的AWGN噪声的信号功率 y1=y(1:Nframe*Nsym); sigPow = sigPow + y1*y1'; continue; end % 加AWGN噪声________________________________________________ snr = EbN0(i)+10*log10(Nbps*(Nused/Nfft)); % SNR vs. Eb/N0 by Eq.(4.28) noise_mag = sqrt((10.^(-snr/10))*sigPow/2); y_GI = y + noise_mag*(randn(size(y))+1i*randn(size(y))); % Rx_____________________________________________________________ kk1=(NgType==2)*Ng+(1:Nsym); kk2=1:Nfft; kk3=1:Nused; kk4=Nused/2+Nvc+1:Nfft; kk5=(Nvc~=0)+(1:Nused/2); if Ch==1 H= fft([h zeros(1,Nfft-Lch)]); % 信道频率响应 H_shift(kk3)= [H(kk4) H(kk5)]; end for k=1:Nframe Y(kk2)= fft(remove_GI(Ng,Nsym,NgType,y_GI(kk1))); Y_shift=[Y(kk4) Y(kk5)]; if Ch==0 Xmod_r(kk3) = Y_shift; else Xmod_r(kk3)=Y_shift./H_shift; % 均衡器 end kk1=kk1+Nsym; kk2=kk2+Nfft; kk3=kk3+Nused; kk4=kk4+Nfft; kk5=kk5+Nfft; end X_r=qamdemod(Xmod_r*norms(Nbps),M,'gray'); Neb=Neb+sum(sum(de2bi(X_r,Nbps)~=de2bi(X,Nbps))); % 判决 Ntb=Ntb+Nused*Nframe*Nbps; %[Ber,Neb,Ntb]=ber(bit_Rx,bit,Nbps); if Neb>Target_neb,break;end end if i==0 sigPow= sigPow/Nsym/Nframe/N_iter; else Ber = Neb/Ntb; fprintf('EbN0=%3d[dB], BER=%4d/%8d =%11.3e\n', EbN0(i), Neb,Ntb,Ber) fprintf(fid, '%d\t%11.3e\n', EbN0(i), Ber); if Ber<1e-6, break; end end end if (fid~=0) fclose(fid); end disp('Simulation is finished') a= load(file_name); semilogy(a(:,1),a(:,2),'g--*'), grid on NgType=2; % 对于CP或ZP，NgType=1或2 file_name=['OFDM_BER_' chType '_' nt '_' 'GL' num2str(Ng) '.dat']; fid=fopen(file_name, 'w+'); for i=0:length(EbN0) randn('state',0);rand('state',0); %Ber2=ber