COMPUTER EXERCISE: POPULATION DYNAMICS IN SPACE September 3, 2013
|
|
- Sharon Quinn
- 5 years ago
- Views:
Transcription
1 COMPUTER EXERCISE: POPULATION DYNAMICS IN SPACE September 3, 2013 Objectives: Introduction to coupled maps lattice as a basis for spatial modeling Solve a spatial Ricker model to investigate how wave speed changes as a function of dispersal rate and recruitment Introduction Preliminary. Cellular automata and coupled map lattices are two approaches for studying population dynamics distributed in space, represented as a grid of adjacent sites. In either case the study involves evaluating a discrete time dynamical system for the growth of the population at each site followed by a dispersal step wherein individuals are redistributed to nearby sites according to some rule. The difference between cellular automata and coupled map lattices is that in cellular automata the state of the sites is binary (occupied/unoccupied) whereas in a coupled map lattice the population size may be any value. Cellular automata and coupled map lattices differ from the reaction-diffusion model in that space is represented as a collection of discrete sites. These models also differ from most metapopulation models in that sites are considered to be adjacent and dispersal is only to sites within a neighborhood of origination. Cellular automata and coupled map lattices provide a good technique for studying population growth within a continuous space where analytical solutions are intractable (e.g., nonlinearity of equations). Boundary conditions. One issue that always comes up in numerical analysis of spatial population dynamics is what to do at the boundaries of the rectangular grid. One possibility is to allow the left boundary to be adjacent to the right boundary and for the top boundary to be adjacent to the bottom boundary. This assumption allows neatly for a continuous space and implies that the spatial geometry of the system is actually a torus, as if the population dynamics were of a microbial population on the surface of a donut. Illustration of a torus. Image:
2 An alternative we will use here is to assume that the grid represents a habitat patch in a matrix of inhospitable territory, in which case individuals that disperse off the grid never return. Population dynamics. In the assigned reading we looked at the spread of a population that was (locally) growing exponentially and dispersing according to the diffusion assumption, giving rise to the reaction diffusion equation. This scenario yielded traveling waves with an asymptotically constant speed proportional to the square root of the product of the intrinsic rate of increase and diffusivity. c*=2 rd In contrast, it is also known that growth equations with leptokurtic (fat-tailed) dispersal kernels we representing long distance dispersal give rise to accelerating traveling waves. But what happens when population dynamics are density-dependent, for instance in fisheries where population growth is often represented by the Ricker stock-recruitment model? n n+ 1 =n t e r (1 n t/ k ) To answer this question (and introduce along the way such techniques as image plots), we will study a coupled map lattice in which the Ricker equation gives the local population dynamics in each each ( larval recruitment ) followed by a redistribution step ( local dispersal ). Exercise 1. The Ricker model defined above will be central to this lab exercise. Since we will often need to update population size, the first step is to write a function called ricker to update population size according to this equation. Spend a little time familiarizing yourself with this model. Perhaps you will want to simulate some trajectories with different values of the parameter r (which is conventionally assumed to represent potential population growth) and k (which governs density dependence, although not in exactly the same way as the parameter with the same name in the logistic model.) Assume carrying capacity is 100. Plot the so-called stock recruitment curve, the line that relates population size at time t to population size at time t+1 for a range of values of r. (Hint: it may be useful to look at values of r along a logarithmic scale. I used values between 2-2 and 2 1/2 ). What are the equilibrium population sizes of the Ricker model? What values of r give rise to population growth? 2. To confirm that your function is working, plot some trajectories for local population growth. As above, assume carrying capacity is 100; initially set r=0.25. How would you describe the dynamics of this system? It is well known that the Ricker model may produce complex dynamics due to overcompensation. (What do you suppose overcompensation is?) Particularly, chaos is approached via the period doubling bifurcation. Find a value of r at which the Ricker model produces a stable two-cycle. What value do you come up with? 3. Now we are in a position to proceed with our coupled map lattice. To start we need to set up an array that we will think of as the space on which population dynamics are occurring. So that your computations don't take too long, I would recommend a lattice no larger than 100,000 sites. (I used 1,000 sites as a trial and that worked just fine.) Use the function matrix to set up a square array of the right dimensions. For later use, initialize the array with the value 0 at all
3 sites. 4. Now, unless you did something unusual in writing the function ricker, it should actually be able to update the entire matrix in one step. That is, if you have an array N that contains the population size at each site, you should be able to compute the population size at the next time step using the code N< ricker(n,r=0.25,k=100). What we lack, however, is a function to update the spatial distribution after dispersal. To be more concrete, our target a computer program executing the following operations: Initialize variables For all times iterate the following i. update population size ii. redistribute iii. store output Inspect results Therefore, the next step is to write a function dispersal that will accomplish the redistribution in step (ii). Assume that this function will take two arguments, the current state of the system (your array N) and a site-to-neighbor dispersal rate m. (We will assume that dispersal is to the von Neumann neighborhood, the four cells orthogonally surrounding a central cell.) A word to the wise: this is probably the trickiest part of the whole exercise. (Hint: Pay close attention to the boundaries. Double hint: let the outermost border of sites be hostile matrix. This simplifies the code considerably). 5. Now you have all the conceptual tools you need to study the Ricker coupled map lattice and do determine its spread rate. Here are a couple of programming tools that you may find helpful. The function image takes a matrix and makes a heat map from it. This is useful for visualizing the state of the system. The operator %in% allows you to test if a value is an element of a set. My code uses the following line to plot the state of the system only every so often rather than at each update step: if(t %in% plot.times) image(n, col=heat.colors(100), zlim=c(0,100)) The function par allows one to reset graphical parameters. Particularly, the argument mfrow lets one set up a multi-panel plot by specifying a grid (composed of a number of rows and a number of columns) of individual plots that are successively populated. My code uses the following line to achieve this. x11(w=12, h=12); par(mfrow=c(5,5)) Starting with the values used above (r=0.25, k=100) and assuming that 1% of individuals move in each of the cardinal directions, iterate the coupled map lattice and watch the invasion play out (this will take a couple hundred time steps). If your analysis generates that following plot, then you have been successful!
4 6. Now you are prepared to answer our question: How does invasion speed change over time in the Ricker coupled map lattice? Answering this question will require (1) developing a method to track the maximum distance the population has spread from the origin, and (2) inspecting how this quantity changes over time. (Hint: A construction nesting two R functions may be helpful here, i.e., min(which(...)) where the ellipsis indicates that I have left out some of the code.) How does spread rate change over time? Does the invasion reach a constant speed? Accelerate? Decelerate? How does spread rate change with different levels of population growth (different values of r)? Different levels of dispersal (different values of m)?
5 Different levels of density dependent (different values of k)? 7. Another useful visualization is the profile of population size along a transect, i.e., if I were to choose a starting point and walk across the space in one direction (a transect) measuring the local abundance of the population as I go, how would abundance change in the course of my trip? Return to your main program and seek do develop such a profile plot. Compare this result with the image plot obtained above. Do these show different pieces of information. 8. If you've made it this far you've come a long way toward understanding how invasion speed is affected by density-dependent population growth. The Ricker model allows us to ask one more question however. As we discovered above and will study in a more thorough way in a few weeks, the Ricker model allows for persistent oscillations. This raises the question about how fast invasions will happen when dynamics are complex as a result of overcompensation? To answer this question, iterate the coupled map lattice using the value of r for which you previously showed there to be a stable two-cycle. Generate the corresponding image plots and profile plots. What is the effect of overcompensation on the invasion process? To demonstrate your work, turn in written answers to these questions and any plot that may help illustrate how you arrived at your answer. Due: September 10, 2013
CELLULAR AUTOMATA IN MATHEMATICAL MODELING JOSH KANTOR. 1. History
CELLULAR AUTOMATA IN MATHEMATICAL MODELING JOSH KANTOR 1. History Cellular automata were initially conceived of in 1948 by John von Neumann who was searching for ways of modeling evolution. He was trying
More informationGPU-based Distributed Behavior Models with CUDA
GPU-based Distributed Behavior Models with CUDA Courtesy: YouTube, ISIS Lab, Universita degli Studi di Salerno Bradly Alicea Introduction Flocking: Reynolds boids algorithm. * models simple local behaviors
More informationMETAPOPULATION DYNAMICS
16 METAPOPULATION DYNAMICS Objectives Determine how extinction and colonization parameters influence metapopulation dynamics. Determine how the number of patches in a system affects the probability of
More informationMapping Distance and Density
Mapping Distance and Density Distance functions allow you to determine the nearest location of something or the least-cost path to a particular destination. Density functions, on the other hand, allow
More informationDevelopment of a Maxwell Equation Solver for Application to Two Fluid Plasma Models. C. Aberle, A. Hakim, and U. Shumlak
Development of a Maxwell Equation Solver for Application to Two Fluid Plasma Models C. Aberle, A. Hakim, and U. Shumlak Aerospace and Astronautics University of Washington, Seattle American Physical Society
More informationLab 3: From Data to Models
Lab 3: From Data to Models One of the goals of mathematics is to explain phenomena represented by data. In the business world, there is an increasing dependence on models. We may want to represent sales
More informationEpidemic spreading on networks
Epidemic spreading on networks Due date: Sunday October 25th, 2015, at 23:59. Always show all the steps which you made to arrive at your solution. Make sure you answer all parts of each question. Always
More informationUNIT 9C Randomness in Computation: Cellular Automata Principles of Computing, Carnegie Mellon University
UNIT 9C Randomness in Computation: Cellular Automata 1 Exam locations: Announcements 2:30 Exam: Sections A, B, C, D, E go to Rashid (GHC 4401) Sections F, G go to PH 125C. 3:30 Exam: All sections go to
More informationCellular Automata. Nicholas Geis. January 22, 2015
Cellular Automata Nicholas Geis January 22, 2015 In Stephen Wolfram s book, A New Kind of Science, he postulates that the world as we know it and all its complexities is just a simple Sequential Dynamical
More information10.2 Diffusion and Cellular Automata
10.2 Diffusion and Cellular Automata Simulating Motion: Cellular Automata If all we have to work with is a grid of cells (spreadsheet), how can we simulate a random walk? Moving a value from one cell to
More informationarxiv: v1 [cond-mat.dis-nn] 30 Dec 2018
A General Deep Learning Framework for Structure and Dynamics Reconstruction from Time Series Data arxiv:1812.11482v1 [cond-mat.dis-nn] 30 Dec 2018 Zhang Zhang, Jing Liu, Shuo Wang, Ruyue Xin, Jiang Zhang
More informationApplication of Finite Volume Method for Structural Analysis
Application of Finite Volume Method for Structural Analysis Saeed-Reza Sabbagh-Yazdi and Milad Bayatlou Associate Professor, Civil Engineering Department of KNToosi University of Technology, PostGraduate
More informationTwo-dimensional Totalistic Code 52
Two-dimensional Totalistic Code 52 Todd Rowland Senior Research Associate, Wolfram Research, Inc. 100 Trade Center Drive, Champaign, IL The totalistic two-dimensional cellular automaton code 52 is capable
More informationL Modeling and Simulating Social Systems with MATLAB
851-0585-04L Modeling and Simulating Social Systems with MATLAB Lecture 4 Cellular Automata Karsten Donnay and Stefano Balietti Chair of Sociology, in particular of Modeling and Simulation ETH Zürich 2011-03-14
More informationContinuum-Microscopic Models
Scientific Computing and Numerical Analysis Seminar October 1, 2010 Outline Heterogeneous Multiscale Method Adaptive Mesh ad Algorithm Refinement Equation-Free Method Incorporates two scales (length, time
More informationLecture VII : Random systems and random walk
Lecture VII : Random systems and random walk I. RANDOM PROCESSES In nature, no processes are truly deterministic. However, while dealing with many physical processes such as calculating trajectories of
More informationHomework # 2 Due: October 6. Programming Multiprocessors: Parallelism, Communication, and Synchronization
ECE669: Parallel Computer Architecture Fall 2 Handout #2 Homework # 2 Due: October 6 Programming Multiprocessors: Parallelism, Communication, and Synchronization 1 Introduction When developing multiprocessor
More information1. Mathematical Modelling
1. describe a given problem with some mathematical formalism in order to get a formal and precise description see fundamental properties due to the abstraction allow a systematic treatment and, thus, solution
More informationConway s Game of Life Wang An Aloysius & Koh Shang Hui
Wang An Aloysius & Koh Shang Hui Winner of Foo Kean Pew Memorial Prize and Gold Award Singapore Mathematics Project Festival 2014 Abstract Conway s Game of Life is a cellular automaton devised by the British
More informationHONORS ACTIVITY #2 EXPONENTIAL GROWTH & DEVELOPING A MODEL
Name HONORS ACTIVITY #2 EXPONENTIAL GROWTH & DEVELOPING A MODEL SECTION I: A SIMPLE MODEL FOR POPULATION GROWTH Goal: This activity introduces the concept of a model using the example of a simple population
More informationMath 1525: Lab 4 Spring 2002
Math 1525: Lab 4 Spring 2 Modeling---Best Fit Function: In this lab we will see how to use Excel to find a "best-fit equation" or model for your data. Example: When a new motion picture comes out, some
More informationAn Evolution of Mathematical Tools
An Evolution of Mathematical Tools From Conceptualization to Formalization Here's what we do when we build a formal model (or do a computation): 0. Identify a collection of objects/events in the real world.
More informationCellular Automata. Cellular Automata contains three modes: 1. One Dimensional, 2. Two Dimensional, and 3. Life
Cellular Automata Cellular Automata is a program that explores the dynamics of cellular automata. As described in Chapter 9 of Peak and Frame, a cellular automaton is determined by four features: The state
More informationParallel Algorithms: Adaptive Mesh Refinement (AMR) method and its implementation
Parallel Algorithms: Adaptive Mesh Refinement (AMR) method and its implementation Massimiliano Guarrasi m.guarrasi@cineca.it Super Computing Applications and Innovation Department AMR - Introduction Solving
More informationLandscape Ecology. Lab 2: Indices of Landscape Pattern
Introduction In this lab exercise we explore some metrics commonly used to summarize landscape pattern. You will begin with a presettlement landscape entirely covered in forest. You will then develop this
More informationExponential and Logarithmic Functions. College Algebra
Exponential and Logarithmic Functions College Algebra Exponential Functions Suppose you inherit $10,000. You decide to invest in in an account paying 3% interest compounded continuously. How can you calculate
More informationSTANDARDS OF LEARNING CONTENT REVIEW NOTES ALGEBRA II. 3 rd Nine Weeks,
STANDARDS OF LEARNING CONTENT REVIEW NOTES ALGEBRA II 3 rd Nine Weeks, 2016-2017 1 OVERVIEW Algebra II Content Review Notes are designed by the High School Mathematics Steering Committee as a resource
More informationCS 1567 Intermediate Programming and System Design Using a Mobile Robot Aibo Lab3 Localization and Path Planning
CS 1567 Intermediate Programming and System Design Using a Mobile Robot Aibo Lab3 Localization and Path Planning In this lab we will create an artificial landscape in the Aibo pen. The landscape has two
More informationPlanar Graphs with Many Perfect Matchings and Forests
Planar Graphs with Many Perfect Matchings and Forests Michael Biro Abstract We determine the number of perfect matchings and forests in a family T r,3 of triangulated prism graphs. These results show that
More informationLosAlamos National Laboratory LosAlamos New Mexico HEXAHEDRON, WEDGE, TETRAHEDRON, AND PYRAMID DIFFUSION OPERATOR DISCRETIZATION
. Alamos National Laboratory is operated by the University of California for the United States Department of Energy under contract W-7405-ENG-36 TITLE: AUTHOR(S): SUBMllTED TO: HEXAHEDRON, WEDGE, TETRAHEDRON,
More informationComparison of different solvers for two-dimensional steady heat conduction equation ME 412 Project 2
Comparison of different solvers for two-dimensional steady heat conduction equation ME 412 Project 2 Jingwei Zhu March 19, 2014 Instructor: Surya Pratap Vanka 1 Project Description The purpose of this
More informationWatershed Sciences 4930 & 6920 GEOGRAPHIC INFORMATION SYSTEMS
HOUSEKEEPING Watershed Sciences 4930 & 6920 GEOGRAPHIC INFORMATION SYSTEMS CONTOURS! Self-Paced Lab Due Friday! WEEK SIX Lecture RASTER ANALYSES Joe Wheaton YOUR EXCERCISE Integer Elevations Rounded up
More informationModelling and Quantitative Methods in Fisheries
SUB Hamburg A/553843 Modelling and Quantitative Methods in Fisheries Second Edition Malcolm Haddon ( r oc) CRC Press \ y* J Taylor & Francis Croup Boca Raton London New York CRC Press is an imprint of
More informationModeling and Simulating Social Systems with MATLAB
Modeling and Simulating Social Systems with MATLAB Lecture 4 Cellular Automata Olivia Woolley, Tobias Kuhn, Dario Biasini, Dirk Helbing Chair of Sociology, in particular of Modeling and Simulation ETH
More information6.001 Notes: Section 4.1
6.001 Notes: Section 4.1 Slide 4.1.1 In this lecture, we are going to take a careful look at the kinds of procedures we can build. We will first go back to look very carefully at the substitution model,
More informationSensitivity Analysis. Nathaniel Osgood. NCSU/UNC Agent-Based Modeling Bootcamp August 4-8, 2014
Sensitivity Analysis Nathaniel Osgood NCSU/UNC Agent-Based Modeling Bootcamp August 4-8, 2014 Types of Sensitivity Analyses Variables involved One-way Multi-way Type of component being varied Parameter
More informationRaster Analysis and Functions. David Tenenbaum EEOS 465 / 627 UMass Boston
Raster Analysis and Functions Local Functions By-cell operations Operated on by individual operators or by coregistered grid cells from other themes Begin with each target cell, manipulate through available
More informationHigh-Performance Computing
Informatik und Angewandte Kognitionswissenschaft Lehrstuhl für Hochleistungsrechnen Rainer Schlönvoigt Thomas Fogal Prof. Dr. Jens Krüger High-Performance Computing http://hpc.uni-duisburg-essen.de/teaching/wt2013/pp-nbody.html
More informationA Source Localization Technique Based on a Ray-Trace Technique with Optimized Resolution and Limited Computational Costs
Proceedings A Source Localization Technique Based on a Ray-Trace Technique with Optimized Resolution and Limited Computational Costs Yoshikazu Kobayashi 1, *, Kenichi Oda 1 and Katsuya Nakamura 2 1 Department
More informationLAB #2: SAMPLING, SAMPLING DISTRIBUTIONS, AND THE CLT
NAVAL POSTGRADUATE SCHOOL LAB #2: SAMPLING, SAMPLING DISTRIBUTIONS, AND THE CLT Statistics (OA3102) Lab #2: Sampling, Sampling Distributions, and the Central Limit Theorem Goal: Use R to demonstrate sampling
More informationData Partitioning. Figure 1-31: Communication Topologies. Regular Partitions
Data In single-program multiple-data (SPMD) parallel programs, global data is partitioned, with a portion of the data assigned to each processing node. Issues relevant to choosing a partitioning strategy
More informationMicroscopic Measurement
Microscopic Measurement Estimating Specimen Size : The area of the slide that you see when you look through a microscope is called the " field of view ". If you know the diameter of your field of view,
More informationTo complete the computer assignments, you ll use the EViews software installed on the lab PCs in WMC 2502 and WMC 2506.
An Introduction to EViews The purpose of the computer assignments in BUEC 333 is to give you some experience using econometric software to analyse real-world data. Along the way, you ll become acquainted
More informationCREATING & RUNNING A VERY SIMPLE MODEL
CREATING & RUNNING A VERY SIMPLE MODEL Goal: This activity introduces a simple analog model of population change using dynamic systems modeling terminology and diagrams. It will serve as a launching point
More informationCUDA. Fluid simulation Lattice Boltzmann Models Cellular Automata
CUDA Fluid simulation Lattice Boltzmann Models Cellular Automata Please excuse my layout of slides for the remaining part of the talk! Fluid Simulation Navier Stokes equations for incompressible fluids
More informationUNIVERSITY OF WATERLOO Faculty of Mathematics
UNIVERSITY OF WATERLOO Faculty of Mathematics Exploring the application of Space Partitioning Methods on river segments S.S. Papadopulos & Associates Bethesda, MD, US Max Ren 20413992 3A Computer Science/BBA
More informationCollege Algebra Exam File - Fall Test #1
College Algebra Exam File - Fall 010 Test #1 1.) For each of the following graphs, indicate (/) whether it is the graph of a function and if so, whether it the graph of one-to one function. Circle your
More informationDraft Notes 1 : Scaling in Ad hoc Routing Protocols
Draft Notes 1 : Scaling in Ad hoc Routing Protocols Timothy X Brown University of Colorado April 2, 2008 2 Introduction What is the best network wireless network routing protocol? This question is a function
More informationGG450 4/5/2010. Today s material comes from p and in the text book. Please read and understand all of this material!
GG450 April 6, 2010 Seismic Reflection I Today s material comes from p. 32-33 and 81-116 in the text book. Please read and understand all of this material! Back to seismic waves Last week we talked about
More informationSpecific Objectives Students will understand that that the family of equation corresponds with the shape of the graph. Students will be able to create a graph of an equation by plotting points. In lesson
More informationRevision of the SolidWorks Variable Pressure Simulation Tutorial J.E. Akin, Rice University, Mechanical Engineering. Introduction
Revision of the SolidWorks Variable Pressure Simulation Tutorial J.E. Akin, Rice University, Mechanical Engineering Introduction A SolidWorks simulation tutorial is just intended to illustrate where to
More informationUnit II Graphing Functions and Data
Unit II Graphing Functions and Data These Materials were developed for use at and neither nor the author, Mark Schneider, assume any responsibility for their suitability or completeness for use elsewhere
More informationCMPSCI 187: Programming With Data Structures. Lecture 5: Analysis of Algorithms Overview 16 September 2011
CMPSCI 187: Programming With Data Structures Lecture 5: Analysis of Algorithms Overview 16 September 2011 Analysis of Algorithms Overview What is Analysis of Algorithms? L&C s Dishwashing Example Being
More informationSupervised vs.unsupervised Learning
Supervised vs.unsupervised Learning In supervised learning we train algorithms with predefined concepts and functions based on labeled data D = { ( x, y ) x X, y {yes,no}. In unsupervised learning we are
More information1 2 (3 + x 3) x 2 = 1 3 (3 + x 1 2x 3 ) 1. 3 ( 1 x 2) (3 + x(0) 3 ) = 1 2 (3 + 0) = 3. 2 (3 + x(0) 1 2x (0) ( ) = 1 ( 1 x(0) 2 ) = 1 3 ) = 1 3
6 Iterative Solvers Lab Objective: Many real-world problems of the form Ax = b have tens of thousands of parameters Solving such systems with Gaussian elimination or matrix factorizations could require
More informationOverview of Sentence Order Reference Document Development Process
Overview of Sentence Order Reference Document Development Process Scott Came Justice Integration Solutions, Inc. September 14, 2004 Purpose The purpose of this document is to outline the process/methodology
More informationDYNAMIC ANALYSIS OF A GENERATOR ON AN ELASTIC FOUNDATION
DYNAMIC ANALYSIS OF A GENERATOR ON AN ELASTIC FOUNDATION 7 DYNAMIC ANALYSIS OF A GENERATOR ON AN ELASTIC FOUNDATION In this tutorial the influence of a vibrating source on its surrounding soil is studied.
More informationAn introduction to plotting data
An introduction to plotting data Eric D. Black California Institute of Technology February 25, 2014 1 Introduction Plotting data is one of the essential skills every scientist must have. We use it on a
More informationAREA Judo Math Inc.
AREA 2013 Judo Math Inc. 6 th grade Problem Solving Discipline: Black Belt Training Order of Mastery: Area 1. Area of triangles by composition 2. Area of quadrilaterals by decomposing 3. Draw polygons
More informationHigh Performance Computing: Tools and Applications
High Performance Computing: Tools and Applications Edmond Chow School of Computational Science and Engineering Georgia Institute of Technology Lecture 15 Numerically solve a 2D boundary value problem Example:
More informationJoint Advanced Student School 2007 Martin Dummer
Sierpiński-Curves Joint Advanced Student School 2007 Martin Dummer Statement of the Problem What is the best way to store a triangle mesh efficiently in memory? The following points are desired : Easy
More informationSection Exponential Functions(Part I Growth)
Section 4.1 - Exponential Functions(Part I Growth) The number of cars in this city is growing exponentially every year. You may have heard quotes such as this. Let s take a look at some other exponential
More informationChapter 23. Linear Motion Motion of a Bug
Chapter 23 Linear Motion The simplest example of a parametrized curve arises when studying the motion of an object along a straight line in the plane We will start by studying this kind of motion when
More informationLesson 2: Analyzing a Data Set
Student Outcomes Students recognize linear, quadratic, and exponential functions when presented as a data set or sequence, and formulate a model based on the data. Lesson Notes This lesson asks students
More informationGIS Data Models. 4/9/ GIS Data Models
GIS Data Models 1 Conceptual models of the real world The real world can be described using two conceptually different models: 1. As discrete objects, possible to represent as points, lines or polygons.
More informationAP Calculus BC Summer Assignment
AP Calculus BC Summer Assignment Name Due Date: First Day of School Welcome to AP Calculus BC! This is an exciting, challenging, fast paced course that is taught at the college level. We have a lot of
More informationChapter 3: Rate Laws Excel Tutorial on Fitting logarithmic data
Chapter 3: Rate Laws Excel Tutorial on Fitting logarithmic data The following table shows the raw data which you need to fit to an appropriate equation k (s -1 ) T (K) 0.00043 312.5 0.00103 318.47 0.0018
More informationMath 2524: Activity 1 (Using Excel) Fall 2002
Math 2524: Activity 1 (Using Excel) Fall 22 Often in a problem situation you will be presented with discrete data rather than a function that gives you the resultant data. You will use Microsoft Excel
More informationComplex Dynamics in Life-like Rules Described with de Bruijn Diagrams: Complex and Chaotic Cellular Automata
Complex Dynamics in Life-like Rules Described with de Bruijn Diagrams: Complex and Chaotic Cellular Automata Paulina A. León Centro de Investigación y de Estudios Avanzados Instituto Politécnico Nacional
More informationUniform Motion Lab. The position equation for an object moving with a constant velocity is:
Uniform Motion Lab INTRODUCTION: In this experiment we will investigate motion without acceleration. Motion without acceleration is uniform (constant velocity) motion, which means it describes the motion
More informationMachine Learning / Jan 27, 2010
Revisiting Logistic Regression & Naïve Bayes Aarti Singh Machine Learning 10-701/15-781 Jan 27, 2010 Generative and Discriminative Classifiers Training classifiers involves learning a mapping f: X -> Y,
More informationCentrality Book. cohesion.
Cohesion The graph-theoretic terms discussed in the previous chapter have very specific and concrete meanings which are highly shared across the field of graph theory and other fields like social network
More informationProportional Relationships: Connections
Proportional Relationships: Connections Henri Picciotto A set of activities using multiple mathematical tools (and some real world applications) to help middle school students make connections between
More information1.2 Numerical Solutions of Flow Problems
1.2 Numerical Solutions of Flow Problems DIFFERENTIAL EQUATIONS OF MOTION FOR A SIMPLIFIED FLOW PROBLEM Continuity equation for incompressible flow: 0 Momentum (Navier-Stokes) equations for a Newtonian
More informationMidterm Exam 2B Answer key
Midterm Exam 2B Answer key 15110 Principles of Computing Fall 2015 April 6, 2015 Name: Andrew ID: Lab section: Instructions Answer each question neatly in the space provided. There are 6 questions totaling
More information7 Fractions. Number Sense and Numeration Measurement Geometry and Spatial Sense Patterning and Algebra Data Management and Probability
7 Fractions GRADE 7 FRACTIONS continue to develop proficiency by using fractions in mental strategies and in selecting and justifying use; develop proficiency in adding and subtracting simple fractions;
More informationHeat Kernels and Diffusion Processes
Heat Kernels and Diffusion Processes Definition: University of Alicante (Spain) Matrix Computing (subject 3168 Degree in Maths) 30 hours (theory)) + 15 hours (practical assignment) Contents 1. Solving
More informationCompressible Flow in a Nozzle
SPC 407 Supersonic & Hypersonic Fluid Dynamics Ansys Fluent Tutorial 1 Compressible Flow in a Nozzle Ahmed M Nagib Elmekawy, PhD, P.E. Problem Specification Consider air flowing at high-speed through a
More informationDirect Variations DIRECT AND INVERSE VARIATIONS 19. Name
DIRECT AND INVERSE VARIATIONS 19 Direct Variations Name Of the many relationships that two variables can have, one category is called a direct variation. Use the description and example of direct variation
More informationCMSC/BIOL 361: Emergence Cellular Automata: Introduction to NetLogo
Disclaimer: To get you oriented to the NetLogo platform, I ve put together an in-depth step-by-step walkthrough of a NetLogo simulation and the development environment in which it is presented. For those
More informationSimi imilar Shapes lar Shapes Nesting Squares Poly lyhedr hedra and E a and Euler ler s Form s Formula ula
TABLE OF CONTENTS Introduction......................................................... 5 Teacher s Notes....................................................... 6 NCTM Standards Alignment Chart......................................
More informationMatLab Project # 1 Due IN TUTORIAL Wednesday October 30
Mathematics 110 University of Victoria Fall 2013 MatLab Project # 1 Due IN TUTORIAL Wednesday October 30 Name ID V00 Section A0 Tutorial T0 Instructions: After completing this project, copy and paste your
More informationDigital Image Processing ERRATA. Wilhelm Burger Mark J. Burge. An algorithmic introduction using Java. Second Edition. Springer
Wilhelm Burger Mark J. Burge Digital Image Processing An algorithmic introduction using Java Second Edition ERRATA Springer Berlin Heidelberg NewYork Hong Kong London Milano Paris Tokyo 12.1 RGB Color
More informationVon Neumann Analysis for Higher Order Methods
1. Introduction Von Neumann Analysis for Higher Order Methods Von Neumann analysis is a widely used method to study how an initial wave is propagated with certain numerical schemes for a linear wave equation
More informationUsing Excel for Graphical Analysis of Data
Using Excel for Graphical Analysis of Data Introduction In several upcoming labs, a primary goal will be to determine the mathematical relationship between two variable physical parameters. Graphs are
More informationExercises Optimal binary search trees root
5.5 Optimal binary search trees 403 e w 5 5 j 4.75 i j 4.00 i 3.75.00 3 3 0.70 0.80 3.5.0. 4 0.55 0.50 0.60 4 0.90 0.70 0.60 0.90 5 0.45 0.35 0. 0.50 5 0 0.45 0.40 0.5 0. 0.50 6 0 0. 0.5 0.5 0.0 0.35 6
More informationModeling Evaporating Liquid Spray
Tutorial 16. Modeling Evaporating Liquid Spray Introduction In this tutorial, FLUENT s air-blast atomizer model is used to predict the behavior of an evaporating methanol spray. Initially, the air flow
More informationCurve fitting. Lab. Formulation. Truncation Error Round-off. Measurement. Good data. Not as good data. Least squares polynomials.
Formulating models We can use information from data to formulate mathematical models These models rely on assumptions about the data or data not collected Different assumptions will lead to different models.
More information6. Parallel Volume Rendering Algorithms
6. Parallel Volume Algorithms This chapter introduces a taxonomy of parallel volume rendering algorithms. In the thesis statement we claim that parallel algorithms may be described by "... how the tasks
More informationA Logistics Model Group Activity 8 STEM Project Week #11. Plot the data on the grid below. Be sure to label the x and y axis and label the window.
A Logistics Model Group Activity 8 STEM Project Week #11 Consider fencing off several thousand acres of land and placing 1000 rabbits on the land. Initially the rabbits would grow at a constant percent
More informationTutorial 1. Introduction to Using FLUENT: Fluid Flow and Heat Transfer in a Mixing Elbow
Tutorial 1. Introduction to Using FLUENT: Fluid Flow and Heat Transfer in a Mixing Elbow Introduction This tutorial illustrates the setup and solution of the two-dimensional turbulent fluid flow and heat
More informationHigh-Performance Computing
Informatik und Angewandte Kognitionswissenschaft Lehrstuhl für Hochleistungsrechnen Thomas Fogal Prof. Dr. Jens Krüger High-Performance Computing http://hpc.uni-due.de/teaching/wt2014/nbody.html Exercise
More informationA numerical microscope for plasma physics
A numerical microscope for plasma physics A new simulation capability developed for heavy-ion inertial fusion energy research will accelerate plasma physics and particle beam modeling, with application
More informationHomework # 1 Due: Feb 23. Multicore Programming: An Introduction
C O N D I T I O N S C O N D I T I O N S Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.86: Parallel Computing Spring 21, Agarwal Handout #5 Homework #
More informationUniversity of Florida CISE department Gator Engineering. Clustering Part 4
Clustering Part 4 Dr. Sanjay Ranka Professor Computer and Information Science and Engineering University of Florida, Gainesville DBSCAN DBSCAN is a density based clustering algorithm Density = number of
More informationIntroduction to Wavelets
Lab 11 Introduction to Wavelets Lab Objective: In the context of Fourier analysis, one seeks to represent a function as a sum of sinusoids. A drawback to this approach is that the Fourier transform only
More informationCSci 1113, Spring 2018 Lab Exercise 3 (Week 4): Repeat, Again and Again
CSci 1113, Spring 2018 Lab Exercise 3 (Week 4): Repeat, Again and Again Iteration Imperative programming languages such as C++ provide high-level constructs that support both conditional selection and
More informationCell based GIS. Introduction to rasters
Week 9 Cell based GIS Introduction to rasters topics of the week Spatial Problems Modeling Raster basics Application functions Analysis environment, the mask Application functions Spatial Analyst in ArcGIS
More informationReasoning and writing about algorithms: some tips
Reasoning and writing about algorithms: some tips Theory of Algorithms Winter 2016, U. Chicago Notes by A. Drucker The suggestions below address common issues arising in student homework submissions. Absorbing
More informationEngineering Effects of Boundary Conditions (Fixtures and Temperatures) J.E. Akin, Rice University, Mechanical Engineering
Engineering Effects of Boundary Conditions (Fixtures and Temperatures) J.E. Akin, Rice University, Mechanical Engineering Here SolidWorks stress simulation tutorials will be re-visited to show how they
More information