import pandas as pd import numpy as np

读取Excel文件

data = pd.read_excel(r'D:\pythonProject3\会员信息\附件二：会员信息数据.xlsx')

获取商品和会员的位置信息

商品经度 = data['商品GPS经度'] 商品纬度 = data['商品GPS纬度'] 会员经度 = data['会员GPS经度'] 会员纬度 = data['会员GPS纬度']

获取商品价格

商品价格 = data['任务标价']

计算距离

距离 = np.sqrt((商品经度 - 会员经度)**2 + (商品纬度 - 会员纬度)**2)

计算引力子

引力子 = 商品价格 / (距离**2)

拟合常量k和n

k, n = np.polyfit(商品价格, 引力子, 1)

输出带常量的公式

print("y = {} * q1 * q2 / r^{}".format(k, n))