Double knife-edge diffraction

example Vavilov S. A., Lytaev M. S. Modeling Equation for Multiple Knife-Edge Diffraction //IEEE Transactions on Antennas and Propagation. – 2020. – Vol. 68. – Iss. 5. – pp. 3869-3877.

import os
os.chdir('../../')
from rwp.kediffraction import *
from rwp.antennas import *
from rwp.environment import *
from rwp.vis import *

Preparing environment

env = Troposphere(flat=True)
env.z_max = 150
env.knife_edges = [KnifeEdge(range=200, height=50), KnifeEdge(range=800, height=50)]
antenna = GaussAntenna(freq_hz=300e6, height=50, beam_width=15, elevation_angle=0, polarz='H')

Starting calculation

kdc = KnifeEdgeDiffractionCalculator(src=antenna, env=env, max_range_m=1000)
field = kdc.calculate()

Visualising results

vis = FieldVisualiser(field, env=env, trans_func=lambda v: 10 * cm.log10(1e-16 + abs(v)), label='ke')
plt = vis.plot2d(min=-40, max=0)
plt.title('The intensity of the field component 10log10|u|')
plt.xlabel('Range (m)')
plt.ylabel('Height (m)')
plt.show()

plt = vis.plot_hor(51)
plt.title('The intensity of the field component 10log10|u| at the height 51 m')
plt.xlabel('Range (m)')
plt.ylabel('10log|u| (dB)')
plt.show()