selfgeometry2tgncorrelation.m

这是一个关于ofdm在802.11协议下的源码

% SelfGeometry2TGnCorrelation
%
% Computation of the correlation matrices of the MIMO radio
% channel model described in IEEE 802 document 11-03/940r0
% for a limited set of hard-coded scenarios.
%
%
% STANDARD DISCLAIMER
%
% The Computer Science Institute of the University of Namur (hereafter
% "FUNDP-INFO") is furnishing this item "as is". FUNDP-INFO does not
% provide any warranty of the item whatsoever, whether express,
% implied, or statutory, including, but not limited to, any warranty
% of merchantability or fitness for a particular purpose or any
% warranty that the contents of the item will be error-free.
%
% In no respect shall FUNDP-INFO incur any liability for any damages,
% including, but not limited to, direct, indirect, special, or
% consequential damages arising out of, resulting from, or any way
% connected to the use of the item, whether or not based upon
% warranty, contract, tort, or otherwise; whether or not injury was
% sustained by persons or property or otherwise; and whether or not
% loss was sustained from, or arose out of, the results of, the
% item, or any services that may be provided by FUNDP-INFO.
%
% (c) Laurent Schumacher, FUNDP-INFO - November 2003

clear all;
close all;
delete TGnCorrelationMatrices.mat;

TGnScenario        = ['A','B','C','D','E','F'];
NumberOfTxAntennas = [1,2,4];
SpacingTx          = [.5,1,4];
NumberOfRxAntennas = [1,2,4];
SpacingRx          = [.5,1,4];

qq = 1;
for (kk=1:size(TGnScenario,2))
    for (ll=1:size(NumberOfTxAntennas,2))
        for (mm=1:size(NumberOfRxAntennas,2))
            for (nn=1:size(SpacingTx,2))
                for (oo=1:size(SpacingRx,2))
                    Scenario = [TGnScenario(kk),'_',num2str(NumberOfTxAntennas(ll)),...
                            'x',num2str(NumberOfRxAntennas(mm)),'_'];
                    if (SpacingTx(nn) < 1)
                        Scenario = [Scenario,'i',num2str(inv(SpacingTx(nn))),'_'];
                    else
                        Scenario=[Scenario,num2str(SpacingTx(nn)),'_'];
                    end;
                    if (SpacingRx(oo) < 1)
                        Scenario = [Scenario,'i',num2str(inv(SpacingRx(oo)))];
                    else
                        Scenario=[Scenario,num2str(SpacingRx(oo))];
                    end;
                    time = clock;
                    disp([num2str(time(4)),':',num2str(time(5)),'.',num2str(floor(time(6))),...
                            ' - Iteration ',num2str(qq),'/',num2str(size(TGnScenario,2)*...
                            size(NumberOfTxAntennas,2)*size(NumberOfRxAntennas,2)*...
                            size(SpacingTx,2)*size(SpacingRx,2)),' - ',Scenario]);

                    %
                    % Description of scenarios
                    %
                    
                    switch TGnScenario(kk)
                        
                    case 'A'
                        % Flat fading with 0 ns rms delay spread (one tap at 0 ns delay model)
                        PDP_dB = [0;   % Average power [dB]
                                  0];  % Relative delay (ns)
                        % Power roll-off coefficients
                        Power_per_angle_dB = 0;
                        % Tx
                        AoD_Tx_deg = 45;
                        AS_Tx_deg  = 40;
                        Type_Tx    = 3;
                        % Rx
                        AoA_Rx_deg = 45;
                        AS_Rx_deg  = 40;
                        Type_Rx    = 3;
                    case 'B'
                        % Typical residential environment, LOS conditions, 15 ns rms delay spread, and 10 dB Ricean K-factor at first delay
                        PDP_dB = [0 -5.4287 -2.5162 -5.8905 -9.1603 -12.5105 -15.6126 -18.7147 -21.8168;  % Average power [dB]
                                  0 10e-9   20e-9   30e-9   40e-9   50e-9    60e-9    70e-9    80e-9   ]; % Relative delay (ns)
                        % Power roll-off coefficients
                        Power_per_angle_dB = [   0   -5.4287  -10.8574  -16.2860  -21.7147      -Inf      -Inf      -Inf      -Inf;
                                              -Inf      -Inf   -3.2042   -6.3063   -9.4084  -12.5105  -15.6126  -18.7147  -21.8168];
                        % Tx
                        AoD_Tx_deg = [225.1084.*ones(1,5)  -Inf.*ones(1,4);
                                      -Inf.*ones(1,2)      106.5545.*ones(1,7)];
                        AS_Tx_deg  = [14.4490.*ones(1,5)  -Inf.*ones(1,4);
                                      -Inf.*ones(1,2)     25.4311.*ones(1,7)];
                        Type_Tx    = 3.*ones(1, size(AoD_Tx_deg, 2));
                        % Rx
                        AoA_Rx_deg = [4.3943.*ones(1,5)  -Inf.*ones(1,4);
                                      -Inf.*ones(1,2)    118.4327.*ones(1,7)];
                        AS_Rx_deg  = [14.4699.*ones(1,5)  -Inf.*ones(1,4);
                                      -Inf.*ones(1,2)     25.2566.*ones(1,7)];
                        Type_Rx    = 3.*ones(1, size(AoA_Rx_deg, 2));
                    case 'C'
                        % Typical residential or small office environment, LOS/NLOS conditions, 30 ns rms delay spread, and 3 dB Ricean K-factor at the first delay
                        PDP_dB = [0 -2.1715 -4.3429 -6.5144 -8.6859 -10.8574 -4.3899 -6.5614 -8.7329 -10.9043 -13.7147 -15.8862 -18.0577 -20.2291; % Average power [dB]
                                  0 10e-9   20e-9   30e-9   40e-9   50e-9    60e-9   70e-9   80e-9   90e-9    110e-9   140e-9   170e-9   200e-9];  % Relative delay (ns)
                        % Power roll-off coefficients
                        Power_per_angle_dB = [   0 -2.1715 -4.3429 -6.5144 -8.6859 -10.8574 -13.0288 -15.2003 -17.3718 -19.5433     -Inf     -Inf     -Inf     -Inf;
                                              -Inf    -Inf    -Inf    -Inf    -Inf     -Inf  -5.0288  -7.2003  -9.3718 -11.5433 -13.7147 -15.8862 -18.0577 -20.2291];
                        % Tx
                        AoD_Tx_deg = [13.5312.*ones(1,10)  -Inf.*ones(1,4);
                                      -Inf.*ones(1,6)      56.4329.*ones(1,8)];
                        AS_Tx_deg  = [24.7897.*ones(1,10)  -Inf.*ones(1,4);
                                      -Inf.*ones(1,6)      22.5729.*ones(1,8)];
                        Type_Tx    = 3.*ones(1, size(AoD_Tx_deg, 2));
                        % Rx
                        AoA_Rx_deg = [290.3715.*ones(1,10)  -Inf.*ones(1,4);
                                      -Inf.*ones(1,6)       332.3754.*ones(1,8)];
                        AS_Rx_deg  = [24.6949.*ones(1,10)  -Inf.*ones(1,4);
                                      -Inf.*ones(1,6)      22.4530.*ones(1,8)];
                        Type_Rx    = 3.*ones(1, size(AoA_Rx_deg, 2));
                    case 'D'
                        % Medbo model A - Typical office environment, NLOS conditions, and 50 ns rms delay spread
                        PDP_dB = [0 -0.9  -1.7  -2.6  -3.5  -4.3  -5.2  -6.1  -6.9  -7.8  -4.7   -7.3   -9.9   -12.5  -13.7  -18    -22.4  -26.7;   % Average power [dB]
                                  0 10e-9 20e-9 30e-9 40e-9 50e-9 60e-9 70e-9 80e-9 90e-9 110e-9 140e-9 170e-9 200e-9 240e-9 290e-9 340e-9 390e-9]; % Relative delay (ns)
                        % Power roll-off coefficients
                        Power_per_angle_dB = [0    -0.9 -1.7 -2.6 -3.5 -4.3 -5.2 -6.1 -6.9 -7.8 -9.0712046 -11.199064  -13.795428 -16.391791 -19.370991 -23.201722 -Inf       -Inf;
                                              -Inf -Inf -Inf -Inf -Inf -Inf -Inf -Inf -Inf -Inf -6.6756386  -9.5728825 -12.175385 -14.777891 -17.435786 -21.992788 -25.580689 -Inf;
                                              -Inf -Inf -Inf -Inf -Inf -Inf -Inf -Inf -Inf -Inf -Inf       -Inf        -Inf       -Inf       -18.843300 -23.238125 -25.246344 -26.7];
                        % Tx
                        AoD_Tx_deg = [332.1027.*ones(1,16)  -Inf.*ones(1,2);
                                      -Inf.*ones(1,10)        49.3840.*ones(1,7)  -Inf;
                                      -Inf.*ones(1,14)       275.9769.*ones(1,4)];
                        AS_Tx_deg  = [27.4412.*ones(1,16)  -Inf.*ones(1,2);
                                      -Inf.*ones(1,10)     32.1430.*ones(1,7)  -Inf;
                                      -Inf.*ones(1,14)     36.8825.*ones(1,4)];
                        Type_Tx    = 3.*ones(1, size(AoD_Tx_deg, 2));
                        % Rx