Single knife-edge diffraction

example from 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=500, height=50)]
antenna = Source(wavelength=1, height_m=50)

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(50)
plt.title('The intensity of the field component 10log10|u| at the height 50 m')
plt.xlabel('Range (m)')
plt.ylabel('10log|u| (dB)')
plt.show()