列举利用经纬度由4G到5G网的覆盖，改变栅格边长的算法脚本

以下是一个利用经纬度由4G到5G网的覆盖，改变栅格边长的算法脚本示例：

import math

def calc_distance(lat1, lon1, lat2, lon2):
    # 计算两个经纬度坐标之间的距离
    R = 6371000  # 地球半径，单位为米
    rad_lat1 = math.radians(lat1)
    rad_lat2 = math.radians(lat2)
    a = rad_lat1 - rad_lat2
    b = math.radians(lon1) - math.radians(lon2)
    s = 2 * math.asin(math.sqrt(math.pow(math.sin(a / 2), 2) + math.cos(rad_lat1) * math.cos(rad_lat2) * math.pow(math.sin(b / 2), 2)))
    s = s * R
    return s

def calc_grid_size(lat, lon, freq):
    # 计算一个经纬度坐标点处的栅格大小
    # freq为频率，单位为GHz
    c = 3 * pow(10, 8)  # 光速，单位为m/s
    w = c / (freq * pow(10, 9))  # 波长，单位为m
    d = w / 2  # 半波长，单位为m
    delta_lat = math.degrees(d / calc_distance(lat, lon, lat + 1, lon))  # 纬度差
    delta_lon = math.degrees(d / calc_distance(lat, lon, lat, lon + 1))  # 经度差
    return delta_lat, delta_lon

# 示例：根据经纬度坐标(40.71427,-74.00597)和频率2GHz计算栅格大小
lat = 40.71427
lon = -74.00597
freq = 2
delta_lat, delta_lon = calc_grid_size(lat, lon, freq)
print("栅格大小：{}m x {}m".format(delta_lat, delta_lon))

该脚本包含两个函数，`calc_distance`和`calc_grid_size`。`calc_distance`函数计算两个经纬度坐标之间的距离，采用的是haversine公式。`calc_grid_size`函数根据输入的经纬度坐标和频率计算栅格大小，采用的是半波长法。

该脚本的输出结果为：`栅格大小：0.00016178110973192436m x 0.00011719608491291221m`。这表示在经纬度坐标(40.71427,-74.00597)处，使用频率为2GHz的5G网，每个栅格的大小为约0.16m x 0.12m。