Physics Book - User contributions [en]

Vectors

2020-04-20T20:39:32Z

Zsanchez6: Fixed a couple of typos

== Claimed by: Zach Sanchez (Spring 2020) ==

This page defines and describes vectors.

==The Main Idea==

In mathematics and physics, a vector is a quantity with both a magnitude and a direction in space. The magnitude of a vector is a scalar value which can represent a variety of characteristics, depending on the situation. The magnitude of the vector (and the vector itself) has units corresponding to the characteristic it represents. The direction refers to the angle from the positive x -axis. Thus, the magnitude refers to the length and the direction refers to the angle. Note that a vector does not necessarily extend from one point in real, physical space to another; unless the magnitude is in units of length, the magnitude of the vector usually represents a property that exists at a single point in real, physical space, or at no position in particular.

===A Mathematical Model===

====Simple Examples of Vector Quantities====

To understand what it means for a vector to have both a magnitude and a direction, and to understand why the magnitude and direction together are often considered a single quantity, it can be helpful to consider an example. The [[Velocity]] of an object is an example of a vector quantity. The speed of the object, often given in meters per second (m/s), is a scalar value describing how quickly the object is moving. Speed is the magnitude of the velocity vector. However, the object's movement happens in a particular spatial direction, which the speed alone does not tell us. The direction of the object's movement is also a part of the velocity vector. Together, speed and direction comprise the velocity vector and give a complete description of an object's motion at a point in time. Another example of a vector is position; the distance of a point from the origin of a coordinate system can be represented as the magnitude of a vector, and this distance together with a direction describes exactly where a point can be found.

====Beginning Concepts====

A vector is typically written as a letter with an arrow over it. For example: <math>\vec{a} = \langle 1, 2, 3 \rangle</math>. The arrow may one-sided to make it easier to write. Alternatively, a vector might simply be written as a boldface letter: a. Which letter is used depends on the context; for example, <math>\vec{v}</math> represents velocity.

A specific component of a vector (see the section titled "forms") is denoted by a subscript: cx. For example, qy represents the y component of some vector <math>\vec{q}</math>.

Surrounding a vector by | symbols denotes its magnitude: <math>|\vec{b}|</math>. Alternatively, there are exists the notation <math>\overline{v}</math> and <math>\lVert \vec{v} \rVert_2</math>, with the latter denoting the Euclidean norm, which is the type of vector norm that you will see most often.

The magnitude of some vector is defined to be <math>\lVert \vec{v} \rVert_2 = \sqrt{v_x^2 + v_y^2 + v_z^2}</math>, where the axes of the 3D-map are mathematically notated as <math>\hat{x}, \hat{y}, \hat{z}</math>.

A unit vector is a vector whose magnitude is equal to one, and normalization is the process of setting a vector's magnitude to 1. A unit vector is denoted by a letter with a ^ symbol (called a "hat") written over it: <math>\hat{d}</math> (read as "d-hat"). Certain letters represent specific unit vectors. For example, <math>\hat{i}</math>, <math>\hat{j}</math>, and <math>\hat{k}</math> (alternatively, <math>\hat{x}</math>, <math>\hat{y}</math>, and <math>\hat{z}</math>) are unit vectors pointing in the +x, +y, and +z directions respectively. Non-Cartesian coordinate systems often have their own unit vectors; for example, 2D polar coordinates make use of the <math>\hat{r}</math> and <math>\hat{\theta}</math> unit vectors.

Unit vectors are defined to be for an arbitrary vector: <math>\hat{v} = \frac{\vec{v}}{\lVert \vec{v} \rVert_2}</math>

====Visually Representing Vectors====

Vectors are visually represented by arrows. The length of the arrow represents the magnitude of the vector, while the direction the arrow points in represents the direction of the vector. If a vector exists at a particular point in space, the "tail" of the arrow (the end without the V shape) should be placed at that point.

This example shows a visual representation of the velocity vector of a ball, which is moving to the right at a speed of 5m/s.

[[File:Vectorvisualrepresentation.png]]

====Important Vector Operations====

It is possible to perform a variety of mathematical operations on vectors, both with other vectors and with scalars. These operations appear in a variety of formulas in physics. To make the operations easier to learn, they are defined below assuming all vectors to be 3-dimensional; the more general n-dimensional definitions look more confusing. If necessary, it is easy to guess how to perform each operation with n-dimensional vectors by extrapolating from the 3-dimensional case.

Addition:

<math>\vec{a} + \vec{b} = \langle (a_x + b_x), (a_y + b_y), (a_z + b_z) \rangle</math>.

In other words, to add two vectors, simply add their like components to form the new components. Visually, if the tail of one vector is placed at the tip of another, their sum will extend from the tail of the second vector to the tip of the first:

[[File:Vectoraddition.png|600px]]

Subtraction:

<math>\vec{a} - \vec{b} = \langle (a_x - b_x), (a_y - b_y), (a_z - b_z)</math>.

In other words, to subtract two vectors, simply subtract their like components to form the new components.

Multiplication by scalar:

<math>k \vec{a} = \langle k \cdot a_x, k \cdot a_y, k \cdot a_z \rangle</math>.

In other words, multiplying a vector by a scalar multiplies each of that vector's components by that scalar. Note that this only affects a vector's magnitude, not its direction, unless the scalar is negative, in which case the direction of the vector is reversed.

[[File:Vectorscalarmultiplication.png|400px]]

Division by scalar:

Division by a scalar behaves exactly like scalar multiplication.

<math>\frac{\vec{a}}{k} = \langle \frac{a_x}{k}, \frac{a_y}{k}, \frac{a_z}{k} \rangle</math>.

Dot product (also called scalar product):

<math>\vec{a}\cdot\vec{b} =(a_x \cdot b_x) + (a_y \cdot b_y) + (a_z \cdot b_z)</math>.

In other words, the dot product of two vectors is the sum of the products of their like components. Note that this is a scalar value.

It is important to note that the dot product of two vectors has a specific value: <math>\vec{a}\cdot\vec{b} = |\vec{a}| \cdot |\vec{b}| \cdot \cos(\theta)</math>, where <math>\theta</math> is the angle between the vectors.

Cross product (also called vector product):

The simplest definition, which is also the one found on the formula sheet, is:

<math>\vec{a} \times \vec{b} = = \langle (a_yb_z - a_zb_y), (a_zb_x - a_xb_z), (a_xb_y - a_yb_x) \rangle = (a_yb_z - a_zb_y)\hat{i} + (a_zb_x - a_xb_z)\hat{j} + (a_xb_y - a_yb_x)\hat{k}</math>.

This is equivalent to the following matrix determinant, which may be easier to remember:

<math>\begin{vmatrix}
\hat{i} & \hat{j} & \hat{k} \\
a_x & a_y & a_z \\
b_x & b_y & b_z
\end{vmatrix}</math>

Note that this is a vector quantity. It is important to note that the magnitude of the cross product of two vectors has a specific value: <math>|\vec{a}\times\vec{b}| = |\vec{a}| \cdot |\vec{b}| \cdot \sin(\theta)</math>, where <math>\theta</math> is the angle between the vectors.

The direction of the cross product of two vectors is perpendicular to the plane in which those vectors lie and is given by the [[Right Hand Rule]]. 2D vectors do not have cross products. While the other operations listed here are commutative, associative (where the associative property is defined), and distributive over addition; cross product multiplication is not associative and is anticommutative (<math>\vec{a}\times\vec{b} = -\vec{b}\times\vec{a}</math>), meaning that if the order of the factors is reversed, their cross product will be reversed in direction.

For more mathematically-advanced students, I will invoke some higher-level linear algebra. Recall that the null space, defined as <math>Null(A) = \{ \vec{v} \in \mathbb{R}^n \mid A\vec{v} = \vec{0} \}</math> for an arbitrary matrix <math>A \in \mathbb{R}^{n \times n}</math>, is orthogonal, or normal, to its row space. That is, for all <math>\vec{v} \in Null(A)</math> and for all columns <math>a_{j:}</math>, we have <math>\sum^n_{i=1} a_{i,j} v_i = 0</math>.

For two vectors <math>\vec{a},\vec{b} \in \mathbb{R}^{3 \times 1}</math> (note that this is only valid for 3x1 vectors), if we think of <math>\vec{a}</math> as corresponding to a row in a matrix and <math>\vec{b}</math> to another, then the null space is equivalent to their cross-product. This is why the right-hand rule works. If we visualize the plane spanned by two vectors with our index and middle fingers, then their null space will be perpendicular to the intersection point of the two vectors, or rather the thumb. Likewise, if we use all of our fingers and our thumb to make an L-shape representing the span of the two vectors, upon curling our fingers, we will have the direction of the cross product, which is orthogonal to the two vectors represented by your thumb and fingers.

====Forms====

The information necessary to describe a specific vector can be presented in several forms.

Magnitude and direction form

In this form, the magnitude and the direction of the vector are explicitly stated. The statement describing direction might be a cardinal direction (ex. "north"), a direction on a graph (ex. "the +x direction"), or an angle (ex. "210<math>^\circ</math> from the x axis counterclockwise"), depending on the situation. Magnitude and direction form is often used in word problems because it is easy for humans to understand.

Component form

In this form, the vector is divided into components, each representing a different coordinate direction. In 2D space, these are the x and y directions. In 3D space, these are the x, y, and z directions. Each component tells how much the vector extends in that particular direction. Often, the three components are written enclosed by angle brackets and separated by commas. For example, the vector <2,0,-3> describes a vector that extends 2 units in the +x direction, 0 units in the y direction, and 3 units in the -z direction. Most vector operations described below can only be performed if the vectors are in component form, so this form may be necessary to do math for certain problems. Furthermore, programming languages store vectors in component form.

Unit vector form

In this form, the vector is expressed as a sum of unit vectors, each corresponding to a different coordinate direction. The symbols <math>\hat{i}</math>, <math>\hat{j}</math>, and <math>\hat{k}</math> OR the symbols <math>\hat{x}</math>, <math>\hat{y}</math>, and <math>\hat{z}</math> are used to represent unit vectors in the x, y, and z directions respectively. Consider the vector <2,0,-3>. It can be expressed in unit vector form as <math>2\hat{i} - 3\hat{k}</math>, meaning 2 times the x direction unit vector minus 3 times the z direction unit vector (see vector operations). While often considered its own form, unit vector form is very similar to component form, as the information describing the vector is stored in the same values. All references to component form in the rest of this page also apply to unit vector form.

[[File:Vectorsdifferentforms.png]]

Note that regardless of which form is used, an n-dimensional vector requires n values to mathematically describe. For example, consider a 3-dimensional vector. Describing this vector in component form requires an x value, a y value, and a z value. Describing it in magnitude and direction form requires one value to give the magnitude of the vector and two to give the its direction (the direction of a 3D vector could be described using, say, its polar angle <math>\theta</math> and its azimuthal angle <math>\phi</math>). A 1-dimensional vector (such as the velocity of a particle whose movement is constrained to the x axis) can be described using only 1 value whose sign indicates the vector's direction.

Converting between forms

It is possible to convert vectors from one form to another using simple trigonometry.

To find the magnitude of a vector in component form, use the Pythagorean theorem: add the squares of the components and take the square root of the result. For a 2D vector, <math>|\vec{a}| = \sqrt{a_x^2 + a_y^2}</math>.

To find the direction of a vector in component form, use inverse trigonometric functions. For a 2D vector, <math>\theta = \tan^{-1}(\frac{a_y}{a_x})</math>, where <math>\theta</math> is the angle vector <math>\vec{a}</math> makes with the x axis in the counterclockwise direction.

To find the components of a vector in magnitude and direction form, use trigonometric functions. For a 2D vector, <math>a_x = |\vec{a}|\cos\theta</math> and <math>a_y = |\vec{a}|\sin\theta</math>, where <math>\theta</math> is the angle vector <math>\vec{a}</math> makes with the x axis in the counterclockwise direction.

===A Computational Model===

In VPython, vector objects are in component form; each one has an x, y, and z component. Recall that in VPython, using the default camera orientation, the +x axis points to the right, the +y axis points upwards, and the +z axis points out of the plane of the screen towards the viewer. The constructor for a vector object is the word "vec" or "vector," and it takes three arguments, which define its x, y and z components respectively. A line to create a vector called "velocity" might look like this:

velocity = vec(3,-1,2)

To access or modify a specific component of a vector object, its name should be followed by a period and an x, y, or z. For example, to change the x component of the above velocity vector from 3 to 5, the following line might be used:

velocity.x = 5

In VPython, vectors have many uses. The position of each object is defined as a vector; the position vector's tail lies at the origin and its head lies at the center of the object in question. Furthermore, the dimensions of a rectangular prism (a "box" object) are defined as a vector; the x component determines its width, the y component its height, and the z component its thickness.

To graphically represent a vector such as an electric field, "arrow" objects should be used. In addition to taking a position vector (which determines the position of the arrow's tail), arrow objects take an "axis" vector, which determines their size and shape. To represent a vector, simply make it the axis of an arrow object.

==Examples==

===Simple===

Vector <math>\vec{a}</math> is <2,4,2>. Vector <math>\vec{b}</math> is <-1,1,3>. What is the magnitude of the vector <math>\vec{a} - 2\vec{b}</math>?

Solution:

<math>\vec{a} - 2\vec{b} = <2,4,2> - 2 * <-1,1,3> </math>

<math> = <2,4,2> - <-2,2,6> </math>

<math> = <4,2,-4></math>

We are asked to find the magnitude of this vector, so let us use the Pythagorean theorem with its components:

<math>|<4,2,-4>| = \sqrt{4^2 + 2^2 + (-4)^2}</math>

<math> = \sqrt{16 + 4 + 16} </math>

<math> = \sqrt{36} </math>

<math> = 6 </math>

===Intermediate===

An airplane is travelling in still air at 240m/s in the direction 35<math>^\circ</math> south of west. A wind begins to blow; the wind has a speed of 80m/s in the direction 15<math>^\circ</math> east of north. What should be the new velocity of the plane relative to the air around it to maintain its original trajectory? You may give your answer in component form (+x is east, +y is north).

[[File:Vectorsplaneproblem.png]]

Solution:

The vector sum of the new velocity of the plane <math>\vec{v_{p,1}}</math> and the velocity of the wind <math>\vec{v_w}</math> should equal the original velocity of the plane <math>\vec{v_{p,0}}</math> (see [[Relative Velocity]]):

[[File:Vectorsplanesolution.png]]

<math> \vec{v_{p,1}} + \vec{v_w} = \vec{v_{p,0}} </math>

<math> \vec{v_{p,1}} = \vec{v_{p,0}} - \vec{v_w} </math>

Let us convert the given vectors to component form for easier subtraction. The +x and +y directions will be east and north respectively.

<math>\vec{v_{p,0}} = <240\cos(215^\circ), 240\sin(215^\circ)></math>m/s (35<math>^\circ</math> south of west is 215<math>^\circ</math> above the x axis)

<math>\vec{v_{p,0}} = <-196.6, -137.7> </math>m/s

<math>\vec{v_w} = <80\cos(75^\circ), 80\sin(75^\circ)></math>m/s (15<math>^\circ</math> east of north is 75<math>^\circ</math> above the x axis)

<math>\vec{v_w} = <20.7, 77.3> </math>m/s

Now let us subtract:

<math> \vec{v_{p,1}} = \vec{v_{p,0}} - \vec{v_w} </math>

<math> \vec{v_{p,1}} = <-196.6, -137.7> - <20.7, 77.3> </math>m/s

<math> \vec{v_{p,1}} = <-217.3, -214.9> </math>m/s

===Difficult===

What is the angle between the vectors <2,5,-2> and <3,-4,-1>?

Solution:

The dot product between two vectors is equal to the product of their magnitudes times the cosine of the angle between them. Let us use this property to find the angle between the given vectors.

<math> <2,5,-2> \cdot <3,-4,-1> = |<2,5,-2>| * |<3,-4,-1>| * \cos\theta </math>

Rearranging this yields

<math> \theta = \cos^{-1}\frac{<2,5,-2> \cdot <3,-4,-1>}{|<2,5,-2>| * |<3,-4,-1>|} </math>

Let us evaluate the dot product and simplify:

<math> \theta = \cos^{-1}\frac{2(3) + 5(-4) + (-2)(-1)}{\sqrt{2^2 + 5^2 + (-2)^2} * \sqrt{3^2 + (-4)^2 + (-1)^2}} </math>

<math> \theta = \cos^{-1}\frac{-12}{\sqrt{33 * 26}} </math>

<math> \theta = 114^\circ </math>

==Connectedness==

Vectors are used in many fields and levels of physics. Kinematics, for example, studies the relationships between the position vector and its time derivatives (velocity, acceleration). Force is also a vector quantity, as are certain system properties such as linear and angular momentum. Vector fields are commonly used in electricity and magnetism, as well as in fluid dynamics.

==History==

It is unknown who first developed the idea of vectors, but the oldest known reference to vectors is in the work Mechanics by Hero of Alexandria (first century AD), which described their addition. At this point, however, the idea of a vector was little more than a line segment with a specific orientation; they had a length extending from one point in physical space to another but were not used to represent anything else.

In the early 19th century, several mathematicians and physicists (including Caspar Wessel (1745-1818), Jean Robert Argand (1768-1822), Carl Friedrich Gauss (1777-1855)), and William Rowan Hamilton (1805-1865) used 2D vectors to represent complex numbers; one component would represent the real value and another would represent the imaginary value. Hamilton would also become the first to use the word "vector." August Ferdinand Möbius (1790-1868) contributed to vector math in his 1827 book The Barycentric Calculus, in which he developed the convention of labeling vectors with letters and defined the multiplication of a vector by a scalar. Hermann Grassmann (1809-1877) wrote in his 1844 book Ausdehnungslehre (German for "The Calculus of Extension") that vectors could exist in space of any number of dimensions and described much of what would become linear algebra, which makes ample use of vectors.

The modern language and conventions surrounding vectors come largely from notes created by J. Willard Gibbs (1839--1903), a professor at Yale University.

==See also==

===External links===
Mathematical Computations on Vectors: [http://ocw.mit.edu/courses/mathematics/18-02sc-multivariable-calculus-fall-2010/1.-vectors-and-matrices/part-a-vectors-determinants-and-planes/session-1-vectors/MIT18_02SC_notes_0.pdf [http://ocw.mit.edu/courses/mathematics/18-02sc-multivariable-calculus-fall-2010/1.-vectors-and-matrices/part-a-vectors-determinants-and-planes/session-1-vectors/MIT18_02SC_notes_0.pdf]

Computational Work with Vectors: [http://vpython.org/contents/docs/vector.html http://vpython.org/contents/docs/vector.html]

Basics of Vectors: [https://www.physics.uoguelph.ca/tutorials/vectors/vectors.html https://www.physics.uoguelph.ca/tutorials/vectors/vectors.html]

===Further Reading===

Vector Analysis by Josiah Willard Gibbs

Introduction to Matrices and Vectors by Jacob T. Schwartz

==References==
[https://www.mathsisfun.com/algebra/vectors.html https://www.mathsisfun.com/algebra/vectors.html]

[http://ocw.mit.edu/courses/mathematics/18-02sc-multivariable-calculus-fall-2010/1.-vectors-and-matrices/part-a-vectors-determinants-and-planes/session-1-vectors/MIT18_02SC_notes_0.pdf http://ocw.mit.edu/courses/mathematics/18-02sc-multivariable-calculus-fall-2010/1.-vectors-and-matrices/part-a-vectors-determinants-and-planes/session-1-vectors/MIT18_02SC_notes_0.pdf ]

[http://mathinsight.org/vector_introduction http://mathinsight.org/vector_introduction]

[http://www.math.mcgill.ca/labute/courses/133f03/VectorHistory.html http://www.math.mcgill.ca/labute/courses/133f03/VectorHistory.html]

[[Category: Geometry]]

Vectors

2020-04-20T20:34:33Z

Zsanchez6: Cleaned up the ordering to make it more logical, added some definitions to help clarify some concepts, added in more advanced math for the cross product to explain the right-hand rule

== Claimed by: Zach Sanchez (Spring 2020) ==

This page defines and describes vectors.

==The Main Idea==

In mathematics and physics, a vector is a quantity with both a magnitude and a direction in space. The magnitude of a vector is a scalar value which can represent a variety of characteristics, depending on the situation. The magnitude of the vector (and the vector itself) has units corresponding to the characteristic it represents. The direction refers to the angle from the positive x -axis. Thus, the magnitude refers to the length and the direction refers to the angle. Note that a vector does not necessarily extend from one point in real, physical space to another; unless the magnitude is in units of length, the magnitude of the vector usually represents a property that exists at a single point in real, physical space, or at no position in particular.

===A Mathematical Model===

====Simple Examples of Vector Quantities====

To understand what it means for a vector to have both a magnitude and a direction, and to understand why the magnitude and direction together are often considered a single quantity, it can be helpful to consider an example. The [[Velocity]] of an object is an example of a vector quantity. The speed of the object, often given in meters per second (m/s), is a scalar value describing how quickly the object is moving. Speed is the magnitude of the velocity vector. However, the object's movement happens in a particular spatial direction, which the speed alone does not tell us. The direction of the object's movement is also a part of the velocity vector. Together, speed and direction comprise the velocity vector and give a complete description of an object's motion at a point in time. Another example of a vector is position; the distance of a point from the origin of a coordinate system can be represented as the magnitude of a vector, and this distance together with a direction describes exactly where a point can be found.

====Unit Vectors====

A unit vector is a vector whose magnitude is 1 by design.

Some unit vectors are normalized regular vectors. For example, <math>\hat{r}</math> is the normalized version of some vector <math>\vec{r}</math>. Normalization is the process of setting a vector's magnitude to 1. This is done to discard information about the vector's magnitude while retaining information about its direction. To normalize a vector, the following formula can be used: <math>\hat{a} = \frac{\vec{a}}{|\vec{a}|}</math>

Certain letters represent specific unit vectors. For example, <math>\hat{i}</math>, <math>\hat{j}</math>, and <math>\hat{k}</math> (alternatively, <math>\hat{x}</math>, <math>\hat{y}</math>, and <math>\hat{z}</math>) are unit vectors pointing in the +x, +y, and +z directions respectively. Non-Cartesian coordinate systems often have their own unit vectors; for example, 2D polar coordinates make use of the <math>\hat{r}</math> and <math>\hat{\theta}</math> unit vectors.

====Notation====
A vector is typically written as a letter with an arrow over it: <math>\vec{a}</math>. The arrow may one-sided to make it easier to write. Alternatively, a vector might simply be written as a boldface letter: a. Which letter is used depends on the context; for example, <math>\vec{v}</math> represents velocity.

Surrounding a vector by | symbols denotes its magnitude: <math>|\vec{b}|</math>. Alternatively, there are exists the notation <math>\overline{v}</math> and <math>\lVert \vec{v} \rVert_2</math>, with the latter denoting the Euclidean norm, which is the type of vector norm that you will see most often.

The magnitude of some vector is defined to be <math>\lVert \vec{v} \rVert_2 = \sqrt{v_x^2 + v_y^2 + v_z^2}</math>, where the axes of the 3D-map are mathematically notated as <math>\hat{x}, \hat{y}, \hat{z}</math>.

A specific component of a vector (see the section titled "forms") is denoted by a subscript: cx. For example, qy represents the y component of some vector <math>\vec{q}</math>.

A unit vector is denoted by a letter with a ^ symbol (called a "hat") written over it: <math>\hat{d}</math> (read as "d-hat").

Unit vectors are defined to be for some vector: <math>\hat{v} = \frac{\vec{v}}{\lVert \vec{v} \rVert_2}

====Visually Representing Vectors====

Vectors are visually represented by arrows. The length of the arrow represents the magnitude of the vector, while the direction the arrow points in represents the direction of the vector. If a vector exists at a particular point in space, the "tail" of the arrow (the end without the V shape) should be placed at that point.

This example shows a visual representation of the velocity vector of a ball, which is moving to the right at a speed of 5m/s.

[[File:Vectorvisualrepresentation.png]]

====Important Vector Operations====

It is possible to perform a variety of mathematical operations on vectors, both with other vectors and with scalars. These operations appear in a variety of formulas in physics. To make the operations easier to learn, they are defined below assuming all vectors to be 3-dimensional; the more general n-dimensional definitions look more confusing. If necessary, it is easy to guess how to perform each operation with n-dimensional vectors by extrapolating from the 3-dimensional case.

Addition:

<math>\vec{a} + \vec{b} = \langle (a_x + b_x), (a_y + b_y), (a_z + b_z) \rangle</math>.

In other words, to add two vectors, simply add their like components to form the new components. Visually, if the tail of one vector is placed at the tip of another, their sum will extend from the tail of the second vector to the tip of the first:

[[File:Vectoraddition.png|600px]]

Subtraction:

<math>\vec{a} - \vec{b} = \langle (a_x - b_x), (a_y - b_y), (a_z - b_z)</math>.

In other words, to subtract two vectors, simply subtract their like components to form the new components.

Multiplication by scalar:

<math>k \vec{a} = \langle k \cdot a_x, k \cdot a_y, k \cdot a_z \rangle</math>.

In other words, multiplying a vector by a scalar multiplies each of that vector's components by that scalar. Note that this only affects a vector's magnitude, not its direction, unless the scalar is negative, in which case the direction of the vector is reversed.

[[File:Vectorscalarmultiplication.png|400px]]

Division by scalar:

Division by a scalar behaves exactly like scalar multiplication.

<math>\frac{\vec{a}}{k} = \langle \frac{a_x}{k}, \frac{a_y}{k}, \frac{a_z} \rangle</math>.

Dot product (also called scalar product):

<math>\vec{a}\cdot\vec{b} =(a_x \cdot b_x) + (a_y \cdot b_y) + (a_z \cdot b_z)</math>.

In other words, the dot product of two vectors is the sum of the products of their like components. Note that this is a scalar value.

It is important to note that the dot product of two vectors has a specific value: <math>\vec{a}\cdot\vec{b} = |\vec{a}| \cdot |\vec{b}| \cdot \cos(\theta)</math>, where <math>\theta</math> is the angle between the vectors.

Cross product (also called vector product):

The simplest definition, which is also the one found on the formula sheet, is:

<math>\vec{a} \times \vec{b} = = \langle (a_yb_z - a_zb_y), (a_zb_x - a_xb_z), (a_xb_y - a_yb_x) \rangle = (a_yb_z - a_zb_y)\hat{i} + (a_zb_x - a_xb_z)\hat{j} + (a_xb_y - a_yb_x)\hat{k}</math>.

This is equivalent to the following matrix determinant, which may be easier to remember:

<math>\begin{vmatrix}
\hat{i} & \hat{j} & \hat{k} \\
a_x & a_y & a_z \\
b_x & b_y & b_z
\end{vmatrix}</math>

Note that this is a vector quantity. It is important to note that the magnitude of the cross product of two vectors has a specific value: <math>|\vec{a}\times\vec{b}| = |\vec{a}| \cdot |\vec{b}| \cdot \sin(\theta)</math>, where <math>\theta</math> is the angle between the vectors.

The direction of the cross product of two vectors is perpendicular to the plane in which those vectors lie and is given by the [[Right Hand Rule]]. 2D vectors do not have cross products. While the other operations listed here are commutative, associative (where the associative property is defined), and distributive over addition; cross product multiplication is not associative and is anticommutative (<math>\vec{a}\times\vec{b} = -\vec{b}\times\vec{a}</math>), meaning that if the order of the factors is reversed, their cross product will be reversed in direction.

For more mathematically-advanced students, I will invoke some higher-level linear algebra. Recall that the null space, defined as <math>Null(A) = \{ \vec{v} \in \mathbb{R}^n \mid A\vec{v} = \vec{0} \}</math> for an arbitrary matrix <math>A \in \mathbb{R}^{n \times n}</math>, is orthogonal, or normal, to its row space. That is, for all <math>\vec{v} \in Null(A)</math> and for all columns <math>a_{j:}</math>, we have <math>\sum^n_{i=1} a_{i,j} v_i = 0</math>.

For two vectors <math>\vec{a},\vec{b} \in \mathbb{R}^{3 \times 1}</math> (note that this is only valid for 3x1 vectors), if we think of <math>\vec{a}</math> as corresponding to a row in a matrix and <math>\vec{b}</math> to another, then the null space is equivalent to their cross-product. This is why the right-hand rule works. If we visualize the plane spanned by two vectors with our index and middle fingers, then their null space will be perpendicular to the intersection point of the two vectors, or rather the thumb. Likewise, if we use all of our fingers and our thumb to make an L-shape representing the span of the two vectors, upon curling our fingers, we will have the direction of the cross product, which is orthogonal to the two vectors represented by your thumb and fingers.

====Forms====

The information necessary to describe a specific vector can be presented in several forms.

Magnitude and direction form

In this form, the magnitude and the direction of the vector are explicitly stated. The statement describing direction might be a cardinal direction (ex. "north"), a direction on a graph (ex. "the +x direction"), or an angle (ex. "210<math>^\circ</math> from the x axis counterclockwise"), depending on the situation. Magnitude and direction form is often used in word problems because it is easy for humans to understand.

Component form

In this form, the vector is divided into components, each representing a different coordinate direction. In 2D space, these are the x and y directions. In 3D space, these are the x, y, and z directions. Each component tells how much the vector extends in that particular direction. Often, the three components are written enclosed by angle brackets and separated by commas. For example, the vector <2,0,-3> describes a vector that extends 2 units in the +x direction, 0 units in the y direction, and 3 units in the -z direction. Most vector operations described below can only be performed if the vectors are in component form, so this form may be necessary to do math for certain problems. Furthermore, programming languages store vectors in component form.

Unit vector form

In this form, the vector is expressed as a sum of unit vectors, each corresponding to a different coordinate direction. The symbols <math>\hat{i}</math>, <math>\hat{j}</math>, and <math>\hat{k}</math> OR the symbols <math>\hat{x}</math>, <math>\hat{y}</math>, and <math>\hat{z}</math> are used to represent unit vectors in the x, y, and z directions respectively. Consider the vector <2,0,-3>. It can be expressed in unit vector form as <math>2\hat{i} - 3\hat{k}</math>, meaning 2 times the x direction unit vector minus 3 times the z direction unit vector (see vector operations). While often considered its own form, unit vector form is very similar to component form, as the information describing the vector is stored in the same values. All references to component form in the rest of this page also apply to unit vector form.

[[File:Vectorsdifferentforms.png]]

Note that regardless of which form is used, an n-dimensional vector requires n values to mathematically describe. For example, consider a 3-dimensional vector. Describing this vector in component form requires an x value, a y value, and a z value. Describing it in magnitude and direction form requires one value to give the magnitude of the vector and two to give the its direction (the direction of a 3D vector could be described using, say, its polar angle <math>\theta</math> and its azimuthal angle <math>\phi</math>). A 1-dimensional vector (such as the velocity of a particle whose movement is constrained to the x axis) can be described using only 1 value whose sign indicates the vector's direction.

Converting between forms

It is possible to convert vectors from one form to another using simple trigonometry.

To find the magnitude of a vector in component form, use the Pythagorean theorem: add the squares of the components and take the square root of the result. For a 2D vector, <math>|\vec{a}| = \sqrt{a_x^2 + a_y^2}</math>.

To find the direction of a vector in component form, use inverse trigonometric functions. For a 2D vector, <math>\theta = \tan^{-1}(\frac{a_y}{a_x})</math>, where <math>\theta</math> is the angle vector <math>\vec{a}</math> makes with the x axis in the counterclockwise direction.

To find the components of a vector in magnitude and direction form, use trigonometric functions. For a 2D vector, <math>a_x = |\vec{a}|\cos\theta</math> and <math>a_y = |\vec{a}|\sin\theta</math>, where <math>\theta</math> is the angle vector <math>\vec{a}</math> makes with the x axis in the counterclockwise direction.

===A Computational Model===

In VPython, vector objects are in component form; each one has an x, y, and z component. Recall that in VPython, using the default camera orientation, the +x axis points to the right, the +y axis points upwards, and the +z axis points out of the plane of the screen towards the viewer. The constructor for a vector object is the word "vec" or "vector," and it takes three arguments, which define its x, y and z components respectively. A line to create a vector called "velocity" might look like this:

velocity = vec(3,-1,2)

To access or modify a specific component of a vector object, its name should be followed by a period and an x, y, or z. For example, to change the x component of the above velocity vector from 3 to 5, the following line might be used:

velocity.x = 5

In VPython, vectors have many uses. The position of each object is defined as a vector; the position vector's tail lies at the origin and its head lies at the center of the object in question. Furthermore, the dimensions of a rectangular prism (a "box" object) are defined as a vector; the x component determines its width, the y component its height, and the z component its thickness.

To graphically represent a vector such as an electric field, "arrow" objects should be used. In addition to taking a position vector (which determines the position of the arrow's tail), arrow objects take an "axis" vector, which determines their size and shape. To represent a vector, simply make it the axis of an arrow object.

==Examples==

===Simple===

Vector <math>\vec{a}</math> is <2,4,2>. Vector <math>\vec{b}</math> is <-1,1,3>. What is the magnitude of the vector <math>\vec{a} - 2\vec{b}</math>?

Solution:

<math>\vec{a} - 2\vec{b} = <2,4,2> - 2 * <-1,1,3> </math>

<math> = <2,4,2> - <-2,2,6> </math>

<math> = <4,2,-4></math>

We are asked to find the magnitude of this vector, so let us use the Pythagorean theorem with its components:

<math>|<4,2,-4>| = \sqrt{4^2 + 2^2 + (-4)^2}</math>

<math> = \sqrt{16 + 4 + 16} </math>

<math> = \sqrt{36} </math>

<math> = 6 </math>

===Intermediate===

An airplane is travelling in still air at 240m/s in the direction 35<math>^\circ</math> south of west. A wind begins to blow; the wind has a speed of 80m/s in the direction 15<math>^\circ</math> east of north. What should be the new velocity of the plane relative to the air around it to maintain its original trajectory? You may give your answer in component form (+x is east, +y is north).

[[File:Vectorsplaneproblem.png]]

Solution:

The vector sum of the new velocity of the plane <math>\vec{v_{p,1}}</math> and the velocity of the wind <math>\vec{v_w}</math> should equal the original velocity of the plane <math>\vec{v_{p,0}}</math> (see [[Relative Velocity]]):

[[File:Vectorsplanesolution.png]]

<math> \vec{v_{p,1}} + \vec{v_w} = \vec{v_{p,0}} </math>

<math> \vec{v_{p,1}} = \vec{v_{p,0}} - \vec{v_w} </math>

Let us convert the given vectors to component form for easier subtraction. The +x and +y directions will be east and north respectively.

<math>\vec{v_{p,0}} = <240\cos(215^\circ), 240\sin(215^\circ)></math>m/s (35<math>^\circ</math> south of west is 215<math>^\circ</math> above the x axis)

<math>\vec{v_{p,0}} = <-196.6, -137.7> </math>m/s

<math>\vec{v_w} = <80\cos(75^\circ), 80\sin(75^\circ)></math>m/s (15<math>^\circ</math> east of north is 75<math>^\circ</math> above the x axis)

<math>\vec{v_w} = <20.7, 77.3> </math>m/s

Now let us subtract:

<math> \vec{v_{p,1}} = \vec{v_{p,0}} - \vec{v_w} </math>

<math> \vec{v_{p,1}} = <-196.6, -137.7> - <20.7, 77.3> </math>m/s

<math> \vec{v_{p,1}} = <-217.3, -214.9> </math>m/s

===Difficult===

What is the angle between the vectors <2,5,-2> and <3,-4,-1>?

Solution:

The dot product between two vectors is equal to the product of their magnitudes times the cosine of the angle between them. Let us use this property to find the angle between the given vectors.

<math> <2,5,-2> \cdot <3,-4,-1> = |<2,5,-2>| * |<3,-4,-1>| * \cos\theta </math>

Rearranging this yields

<math> \theta = \cos^{-1}\frac{<2,5,-2> \cdot <3,-4,-1>}{|<2,5,-2>| * |<3,-4,-1>|} </math>

Let us evaluate the dot product and simplify:

<math> \theta = \cos^{-1}\frac{2(3) + 5(-4) + (-2)(-1)}{\sqrt{2^2 + 5^2 + (-2)^2} * \sqrt{3^2 + (-4)^2 + (-1)^2}} </math>

<math> \theta = \cos^{-1}\frac{-12}{\sqrt{33 * 26}} </math>

<math> \theta = 114^\circ </math>

==Connectedness==

Vectors are used in many fields and levels of physics. Kinematics, for example, studies the relationships between the position vector and its time derivatives (velocity, acceleration). Force is also a vector quantity, as are certain system properties such as linear and angular momentum. Vector fields are commonly used in electricity and magnetism, as well as in fluid dynamics.

==History==

It is unknown who first developed the idea of vectors, but the oldest known reference to vectors is in the work Mechanics by Hero of Alexandria (first century AD), which described their addition. At this point, however, the idea of a vector was little more than a line segment with a specific orientation; they had a length extending from one point in physical space to another but were not used to represent anything else.

In the early 19th century, several mathematicians and physicists (including Caspar Wessel (1745-1818), Jean Robert Argand (1768-1822), Carl Friedrich Gauss (1777-1855)), and William Rowan Hamilton (1805-1865) used 2D vectors to represent complex numbers; one component would represent the real value and another would represent the imaginary value. Hamilton would also become the first to use the word "vector." August Ferdinand Möbius (1790-1868) contributed to vector math in his 1827 book The Barycentric Calculus, in which he developed the convention of labeling vectors with letters and defined the multiplication of a vector by a scalar. Hermann Grassmann (1809-1877) wrote in his 1844 book Ausdehnungslehre (German for "The Calculus of Extension") that vectors could exist in space of any number of dimensions and described much of what would become linear algebra, which makes ample use of vectors.

The modern language and conventions surrounding vectors come largely from notes created by J. Willard Gibbs (1839--1903), a professor at Yale University.

==See also==

===External links===
Mathematical Computations on Vectors: [http://ocw.mit.edu/courses/mathematics/18-02sc-multivariable-calculus-fall-2010/1.-vectors-and-matrices/part-a-vectors-determinants-and-planes/session-1-vectors/MIT18_02SC_notes_0.pdf [http://ocw.mit.edu/courses/mathematics/18-02sc-multivariable-calculus-fall-2010/1.-vectors-and-matrices/part-a-vectors-determinants-and-planes/session-1-vectors/MIT18_02SC_notes_0.pdf]

Computational Work with Vectors: [http://vpython.org/contents/docs/vector.html http://vpython.org/contents/docs/vector.html]

Basics of Vectors: [https://www.physics.uoguelph.ca/tutorials/vectors/vectors.html https://www.physics.uoguelph.ca/tutorials/vectors/vectors.html]

===Further Reading===

Vector Analysis by Josiah Willard Gibbs

Introduction to Matrices and Vectors by Jacob T. Schwartz

==References==
[https://www.mathsisfun.com/algebra/vectors.html https://www.mathsisfun.com/algebra/vectors.html]

[http://ocw.mit.edu/courses/mathematics/18-02sc-multivariable-calculus-fall-2010/1.-vectors-and-matrices/part-a-vectors-determinants-and-planes/session-1-vectors/MIT18_02SC_notes_0.pdf http://ocw.mit.edu/courses/mathematics/18-02sc-multivariable-calculus-fall-2010/1.-vectors-and-matrices/part-a-vectors-determinants-and-planes/session-1-vectors/MIT18_02SC_notes_0.pdf ]

[http://mathinsight.org/vector_introduction http://mathinsight.org/vector_introduction]

[http://www.math.mcgill.ca/labute/courses/133f03/VectorHistory.html http://www.math.mcgill.ca/labute/courses/133f03/VectorHistory.html]

[[Category: Geometry]]

Python Syntax

2017-11-29T18:18:45Z

Zsanchez6:

'''Edited by Zach Sanchez - Fall 2017'''

==The Main Idea==
VPython is an extension of the Python programming language that contains a 3D graphics module called Visual, which allows users to create simple simulations of physics problems. Its primary use is for educational purposes, although it has been used in research before in the past. We will use VPython on occasion to aid us in visualizing and experimenting with the ideas that we learn in this course.

==Downloading vPython==
In all cases you will need to have Python 2.7.9 installed. Future versions may also work, but be warned that it may cause VPython to behave erratically.

[http://vpython.org/contents/download_windows.html For Windows]

[http://vpython.org/contents/download_mac.html For Mac]

[http://vpython.org/contents/download_linux.html For Linux]

Note for Linux: You can only install Classic VPython, which is no longer supported and may be missing some features and bug fixes that are present in the current version of VPython.

Alternatively, you can use GlowScript, a virtual environment that allows you to execute VPython code in your browser. You can create an account at [http://www.glowscript.org/ http://www.glowscript.org/].

==Python Basics==
It's always helpful to have a background knowledge of Python, however basic, before proceeding with the applications more relevant to this course.

# Comments start with a number symbol. Use these liberally in order to keep track of the changes you make/the code you write.
# We're always going to start our programs with the following headers to ensure that math works correctly.
from__future__ import division
from visual import*

# Math works as you would expect.
1 + 1 # -> 2
2 - 1 # -> 1
3 * 2 # -> 6
8 / 4 # -> 4
(1 + 3) * 2 # -> 8
1 + (3 * 2) # -> 7

# That is, until you need to do an exponent. But you shouldn't encounter anything much weirder than that math-wise.
2 ** 3 # -> 8

# There's also various math functions you can use, such as for finding the square root.
sqrt(4) # -> 2

# However, these math operations are useless if we don't store the values anywhere.
x = 1 + 1
# Now if you try to access x in the future, it will have the value of '2'.

# You can print your variables like this:
print(x)
# You can also add in some text with them like this:
print("The value of x is ", x)

# Sometimes we only want to execute something based on a certain condition.
if x == 2:
y = 3
# The '==' means that you're directly comparing two values. If they're equal, whatever is on the inside of the if-statement will be executed.
# Now is a good time to note that whitespace is very important in Python.
# 'y = 3' will always execute if it isn't tabbed, so even if x is equal to 7, y will still become 3.
# If you ever see me use whitespace in this tutorial, make sure you do as well.
# Some other useful comparators: != (not equal), >, >= (greater than or equal to), <, <= (less than or equal to), 'and', 'or', 'not'
# You can add onto the above if-statement by tacking on:
elif x == 3:
y = 4
else:
y = 5
# With elif meaning "else if" for a second possible condition, and "else" for something that you want to execute if none of the earlier conditions are met.

# You can also have a loop that executes code for a certain amount of time or based on a certain condition.
while x < 5:
print("x is ", x)
x = x + 1
# This will output the value of x up until it reaches 4.

Thus concludes our very, very brief foray into Python. That was simply the bare minimum required to understand what will be needed for this class. More concepts will be introduced as needed, whether that be in this tutorial or in a lab. For a more complete understanding of Python, you can look up other Python tutorials elsewhere on the Internet.

==VPython Basics==

In this section we will analyze some of the objects utilized by VPython. Each object is more or less self-explanatory with equally clear fields within.

# The scene has a width and a height; you can set them like so.
scene.width = 1024
scene.height = 760

# Beyond that, we have vectors.
# You can access each component with pos.x, pos.y, or pos.z depending on which component you want.
# Now would be a good time to mention that you access an object's field with a period following the variable name.
pos = vector(0, 0, 0)
xComponent = pos.x
yComponent = pos.y
zComponent = pos.z

# You can add two or more vectors together and you can multiply a vector by a scalar,
# but you can't add a scalar to a vector.
# To get the magnitude of or to normalize a vector, use the appropriate function.
initialPos = vector(10, 3, 0)
finalPos = vector(30, 15, 0)
deltaR = finalPos - initialPos # -> vector(20, 12, 0)
rMag = mag(finalPos)
rHat = norm(finalPos)
errorDemo = magR + finalPos # -> error; causes your program to crash

# If you were to make, say, a ball, you would want to draw a sphere.
# It has a position field, a radius field, and a color field.
# When creating a new instance of an object, each field is comma separated.
# You access each field the same way as we did with the position vector.
ball = sphere(pos = vector(10, 10, 10), radius = 9e2, color = color.red)
ballColor = ball.color
ballPos = ball.pos
ballRadius = ball.radius

# You can draw arrows.
# These also have a color and position field, but instead of a radius, they have an axis that determines their length.
velocityArrow = arrow(color = color.blue, pos = vector(-10, -10, -10), axis = velocity * vScale)

# There are also boxes.
# They have a position vector and a size vector.
myBox = box(pos = vector(50, 50, 50), size = vector(20, 15, 12))

# You can draw helixes, which are useful for depicting springs.
# They have a position field, a color field, a thickness field, a radius field, and a field for the number of coils.
spring = helix(pos=ceiling.pos, color=color.cyan, thickness=.003, coils=40, radius=0.015)

# If you want to trace a moving object, you can have a curve follow it by adding to it every time step.
posTrace = curve(color = color.yellow)
trail.append(ball.pos)

# You also have the option of making graphs as your program progresses; this is particularly useful for potential vs kinetic energy.
Kgraph = gcurve(color = color.cyan)
Ugraph = gcurve(color = color.yellow)
KplusUgraph = gcurve(color = color.red)
# For each time step...
Kgraph.plot(pos = (t, Kenergy))
Ugraph.plot(pos = (t, Uenergy))
KplusUgraph.plot(pos = (t, Kenergy + Uenergy))

By updating the fields of whatever objects you decide to use and by using any intermediate variables necessary, you're able to draw and compute what you'll need for this course.

Python Syntax

2017-11-29T18:15:07Z

Zsanchez6: /* VPython Basics */

Edited by Zach Sanchez - Fall 2017 (in progress)

==The Main Idea==
VPython is an extension of the Python programming language that contains a 3D graphics module called Visual, which allows users to create simple simulations of physics problems. Its primary use is for educational purposes, although it has been used in research before in the past. We will use VPython on occasion to aid us in visualizing and experimenting with the ideas that we learn in this course.

==Downloading vPython==
In all cases you will need to have Python 2.7.9 installed. Future versions may also work, but be warned that it may cause VPython to behave erratically.

[http://vpython.org/contents/download_windows.html For Windows]

[http://vpython.org/contents/download_mac.html For Mac]

[http://vpython.org/contents/download_linux.html For Linux]

Note for Linux: You can only install Classic VPython, which is no longer supported and may be missing some features and bug fixes that are present in the current version of VPython.

Alternatively, you can use GlowScript, a virtual environment that allows you to execute VPython code in your browser. You can create an account at [http://www.glowscript.org/ http://www.glowscript.org/].

==Python Basics==
It's always helpful to have a background knowledge of Python, however basic, before proceeding with the applications more relevant to this course.

# Comments start with a number symbol. Use these liberally in order to keep track of the changes you make/the code you write.
# We're always going to start our programs with the following headers to ensure that math works correctly.
from__future__ import division
from visual import*

# Math works as you would expect.
1 + 1 # -> 2
2 - 1 # -> 1
3 * 2 # -> 6
8 / 4 # -> 4
(1 + 3) * 2 # -> 8
1 + (3 * 2) # -> 7

# That is, until you need to do an exponent. But you shouldn't encounter anything much weirder than that math-wise.
2 ** 3 # -> 8

# There's also various math functions you can use, such as for finding the square root.
sqrt(4) # -> 2

# However, these math operations are useless if we don't store the values anywhere.
x = 1 + 1
# Now if you try to access x in the future, it will have the value of '2'.

# You can print your variables like this:
print(x)
# You can also add in some text with them like this:
print("The value of x is ", x)

# Sometimes we only want to execute something based on a certain condition.
if x == 2:
y = 3
# The '==' means that you're directly comparing two values. If they're equal, whatever is on the inside of the if-statement will be executed.
# Now is a good time to note that whitespace is very important in Python.
# 'y = 3' will always execute if it isn't tabbed, so even if x is equal to 7, y will still become 3.
# If you ever see me use whitespace in this tutorial, make sure you do as well.
# Some other useful comparators: != (not equal), >, >= (greater than or equal to), <, <= (less than or equal to), 'and', 'or', 'not'
# You can add onto the above if-statement by tacking on:
elif x == 3:
y = 4
else:
y = 5
# With elif meaning "else if" for a second possible condition, and "else" for something that you want to execute if none of the earlier conditions are met.

# You can also have a loop that executes code for a certain amount of time or based on a certain condition.
while x < 5:
print("x is ", x)
x = x + 1
# This will output the value of x up until it reaches 4.

Thus concludes our very, very brief foray into Python. That was simply the bare minimum required to understand what will be needed for this class. More concepts will be introduced as needed, whether that be in this tutorial or in a lab. For a more complete understanding of Python, you can look up other Python tutorials elsewhere on the Internet.

==VPython Basics==

In this section we will analyze some of the objects utilized by VPython. Each object is more or less self-explanatory with equally clear fields within.

# The scene has a width and a height; you can set them like so.
scene.width = 1024
scene.height = 760

# Beyond that, we have vectors.
# You can access each component with pos.x, pos.y, or pos.z depending on which component you want.
# Now would be a good time to mention that you access an object's field with a period following the variable name.
pos = vector(0, 0, 0)
xComponent = pos.x
yComponent = pos.y
zComponent = pos.z

# You can add two or more vectors together and you can multiply a vector by a scalar,
# but you can't add a scalar to a vector.
# To get the magnitude of or to normalize a vector, use the appropriate function.
initialPos = vector(10, 3, 0)
finalPos = vector(30, 15, 0)
deltaR = finalPos - initialPos # -> vector(20, 12, 0)
rMag = mag(finalPos)
rHat = norm(finalPos)
errorDemo = magR + finalPos # -> error; causes your program to crash

# If you were to make, say, a ball, you would want to draw a sphere.
# It has a position field, a radius field, and a color field.
# When creating a new instance of an object, each field is comma separated.
# You access each field the same way as we did with the position vector.
ball = sphere(pos = vector(10, 10, 10), radius = 9e2, color = color.red)
ballColor = ball.color
ballPos = ball.pos
ballRadius = ball.radius

# You can draw arrows.
# These also have a color and position field, but instead of a radius, they have an axis that determines their length.
velocityArrow = arrow(color = color.blue, pos = vector(-10, -10, -10), axis = velocity * vScale)

# There are also boxes.
# They have a position vector and a size vector.
myBox = box(pos = vector(50, 50, 50), size = vector(20, 15, 12))

# You can draw helixes, which are useful for depicting springs.
# They have a position field, a color field, a thickness field, a radius field, and a field for the number of coils.
spring = helix(pos=ceiling.pos, color=color.cyan, thickness=.003, coils=40, radius=0.015)

# If you want to trace a moving object, you can have a curve follow it by adding to it every time step.
posTrace = curve(color = color.yellow)
trail.append(ball.pos)

# You also have the option of making graphs as your program progresses; this is particularly useful for potential vs kinetic energy.
Kgraph = gcurve(color = color.cyan)
Ugraph = gcurve(color = color.yellow)
KplusUgraph = gcurve(color = color.red)
# For each time step...
Kgraph.plot(pos = (t, Kenergy))
Ugraph.plot(pos = (t, Uenergy))
KplusUgraph.plot(pos = (t, Kenergy + Uenergy))

By updating the fields of whatever objects you decide to use and by using any intermediate variables necessary, you're able to draw and compute what you'll need for this course.

Python Syntax

2017-11-29T18:10:09Z

Zsanchez6:

Edited by Zach Sanchez - Fall 2017 (in progress)

==The Main Idea==
VPython is an extension of the Python programming language that contains a 3D graphics module called Visual, which allows users to create simple simulations of physics problems. Its primary use is for educational purposes, although it has been used in research before in the past. We will use VPython on occasion to aid us in visualizing and experimenting with the ideas that we learn in this course.

==Downloading vPython==
In all cases you will need to have Python 2.7.9 installed. Future versions may also work, but be warned that it may cause VPython to behave erratically.

[http://vpython.org/contents/download_windows.html For Windows]

[http://vpython.org/contents/download_mac.html For Mac]

[http://vpython.org/contents/download_linux.html For Linux]

Note for Linux: You can only install Classic VPython, which is no longer supported and may be missing some features and bug fixes that are present in the current version of VPython.

Alternatively, you can use GlowScript, a virtual environment that allows you to execute VPython code in your browser. You can create an account at [http://www.glowscript.org/ http://www.glowscript.org/].

==Python Basics==
It's always helpful to have a background knowledge of Python, however basic, before proceeding with the applications more relevant to this course.

# Comments start with a number symbol. Use these liberally in order to keep track of the changes you make/the code you write.
# We're always going to start our programs with the following headers to ensure that math works correctly.
from__future__ import division
from visual import*

# Math works as you would expect.
1 + 1 # -> 2
2 - 1 # -> 1
3 * 2 # -> 6
8 / 4 # -> 4
(1 + 3) * 2 # -> 8
1 + (3 * 2) # -> 7

# That is, until you need to do an exponent. But you shouldn't encounter anything much weirder than that math-wise.
2 ** 3 # -> 8

# There's also various math functions you can use, such as for finding the square root.
sqrt(4) # -> 2

# However, these math operations are useless if we don't store the values anywhere.
x = 1 + 1
# Now if you try to access x in the future, it will have the value of '2'.

# You can print your variables like this:
print(x)
# You can also add in some text with them like this:
print("The value of x is ", x)

# Sometimes we only want to execute something based on a certain condition.
if x == 2:
y = 3
# The '==' means that you're directly comparing two values. If they're equal, whatever is on the inside of the if-statement will be executed.
# Now is a good time to note that whitespace is very important in Python.
# 'y = 3' will always execute if it isn't tabbed, so even if x is equal to 7, y will still become 3.
# If you ever see me use whitespace in this tutorial, make sure you do as well.
# Some other useful comparators: != (not equal), >, >= (greater than or equal to), <, <= (less than or equal to), 'and', 'or', 'not'
# You can add onto the above if-statement by tacking on:
elif x == 3:
y = 4
else:
y = 5
# With elif meaning "else if" for a second possible condition, and "else" for something that you want to execute if none of the earlier conditions are met.

# You can also have a loop that executes code for a certain amount of time or based on a certain condition.
while x < 5:
print("x is ", x)
x = x + 1
# This will output the value of x up until it reaches 4.

Thus concludes our very, very brief foray into Python. That was simply the bare minimum required to understand what will be needed for this class. More concepts will be introduced as needed, whether that be in this tutorial or in a lab. For a more complete understanding of Python, you can look up other Python tutorials elsewhere on the Internet.

==VPython Basics==

In this section we will analyze some of the objects utilized by VPython. Each object is more or less self-explanatory with equally clear fields within.

# The scene has a width and a height; you can set them like so.
scene.width = 1024
scene.height = 760

# Beyond that, we have vectors.
# You can access each component with pos.x, pos.y, or pos.z depending on which component you want.
# Now would be a good time to mention that you access an object's field with a period following the variable name.
pos = vector(0, 0, 0)
xComponent = pos.x
yComponent = pos.y
zComponent = pos.z

# You can add two or more vectors together and you can multiply a vector by a scalar,
# but you can't add a scalar to a vector.
# To get the magnitude of or to normalize a vector, use the appropriate function.
initialPos = vector(10, 3, 0)
finalPos = vector(30, 15, 0)
deltaR = finalPos - initialPos # -> vector(20, 12, 0)
rMag = mag(finalPos)
rHat = norm(finalPos)
errorDemo = magR + finalPos # -> error; causes your program to crash

# If you were to make, say, a ball, you would want to draw a sphere.
# It has a position field, a radius field, and a color field.
# When creating a new instance of an object, each field is comma separated.
# You access each field the same way as we did with the position vector.
ball = sphere(pos = vector(10, 10, 10), radius = 9e2, color = color.red)
ballColor = ball.color
ballPos = ball.pos
ballRadius = ball.radius

# You can draw arrows.
# These also have a color and position field, but instead of a radius, they have an axis that determines their length.
velocityArrow = arrow(color = color.blue, pos = vector(-10, -10, -10), axis = velocity * vScale)

# There are also boxes.
# They have a position vector and a size vector.
myBox = box(pos = vector(50, 50, 50), size = vector(20, 15, 12))

# You can draw helixes, which are useful for depicting springs.
# They have a position field, a color field, a thickness field, a radius field, and a field for the number of coils.
spring = helix(pos=ceiling.pos, color=color.cyan, thickness=.003, coils=40, radius=0.015)

# If you want to trace a moving object, you can have a curve follow it by adding to it every time step.
posTrace = curve(color = color.yellow)
trail.append(ball.pos)

# You also have the option of making graphs as your program progresses; this is particularly useful for potential vs kinetic energy.
Kgraph = gcurve(color = color.cyan)
Ugraph = gcurve(color = color.yellow)
KplusUgraph = gcurve(color = color.red)
# For each time step...
Kgraph.plot(pos = (t, Kenergy))
Ugraph.plot(pos = (t, Uenergy))
KplusUgraph.plot(pos = (t, Kenergy + Uenergy))

Using all of these concepts, you're able to draw and compute what you'll need for this course.

Python Syntax

2017-11-01T19:49:25Z

Zsanchez6: Added a basic Python tutorial, cleaned up the page some. Will be back to clean up some more and improve upon the VPython tutorial itself.

Edited by Zach Sanchez - Fall 2017 (in progress)

==The Main Idea==
VPython is an extension of the Python programming language that contains a 3D graphics module called Visual, which allows users to create simple simulations of physics problems. Its primary use is for educational purposes, although it has been used in research before in the past. We will use VPython on occasion to aid us in visualizing and experimenting with the ideas that we learn in this course.

==Downloading vPython==
In all cases you will need to have Python 2.7.9 installed. Future versions may also work, but be warned that it may cause VPython to behave erratically.

[http://vpython.org/contents/download_windows.html For Windows]

[http://vpython.org/contents/download_mac.html For Mac]

[http://vpython.org/contents/download_linux.html For Linux]

Note for Linux: You can only install Classic VPython, which is no longer supported and may be missing some features and bug fixes that are present in the current version of VPython.

Alternatively, you can use GlowScript, a virtual environment that allows you to execute VPython code in your browser. You can create an account at [http://www.glowscript.org/ http://www.glowscript.org/].

==Python Basics==
It's always helpful to have a background knowledge of Python, however basic, before proceeding with the applications more relevant to this course.

# Comments start with a number symbol. Use these liberally in order to keep track of the changes you make/the code you write.
# We're always going to start our programs with the following headers to ensure that math works correctly.
from__future__ import division
from visual import*

# Math works as you would expect.
1 + 1 # -> 2
2 - 1 # -> 1
3 * 2 # -> 6
8 / 4 # -> 4
(1 + 3) * 2 # -> 8
1 + (3 * 2) # -> 7

# That is, until you need to do an exponent. But you shouldn't encounter anything much weirder than that math-wise.
2 ** 3 # -> 8

# However, these math operations are useless if we don't store the values anywhere.
x = 1 + 1
# Now if you try to access x in the future, it will have the value of '2'.

# You can print your variables like this:
print(x)
# You can also add in some text with them like this:
print("The value of x is ", x)

# Sometimes we only want to execute something based on a certain condition.
if x == 2:
y = 3
# The '==' means that you're directly comparing two values. If they're equal, whatever is on the inside of the if-statement will be executed.
# Now is a good time to note that whitespace is very important in Python. 'y = 3' will always execute if it isn't tabbed, so even if x is equal to 7, y will still become 3.
# If you ever see me use whitespace in this tutorial, make sure you do as well.
# Some other useful comparators: != (not equal), >, >= (greater than or equal to), <, <= (less than or equal to), 'and', 'or', 'not'
# You can add onto the above if-statement by tacking on:
elif x == 3:
y = 4
else:
y = 5
# With elif meaning "else if" for a second possible condition, and "else" for something that you want to execute if none of the earlier conditions are met.

# You can also have a loop that executes code for a certain amount of time or based on a certain condition.
while x < 5:
print("x is ", x)
x = x + 1
# This will output the value of x up until it reaches 4.

Thus concludes our very, very brief foray into Python. That was simply the bare minimum required to understand what will be needed for this class. More concepts will be introduced as needed, whether that be in this tutorial or in a lab. For a more complete understanding of Python, you can look up other Python tutorials elsewhere on the Internet.

==Mathematical Model==
Vpython can compute any equation, but some that may be most helpful and most useful for Physics can be found below. (Keep in mind you can alter these numbers to be whatever you need, these are just to provide an example):

To update momentum:

pf = pi + Fnet*deltat

To update position:

objectf.pos = objecti.pos + (pcart/mcart)*deltat

To create a vector:

vector(0,0,0) -- fill in with whatever numbers the vector should be

Gravitational Force:

# Constants
G = 6.7e-11
mEarth = 6e24
mcraft = 15e3
deltat = 60
t = 0

Finds the change in position:

r=craft.pos-Earth.pos
m=mcraft

To find the magnitude of the change in position:

rmag= mag(r)

To calculate the new magnitude of gravitational force:

Fmag=(G*mcraft*mEarth)/(rmag**2)

To calculate the direction of the change in position:

rhat=r/rmag

To calculate net force:

Fnet=-Fmag*rhat

To calculate spring force:

L0 = 0.3
Lvec = ball.pos - ceiling.pos
Lhat = norm(Lvec)
Lmag = mag(Lvec)
Fspr = (-ks)*(Lmag - L0)*(Lhat)

To calculate kinetic energy:

Kinetic = (1/2)*(mball*(vel**2))

To create a graph:

gdisplay(width=500, height=250, x=600, y=1)
ygraph= gcurve(color=color.cyan)

To add plots to the graph:

ygraph.plot(pos=(t, fnet.y))

==Computational Model==

VPython is used to create computational, 3D models of various real world situations in order to better visualize how different equations can manipulate different scenarios. This is very valuable since many of the equations and situations that are coded in vPython are extremely difficult to make a functioning model of in real life.

For instance, the picture below is an example of program that was programmed to show the orbit of a craft around Earth.

[[File:Earthorbits.jpg]]

This picture shows a graph that tracks.
[[File:PositionGraph.jpg]]

==Examples==

Simple:

Creating Shapes:

Sphere:

sphere= sphere(pos=vector(-4,-2,5), radius=.4, color=color.red)

Arrow:

bt=arrow(pos=sphere.pos, axis=sphere2.pos-sphere.pos, color=color.cyan)

Vector:

vector=vector(0, 0, 0)

Trail:

trail = curve(color=sphere.color)
trail.append(pos=sphere.pos)

Setting Scene Range:

scene.range=11*sphere.radius

Helix:

spring = helix(pos=ceiling.pos, color=color.cyan, thickness=.003, coils=40, radius=0.015)

Intermediate:

Graphs:

Setup graphing windows:

gdisplay(width=500, height=250, x=600, y=1)
ygraph = gcurve(color=color.yellow)
gdisplay(width=500, height=250, x=600, y=300)

Plotting:

pgraph = gcurve(color=color.blue)
ygraph.plot(pos=(time, Fnet.y))
pgraph.plot(pos=(time, sphere.y))

Difficult:

Using Loops to update Equations:

Constants:

G = ?

mEarth = ?

mmoon = ?

mcraft = ?

deltat = ?

t = ?

Objects and initial values:

Earth = sphere(pos=vector(0,0,0), radius=6.4e6, color=color.cyan)
scene.range=11*Earth.radius
Moon = sphere(pos=(4e8, 0, 0), radius=1.75e6, color=color.white)

Add a radius for the spacecraft. It should be BIG, so it can be seen:

craft = sphere(pos=vector(-6.656e7,-3.648e6,0), radius= 10000, color=color.yellow)
vcraft = vector(206, 2645,0)
pcraft = mcraft*vcraft
pArrow=arrow(color=color.green)
fArrow=arrow(color=color.cyan)
dpArrow=arrow(color=color.red)
Fnet_tangent_arrow = arrow(color=color.yellow)
Fnet_perp_arrow= arrow(color=color.magenta)

This creates a trail for the spacecraft:

trail = curve(color=craft.color)

And this prevents zooming in or out:

scene.autoscale = 0
pscale=Earth.radius/mag(pcraft)
fscale=Earth.radius/((G*mEarth*mcraft)*mag(craft.pos-Earth.pos)**2)
dpscale=500*Earth.radius/mag(pcraft)
print("p=", pcraft)

Calculations:

Sets time for loop to run:

while t < 165240:

This slows down the animation (runs faster with bigger number):

rate(10000)

Add statements here for the iterative update of gravitational
force, momentum, and position.

r = craft.pos-Earth.pos
rmag = sqrt(r.x**(2)+r.y**(2)+r.z**(2))
Fmag= G*mEarth*mcraft/(rmag**2)
rhat= r/rmag
rmoon= craft.pos - Moon.pos
rmoonmag= mag(rmoon)
rmoonhat= norm(rmoon)
Fmoonmag= G*mmoon*mcraft/(rmoonmag**2)
Fmoon= -Fmoonmag*rmoonhat
p_init= mag(pcraft)
pcraft_i=pcraft+vector(0,0,0)
Fearth= -Fmag*rhat
Fnet= Fearth + Fmoon
pcraft=Fnet*deltat+pcraft
p_final=mag(pcraft)

Fnet_tangent = (p_final-p_init)*norm(pcraft)/deltat
Fnet_tangent_arrow.pos=craft.pos
Fnet_tangent_arrow.axis=Fnet_tangent*fscale
Fnet_perp = Fnet-Fnet_tangent
Fnet_perp_arrow.pos=craft.pos
Fnet_perp_arrow.axis=Fnet_perp*fscale
vcraft=pcraft/mcraft
craft.pos=vcraft*deltat+craft.pos
pArrow.pos=craft.pos
pArrow.axis=pcraft*pscale
fArrow.pos=craft.pos
fArrow.axis=Fnet*fscale
deltap= pcraft-pcraft_i
dpArrow.pos=craft.pos
dpArrow.axis=deltap*dpscale
scene.center=craft.pos
scene.range=craft.radius*600

Uncomment these two lines to exit the loop if
the spacecraft crashes onto the Earth.

if rmag < Earth.radius:
break

trail.append(pos=craft.pos)
t = t+deltat