足球的最佳射门位置

　　足球是世界上最受欢迎的体育运动，每四年举行一次的国际足联俱乐部（FIFA）世界杯足球赛更是吸引了全世界无数球迷观看。一场足球比赛中最激动人心的时刻莫过于射门和进球，而射门和进球也是一支球队能否获胜的关键。一记射门能否进球得分取决于多种因素，这些因素不仅包括运动员本身的球技，也包括射门所在位置。
　　射门位置离球门越远，越不容易进球。因此，一记超远吊射也常常为人津津乐道。2022年波兰足球联赛扎布热对战坡根的比赛中，足球老将波多尔斯在中场线附近以一记超远吊射破门得分（图1）。著名球星贝克汉姆和鲁尼也用中场吊射的方式进过球。
　　此外，射门位置的球门视角也极大影响着进球的难易程度，正对着球门方向的射门进球的几率要大于两边。虽然小视角射门得分的难度很大，但在足球比赛中也是出现过的。例如，在2010年世界杯巴西对战朝鲜的小组赛中，麦孔在底线附近接近零度视角的情况下射门得分（图2）。
　　甚至球门视角为零的角球直接射门得分的情况也是存在的。例如，在201920赛季西班牙超级杯半决赛皇家马德里对战克瓦伦西亚的比赛中，克罗斯就用角球直接射门得分（图3）。著名球星贝克汉姆和亨利也完成过角球得分的名场面。
　　相同球门距离的位置可能有着不同的球门视角，而相同球门视角的位置也可能有着不同的球门距离（图4）。不同的射门位置对进球的影响有多大？或者说如果在某个特定的位置射门，进球的预期概率是多少？本文将从理论和统计两个角度分别建立模型回答这个问题。
　　模型
　　理论模型
　　射门位置到球门距离和射门位置的球门视角是影响射门能否得分的两个重要因素。为研究球门距离和球门视角对能否进球的影响，本文建立了如图5所示的坐标系。球门宽度为7。32m，射门位置（，）到球门（线中点）的距离为，到两根球门柱的距离分别为和，球门视角为。
　　射门位置到球门的距离可以很容易地由射门位置的坐标（，）表示射门位置的球门视角也可以由射门位置的坐标（，）表示。由正弦定理可知由上式及余弦定理可得如果假定，则上式中的和均可表示为和的函数因此有根据三角函数关系容易解得接下来，本文定义得分势函数来表示一记射门得分的难易程度。如上所述，得分势函数中的主要参数应该是球门距离和球门视角。直观上很容易理解，射门位置离球门越近，进球就越容易。因为球在飞行的过程中速度会降低，球在空中飞行的距离越长，球到达球门的速度和力量就越小，就越容易被守门员守住。此外，射门位置离球门越远，球在空中飞行的时间就越长，留给守门员的反应时间也就越长。另一方面，球门视角越大也就越容易得分。因为球门视角越小，越不容易瞄准。因此，得分势函数可以定义为与球门距离成反比，与球门视角成正比：根据上式可以绘制出图6所示的得分势函数及其等势线。得分势函数反映了在不同位置射门时进球的相对难易程度。
　　在同一条等势线上的位置射门进球的难易程度相当，而势越大的位置射门进球越容易。例如，在图6中的绿线上各位置射门进球的难易程度相当，在绿色等势线上各位置射门进球的难度要高于红色等势线。
　　统计模型
　　需要注意的是，理论模型给出的得分势函数反映的是得分相对难易程度，并不是一记射门能否得分的直接概率。为了验证理论模型的合理性，以及计算球场上不同位置射门的得分概率，需要统计大量真实的足球比赛数据。本文接下来的分析使用了由Wyscout收集的足球比赛事件数据集。该数据集记录了七个赛事（西甲、意甲、德甲、英超、法甲、2018年FIFA世界杯、2016年欧洲杯）整个赛季的所有比赛事件（传球、射门、犯规等）的时空数据（包含事件的有关位置、时间、结果、球员和特征等信息）。数据集中各事件的位置都由一对坐标（，）表示，和是范围内的整数。坐标表示事件发生点到对方球门底线边界的距离（以百分比表示），而坐标表示事件发生点到场地右侧边界的距离（以百分比表示）。
　　由于本文主要关注于射门的位置和射门的结果（是否得分），因此将球场在长度和宽度方向等间距地划分为5050个格子，并统计了每个格子中发生射门事件的次数，结果如图7。然而，一个格子射门的次数多并不意味着进球的次数多。为此，本文还统计了每个格子中发生射门事件且得分的次数，结果如图8。
　　由此，可以计算出每个格子中射门的得分频率结果如图9，如果统计的射门事件足够多，那么每个格子中射门的得分频率将趋于得分概率。显然，由于有限的数据，本文获得的得分频率图并不光滑连续。
　　但从空间得分频率图不难发现，离球门越近，越正对着球门的位置越容易得分。需要注意的是，图9中靠近边界的少量格子存在一些高频率的异常格子。这实际上是由于发生在这些格子中的射门次数很少，不具有统计意义。例如某个格子中只发生一次射门，并且这次射门进球了，那么该格子的进球频率就会为1。
　　为了估计空间上每个位置射门的得分概率，需要对得分频率进行拟合。根据前文理论模型的分析，某个位置射门得分的难易程度与该位置到球门距离和球门视角有关。为此，本文由各格子的坐标计算出了各格子的球门距离和球门视角。并在此基础上拟合进球概率对和的依赖关系。这个问题中，进球与否是离散的二分类因变量，而球门距离和球门视角则是连续的自变量。对于因变量为二分类或多分类时，逻辑回归是非常重要的模型。逻辑回归模型的表达形式为逻辑回归是通过极大似然估计求解常数项和偏回归系数（1，2，，），基本思想是当从总体中随机抽取个样本后，最合理的参数估计量应该使得这个样本观测值的概率最大。逻辑回归模型的概率表达式为逻辑回归可以很容易由MATLAB函数glmfit实现。因此，本文选用逻辑回归模型对进球概率与球门视角和球门距离进行分析。首先，本文单独考虑进球概率对两个因素的各自依赖关系。利用逻辑回归，可以得到进球概率与球门视角的依赖关系：相应的拟合结果如图10所示。从图中不难看出，当球门视角逐渐由0增大到120时，进球概率也从几乎为0增加到接近1。
　　类似地，可以得到进球概率与球门距离的依赖关系：相应的拟合结果如图11所示。从图中不难看出，进球概率随着球门距离的增大，逐渐趋于0。距离球门超过20米的位置射门进球就比较困难了。
　　接下来，本文考虑进球概率同时对球门视角和球门距离的依赖关系。利用逻辑回归，可以得到：相应的空间概率分布如图12所示。从概率图中不难看出，球门距离越近、球门视角越大的位置越容易射门进球。对比图12和图6不难发现，逻辑回归得到的空间概率分布与理论模型给出的得分势函数具有很大的相似性。图12中30等概率曲线与图6中1等势线与几乎一致。但两图也有一些区别，图6中的等势线都通过两个球门柱，但图12中较小概率（7）的等概率曲线与底线的交点并不在门柱位置。
　　在此基础上，本文还综合考虑了射门位置到底线的距离、到侧边的距离、球门距离与球门视角的积等多种因素对得分概率的影响。通过对各拟合结果的比较分析，本文最终得到并选择如下模型作为射门得分概率的最优模型：相应的空间概率分布如图13所示。从图中不难看出，最优模型的等概率曲线都经过门柱。这与角球直接射门得分的事实并不矛盾，因为本文所统计的射门都是发生在正常进攻过程中的，并不包括角球射门的情况。实际上，动态攻守过程中，在底线附近直接射门进球是极为困难的。
　　通过等概率曲线可以看出，把球盘带到绿线以内射门，进球的概率就能超过7。而把球盘带到红线以内再射门，进球的概率就能超过15。
　　结论
　　不同的射门位置有着不同的球门距离和球门视角，射门位置很大程度地影响着能否进球。基于射门位置的球门距离和球门视角，本文分别从理论和统计两种角度建立模型研究了在不同位置射门的进球预期概率。在理论模型中，本文构造了得分势函数。得分势函数能给出不同位置射门时进球的相对难易程度。为了得到不同位置射门得分的概率，本文对足球比赛事件数据进行了统计分析，应用逻辑回归模型对进球概率与球门视角、球门距离等因素的依赖关系进行了拟合，并最终给出了空间概率图和等概率曲线。
　　本文的结果给出了球场各位置射门得分的预期概率，这可以告诉球员更应该在哪些位置射门。总体而言，与我们的直觉是一致的：射门位置到球门距离越近、越正对着球门，射门进球的机会越大。当然，这个概率并不适用于所有球员，但得分概率图可以让所有球员了解到不同位置射门得分的难易程度，以便球员权衡自身能力和对方防守情况做出合理的决策。此外，得分概率图的梯度方向还为球员指明了进攻方向。例如，图13中绿线附近射门得分的机会是7，而沿梯度方向跑几步到达红线附近再射门就可以使得分机会翻倍。
　　参考资料
　　〔1〕
　　CsarAMorales。Amathematicsbasednewpenaltyareainfootball：tacklingping。Journalofsportssciences，34（24）：22332237，2016。
　　〔2〕
　　LucaPappalardoandetal。Apublicdatasetofspatiotemporalmatcheventsinsoccercompetitions。Scientificdata，6（1）：115，2019：
　　〔3〕
　　DavidSumpter。Mathematicalmodellingoffootball，2020：
　　〔4〕
　　DavidSumpter。Soccermatics：MathematicalAdventuresintheBeautifulGameProEdition。BloomsburyPublishing，2017。238241。〔5〕
　　Wikipediacontributors。LogisticregressionWikipedia，thefreeencyclopedia，2022：〔6〕
　　MathWorks。glmfit：Fitgeneralizedlinearregressionmodel。2022：
　　转载内容仅代表作者观点
　　不代表中科院物理所立场
　　如需转载请联系原公众号
　　来源：数学模型
　　编辑：QuantumBard