aeronauticalchannel.m

国外发布的航空信道matlab模型

% Disclaimer: The channel simulation software "aeronautical channel model" has been
% developed by DLR under ESA-Contract.  The use of this software is free of charge. 
% Anyone using this software, will use it exclusively at his/her own risk and responsibility.
% Neither ESA nor DLR nor their partners, suppliers, affiliates nor any other party involved
% in creating, producing, or delivering the software is liable for any direct, indirect,
% incidental, special, consequential, punitive or other damages whatsoever including business
% interruption, loss of use, data, information, profits (regardless of the form of action,
% including but not limited to contract, negligence or other tortious act) arising out of or
% in connection with the use of the software. Without limiting the foregoing, everything in
% the software is provided "as it is" without representation, warranty or condition of any
% kind, either express, implied, or statutory including, but not limited to, implied
% representations, warranties or conditions of merchantabiltiy, fitness for a particular
% purpose, durability, title, non-infringement of intellectual property rights or
% inter-operability of products or services.

clear all
close all
clc

SampFreq = 25.4; %Hz
StartAltitude = 800; %m
DescendingRate = 3; %m/s
AircraftType = 'VFW614';
MarcovProbs = {[0.9866 0.0087 0.0047 0     ; ...
                0.6087 0.3043 0.0870 0     ; ...
                0.2143 0.3571 0.4286 0     ; ...
                0.3333 0.3333 0.3334 0     ] ...
                ,...
               [0.9842 0.0130 0.0028 0     ; ...
                0.6667 0.2222 0.0889 0.0222; ...
                0.0667 0.1167 0.5000 0.3166; ...
                0      0      0.3279 0.6721] ...
                ,...
               [0.9645 0.0310 0.0045 0     ; ...
                0.7308 0.1538 0.1154 0     ; ...
                0.6250 0.1250 0.2500 0     ; ...
                0.3333 0.3333 0.3334 0     ] ...
                ,...
               [1      0      0      0     ; ...
                1      0      0      0     ; ...
                1      0      0      0     ; ...
                1      0      0      0     ]};   
MarkovAtten = {[-50 -23 -19 -15],... %state 1
               [-50 -23 -19 -15],... %state 2
               [-50 -23 -19 -15],... %state 3
               [-50 -23 -19 -15]};   %state 4 in dB
Altitudevector = [1500 400 100 10];%m
Elevation   = 10;    %degrees
Azimuth     = 120;   %degrees

state=Init_AeronauticalChannel(SampFreq,StartAltitude,DescendingRate,AircraftType,MarcovProbs,MarkovAtten,Altitudevector,Elevation,Azimuth);

disp('generating samples')

[state,Output1]=GenerateAeronauticalChannel(state,10000);
[state,Output2]=GenerateAeronauticalChannel(state,100);