plot_week1.asv

ＭＩＭＯ－ＯＦＤＭ系统仿真

clear all;
close all;
clc;

% plot LS channel estimation
figure;

BER_AWGN0 = [0.183034179687500,0.0907099609375000,0.0264218750000000,0.00273144531250000,4.39453125000000e-05];
semilogy(0:3:12, BER_AWGN0, 'b-','linewidth',2);
hold on;
grid on;

BER_AWGN1 = [0.338705078125000,0.241589843750000,0.135983398437500,0.0520644531250000,0.00982128906250000,0.000444726562500000,1.75781250000000e-06];
semilogy(0:3:18,BER_AWGN1,'r--','linewidth',2);

SNR = 0:3:21;
BER_Flat0 = [0.264300781250000,0.200087890625000,0.136792968750000,0.0879238281250000,0.0482089843750000,0.0268691406250000,0.0129140625000000,0.00705664062500000;];
semilogy(SNR,BER_Flat0,'^b-','linewidth',2);
BER_Flat1 = [0.371609375000000,0.312480468750000,0.242427734375000,0.169503906250000,0.109542968750000,0.0649707031250000,0.0379824218750000,0.0218398437500000];
semilogy(SNR,BER_Flat1,'^r-.','linewidth',2);

BER_Fre0 = [0.266335937500000,0.198181640625000,0.136345703125000,0.0851347656250000,0.0497343750000000,0.0266542968750000,0.0134453125000000,0.00689648437500000;];
semilogy(SNR,BER_Fre0,'*g-','linewidth',2);
BER_Fre1 = [0.369308593750000,0.308828125000000,0.239726562500000,0.171539062500000,0.110027343750000,0.0668281250000000,0.0370292968750000,0.0189335937500000;];
semilogy(SNR,BER_Fre1,'*g-.','linewidth',2);

xlabel('\fontsize{12}SNR(dB)');
ylabel('\fontsize{12}BER');
title('\fontsize{12} 2\times2MIMO i.i.d channel');
legend('\fontsize{12}AWGN Ideal Channel','\fontsize{12}AWGN Channel Estimation','\fontsize{12}AWGN Channel Estimation','\fontsize{12}AWGN Channel Estimation');

% plot fine channel estimation

SNR = 0:3:21;
BER_AWGN1 = [0.187267578125000,0.0940488281250000,0.0280292968750000,0.00310937500000000,4.88281250000000e-05];
MSE_AWGN1 = [0.0314375728454815,0.0158490101821601,0.00811370586115332,0.00398979606619284,0.00199653372933621,0.000994418860939845,0.000501146493436035,0.000249815222134386;]
BER_Flat1 = [0.269353515625000,0.206185546875000,0.142220703125000,0.0839433593750000,0.0484726562500000,0.0255859375000000,0.0148125000000000,0.00885351562500000];
MSE_Flat1 = [0.0301793305937498,0.0156625026479672,0.00801910045891135,0.00395120743050976,0.00202926227680630,0.00100306158520585,0.000497570890548635,0.000255658315766399;];
BER_Fre1 = [0.282744140625000,0.214691406250000,0.150863281250000,0.0959570312500000,0.0547871093750000,0.0311953125000000,0.0160390625000000,0.00792578125000000;];
MSE_Fre1 = [0.187962517905302,0.0932115053015102,0.0470585630893282,0.0234113773809637,0.0118417130094452,0.00593816784103308,0.00299766930468812,0.00148704747787788;];

figure;
semilogy(0:3:12, BER_AWGN0, 'b-','linewidth',2);
hold on;
grid on;
semilogy(0:3:12,BER_AWGN1,'*r-','linewidth',2);
semilogy(SNR,BER_Flat0,'^b-','linewidth',2);
semilogy(SNR,BER_Flat1,'^r-.','linewidth',2);
semilogy(SNR,BER_Fre0,'*g-','linewidth',2);
semilogy(SNR,BER_Fre1,'*g-.','linewidth',2);

xlabel('\fontsize{12}SNR(dB)');
ylabel('\fontsize{12}BER');
title('\fontsize{12} 2\times2MIMO i.i.d channel');
legend('\fontsize{12}AWGN Ideal Channel','\fontsize{12}AWGN Channel Estimation','\fontsize{12}AWGN Channel Estimation','\fontsize{12}AWGN Channel Estimation');


% plot correlation channel,用精估计，发射rou=0.5
BER_Flat1 = [0.312591796875000,0.247984375000000,0.180445312500000,0.113525390625000,0.0678105468750000,0.0368417968750000,0.0212500000000000,0.0129101562500000;];
MSE_Flat1 = [0.0301793305937498,0.0156625026479672,0.00801910045891135,0.00395120743050976,0.00202926227680630,0.00100306158520585,0.000497570890548635,0.000255658315766399;]