In deze scriptie wordt onderzoek gedaan naar quasi-stationaire verdelingen, in dit geval in een Markovmodel voor virusverspreiding. Het onderzoek is gebaseerd op eerder onderzoek van J. Gani, waar...Show moreIn deze scriptie wordt onderzoek gedaan naar quasi-stationaire verdelingen, in dit geval in een Markovmodel voor virusverspreiding. Het onderzoek is gebaseerd op eerder onderzoek van J. Gani, waar een model is gegeven dat op kortetermijn een hele andere verdeling laat zien dan op de langetermijn. Dit zou verklaard kunnen worden aan de hand van quasi-stationariteit. Quasi-stationaire verdelingen zijn schijnbaar stabiele verdelingen die optreden voordat een Markovketen naar de stationaire verdeling convergeert. Dit soort verdelingen kunnen optreden als de convergentiesnelheid naar de stationaire verdeling erg traag is. In deze scriptie wordt naar twee zaken onderzoek gedaan. Ten eerste, wordt onderzocht hoe we dit gedrag in het Gani-model kunnen verklaren. Hiertoe introduceren we een uitdrukking waarmee we het model daadwerkelijk kunnen modelleren en onderzoeken we een aantal eigenschappen van dit model gebaseerd op theorie uit de Markovketens. Deze eigenschappen, gepaard met de definitie van quasi-stationariteit zal tot antwoord op deze vraag leiden. Ten tweede, wordt onderzocht of hetzelfde gedrag in een realistischer model ook optreedt. Het realistische model zal gebaseerd zijn op het Gani-model met enkele aanpassingen. Ook is onderzocht hoe dit gedrag in modellen te induceren is en welke factoren dus precies invloed hebben op dit gedrag.Show less
This thesis gives an overview of the Shannon Switching Game, Lehman’s necessary and sufficient conditions for the game to be short, cut or neutral and the corresponding strategies, and an...Show moreThis thesis gives an overview of the Shannon Switching Game, Lehman’s necessary and sufficient conditions for the game to be short, cut or neutral and the corresponding strategies, and an elaboration of Bruno and Weinberg’s application of Kishi and Kajitani’s results to finding a pair of disjoint cospanning trees in a graph if the game is short. The main contributions of this thesis are providing Lehman’s results and their proofs with graphical examples to make them easier to read, and the elaboration of Bruno and Weinberg’s algorithm mentioned above.Show less
In Computed Tomography, it is often standard to scan an object with the projection angles spread equidistantly over the full rotation of the object. If the number of projection angles is relatively...Show moreIn Computed Tomography, it is often standard to scan an object with the projection angles spread equidistantly over the full rotation of the object. If the number of projection angles is relatively small (e.g. because of radiation damage), the choice of the distribution of these angles is important. Especially if the object has certain main directions, in which the most important features are aligned, the reconstruction can be improved by selecting angles in these directions. In this thesis, possibilities for improving the choice of projection angles are investigated. Three angle selection algorithms are discussed: a greedy algorithm, a coordinate descent algorithm and an adaptation of the ensemble Kalman Filter algorithm. The algorithms try to find the minimum of a cost function, which is based on the L 2 -norm. The performance of these algorithms is shown on three computer generated phantoms and one real-world image, generated from scanning a wooden tree stem in the FleX-Ray scanner at CWI. The results of choosing angles with the algorithms are compared with choosing equidistant projection angles. The results show that especially the coordinate descent method is capable of finding projection angles that lower the cost function. In real life situations the true image is not available. Therefore, the possibility of using training data to estimate the cost function in that case, is investigated. We assume that these training samples come from the same probability distribution as the true image. Then, averaging over the cost function of these training samples improves the choice of projection angles with respect to equidistantly chosen angles.Show less
In this thesis, we will consider two versions of a discount single server system with exponentially distributed arrivals and departures on a discrete time scale. In the first model, to be discussed...Show moreIn this thesis, we will consider two versions of a discount single server system with exponentially distributed arrivals and departures on a discrete time scale. In the first model, to be discussed in Chapter 1, we add a reward for every customer who has been served, and a fee for each time step a customer is in service. We also add arrival control. The model is described in more detail in Section 1.1. The goal is to find an optimal strategy that minimizes the expected total discounted cost (Section 1.2). To find this strategy, we use a direct, explicit algorithm called value iteration. This strategy turns out to be a threshold strategy: a given amount of customers is allowed to enter the system, and once this amount of customers is present, any newly arriving customer is refused. Furthermore, we prove monotonicity properties of the model in Section 1.3. Lastly, the n-horizon cost function of value iteration is decomposed into a subsequent application of different operators (Section 1.4). In the second chapter, the model of Chapter 1 is used, with the addition of departure control. We add an extra, faster server, so we have a slower Server 1 and a faster Server 2, and we can choose at any transition moment which server to use. Using Server 2 involves extra costs per time unit. A detailed description of the model can be found in Section 2.1. In this model, we again want to find an optimal policy (Section 2.3), but now we first rewrite the explicit algorithm as a consecutive application of operators (Section 2.2). It turns out that we can split these operators into operators deciding whether or not an incoming customer should be accepted, and operators deciding which server to use. This was a surprising result, as we did not know in advance that the two decisions can be made independently from each other. Using these operators, we find that the optimal policy is a two-dimensional threshold strategy. In Section 2.4, we give a theorem on the relationship between the two thresholds. Furthermore, we prove the convergence of the thresholds, and thus the strategy, and of the nhorizon relative cost function, by using two initial functions such that one approaches the optimal value function from below and the other from above (Section 2.5). This shows that the choice of the starting functions affects the results, which was unknown so far. Finally, we prove that such functions exist and give a numeric example to demonstrate the results graphically (Section 2.6).Show less
To determine an optimal strategy for a multi-class server, with impatient customers, can be difficult. There are heuristics, such as the Whittle index, that perform well under restrictive...Show moreTo determine an optimal strategy for a multi-class server, with impatient customers, can be difficult. There are heuristics, such as the Whittle index, that perform well under restrictive conditions. We discuss a few of these heuristics and their merits and drawbacks. To overcome these drawbacks, we experiment with a data-based approach, by generating data using value iteration and predicting a strategy with a neural network. Furthermore, we explore the effect of different parameters on the performance of the neural network, for a particular case. Moreover, we look into more extensive methods such as combining the prediction of multiple neural networks and using neural networks to choose the training data.Show less
In his thesis, we will first outline the theory of stochastic cooperative game theory without transfer payments, as developed by Borm, Tijs and Timmer in [6]. We will compare and contrast it with...Show moreIn his thesis, we will first outline the theory of stochastic cooperative game theory without transfer payments, as developed by Borm, Tijs and Timmer in [6]. We will compare and contrast it with classical cooperative game theory. Subsequently, we will use this theory to model a two-player game in which owners of intermittent renewable energy sources (such as wind mills and solar panels) could cooperate in order to minimize their fine due to prediction errors of their production. We will describe two different kinds of cooperation and compare and contrast them in terms of their respective beneficiality to the players involved by computing the core. Furthermore, we will consider the case in which the prediction errors of the players are correlated. Subsequently, we compute the Shapley values of the two-player game, after which we generalize our results to the case of three participating players. Throughout this thesis, we stress the practical relevance of this work to real-world examples by pointing towards and incorporating previous experimental work on the actual prediction error distributions of renewable energy production devices.Show less
The supply chain network of companies that produce goods and deliver them to customers, can be captured in a mathematical model. ORTEC Supply Chain Design (OSCD) uses a mixed integer mathematical...Show moreThe supply chain network of companies that produce goods and deliver them to customers, can be captured in a mathematical model. ORTEC Supply Chain Design (OSCD) uses a mixed integer mathematical program to design a supply chain network while optimizing a certain output (e.g. minimize cost). The input parameters of such a program, like transport cost and customer demand, not only have to be estimated, but they also have ranges of uncertainty. Consequently, the output (the network design) is uncertain as well. In this work the impact of that uncertainty on the output, in the mathematical model of OSCD, is investigated. Generic effects that hold for several case studies of OSCD are presented. In addition a tool for precise estimation in unique cases is proposed. This work could be used to further extend the relationship between input and output of more complex case studies of OSCD.Show less
A formalization for the Paradigm coordination language is introduced, in terms of graph products and Markov processes. The n-arm robot problem is presented, where a robot composed of several...Show moreA formalization for the Paradigm coordination language is introduced, in terms of graph products and Markov processes. The n-arm robot problem is presented, where a robot composed of several jointed arms tries to plan and to move. We find that the n-arm robot problem reduces to a concurrent competition for the space in which the arms evolve. This competition for space reduces to the classical reader/writers problem. Moving in space is equivalent to writing and sensing in space is equivalent to reading. Several coordination schemes for the arms are developed in Paradigm. Multiple arms can plan paths concurrently, but reserving space to move has to be done sequentially, and moving can be done concurrently. The solutions are simulated in Matlab, and the code is provided. Thereupon, the solutions are compared, and it is shown that not all coordination schemes are equivalent. Finally an outline of future work is provided.Show less
In this thesis, we will give a way to build mechanical metamaterials. We will do this by using a triangular tiling, in which we put spins on the edges of the tiles. These spins have to point either...Show moreIn this thesis, we will give a way to build mechanical metamaterials. We will do this by using a triangular tiling, in which we put spins on the edges of the tiles. These spins have to point either into or out of the triangles and have to satisfy the rule that for every triangle two spins have to point out, and one in, or two spins point in and one out. If we can construct a tiling that is completely filled with these triangles in such a way that all spins on sides of adjacent triangles are pointing in the same direction, we will call this a feasible configuration. Firstly, we derive the number of feasible configurations for the tiling and consider a way to estimate these values. Secondly, we derive the distribution for the number of configurations when there are i spins on the boundary pointing in. Finally, we consider the number of spins in a periodic tiling that can be reversed, independent of all other spins, and derive an upper value and a lower boundary for this.Show less
In this thesis, we will give a way to build mechanical metamaterials. We will do this by using a triangular tiling, in which we put spins on the edges of the tiles. These spins have to point either...Show moreIn this thesis, we will give a way to build mechanical metamaterials. We will do this by using a triangular tiling, in which we put spins on the edges of the tiles. These spins have to point either into or out of the triangles and have to satisfy the rule that for every triangle two spins have to point out, and one in, or two spins point in and one out. If we can construct a tiling that is completely filled with these triangles in such a way that all spins on sides of adjacent triangles are pointing in the same direction, we will call this a feasible configuration. Firstly, we derive the number of feasible configurations for the tiling and consider a way to estimate these values. Secondly, we derive the distribution for the number of configurations when there are i spins on the boundary pointing in. Finally, we consider the number of spins in a periodic tiling that can be reversed, independent of all other spins, and derive an upper value and a lower boundary for this.Show less