max_distance22.m

this code is for coverage distance

%in the name of God
clc
clear all
close all

%********Concept**************************

beta       =60*pi/180;
Ht_min     =[];
Ht_max     =[];

%********calculation**********************

site_distance      = input('node_distance(Km)=');
node_distance      = node_distance*1000;

azimuth_angle      = input('Azimuth_angle(deg)=');
azimuth_angle      = azimuth_angle*pi/180;

a        = node_distance*cos(beta);
H_max1   =  3e5-a;

half_node_distance      = node_distance*sin(beta);

gama_max      =(atan(H_max1/half_node_distance));

Teta_max      = gama_max-(azimuth_angle/2);

for teta  = .01:.001:Teta_max,
    
    H_min     = half_node_distance*tan(teta-(azimuth_angle/2));
    H1_min    =H_min+a;
    Ht_min    =[Ht_min H1_min];

    
    gama      = (azimuth_angle/2)+teta;
    H_max     = half_node_distance*tan(gama);
    H1_max    =H_max+a;
    Ht_max    =[Ht_max H1_max];
    
     Teta      = teta*180/pi;
    
    semilogy(Teta,H1_min,'-.r')
    grid on
    title('Max&Min distance of coverage')
    xlabel('Rotation angle (Teta)-(deg)')
    ylabel('Distance(m)')
    hold on
    
    semilogy(Teta,H1_max,'-.g')
    hold on
end
H2_max  = half_node_distance*tan(gama_max)+a;

H2_min  = half_node_distance*tan(Teta_max)+a;

% disp('Ht_max=');
% disp(Ht_max);
% 
% disp('Ht_min=');
% disp(Ht_min);

Teta_max   =Teta_max*180/pi;

% semilogy(Teta_max,H2_min,'--r')
% hold on
    
semilogy(Teta_max,H2_max,'-.g')
hold off