mimo.m

此程序完成MIMO各个自程序的功能

%MIMO Run MIMO system simulation.
%   MIMO is the main simulation script where all the blocks are glued
%   together. Refer to the PDF documentation. This is the body of the
%   program which allows you to set various simulation conditions and
%   change wide range of parameters. Modular architecture let's you
%   easy focus your effort on the particular block development while
%   the rest of the system is not affected by this module. Other modules
%   can be integrated into the body very easily. You just modify the
%   simulation parameters in the script header. 
%
%   MIMO script is also suited for a code performance measurement. In 
%   such a case you should disable function ECHOes and enable PROG_INFO
%   instead. This keeps you informed about the simulation run, which
%   usually takes a long time. You may also do some post-processing by
%   calling a POST function or archive the simulation report into a
%   log-file.
%
%   See also POST, LOGREPORT, MESG and the script body to learn more.

%   Copyright 2001-2003 Kamil Anis, [e-mail]
%   Dept. of Radioelectronics, 
%   Faculty of Electrical Engineering
%   Czech Technical University in Prague
%   $Revision: 2.2 $  $Date: 2003/5/23 22:01:58 $
%   --
%   <additional stuff goes here>

clear all; close all;
tic;

% system-wide variables
global ECHO SYMB_ENERGY SYMB_SAMPLES SYMB_TIME CH_CONF SMPLPER GAIN...
	SIG_PWR	PULSE INDENT TIMING VERSION

VERSION = '2.2.0';                          % toolkit version
ECHO = 'on';															% enable function echoes
% ECHO = 'off';                               % disable functions echoes
% prog_info = 1;                              % enable progress info
prog_info = 0;                            % disable progress info
% INDENT = 12;                              % set function indenting
fashion;																		% graphics settings

md = 4;                                  		% # of constellation signals
s = 8;                                   		% # of code trellis states

fr_length = 10;                           	% frame length [symbols]
frames = 20;                             	  % # of frames
zf = 3;                                 		% length of zero appendix
src_mode = 'rand';                     		  % source operation mode
% src_mode = 'ramp';
% src_mode = 'const';
src_const = 1;

enc_scheme = 'att';                         % encoding scheme

pulse_shape = 'rrc';                        % modulation impulse shape
% pulse_shape = 'rect';
pulse_cutoff = 8;                       		% mod. impulse shortening
pulse_rolloff = 0.4;                     		% roll-off factor
SYMB_TIME = 0.0001;                      		% symbol duration
SYMB_SAMPLES = 4;                        		% samples per symbol
SYMB_ENERGY = 1;                            % symbol average energy

CH_CONF = [2 2];                         		% MIMO configuration
ch_snr_pts = 10;                               % # of points to be measured
ch_snr_lo = 1;
ch_snr_hi = 10;
ch_snr = 4;                              	% channel SNR [dB]
% ch_snr = linspace(ch_snr_lo,...
%   ch_snr_hi,ch_snr_pts);
% ch_fading = 'none';                      	% MIMO channel fading
% ch_fading = 'awgn';
ch_fading = 'rayleigh';

est_csi = 'perfect';                     	  % CSI mode

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% INIT %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% pre-requisities
SMPLPER = SYMB_TIME / SYMB_SAMPLES; 		% symbol sample period
SIG_PWR = SYMB_ENERGY / SYMB_TIME;      % modulation (constell) power
GAIN = rescale(md,SYMB_ENERGY);         % scale factor
TIMING = 0; % disable DETECT timing when running performance measurement
timing = 0;

name = 'MIMO';
[idt,tag] = iecho(name);

[dlt,slt] = ltable(md,s,'Method',enc_scheme);

PULSE = makepulse('Shape',pulse_shape,'Shortening',pulse_cutoff,...
	'RollOff',pulse_rolloff);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% RUN %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

if prog_info
  clc;
  disp(' ');
  disp([tag,'Running MIMO system performance measurement. Please wait...']);
  disp(' ');
end

for i = 1:length(ch_snr)
	% running modules
	data_src = source(fr_length,frames,md,zf,'Pattern',src_mode);
	
	data_enc = space(data_src,dlt,slt);
	
	sig_mod = modul(data_enc,md);
	
	[sig_corr,ch_coefs] = channel(sig_mod,'Fading',ch_fading,'SNR',ch_snr(i));
	
	% [est_coefs,mse]=est(sig_corr,est_csi,ch_coefs);
	% this block has been temporarily removed. 
	
	sig_mf = mfilter(sig_corr);
	
	[data_est,state_est] = detect(sig_mf,dlt,slt,ch_coefs);
	
	% counting the errors
	[err_i,ser_i] = count(data_src,data_est);
	
	ser(i) = ser_i; err(i) = err_i;
  
  if prog_info
    mesg;
  end
end

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%% VISUAL LAYER & POST-PROCESSING %%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% check pulse shape
% stem(PULSE);

% display design
figure(1); dispdes(dlt,slt);
 
% display decoding process in a graphical form
figure (2); disptrell(dlt,slt,data_est,state_est);

% display code performance
if length(ch_snr) > 1 figure(3); post(ch_snr,ser); end

% save simulation report into a file
logreport;