MIMO-OFDM系统仿真.docVIP

clear all; clc; frmLen = 300; % frame length帧长 numPackets = 10; % number of packets帧数 Nsamp=frmLen*numPackets; EbNo = 2:2:20; % Eb/No varying to 20 dB Nt = 4; % maximum number of Tx antennas Nr = 2; % maximum number of Rx antennas N_FFT=400; Nc=N_FFT; hStr = RandStream(mt19937ar, Seed, 55408); % Create QPSK mod-demod objects P = 4; % modulation order qpskmod = modem.pskmod(M, P, SymbolOrder, Gray); qpskdemod = modem.pskdemod(qpskmod); %% Create comm.OSTBCEncoder and comm.OSTBCCombiner System objects hOSTBCEnc = comm.OSTBCEncoder(... NumTransmitAntennas, Nt); hOSTBCComb = comm.OSTBCCombiner(... NumTransmitAntennas, Nt,... NumReceiveAntennas, Nr); %% rho = 0.7; % Correlation coefficient = antenna correlation Rsym = 10e3; % Input symbol rate Rbit = Rsym * log2(P); % Input bit rate Nos = 4; % Oversampling factor ts = (1/Rbit) / Nos; % Input sample period 12.5us fs=Rbit*Nos; % tau = [0 0.4 0.9]*1e-6; % Path delays, in seconds % pdb = [0 -15 -20]; % Average path gains, in dB % dop = doppler.rounded; % Doppler spectrum, with default parameters fd = 0.5; % h = mimochan(Nt, Nr, ts, fd, tau, pdb); % MIMO channel object % h.KFactor = 4; % Rician K-factor on first path % h.DopplerSpectrum = dop; % Doppler spectrum of MIMO object TxCorrelationMatrix = toeplitz([1 rho rho^2 rho^3]); % Transmit correlation matrix RxCorrelationMatrix =toeplitz([1 rho]); delays=[0, 5e-5, 10e-5]; delay_samp=round(delays./ts); %% %% Create MIMO channel System object hChan = comm.MIMOChannel(... SampleRate, fs,... MaximumDopplerShift, fd,... NumTransmitAntennas, Nt,... NumReceiveAntennas, Nr,... TransmitCorrelationMatrix, TxCorrelationMatrix,... ReceiveCorrelationMatrix, RxCorrel