fig.asv

该程序源代码实现了OFDM通信系统性能的仿真

% ofdm.m ber performance
clear
clc
A=[0 1 2 3 4 5 6 7 8 9 10 ];
B=[1.022786e-001 7.727214e-002	5.581380e-002	3.655599e-002	2.252604e-002 1.180339e-002 5.625000e-003 2.180990e-003 7.421875e-004 2.148438e-004 3.255208e-005 ];
C=[1.000000e+000 1.000000e+000 1.000000e+000 1.000000e+000 1.000000e+000 1.000000e+000 1.000000e+000 9.500000e-001 6.600000e-001 2.900000e-001 5.000000e-002 ];
semilogy(A,B,'bd:','Markersize',7)
grid on
hold on;
% 11  6.510417e-006  1.000000e-002
% bpsk theory
% x=0:.1:10;
% y=erfc(sqrt(x));
% y=y/2;
% semilogy(x,y,'k');
%qpsk one-path rayleigh fading theory
% w=0:.5:40;
% z=1-1./sqrt(1+1./w);
% z=z/2;
% semilogy(w,z,'m');
% %ofdm_fading
%flat==1
% ber=[1.539362e-001 9.473698e-002 6.448307e-002 3.441667e-002 2.178125e-002 1.055599e-002 6.533854e-003...
%     1.739583e-003 5.598958e-004 1.263021e-004  5.468750e-005  2.343750e-005 9.114583e-006 5.208333e-006 3.906250e-006 0.000000e+000];
% per=[1.000000e+000 9.720000e-001 8.660000e-001 6.340000e-001 4.600000e-001 2.480000e-001 1.600000e-001  5.000000e-002 1.400000e-002 4.000000e-003 2.000000e-003 2.000000e-003  2.000000e-003 2.000000e-003 2.000000e-003 0.000000e+000];
% F=[0 3 5 8 10 13 15 20 25 30 32 33 34 35 36 37];
% semilogy(F,ber,'kp-','Markersize',5);
%flat==0
% ber1=[4.926185e-001 4.906810e-001 4.909518e-001 4.909427e-001 4.917318e-001 4.925247e-001 4.931445e-001 4.939753e-001 4.940352e-001 4.936979e-001 4.936758e-001 ];
% per1=[1.000000e+000 1.000000e+000 9.920000e-001 9.780000e-001 9.560000e-001 8.780000e-001 8.320000e-001 8.000000e-001 7.720000e-001 7.660000e-001 7.640000e-001];
% G=[0 3 5 8 10 15 20 25 30 32 33 ];
% semilogy(G,ber1,'bp:','Markersize',5);
%ofdmda.m
%The default mode of "ofdmda.m" is BER and PER under AWGN environment
t3=[0 3 5 6 7 8 10 11 12];
ber3=[1.029728e-001 3.710256e-002 1.237821e-002 5.756410e-003 2.299679e-003 7.003205e-004 3.525641e-005 1.602564e-006 0.000000e+000];
per3=[1.000000e+000 1.000000e+000 1.000000e+000 9.670000e-001 7.670000e-001 3.540000e-001 2.000000e-002 1.000000e-003 0.000000e+000];
semilogy(t3,ber3);
%BER and PER performances under fading environment. 
%no compensation
t4=[0 5 10 15 20 25 30 35 40 45 50];
ber4=[4.958317e-001 4.961763e-001 4.984439e-001 4.998189e-001 4.997500e-001 4.997083e-001 4.989696e-001 4.988862e-001 4.990224e-001 4.991683e-001 4.991619e-001];
per4=[1.000000e+000 9.920000e-001 9.330000e-001 8.700000e-001 8.310000e-001 7.990000e-001 7.740000e-001 7.600000e-001 7.550000e-001 7.540000e-001 7.520000e-001];
semilogy(t4,ber4)
%the effect of perfect compensation under fading environment. 
t5=[0 5 10 15 20 25 30 33 34 35];
ber5=[1.550192e-001 6.521795e-002 2.155449e-002 6.548077e-003 2.057692e-003 6.586538e-004 9.935897e-005 1.923077e-005 1.602564e-006 0.000000e+000];
per5=[9.980000e-001 8.430000e-001 4.180000e-001 1.310000e-001 4.400000e-002 1.200000e-002 4.000000e-003 4.000000e-003 1.000000e-003 0.000000e+000];
semilogy(t5,ber5);
%BER & PER performances under AWGN and fading environments
% The default mode of "ofdmda.m" is BER and PER under AWGN environment
t6=[0 3 5 8 9 10 12 13];
ber6=[1.205785e-001 4.843429e-002 1.852724e-002 1.623397e-003 5.256410e-004 9.455128e-005 1.602564e-006 0.000000e+000];
per6=[1.000000e+000 1.000000e+000 1.000000e+000 6.390000e-001 2.850000e-001 5.500000e-002 1.000000e-003 0.000000e+000];

% BER and PER performances under fading environment

t7=[0 5 10 15 20 25 30 35 40 45 50];
ber7=[4.962869e-001 4.959135e-001 4.981234e-001 4.997131e-001 4.996394e-001 4.993349e-001 4.993029e-001 4.988494e-001 4.990913e-001 4.993734e-001 4.993670e-001];
per7=[1.000000e+000 9.930000e-001 9.390000e-001 8.720000e-001 8.360000e-001 8.010000e-001 7.750000e-001 7.620000e-001 7.550000e-001 7.540000e-001 7.520000e-001];
%  perfect compensation under fading environment
t8=[0 10 20 30 33 35 36 37];
ber8=[1.701731e-001 2.480128e-002 2.325321e-003 1.554487e-004 2.564103e-005 6.410256e-006 3.205128e-006 0.000000e+000];
per8=[1.000000e+000 4.580000e-001 5.100000e-002 6.000000e-003 4.000000e-003 3.000000e-003 1.000000e-003 0.000000e+000];
%  the effect of CE-symbol under fading environment
t9=[0 10 20 30 40 45 46];
ber9=[2.657131e-001 5.004647e-002 4.653846e-003 4.839744e-004 1.282051e-005 4.807692e-006 0.000000e+000];
per9=[1.000000e+000 6.750000e-001 8.100000e-002 8.000000e-003 2.000000e-003 1.000000e-003 0.000000e+000];
% Simulation of "ofdmci.m"
%BER & PER when interference wave exists under fading environmentCE-based fading compensation
0.000000		1.000000e+000	1000	150
10.000000		8.610000e-001	1000	150
20.000000		4.890000e-001	1000	150
30.000000		3.700000e-001	1000	150
40.000000		3.310000e-001	1000	150
100.000000		3.190000e-001	1000	150
200.000000		3.210000e-001	1000	150
300.000000		3.210000e-001	1000	150
1000.000000	3.180000e-001	1000	150
2000.000000	6.493910e-002	3.220000e-001	1000	150
t10=[0 10 20 30 40 100 200 300 1000 2000];
ber10=[2.831474e-001 1.048221e-001 6.938942e-002 6.692308e-002 6.389583e-002
6.467468e-002 6.641827e-002 6.701442e-002 6.557532e-002	];
per10=[];
%perfect fading compensation
% 0.000000	1.822869e-001	1.000000e+000	1000	150
% 5.000000	1.003526e-001	9.450000e-001	1000	150
% 8.000000	6.956731e-002	7.910000e-001	1000	150
% 10.000000	5.686699e-002	6.630000e-001	1000	150
% 20.000000	3.532372e-002	2.720000e-001	1000	150
% 30.000000	3.018429e-002	1.480000e-001	1000	150
% 50.000000	2.916186e-002	1.040000e-001	1000	150
% 100.000000	2.970513e-002	1.000000e-001	1000	150
% 200.000000	3.291667e-002	1.130000e-001	1000	150
% 500.000000	2.916827e-002	1.010000e-001	1000	150
% 1000.000000	2.937981e-002	1.020000e-001	1000	150
% 2000.000000	2.951122e-002	1.010000e-001	1000	150
% 10000.000000	2.913942e-002	1.010000e-001	1000	150