使用python给出‘费涅尔绕射的刃形绕射增益’的代码

这是一个Python代码例子，用于计算费涅尔绕射的刃形绕射增益：import numpy as np

def far_field_gain(r, theta):
"""
Calculates the far field gain at a given distance and
angle from a uniformly-illuminated circular aperture.

Parameters
----------
r : float
Distance from aperture center in meters.
theta : float
Angle from aperture center in radians.

Returns
-------
float
Gain in dB.
"""
k = 2 * np.pi # wave number
G = 2 * np.pi * r**2 / (k * np.sin(theta)) * np.abs(np.cos(theta))
return 20 * np.log10(G)

# Compute the far field gain at a distance of 1 meter and an angle
# of 30 degrees
print(far_field_gain(1.0, np.radians(30.0)))

相关问题

使用python给出‘费涅尔绕射的刃形绕射增益’的代码，并对代码进行合理缩进

示例代码：def feynman_diffraction_gain(wavelength, width, spacing):
"""
计算费涅尔绕射的刃形绕射增益

参数：
wavelength: 波长
width: 刃宽
spacing: 刃间距

返回：
费涅尔绕射的刃形绕射增益
"""
k = 2 * np.pi / wavelength # 波矢
G = (np.sin(k * spacing) / (k * spacing)) ** 2 # 刃形绕射增益
return G * np.sin(np.pi * width / wavelength) ** 2