%M=2,h=1/2,GMSK，N=512,g=[17,15],Ns=16,2倍采樣,L=2,H with 4 taps
clc
clear
EbN0=[0:2:10];
ber(1,:)=[1.0499e-001    7.3060e-002    6.2500e-002    6.8639e-002    6.1395e-002];
ber(2,:)=[5.0913e-002    2.5248e-002    2.3241e-002    2.3785e-002    2.3032e-002];
ber(3,:)=[2.2522e-002    7.6403e-003    6.5735e-003    6.5020e-003    6.5224e-003];
ber(4,:)=[8.4236e-003    2.6688e-003    2.2298e-003    2.0907e-003    2.2080e-003];
ber(5,:)=[3.4642e-003    7.1992e-004    4.2287e-004    4.4517e-004    4.0615e-004];
ber(6,:)=[1.0876e-003    1.5423e-004    6.3998e-005    4.9645e-005    4.8469e-005];
ber2 = [5.4295e-002  1.6745e-002  4.5721e-003  1.6621e-003  2.7238e-004  1.8468e-005]; % after decoding stbc MMSE Eq.
semilogy(EbN0,ber(:,1),'k-o',EbN0,ber(:,2),'k-s',EbN0,ber(:,3),'k->',EbN0,ber(:,4),'k-p',EbN0,ber(:,5),'k-<',EbN0,ber2,'k-v','LineWidth',2);
grid
xlabel('Eb/N0(dB)');
ylabel('BER');
legend('1st Iter','2nd Iter','3rd Iter','4th Iter','5th Iter','ML');