MIMO OFDM MATLAB仿真程序续一.doc

modulate.m: Modulate the random input sequence according to the bit allocations for each channel. function y = modulate(x,b,e, s2,s4,s16,s64,s256) % function y = modulate(x,b,e, s2,s4,s16,s64,s256) % % y - modulated output, in the form of a row vector % x - a vector of input bits, for all the subcarriers (row vector) % b - subcarrier bit allocation (64 elements in this matrix, each one % corresponding to the number of bits to be allocated to the subcarrier % having the same index) % e - subcarrier energy allocation (64 elements in this matrix, each one % corresponding to the energy to be allocated to the subcarrier % having the same index) % s_ - the encoder for a given constellation size y=[]; b2 = zeros(1,length(b)); b2(1) = 1; for i = 1:length(b) b2(i+1) = b(i) + b2(i); end for i = 1:length(b) switch b(i) case {1} y = [y s2(x(b2(i))+1)*sqrt(e(i))]; case {2} y = [y s4(x(b2(i):(b2(i+1)-1))*[2;1]+1)*sqrt(e(i))]; case {4} y = [y s16(x(b2(i):(b2(i+1)-1))*[8;4;2;1]+1)*sqrt(e(i))]; case {6} y = [y s64(x(b2(i):(b2(i+1)-1))*[32;16;8;4;2;1]+1)*sqrt(e(i))]; case {8} y = [y s256(x(b2(i):(b2(i+1)-1))*[128;64;32;16;8;4;2;1]+1)*sqrt(e(i))]; otherwise % must be zero bits allocated y = [y 0]; end end ? ENC2.mat: BPSK Modulator ? ENC4.mat: 4-QAM Modulator (Gray coded) ? ENC16.mat: 16-QAM Modulator (Gray coded) ? ENC64.mat: 64-QAM Modulator (Gray coded) ? ENC256.mat: 256-QAM Modulator (Gray coded) precode.m: Precode the transmitted vector at each time instance by filtering the modulated vector with the right-inverse of the channels right singluar matrix. function [x_pre] = precode(Mt, x_mod, V, N); % [x_pre] = precode(Mt, x_mod, V, N); % % Function precodes the modulated symbols prior to sending through the IFFT module. The precoder % left-multiplies the channel at each tone, to eliminate the right unitary matrix in the SVD % decomp