outagescenarios.m

空间无线信道建模程序matlab

% File OutageScenarios
%
% This is a set of scenarios for the Coded Modulation Library.
%
% The simulations specified in this file are for the information outage
% probability - independent Rayleigh-faded blocks
%
% NOTE: - In order to run OutageScenarios under modulation contraints, 
%       the desired modulation should have been simulated in an AWGN 
%       channel using the CML Library, (i.e. run
%       [sim_param, sim_state] = CmlSimulate( 'CapacityScenarios', record_no)
%       where the record_no is for the desired modulation )
%       - Use the filename that was entered in the previous record as the 
%         input_filename for the 'OutageScenarios' simulation record.
%
% Last updated on July 4, 2006

% determine where your root directory is
load( 'CmlHome.mat' );

% determine where to store your files
base_name = 'outage';
if ispc
    data_directory = strcat( cml_home, '\output\', base_name, '\' );
else
    data_directory = strcat( cml_home, '/output/', base_name, '/' );
end
if ~exist( data_directory, 'dir' )
    mkdir(data_directory);
end

% where are the AWGN capacity curves?
data_directory_input = strcat( cml_home, '\output\capacity\');

% 2-D UGI, R=2, L=1
record = 1;
sim_param(record).comment = '2-D UGI, CC, r = 2, L = 1';
sim_param(record).sim_type = 'outage';
sim_param(record).blocks_per_frame = 1;
sim_param(record).rate = 2;
sim_param(record).filename = strcat( data_directory, 'UGI2DCCL1r2.mat');
sim_param(record).linetype = 'r:';
sim_param(record).combining_type = 'code';
sim_param(record).SNR = [0:1:50];
sim_param(record).SNR_type = 'Es/No in dB';
sim_param(record).mod_order = 0;
sim_param(record).legend = sim_param(record).comment;
sim_param(record).reset = 0;
sim_param(record).max_trials = 1e9*ones( size(sim_param(record).SNR) );
sim_param(record).save_rate = 1e6;
sim_param(record).max_frame_errors = 3e3*ones( 1, length(sim_param(record).SNR) );
sim_param(record).minFER = 2e-5;
sim_param(record).trial_size = 1e4;

% 2-D UGI, R=2, L=2
record = 2;
sim_param(record).comment = '2-D UGI, CC, r = 2, L = 2';
sim_param(record).sim_type = 'outage';
sim_param(record).blocks_per_frame = 2;
sim_param(record).rate = 2;
sim_param(record).filename = strcat( data_directory, 'UGI2DCCL2r2.mat');
sim_param(record).linetype = 'b:';
sim_param(record).combining_type = 'code';
sim_param(record).SNR = [0:1:50];
sim_param(record).SNR_type = 'Es/No in dB';
sim_param(record).mod_order = 0;
sim_param(record).legend = sim_param(record).comment;
sim_param(record).reset = 0;
sim_param(record).max_trials = 1e9*ones( size(sim_param(record).SNR) );
sim_param(record).save_rate = 1e6;
sim_param(record).max_frame_errors = 3e3*ones( 1, length(sim_param(record).SNR) );
sim_param(record).minFER = 2e-5;
sim_param(record).trial_size = 1e4;

% 2-D UGI, R=2, L=3
record = 3;
sim_param(record).comment = '2-D UGI, CC, r = 2, L = 3';
sim_param(record).sim_type = 'outage';
sim_param(record).blocks_per_frame = 3;
sim_param(record).rate = 2;
sim_param(record).filename = strcat( data_directory, 'UGI2DCCL3r2.mat');
sim_param(record).linetype = 'k:';
sim_param(record).combining_type = 'code';
sim_param(record).SNR = [0:1:50];
sim_param(record).SNR_type = 'Es/No in dB';
sim_param(record).mod_order = 0;
sim_param(record).legend = sim_param(record).comment;
sim_param(record).reset = 0;
sim_param(record).max_trials = 1e9*ones( size(sim_param(record).SNR) );
sim_param(record).save_rate = 1e6;
sim_param(record).max_frame_errors = 3e3*ones( 1, length(sim_param(record).SNR) );
sim_param(record).minFER = 2e-5;
sim_param(record).trial_size = 1e4;

% 2-D UGI, R=2, L=4
record = 4;
sim_param(record).comment = '2-D UGI, CC, r = 2, L = 4';
sim_param(record).sim_type = 'outage';
sim_param(record).blocks_per_frame = 4;
sim_param(record).rate = 2;
sim_param(record).filename = strcat( data_directory, 'UGI2DCCL4r2.mat');
sim_param(record).linetype = 'c:';
sim_param(record).combining_type = 'code';
sim_param(record).SNR = [0:1:50];
sim_param(record).SNR_type = 'Es/No in dB';
sim_param(record).mod_order = 0;
sim_param(record).legend = sim_param(record).comment;
sim_param(record).reset = 0;
sim_param(record).max_trials = 1e9*ones( size(sim_param(record).SNR) );
sim_param(record).save_rate = 1e6;
sim_param(record).max_frame_errors = 3e3*ones( 1, length(sim_param(record).SNR) );
sim_param(record).minFER = 2e-5;
sim_param(record).trial_size = 1e4;

% 2-D UGI, R=2, L=10
record = 5;
sim_param(record).comment = '2-D UGI, CC, r = 2, L = 10';
sim_param(record).sim_type = 'outage';
sim_param(record).blocks_per_frame = 10;
sim_param(record).rate = 2;
sim_param(record).filename = strcat( data_directory, 'UGI2DCCL10r2.mat');
sim_param(record).linetype = 'm:';
sim_param(record).combining_type = 'code';
sim_param(record).SNR = [0:1:50];
sim_param(record).SNR_type = 'Es/No in dB';
sim_param(record).mod_order = 0;
sim_param(record).legend = sim_param(record).comment;
sim_param(record).reset = 0;
sim_param(record).max_trials = 1e9*ones( size(sim_param(record).SNR) );
sim_param(record).save_rate = 1e6;
sim_param(record).max_frame_errors = 3e3*ones( 1, length(sim_param(record).SNR) );
sim_param(record).minFER = 2e-5;
sim_param(record).trial_size = 1e4;

% 2-D gray-coded 16QAM, R=2, L=1
record = 6;
sim_param(record).comment = 'Gray coded 16-QAM, CC, r = 2, L = 1';
sim_param(record).sim_type = 'outage';
sim_param(record).blocks_per_frame = 1;
sim_param(record).rate = 2;
sim_param(record).filename = strcat( data_directory, 'QAM2DCCL1r2.mat');
sim_param(record).linetype = 'r-.+';
sim_param(record).combining_type = 'code';
sim_param(record).SNR = [0:1:50];
sim_param(record).SNR_type = 'Es/No in dB';
sim_param(record).mod_order = 16;
sim_param(record).input_filename = strcat( data_directory_input, 'QAM16BICMAWGN' );
sim_param(record).legend = sim_param(record).comment;
sim_param(record).reset = 0;
sim_param(record).max_trials = 1e9*ones( size(sim_param(record).SNR) );
sim_param(record).save_rate = 1e6;
sim_param(record).max_frame_errors = 3e3*ones( 1, length(sim_param(record).SNR) );
sim_param(record).minFER = 2e-5;
sim_param(record).trial_size = 1e4;

% 2-D gray-coded 16QAM, R=2, L=2
record = 7;
sim_param(record).comment = 'Gray coded 16-QAM, CC, r = 2, L = 2';
sim_param(record).sim_type = 'outage';
sim_param(record).blocks_per_frame = 2;
sim_param(record).rate = 2;
sim_param(record).filename = strcat( data_directory, 'QAM2DCCL2r2.mat');
sim_param(record).linetype = 'b-.+';
sim_param(record).combining_type = 'code';
sim_param(record).SNR = [0:1:50];
sim_param(record).SNR_type = 'Es/No in dB';
sim_param(record).mod_order = 16;
sim_param(record).input_filename = strcat( data_directory_input, 'QAM16BICMAWGN' );
sim_param(record).legend = sim_param(record).comment;
sim_param(record).reset = 0;
sim_param(record).max_trials = 1e9*ones( size(sim_param(record).SNR) );
sim_param(record).save_rate = 1e6;
sim_param(record).max_frame_errors = 3e3*ones( 1, length(sim_param(record).SNR) );
sim_param(record).minFER = 2e-5;
sim_param(record).trial_size = 1e4;

% 2-D gray-coded 16QAM, R=2, L=3
record = 8;
sim_param(record).comment = 'Gray coded 16-QAM, CC, r = 2, L = 3';
sim_param(record).sim_type = 'outage';
sim_param(record).blocks_per_frame = 3;
sim_param(record).rate = 2;
sim_param(record).filename = strcat( data_directory, 'QAM2DCCL3r2.mat');
sim_param(record).linetype = 'k-.+';
sim_param(record).combining_type = 'code';
sim_param(record).SNR = [0:1:50];
sim_param(record).SNR_type = 'Es/No in dB';
sim_param(record).mod_order = 16;
sim_param(record).input_filename = strcat( data_directory_input, 'QAM16BICMAWGN' );
sim_param(record).legend = sim_param(record).comment;
sim_param(record).reset = 0;
sim_param(record).max_trials = 1e9*ones( size(sim_param(record).SNR) );
sim_param(record).save_rate = 1e6;
sim_param(record).max_frame_errors = 3e3*ones( 1, length(sim_param(record).SNR) );
sim_param(record).minFER = 2e-5;
sim_param(record).trial_size = 1e4;

% 2-D gray-coded 16QAM, R=2, L=4
record = 9;
sim_param(record).comment = 'Gray coded 16-QAM, CC, r = 2, L = 4';
sim_param(record).sim_type = 'outage';
sim_param(record).blocks_per_frame = 4;
sim_param(record).rate = 2;
sim_param(record).filename = strcat( data_directory, 'QAM2DCCL4r2.mat');
sim_param(record).linetype = 'c-.+';
sim_param(record).combining_type = 'code';
sim_param(record).SNR = [0:1:50];
sim_param(record).SNR_type = 'Es/No in dB';
sim_param(record).mod_order = 16;
sim_param(record).input_filename = strcat( data_directory_input, 'QAM16BICMAWGN' );
sim_param(record).legend = sim_param(record).comment;
sim_param(record).reset = 0;
sim_param(record).max_trials = 1e9*ones( size(sim_param(record).SNR) );
sim_param(record).save_rate = 1e6;
sim_param(record).max_frame_errors = 3e3*ones( 1, length(sim_param(record).SNR) );
sim_param(record).minFER = 2e-5;
sim_param(record).trial_size = 1e4;

% 2-D gray-coded 16QAM, R=2, L=10
record = 10;
sim_param(record).comment = 'Gray coded 16-QAM, CC, r = 2, L = 10';
sim_param(record).sim_type = 'outage';
sim_param(record).blocks_per_frame = 10;
sim_param(record).rate = 2;
sim_param(record).filename = strcat( data_directory, 'QAM2DCCL10r2.mat');
sim_param(record).linetype = 'm-.+';
sim_param(record).combining_type = 'code';
sim_param(record).SNR = [0:1:50];
sim_param(record).SNR_type = 'Es/No in dB';
sim_param(record).mod_order = 16;
sim_param(record).input_filename = strcat( data_directory_input, 'QAM16BICMAWGN' );
sim_param(record).legend = sim_param(record).comment;
sim_param(record).reset = 0;
sim_param(record).max_trials = 1e9*ones( size(sim_param(record).SNR) );
sim_param(record).save_rate = 1e6;
sim_param(record).max_frame_errors = 3e3*ones( 1, length(sim_param(record).SNR) );
sim_param(record).minFER = 2e-5;
sim_param(record).trial_size = 1e4;

% 2-D UGI, R=3.5, L=1
record = 11;
sim_param(record).comment = '2-D UGI, CC, r = 3.5, L = 1';
sim_param(record).sim_type = 'outage';
sim_param(record).blocks_per_frame = 1;
sim_param(record).rate = 7/2;
sim_param(record).filename = strcat( data_directory, 'UGI2DCCL1r7by2.mat');