瑞利信道下，2发2收，最大似然译码，空时编码，MATLAB仿真实现？？

function [FER,FER_uncoded,SER,SER_uncoded, BER, BER_uncoded]=stbc22(channel_model,K,Num,no_tx_antennas,no_rx_antennas,modulation)

%SNR upto 20 dBs
EbNo=[0:2:20];

%N, M: number of transmit and receive antennas

N=no_tx_antennas;
M=no_rx_antennas;

%initialize count
idx = 1;

h=waitbar(0,'Percentage Completed');
set(h,'Position',[230 60 275.25 56.25]);
set(h,'name','Please wait...');
wb=9.09;

for SNR=EbNo

sigma=0.5/(10^(SNR/10));

% Num -> number of packets
for packet_count=1:Num

% we are interested in transmitting 'K' SYMBOLS not bits. Hence, K*2 for QPSK
% etc.
switch (modulation)
case 'BPSK '
data=randint(K,N);
BIT=1;
case 'QPSK '
data=randint(K*2,N);
BIT=2;
case '8PSK '
data=randint(K*3,N);
BIT=3;
case '16QAM'
data=randint(K*4,N);
BIT=4;
otherwise
disp('No Modulation')
end
tx_bits=data.';
temp=[];
temp1=[];
for i=1:N
[temp1 s P]=tx_modulate(tx_bits(i,:),modulation);
temp=[temp; temp1];
temp1=0;
end

%ready to transmit symbols of length 'K'
X=temp.';

fr_length=length(X);

% block coding-Alamouti
x0=X(:,1);% required to verify a 1x1 system
x1=X;
x2(:,1)=-conj(X(:,2));
x2(:,2)=conj(X(:,1));

% form the channel matrix
for n=1:N
if channel_model=='AWGN '
Hr(n,:,:)=ones(fr_length,N);
else
Hr(n,:,:)=(randn(fr_length,N)+j*randn(fr_length,N))/sqrt(2);
end
end

for n=1:M
%transmission matrix
H=reshape(Hr(n,:,:),fr_length,N);
Habs(:,n)=sum(abs(H).^2,2);

%received signal per receiver antenna
r1(:,n)=sum(H.*x1,2)/sqrt(N)+sqrt(sigma)*(randn(fr_length,1)+j*randn(fr_length,1));
r2(:,n)=sum(H.*x2,2)/sqrt(N)+sqrt(sigma)*(randn(fr_length,1)+j*randn(fr_length,1));

% demodulate the received signals
z1(:,n)=r1(:,n).*conj(H(:,1))+conj(r2(:,n)).*H(:,2);
z2(:,n)=r1(:,n).*conj(H(:,2))-conj(r2(:,n)).*H(:,1);
end

%uncoded(1,1)
r01=H(:,1).*x0+sqrt(sigma)*(randn(fr_length,1)+j*randn(fr_length,1));

%form estimates
for m=1:P
d01(:,m)=abs(r01-H(:,1)*s(m)).^2; % uncoded signal
%coded signals
d1(:,m)=abs(sum(z1,2)-s(m)).^2+(-1+sum(Habs,2))*abs(s(m))^2;
d2(:,m)=abs(sum(z2,2)-s(m)).^2+(-1+sum(Habs,2))*abs(s(m))^2;
end

% determine the minimum of estimates

%decision for detecting uncoded
[y0,i0]=min((d01),[],2);
s0d=s(i0).';
clear d01

%decision for detecting s1
[y1,i1]=min(d1,[],2);
s1d=s(i1).';
clear d1
%decision for detecting s2
[y2,i2]=min(d2,[],2);
s2d=s(i2).';
clear d2
% form received symbols
Xd=[s1d s2d];

%determine symbol errors
error_un(packet_count)=sum(X(:,1)~=s0d);% for uncoded
temp1=X>0;
temp2=Xd>0;
error(packet_count)=sum(sum(temp1~=temp2));% for coded
end % end of FOR loop for "packet_count"
%calculate FER, SER and BER for current idx

%for uncoded signal
SER_uncoded(idx)=sum(error_un)/(Num*K);
BER_uncoded(idx)=SER_uncoded(idx)/BIT;
FER_uncoded(idx)=SER_uncoded(idx)*K;

%for coded signal
SER(idx)=sum(error)/(Num*K);
BER(idx)=SER(idx)/BIT;
FER(idx)=SER(idx)*K;

%increment idx
idx=idx + 1;

str_bar=[num2str(wb) '% Completed'];
waitbar(wb/100,h,str_bar);
wb=wb+9.09;

end % end of FOR loop for SNR

close(h);

%%%%%%%以下是tx_modulate的function.m

function [mod_symbols,table, P] = tx_modulate(bits_in, modulation)

%**************************************************************************
%This program modulates the input binary data.The inputs are:
% bits_in -> the binary input bits
%modulation -> we choose one of 'BPSK','QPSK',8PSK',16QAM'
%The outputs are:
%mod_symbols -> modulated output
%table -> the complete set of symbols in a chosen constellation. This is
%required by the demodulator.
%P -> the number of points in a chosen constellation. This is required by
%the demodulator
%**************************************************************************

full_len = length(bits_in);

% BPSK modulation
if ~isempty(findstr(modulation, 'BPSK'))
% Angle [pi/4 -3*pi/4] corresponds to
% Gray code vector [0 1], respectively.
table=exp(j*[0 -pi]); % generates BPSK symbols
table=table([1 0]+1); % Gray code mapping pattern for BPSK symbols
inp=bits_in;
mod_symbols=table(inp+1); % maps transmitted bits into BPSK symbols
P=2;% 2 constellation points
% QPSK modulation
elseif ~isempty(findstr(modulation, 'QPSK'))
% Angle [pi/4 3*pi/4 -3*pi/4 -pi/4] corresponds to
% Gray code vector [00 10 11 01], respectively.
table=exp(j*[-3/4*pi 3/4*pi 1/4*pi -1/4*pi]); % generates QPSK symbols
table=table([0 1 3 2]+1); % Gray code mapping pattern for QPSK symbols
inp=reshape(bits_in,2,full_len/2);
mod_symbols=table([2 1]*inp+1); % maps transmitted bits into QPSK symbols
P=4;% 4 constellation points
elseif ~isempty(findstr(modulation, '8PSK'))
% generate constellations
% Angle [0 pi/4 pi/2 3*pi/4 pi 5*pi/4 3*pi/2 7*pi/4] corresponds to
% Gray code vector [000 001 011 010 110 111 101 100], respectively.
table=exp(j*[0 pi/4 pi/2 3*pi/4 pi 5*pi/4 3*pi/2 7*pi/4]); % generates 8PSK symbols
table=table([0 1 3 2 6 7 5 4]+1); % Gray code mapping pattern for 8PSK symbols
inp=reshape(bits_in,3,full_len/3);
mod_symbols=table([4 2 1]*inp+1); % maps transmitted bits into 8PSK symbols
P=8;%8 constellation points
% 16-QAM modulation
elseif ~isempty(findstr(modulation, '16QAM'))
% generates 16QAM symbols
m=1;
for k=-3:2:3
for l=-3:2:3
table(m) = (k+j*l)/sqrt(10); % power normalization
m=m+1;
end;
end;
table=table([0 1 3 2 4 5 7 6 12 13 15 14 8 9 11 10]+1); % Gray code mapping pattern for 8-PSK symbols
inp=reshape(bits_in,4,full_len/4);
mod_symbols=table([8 4 2 1]*inp+1); % maps transmitted bits into 16QAM symbols
P=16; %16 constellation points

else
error('Unimplemented modulation');
end
希望能有帮助！

T2航站楼办理