基于matlab开发的TD-LTE 物理层解传输分集预编码matlaB仿真程序.rar

clc; clear all; close all; %本次仿真带宽为10M,子载波个数为600个，即频域上展50个资源块，子载波间隔deltaf=15kHz。 %采用QPSK调制 %下面主要是主要是数据经过调制，层映射，预编码后通过信道，下面主要是接传输分集的预编码 nRx=2; Eb_N0_dB = [-15:15]; nErr1(26)=0; M=8400; %产生数据的长度 N=4200; %经过调制后数据的长度 bit_sourse = randint(8400,1); %随机产生8400个符号，即8400个比特 for Y=1:10 for ii = 1:length(Eb_N0_dB) symbol_out = mod_qpsk(bit_sourse); %QPSK调制,4200*1的矩阵 symbol_layer = layermap(symbol_out); %层映射 [sym_precode1,sym_precode2] = precode(symbol_layer); %预编码 h = 1/sqrt(2)*[randn(nRx,N) + j*randn(nRx,N)]; % Rayleigh channel 假设是瑞利信道 n = 1/sqrt(2)*[randn(nRx,N) + j*randn(nRx,N)]; % white gaussian noise, 0dB variance sym_precode11=sym_precode1'; sym_precode22=sym_precode2'; sCode2(1,:)=sym_precode11; %送入天线0的数据 sCode2(2,:)=sym_precode22; %送入天线1的数据 %************************************是针对两个天线发送，两个天线接收********************** %在经过传输预编码后，两路数据送入天线端口，在这里必须考虑是在两个子载波的信道上的冲激相应是相同的，在一个子载波上 %接收端的一根天线接收的数据在这个子载波上要接收来自发送端的两根天线在这个子载波上的数据。 y = zeros(nRx,N); yMod = zeros(nRx*2,N); hMod = zeros(nRx*2,N); for kk = 1:nRx hMod = kron(reshape(h(kk,:),2,N/2),ones(1,2));%针对一根接收天线，发送端的两根天线到一根接收天线的信道冲激相应，并且相应的两个子载波冲激相应相同 % hMod=kron(reshape(h(kk,:),2,N/2)),ones(1,2)); temp=hMod; %hMod(2,[1:2:end]) = -(temp(2,[1:2:end])); %hMod(1,[2:2:end]) = conj(temp(2,[2:2:end])); %hMod(2,[2:2:end]) = conj(temp(1,[2:2:end])); % Channel and noise Noise addition y(kk,:) = sum(hMod.*sCode2,1) + 10^(-Eb_N0_dB(ii)/20)*n(kk,:); %接收的信号加上噪声 % Receiver %根据接收的算法，把接收端的每根天线的第二路数据取共轭 yMod([2*kk-1:2*kk],:) = kron(reshape(y(kk,:),2,N/2),ones(1,2)); yMod(2*kk,:)=conj(yMod(2*kk,:)); % forming the equalization matrix %接收端要解出相应的数据，把信道冲激相应做相应的调整 hEq([2*kk-1:2*kk],:) = hMod; hEq(2*kk-1,[1:2:end]) = conj(hEq(2*kk-1,[1:2:end])); c = (hEq(2*kk-1,[2:2:end])); hEq(2*kk-1,[2:2:end]) = -conj(hEq(2*kk,[2:2:end])); hEq(2*kk,[2:2:end]) = c; end hEqPower = sum(hEq.*conj(hEq),1);%计算信道冲激相应的模 Re_Mod = sum(hEq.*yMod,1)./hEqPower; % [h1*y1 + h2y2*, h2*y1 -h1y2*, ... ]把两根解出的每一个子载波数据直接相加 Re_Mod(2:2:end) = conj(Re_Mod(2:2:end));%把解出的偶数位数据取共轭即解出相应的数据,解完层映射后的数据 bit_receive = demod_qpsk(Re_Mod);%解调后的数据 sym_precode2=sym_precode11; bit_sourse2 = demod_qpsk(sym_precode2); nErr(ii) = size(find([ bit_receive- bit_sourse2]),2);%计算误比特率 end nErr1= nErr +nErr1; end simBer = nErr/(200*N); % simulated ber figure semilogy(Eb_N0_dB,simBer,'mo-','LineWidth',2); axis([-15 15 10^-6 0.5]) grid on xlabel('Eb/No, dB'); ylabel('Bit Error Rate'); title('BER for QPSK modulation with 2Tx, 2Rx Alamouti SFBC (Rayleigh channel)');