Instructions:

❼ Read this script carefully before attempting the experiment.

❼ Review MATLAB before attempting the experiment and use it for all the re- quired coding and graphs. Check CANVAS for MATLAB resources (Learning Resources→Supporting Material folder). See online material and resources as well.

❼ Keep a record of all code, graphs and results.

❼ The code must be well-structured and organised (to get a better mark). Use the concept of functions for better code re-usability and managment.

❼ The code must be provided for every requirement with appropriate explana- tion. Use MATLAB commenting (%) to emphasise on and explain relevant code lines.

❼ For every requirement, the code must run and screenshots of the results must be provided.

❼ Any change to the code needs to be reinserted every time.

❼ Be sure to reference any resource you have used in writing your report.

❼ Use your time wisely. Finish as many tasks as you can during the lab (with demonstrators support). If you can’t finish all tasks today, you can complete your work at home later on.

❼ If you have any feedback on your laboratory experience today, please write it down on the last page of the script.

1 Learning outcomes

The purpose of this experiment is to develop, explore and test Monte Carlo techniques in simulating and finding solutions to real-life random processes.  MATLAB will be used as the tool to do the tests of the experiment,  but it is not the main learning outcome  (i.e.   the experiment is not about MATLAB).

2 Introduction

The Monte Carlo method is a numerical method of solving mathematical problems by the sim- ulation of random variables.  The name Monte Carlo was applied to a class of mathematical methods first by scientists working on the development of nuclear weapons in Los Alamos in the 1940s. The essence of the method is the invention of games of chance whose behaviour and outcome can be used to study some interesting phenomena. While there is no essential link to computers, the effectiveness of numerical or simulated gambling as a serious scientific pursuit is enormously enhanced by the availability of modern digital computers [1].

The term “Monte Carlo” refers to procedures in which quantities of interest are approximated by generating many random realisations of a stochastic process and averaging them in some way. In statistics, the quantities of interest are the distributions of estimators and test statis- tics, the size of a test statistic under the null hypothesis, or the power of a test statistic under some specified alternative hypothesis [2].

How can we use Monte Carlo techniques to find the sampling distribution of an estimator? In the real world, we usually observe just one sample of a certain size N, which will give us just one estimate. The Monte Carlo experiment is a lab situation, where we replicate the real world study many (R) times.  Every time, we draw a different sample of size N from the original population.  Thus, we can calculate the estimate many times and any estimate will be a bit different.   The empirical distribution of these many estimates approximates the true of the estimator. A Monte Carlo experiment involves the following steps [3]:

(1) Draw a (pseudo) random sample of size N for the stochastic elements of the stochastic model from their respective probability distribution functions.

(2) Assume values for the parts of the model or draw them from their respective distribution function.

(3) Calculate the parts of the statistical model.

(4) Calculate the value (e.g. the estimate) you are interested in.

(5) Replicate step (1) to (4) R times.

(6) Examine the empirical distribution of the R values.

The Monte Carlo approach is relevant to different scientific disciplines and problems including (but not limited to) the following areas [4]:

❼ Physical sciences: computational physics, physical chemistry, quantum chromodynam- ics, statistical physics, molecular modelling, particle physics and galaxy modelling.

❼ Designs and visuals/Computer graphics: global illumination, photorealistic images of virtual 3D models, video games architecture and design, computer generated films and special effects in cinema.

❼ Finance and business/Operations research: evaluating investments in projects at a business unit, evaluating financial derivatives, construction of stochastic or probabilistic financial models and in enhancing the treatment of uncertainty in the calculation.

❼ Telecommunications: planning a wireless network, generating user patterns and their states, testing the probability of losing information in a network whether it is below a certain threshold.

❼ Games: game playing related artificial intelligence theory.

3 The Practical Work

Penalty kicks are a critical time of decision-making for both the goalkeeper and the penalty taker in football matches. Given that, for most professional games, the average number of goals scored is around 2.5, a penalty kick can have a major influence on the outcome of a match. Penalty kicks may reach speeds near 125 mph and is usually over within a quarter of a second. Thus, the goalkeeper must make a decision on how to stop the shot before the ball is struck. Statistics show that goalkeepers will most often jump to the left or right, hoping to guess cor- rectly the position to block the kick [5,6].

Consider the situation of a football goal and a blindfolded person trying to shoot the ball from the penalty spot and score a goal.  Let’s assume that the goal has dimensions L and W as shown in Figure 2, and there is an imaginary circle that circumscribes the goal. Two cases will be considered: first when there is no goalkeeper and second when there is a goalkeeper saving the ball.

Figure 2: The goal arrangement.

3.1 Part I: No Goalkeeper Tests (40 Marks)

In this case, there is no goalkeeper, and it is just the penalty taker against the goal. You need to model each shot by treating the co-ordinates of the ball in the goal plane as random variables (i.e. ignore the trajectory of the ball).

❼ Task-1. If a large number of shots is attempted, derive a numerical value for the fraction of balls entering the goal to the total number of balls in the circular area. Assume the penalty taker is blindfolded (i.e. the shots are uniformly distributed within the circle). [5 Marks]

❼ Task-2. Design and write a computer programme to find the probability of scoring by simulating N random penalty shots and repeating this experiment R times and taking the mean of the attempts.  Let N and R be inputs to your code.  Use a uniform random number generator in the simulation. [8 Marks]

❼ Task-3. Produce an appropriate scatter plot illustrating your experiment for N = 1,000 and R = 1, using red crosses to indicate score (i.e.  balls on target) and blue circles to indicate miss (i.e. balls off target). Insert an appropriate legend. [4 Marks]

❼ Task-4. For R = 5, find the probability of scoring for N = 100, N = 1,000, N = 10,000 and N = 100,000. Plot the probability against the value of N . Comment on the shape of the plot, making reference to the theoretical probability calculated in Task-1. Remember to label the axes and to insert an appropriate caption in your report. [7 Marks]

❼ Task-5. For N = 1,000, find the probability of scoring for R = 5 times, R = 10 times, R = 15 times and R = 20 times. Plot this probability against the value of R.  Comment on the shape of the plot. [4 Marks]

❼ Task-6. Compare with appropriate explanation the two cases of Task-4 and Task-5 based on the obtained probability plots. [4 Marks]

❼ Task-7. Repeat Task-2 to Task-6 using a normal (Gaussian) random number gener- ator. Assume the distribution to be centred at the centre of the circle and with standard deviation equal to the radius. Comment (with appropriate explanation) on the differences between the results of the two cases. [8 Marks]

3.2 Part II: With Goalkeeper Tests (30 Marks)

Consider now the above case but with a goalkeeper.  The goalkeeper can assume one of five possible actions (see Figure 3): stays in the middle, jumps to the upper left corner, jumps to the upper right corner, jumps to the lower left corner or jumps to the lower right corner.  A ball is saved if the goalkeeper guesses the correct ball position.  The goal area can be divided into five corresponding regions as shown in the figure.

Figure 3: Five possibilities of a goalkeeper action to a penalty shoot-out.

❼ Task-8. Assuming that the goalkeeper action is modelled as a uniform random process, what is the probability of scoring a goal if the penalty taker kicks 100 balls with uniform random distribution within the circle, as before. Increase the kicks to 1000. Compare the probability values with the case where no goalkeeper was in the goal (Task-1 and Task-3 above). [15 Marks]

❼ Task-9. Repeat Task-8 if the balls are kicked with a Gaussian random distribution (as in Task-7). Compare your results with those obtained in Task-7 and Task-8. [5 Marks]

❼ Task-10. Given the fact that statistically 90% of the time goalkeepers tend to jump to the lower two corners of the goal, what is the probability of scoring in this case after randomly kicking 100 balls?  1,000 balls?  (Compare both uniform and Gaussian distribu- tions) [10 Marks]

Note: For Tasks 8-10 you need to provide code, plots, explanations & comments as in Part I.

4 Review Questions (30 Marks)

(Include these in your Conclusions/Discussion section of your report)

Q1. In terms of what you’ve done in this experiment, comment on the advantages and disad- vantages (or drawbacks) of the Monte Carlo experiment. [5 Marks]

Q2. Discuss the ways in which the above model could be made more accurate and realis- tic. [7 Marks]

Q3. With reference to Task-7 and Task-9, discuss the effect of changing the standard devi- ation of the Guassian distribution on both the accuracy and precision of the penalty shots. [5 Marks]

Q4. If a large number of balls are kicked on the goal (i.e. if N is sufficiently large), the value of T can be estimated using (some function of) the ratio of the number of scores to the total number of the shots.  Hence, find the relation that estimates the value of T .  Verify this using your results for both uniform and Gaussian distributions. [8 Marks]

Q5. From your observation and results of Part II, what is the best strategy that should be adopted by the penalty taker?  What is the best strategy that should be adopted by the goalkeeper? [5 Marks]

References

[1] G Rubino and B Tuffin,  “Rare Event Simulation using Monte Carlo Methods”, Wiley,

2009.

[2] C Lemieux, “Monte Carlo and Quasi-Monte Carlo Sampling”, Springer, 2009. [3] M Kalos and P Whitlock, “Monte Carlo Methods”, Wiley, 2004.

[4] G Fishman, “Monte Carlo Concepts, Algorithms and Applications”, Springer, 1996.

[5] M Bar-Eli and O Azar, “Penalty kicks in soccer: an empirical analysis of shooting strategies and goalkeepers preferences”, Soccer & Society, 2009.

[6] Just     About      Football,       “Penalty      Kick      Statistics      and      Success     rates”, http://justaboutfootball.blogspot.com/2009/02/penalty-kick-statistics.html, 2009.

[7] J Billsberry, “Shootouts Alternatives”, http://www.penaltyshootouts.co.uk/alternatives.html,

2010.

Version history