Le portail francophone de la géomatique

include"mapbasic.def"
include"menu.def"
include"icons.def"
declare sub calcul

diaLog title "simulation"
control button title "simuler"
calling calcul

sub calcul
dim alpha,j as integer
dim PA,k,rcouv,tl,Az,x,y,hi,f,hr,pc,ls,pind,pf,s,m,di,eval,rad,gn as float
dim aa,bb,cc,pire,Aaa,gama,xf,xh,PL1,PR,a,b as float
dim purple as integer
purple = RGB (210,25,255)

define gr 0
open table "c:\BS.TAB" 
open table "c:\orientation.TAB" as o
open table "c:\Rxleve1.TAB" as r
set map layer 1 editable on
k = frontwindow()

fetch first from BS 
  while not eot (BS)
     PA= BS.PA_Power
     rcouv= BS.Rayon_couverture
     tl= BS.Tilt
     Az= BS.Azimuth
     x= BS.X_coord
     y= BS.Y_coord
     hi= BS.Hauteur_Bs
     f=  BS.Fréquence
     hr= BS.Hauteur_R
     pc=BS.Affaiblissement
     ls= BS.P_masque
     pind = BS.P_indoor
     pf= BS.P_fading
     s= BS.Shadowing
     m= BS.Marge
     di= BS.distance_incrément
     eval=  BS.Evaluation
     rad= BS.Radials     
                                                  
    Az= Az-180

  
fetch first from o
  for alpha =0 to 17 step 1 
    gn= o.gain
      [url][/url][url][/url]
      for j= 1 to (rcouv*1000) step 50
      j=j*1000
        Pire= PA + gn - pc
Aaa = 32.4 + 20* (Log(100)/Log(10)) + 20*(Log(f)/Log(10))
gama = aa - bb * hi + cc/hi 
xf =6 * (Log (f/1900)/Log(10))
xh =-10.7 * (Log (hr/2)/Log(10))
aa= 4.6
bb= 0.0075
cc=12.6

PL1= 32.4 + 20 * (Log(j)/Log(10)) +20*(Log(f)/Log(10))
PR = Pire + Gr- PL1- ls- pf -m -pind

a=j * cos((alpha + Az)* deg_2_rad) + x
b=j * sin((alpha + Az)* deg_2_rad) + y





insert into r ( X_coord,Y_coord, Rxlev)values (a,b,PR)
            if PR>=-55 then 
                 create rect
                 into window K
                 (a-25,b+25) (a+25,b+25)
                  pen makepen(0,1,purple)
                  brush makebrush (02,red,red)

            elseif PR>=-65 then 
                 create rect
                 into window K
                 (a-25,b+25) (a+25,b+25)
                  pen makepen(0,1,purple)
                  brush makebrush (02,yellow,yellow)

           elseif PR>=-75 then 
                  create rect
                  into window K
                  (a-25,b+25) (a+25,b+25)
                   pen makepen(0,1,purple)
                   brush makebrush (02,blue,blue)

          elseif PR>=-85 then 
                   create rect
                   into window K
                   (a-25,b+25) (a+25,b+25)
                   pen makepen(0,1,purple)
                   brush makebrush (02,green,green)

           end if
    j=j/1000
   fetch next from  o
           next   
   fetch next from BS 
      next 
   wend
end sub