【正文】 ,39。 end % average over channel realizations Psrqostbc(SNRidx) = mean(PPsrqostbc) Psqostbc(SNRidx) = mean(PPsqostbc) Pstqostbc(SNRidx) = mean(PPstqostbc) end figure(1)semilogy(SNRdBvalues, Pstqostbc,39。 SNRstqostbc = SNR/4/R1(1,1)。 R1 = inv(H39。]。 h239。 h439。 h339。 h139。 %%%% QuasiOrthogonal STBC % equivalent channel HH_qo = [h1 h2 h3 h4。S(1) 1]。S(1) S(1) S(1)1]。 S(1) 1 S(1) S(1)。 % BER for uncoded QPSK input PPsqostbc(kk) = *erfc(sqrt(*SNRsqostbc))。*H )。 % putation of postprocessing SNR for stream 1 H = S_r^(.5)*HH_qo*S_t^(.5)。 h139。 h339。 h3 h4 h1 h2。 h439。 h239。S(1) S(1) S(1) 1]。 S(1) 1 S(1) S(1)。S(1) 1]。S(1) S(1) S(1) 1]。 S(1) 1 S(1) S(1)。 % BER for uncoded QPSK input PPsrqostbc(kk) = *erfc(sqrt(*SNRsrqostbc))。*H )。 % putation of postprocessing SNR for stream 1 H = S_r^(.5)*HH_qo*S_t^(.5)。h139。 h339。 h3 h4 h1 h2。 h439。 h239。S(1) S(1) S(1)1]。 S(1) 1 S(1) S(1)。S(1) 1]。 % subchannels S_t = eye(4)。 h3 = H(1,3)。 % 1x4 MISO channel h1 = H(1,1)。 % SNR in linear scale S = []。for SNRdB=SNRdBvalues SNRdB SNRidx = SNRidx + 1。PPstqostbc = zeros(1,It)。%gai% variable initialization: BER for each channel iterationPPsrqostbc = zeros(1,It)。%gai % SNR range in dB%SNRdBvalues = [0:2:30]。,3)gridhold off% 圖44和圖45 clear all % number of channel realizations. This number needs to be increased when% estimating small probabilities%It = 50000。,39。,39。)legend(39。)ylabel(39。)xlabel(39。,1)axis([SNRdBvalues(1) SNRdBvalues(end) 10e8 1])title(39。,39。,1)semilogy(SNRdBvalues, Pqostbc,39。,39。,1)hold onsemilogy(SNRdBvalues, Psqostbc,39。,39。 end % average over channel realizations Pqostbc(SNRidx) = mean(PPqostbc) Psqostbc(SNRidx) = mean(PPsqostbc) Psiso(SNRidx) = mean(PPsiso) end figure(1)semilogy(SNRdBvalues, Psiso,39。 SNRsiso = SNR/4/R1(1,1)。 R1 = inv(H39。]。 h239。 h439。 h339。 h139。 %%%% QuasiOrthogonal STBC % equivalent channel HH_qo = [h1 h2 h3 h4。S(3) S(3) 1 S(3)。 S_r = [1 S(3) S(3) S(3)。 %S_r = [1 S(1)。S(3) S(3) 1 S(3)。 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%% SISO S_t = [1 S(3) S(3) S(3)。 SNRsqostbc = SNR/4/R1(1,1)。 R1 = inv(H39。]。 h239。 h439。 h339。 h139。 %%%% QuasiOrthogonal STBC % equivalent channel HH_qo = [h1 h2 h3 h4。S(2) S(2) 1 S(2)。 S_r = [1 S(2) S(2) S(2)。 % S_r = [1 S(1)。S(2) S(2) 1 S(2)。 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% S_t = [1 S(2) S(2) S(2)。 SNRqostbc = SNR/4/R1(1,1)。 R1 = inv(H39。]。 h239。 h439。 h339。 h139。 %%%% QuasiOrthogonal STBC % equivalent channel HH_qo = [h1 h2 h3 h4。S(1) S(1) 1 S(1)。 S_r = [1 S(1) S(1) S(1)。 % S_r = [1 S(1)。S(1) S(1) 1 S(1)。 % subchannels % generate channels S_t = [1 S(1) S(1) S(1)。 h3 = H(1,3)。 % 1x4 MISO channel h1 = H(1,1)。 % SNR in linear scale S = [ ]。for SNRdB=SNRdBvalues SNRdB SNRidx = SNRidx + 1。PPsiso = zeros(1,It)。%gai% variable initialization: BER for each channel iterationPPqostbc = zeros(1,It)。%gai % SNR range in dB%SNRdBvalues = [0:2:30]。,3)gridhold off% 圖43 clear all % number of channel realizations. This number needs to be increased when% estimating small probabilities%It = 50000。,39。,39。)legend(39。)ylabel(39。)xlabel(39。,1)axis([SNRdBvalues(1) SNRdBvalues(end) 10e8 1])title(39。,39。,1)semilogy(SNRdBvalues, Pqostbc,39。,39。,1)hold onsemilogy(SNRdBvalues, Psqostbc,39。,39。 end % average over channel realizations Pqostbc(SNRidx) = mean(PPqostbc) Psqostbc(SNRidx) = mean(PPsqostbc) Psiso(SNRidx) = mean(PPsiso) end figure(1)semilogy(SNRdBvalues, Psiso,39。 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%% SISO SNRsiso = SNR*norm(h1)^2。 SNRsqostbc = SNR/4/R1(1,1)。 R1 = inv(H39。]。 h239。 h439。 h339。 h139。 %%%% QuasiOrthogonal STBC % equivalent channel HH_qo = [h1 h2 h3 h4。S(1) S(1) 1 S(1)。 S_r = [1 S(1) S(1) S(1)。 % S_r = [1 S(1)。S(1) S(1) 1 S(1)。 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% S_t = [1 S(1) S(1) S(1)。 SNRqostbc = SNR/4/R(1,1)。 % putation of postprocessing SNR for stream 1 R = inv(HH_o39。 h139。 h339。 h3 h4 h1 h2。 h439。 h239。 h4= H(1,4)。 h2 = H(1,2)。 % iterations on channel realizations for kk=1:It % generate channels H = (randn(1,4) + j * randn(1,4)) / sqrt(2)。 % SNR index SNR = 10^(SNRdB/10)。SNRidx = 0。PPsqostbc = zeros(1,It)。 SNRdBvalues = [0:2:16]。It = 1000。正交空時(shí)碼39。TBH準(zhǔn)正交空時(shí)碼39。單入單出39。BER39。Average SNR \rho (dB)39。Bit error rate for a QPSK input . Rayleigh slow fading channel39。linewidth39。g*39。linewidth39。bs39。linewidth39。r+39。 % BER for uncoded QPSK input PPsiso(kk) = *erfc(sqrt(*SNRsiso))。 % BER for uncoded QPSK input PPqostbc(kk) = *erfc(sqrt(*SNRqostbc))。*HH_qo )。]。 h239。 h439。 h339。 h139。 %%%% QuasiOrthogonal STBC % equivalent channel HH_qo = [h1 h2 h3 h4。 SNRostbc = SNR/4/R(