ITU-R p.1546 calculator: Test_step81.m - File Exchange

Code covered by the BSD License

Highlights from
ITU-R p.1546 calculator

ITUP1546(varargin) ITUP1546 M-file for ITUP1546.fig
D06(f,h1,h2) 17 An approximation to the 0.6 Fresnel clearance path length
Emax(d,t,path) 2 Maximum field-strength values
FindDNominals(d)
P1546FieldStr(d,f,t,heff,vara... P1546FieldStr ITU 1546-3 Field strength calculator
P1546FieldStrbeta(d,f,t,heff,... P1546FieldStr ITU 1546-3 Field strength calculator
Qi(x) 15 An approximation to the inverse complementary cumulative normal distribution
Step_12(f,tca) Step 12: If information on the terrain clearance angle at a
Step_13(d,f,t,eff,eff1) Step 13: Calculate the estimated field strength due to tropospheric
Step_14(h2,f,path,varargin) Step 14: Correct the field strength for receiving/mobile antenna height
Step_15(d,h1,ha,f,env) Step 15: If applicable, reduce the field strength by adding the
Step_16() Step 16: If the field strength at a receiving/mobile antenna adjacent to
Step_17(E,Efs,Ese,ds,dtotal) Step 17: If necessary, limit the resulting field strength to the maximum
Step_18(fs,f) Step 18: If required, convert field strength to equivalent basic transmission loss for the path using
V(Kv,h1)
d =FS_find_dnew(FS,f,heff,T,p... FS_FIND_D ITU 1546 Field strength calculator: Find d for known Fieldstrength
find_d(FS,f,heff,T,path,h2) Peverted code the code find the closest value d value that equates the
h1Calc(d,heff,ha,hb,path) 3 Determination of transmitting/base antenna height, h1
h2corr(h2,f,varargin)
searchclosest(x,v) The following code tidbit is by Dr. Murtaza Khan, modified to return
step6_10(figure,h1,dinf,dsup,... function E = step6_10(figure,h1,dinf,dsup,d,path,f,Emaxvalue,t)
step7_normal(figureStep6,h1,d...
step81(tabulatedValues,h1,din... Step 8: Obtain the field strength exceeded at 50% locations for a
step814_815(tabulatedValues,h... Step 8.1.4: Obtain the field strength exceeded at 50% locations
step82(tabulatedValues,h1,din... Step 8.2: For a transmitting/base antenna height h1 less than 10 m
P1546validation.m
Test_D06.m
Test_Emax.m
Test_FindDNominals.m
Test_P1546FieldStr.m
Test_Qi.m
Test_V.m
Test_h1Calc.m
Test_searchclosest.m
Test_step6_10.m
Test_step7_normal.m
Test_step81.m
Test_step814_815.m
Test_step82.m
View all files

from ITU-R p.1546 calculator by Jef Statham
This an implementation of the ITU-R p.1546 model. Mixed land a sea paths not implemented.

Test_step81.m

clc, close all,
%Test function E = step81(tabulatedValues,h1,dinf,dsup,d,Emax)
%test init
%     time,path,freq,figure 
figure = [50,1,100,1;
          10,1,100,2;
           1,1,100,3;
          50,2,100,4;
          10,3,100,5;
           1,3,100,6;
          10,4,100,7;
           1,4,100,8;
          50,1,600,9;
          10,1,600,10;
           1,1,600,11;
          50,2,600,12;
          10,3,600,13;
           1,3,600,14;
          10,4,600,15;
           1,4,600,16;
          50,1,2000,17;
          10,1,2000,18;
           1,1,2000,19;
          50,2,2000,20;
          10,3,2000,21;
           1,3,2000,22;
          10,4,2000,23;
           1,4,2000,24;];
        %p = 50%          %path = Sea    %frequency = 600 MHz
idx = (figure(:,1)==50) & (figure(:,2)==2) & (figure(:,3)==600);
figure = figure(idx,:);
tabulatedValues = xlsread('Rec_P_1546_2_Tab_values.xls',figure(1,4), 'B6:K84');

%Test 1 begin
h1 = 300;
d = 14.5;
dsup = 15;
dinf = 14;
Emax = 10;
Test = 'Test 1'
a = Emax;
isequal(a,step81(tabulatedValues,h1,dinf,dsup,d,Emax)) 
%Test 1 end
%Test 2 begin
Emax = 1000;
Test = 'Test 2'
a = 83.959977+(83.353634-83.959977)*log10(d/dinf)/log10(dsup/dinf);
isequal(a,step81(tabulatedValues,h1,dinf,dsup,d,Emax)) 
%Test 2 end
%Test 3 begin
h1 = 1400;
Test = 'Test 3'
a = 83.976265+(83.37648-83.976265)*log10(d/dinf)/log10(dsup/dinf);
% commented out lines were for tracing an error with interpolating
    %h1 = 600;
    %step814_815(tabulatedValues,h1,dinf,dsup,d)
    %h1=1200;
b = 83.977367+(83.378068-83.977367)*log10(d/dinf)/log10(dsup/dinf);
    %step814_815(tabulatedValues,h1,dinf,dsup,d)
    %h1 = 1400;
c = a+(b-a)*log10(h1/600)/log10(1200/600);
    %step81(tabulatedValues,h1,dinf,dsup,d,Emax)
isequal(c,step81(tabulatedValues,h1,dinf,dsup,d,Emax)) 
%Test 3 end
%Test 4 
tabulatedValues = xlsread('Rec_P_1546_2_Tab_values.xls',19, 'B6:K84');
h1 = 191;
dinf = 975;
dsup = 1000;
d = 984;
Line37 = step81(tabulatedValues,h1,dinf,dsup,d,Emax)
%Test 4 end
%Test 5
tabulatedValues = xlsread('Rec_P_1546_2_Tab_values.xls',22, 'B6:K84');
h1 = 448;
dinf = 975;
dsup = 1000;
d = 988;
Line47 = step81(tabulatedValues,h1,dinf,dsup,d,Emax)
%test 5 end
%test 6
tabulatedValues = xlsread('Rec_P_1546_2_Tab_values.xls',21, 'B6:K84');
h1 = 660;
dinf = 650;
dsup = 675;
d = 672;
Line45 = step81(tabulatedValues,h1,dinf,dsup,d,Emax)
%test 6 end
tabulatedValues = xlsread('Rec_P_1546_2_Tab_values.xls',10, 'B6:K84');
h1 = 341;
dinf = 40;
dsup = 45;
d = 44;
Line19 = step81(tabulatedValues,h1,dinf,dsup,d,Emax)

Highlights from ITU-R p.1546 calculator

Highlights from
ITU-R p.1546 calculator