Scribd
Upload
47 views

New Alamouti

The document describes a simulation that compares the bit error rate (BER) performance of an Alamouti space-time block coded system with 2 transmit antennas and 1 receive antenna to an uncoded system with 1 transmit and 1 receive antenna. It generates random data, encodes it using Alamouti coding, transmits over an AWGN channel, decodes the received data, and calculates the BER for each configuration over a range of Eb/No values. The BER results are plotted on a graph to compare the performance of the different systems.

Uploaded by

Preetam Polaki
Copyright
© © All Rights Reserved
Available Formats
Download as TXT, PDF, TXT or read online on Scribd
47 views

New Alamouti

The document describes a simulation that compares the bit error rate (BER) performance of an Alamouti space-time block coded system with 2 transmit antennas and 1 receive antenna to an uncoded system with 1 transmit and 1 receive antenna. It generates random data, encodes it using Alamouti coding, transmits over an AWGN channel, decodes the received data, and calculates the BER for each configuration over a range of Eb/No values. The BER results are plotted on a graph to compare the performance of the different systems.

Uploaded by

Preetam Polaki
Copyright
© © All Rights Reserved
Available Formats
Download as TXT, PDF, TXT or read online on Scribd
You are on page 1/ 2

frmLen = 100; % frame length

numPackets = 1000; % number of packets


EbNo = 0:2:20; % Eb/No varying to 20 dB
N = 2; % maximum number of Tx antennas
M = 2; % maximum number of Rx antennas
% Create comm.BPSKModulator and comm.BPSKDemodulator System objects(TM)
P = 2; % modulation order
bpskMod = comm.BPSKModulator;
bpskDemod = comm.BPSKDemodulator('OutputDataType','double');

% Create comm.OSTBCEncoder and comm.OSTBCCombiner System objects


ostbcEnc = comm.OSTBCEncoder;
ostbcComb = comm.OSTBCCombiner;

% Create two comm.AWGNChannel System objects for one and two receive
% antennas respectively. Set the NoiseMethod property of the channel to
% 'Signal to noise ratio (Eb/No)' to specify the noise level using the
% energy per bit to noise power spectral density ratio (Eb/No). The output
% of the BPSK modulator generates unit power signals; set the SignalPower
% property to 1 Watt.
awgn1Rx = comm.AWGNChannel('NoiseMethod', 'Signal to noise ratio
(Eb/No)','SignalPower', 1);
awgn2Rx = clone(awgn1Rx);

% Create comm.ErrorRate calculator System objects to evaluate BER.


errorCalc2 = comm.ErrorRate;

% Since the comm.AWGNChannel System objects as well as the RANDI function


% use the default random stream, the following commands are executed so
% that the results will be repeatable, i.e., same results will be obtained
% for every run of the example. The default stream will be restored at the
% end of the example.
s = rng(55408);

% Pre-allocate variables for speed


H = zeros(frmLen, N, M);
ber_Alamouti = zeros(3,length(EbNo));
% Set up a figure for visualizing BER results
fig = figure;
grid on;
ax = fig.CurrentAxes;
hold(ax,'on');

ax.YScale = 'log';
xlim(ax,[EbNo(1), EbNo(end)]);
ylim(ax,[1e-4 1]);
xlabel(ax,'Eb/No (dB)');
ylabel(ax,'BER');
fig.NumberTitle = 'off';
fig.Renderer = 'zbuffer';
fig.Name = 'Transmit vs. Receive Diversity';
title(ax,'Transmit vs. Receive Diversity');
fig.Position = figposition([15 50 25 30]);

% Loop over several EbNo points


for idx = 1:length(EbNo)
reset(errorCalc2);
% Set the EbNo property of the AWGNChannel System objects
awgn1Rx.EbNo = EbNo(idx);
awgn2Rx.EbNo = EbNo(idx);
% Loop over the number of packets
for packetIdx = 1:numPackets
% Generate data vector per frame
data = randi([0 P-1], frmLen, 1);

% Modulate data
modData = step(bpskMod,data);

% Alamouti Space-Time Block Encoder


encData = step(ostbcEnc,modData);

% Create the Rayleigh distributed channel response matrix


% for two transmit and two receive antennas
H(1:N:end, :, :) = (randn(frmLen/2, N, M) + 1i*randn(frmLen/2, N,
M))/sqrt(2);
% assume held constant for 2 symbol periods
H(2:N:end, :, :) = H(1:N:end, :, :);

% Extract part of H to represent the 1x1, 2x1 and 1x2 channels


H21 = H(:,:,1)/sqrt(2);

% Pass through the channels


chanOut21 = sum(H21.* encData, 2);

% Add AWGN
rxSig21 = step(awgn1Rx,chanOut21);

% Alamouti Space-Time Block Combiner


decData = step(ostbcComb,rxSig21, H21);

% ML Detector (minimum Euclidean distance)


demod21 = step(bpskDemod,decData);

% Calculate and update BER for current EbNo value


% for uncoded 1x1 system
% for Alamouti coded 2x1 system
ber_Alamouti(:,idx) = step(errorCalc2,data, demod21);
% for Maximal-ratio combined 1x2 system

end % end of FOR loop for numPackets

% Plot results
semilogy(ax,EbNo(1:idx), ber_Alamouti(1,1:idx), 'go');
legend(ax,'Alamouti (2Tx, 1Rx)');

drawnow;
end % end of for loop for EbNo

% Perform curve fitting and replot the results


fitBER21 = berfit(EbNo, ber_Alamouti(1,:));
plot(ax,EbNo, fitBER21);
hold(ax,'off');

% Restore default stream


rng(s);

You might also like