Milankowitch Cycles. 1 Coordinates. 2 Rotation matrix of Eulerian angles. (%i1) kill(all); (%o0) done
|
|
- Muriel Lyons
- 5 years ago
- Views:
Transcription
1 369(9-rot.wxm 1 / 1 Milankowitch Cycles (%i1 kill(all; (%o one 1 Coorinates (%i1 epens([x,y,theta,phi,psi,theta_1,t,u,l,r,omega],t; (%o1 [x( t,y( t,( t,'( t,( t, 1 ( t,t( t,u( t,l( t,r( t,!( t ] (%i e_r: [sin(theta*cos(phi, sin(theta*sin(phi, cos(theta]; (%o [cos(' sin(,sin(' sin(,cos( ] Rotation matrix of Eulerian angles (%i3 Rz: matrix([cos(phi, sin(phi, ], [-sin(phi, cos(phi, ], [,, 1]; (%o3 cos( ' sin( ' cos( ' 1 (%i4 Ry: matrix([cos(theta,, sin(theta], [, 1, ], [-sin(theta,, cos(theta]; (%o4 1 (%i5 Rz1: matrix([cos(psi, sin(psi, ], [-sin(psi, cos(psi, ], [,, 1]; (%o5 cos( sin( sin( cos( 1 (%i6 R3: Rz1.Ry.Rz; (%o6 cos( ' cos( sin( cos( ' sin( cos( cos( ' sin( ' cos( cos( +cos( ' sin( cos( ' cos( sin( cos( sin(
2 369(9-rot.wxm / 1 3 bol r_1 (%i7 r_1: transpose([r11,r1,r13].r3; (%o7 cos( ' r13 sin( +r1 ( cos( ' sin( cos( +r11 (cos(' cos( r13 sin( +r1 ( cos( ' cos( sin( cos( +r11 (sin(' cos( cos( +cos( ' sin( r1 sin( +cos( r11 sin( +r13 (%i1 X_1: first(r_1[1]; Y_1: first(r_1[]; Z_1: first(r_1[3]; (%o8 cos( ' r13 sin( +r1 ( cos( ' cos( +r11 ( cos( ' (%o9 r13 sin( +r1 ( cos( ' cos( +r11 ( sin( ' cos( +cos( ' (%o1 r1 sin( + r11 sin( +r13 Check (%i13 r_1x: r11*r3[1,1] +r1*r3[,1] +r13*r3[3,1]; r_1y: r11*r3[1,] +r1*r3[,] +r13*r3[3,]; r_1z: r11*r3[1,3] +r1*r3[,3] +r13*r3[3,3]; (%o11 cos( ' r13 sin( +r1 ( cos( ' cos( +r11 ( cos( ' (%o1 r13 sin( +r1 ( cos( ' cos( +r11 ( sin( ' cos( +cos( ' (%o13 r1 sin( + r11 sin( +r13 (%i14 r_1-[r_1x, r_1y, r_1z]; sin sin (%o14 4 bol r (%i15 [X, Y, Z]: [r*cos(theta_1, r*sin(theta_1, ]; (%o15 [r cos( 1,r sin( 1,] (%i16 [X, Y, Z]: [,, ]; (%o16 [,,] 5 Kinetic energy
3 369(9-rot.wxm 3 / 1 (%i19 T[1]: m/*(iff(x+x_1,t^; T[]: m/*(iff(y+y_1,t^; T[3]: m/*(iff(z+z_1,t^; (%o17 (m (r1 (cos(' t cos( ' t cos( + sin( ' t ' cos( +sin( ' t cos( ' t ' cos( + r11 ( cos(' t ' r13 t + t cos( ' ' t ' r13 sin( ^ / t t ' t ' sin( cos(' (%o18 (m (r1 (sin(' t t cos( ' t ' cos( ' t t ' cos( + r11 ( sin(' t ' cos( +cos( ' r13 t t ' t ' r13 sin( ^ / t t cos( +cos( ' ' t ' sin( sin(' (%o19 (m ( r13 t sin( r1 t +cos( r11 cos( t r1 t r11 sin( ^ / t (%i T: (T[1]+T[]+T[3]$ 6 Potential energy (%i1 U: -m*m*g/r; (%o1 G M m r (%i U: ; (%o 7 Lagrange function
4 369(9-rot.wxm 4 / Simplify r_1: r1=r13=, a point on the equator (%i3 r1: r13: ; (%o3 (%i4 L: ev(t - U; (%o4 (m r11 ( sin(' cos( ' t ' cos( + / +(m r11 ( cos(' sin( ' m / + ' t ' cos( r11 t 8 Lagrange equations II 8.1 phi equation (%i5 D1: iff(l, iff(phi,t; t t t cos( ' t r11 t (%o5 m r11 ( cos( ' sin( ' sin( t t cos( +cos( ' ' t t ' sin( +m r11 ( cos( ' t cos( ' t ' ' t cos( + t ' sin( ^ t t ' sin( ^ ' ( sin(' t ' cos( +cos( ' ( cos(' t ' t cos( ' t
5 369(9-rot.wxm 5 / 1 (%i6 E1: (trigsimp(ratsimp(iff(d1,t - iff(l,phi = ; (%o6 m r11 t + ( m m +m t +m cos( t r11 r11 m t ' r11 t + m t ' t m t ' t r11 cos( +m t ' r11 = (%i7 E11: (solve(e1, iff(phi,t,; (%o7 [ '=( t + ( t t + 8. theta equation + t t ' (%i8 D: iff(l, iff(theta,t; t cos( /(cos( 1] t r11 t + t ' (%o8 m sin( ' r11 sin( ( sin(' t sin( ' t cos( +cos( ' t ' cos( +cos( ' t sin( ' t ' sin( m cos( ' r11 sin( ( cos(' t cos( ' t r11 t ' ' t ' sin( +m r11 cos( t t ' r11 t (%i9 E: trigsimp(ratsimp(iff(d,t - iff(l,theta = ; (%o9 m r11 t ' r11 cos( + m t ' m t t +m t r11 t +(m ' t r11 sin( =
6 369(9-rot.wxm 6 / 1 (%i3 E1: solve(e, iff(theta,t,; (%o3 [ =( sin( t t ' / cos( ] 8.3 psi equation t t + t ' (%i31 D3: iff(l, iff(psi,t; t (%o31 m r11 ( cos( ' sin( t t cos( +cos( ' ' t t ' sin( +m r11 ( cos( ' t cos( ' ' sin( t ( sin(' t ' cos( +cos( ' ( cos(' t ' t ' sin( m r11 sin( r11 t r11 t (%i3 E3: trigsimp(ratsimp(iff(d3,t - iff(l,psi = ; (%o3 m m r11 t m t ' r11 +m r11 = (%i33 E31: solve(e3, iff(psi,t,; (%o33 [ = t t ' 8.4 theta_1 equation t + t cos( ' t t ' r11 ' cos( ] t (%i34 D4: iff(l, iff(theta_1,t; (%o34 ' t r11 cos( +m t ' t t t +
7 369(9-rot.wxm 7 / 1 (%i35 E4: trigsimp(ratsimp(iff(d4,t - iff(l,theta_1 = ; (%o35 = (%i36 E41: solve(e4, iff(theta_1,t,; (%o36 all 8.5 r equation (%i37 D5: iff(l, iff(r,t; (%o37 (%i38 E5: trigsimp(ratsimp(iff(d5,t - iff(l,r = ; (%o38 = (%i39 E51: solve(e5, iff(r,t,; (%o39 all 9 Prepare solution (%i4 GC: first(solve([first(e1,first(e,first(e3,first(e4,first(e5], [iff(phi,t,, iff(theta,t,, iff(psi,t,, iff(theta_1,t, solve: epenent equations eliminate: (4 5 (%i41 transform(eq := (A: ratsubst(phi_, iff(phi,t, Eq, A: ratsubst(phi_, iff(phi,t,, A, A: ratsubst(theta_, iff(theta,t, A, A: ratsubst(theta_, iff(theta,t,, A, A: ratsubst(psi_, iff(psi,t, A, A: ratsubst(psi_, iff(psi,t,, A, A: ratsubst(theta_1, iff(theta_1,t, A, A: ratsubst(theta_1, iff(theta_1,t,, A, A: ratsubst(r_, iff(r,t, A, ratsubst(r_, iff(r,t,, A ; (%o41 transform( Eq :=(A:ratsubst ', ',Eq,A: t ratsubst phi_, ',A,A:ratsubst t,,a,a: t ratsubst theta_,,a,a:ratsubst t,,a,a: t ratsubst psi_,,a,a:ratsubst theta_1, t t 1,A,A: ratsubst ratsubst theta_1, r, t r,a t 1,A,A:ratsubst r, t r,a,
8 369(9-rot.wxm 8 / 1 (%i46 G1: transform(gc[1]$ G: transform(gc[]$ G3: transform(gc[3]$ G4: transform(gc[4]$ G5: transform(gc[5]$ 1 Solution (%i47 str: [r11=1, m=1, G=, M=]; (%o47 [r11=1,m=1,g =,M =] (%i57 Eq1: gfactor(ev(rhs(g1, str, eval$ Eq: gfactor(ev(rhs(g, str, eval$ Eq3: gfactor(ev(rhs(g3, str, eval$ Eq4: gfactor(ev(rhs(g4, str, eval$ Eq5: gfactor(ev(rhs(g5, str, eval$ Eq6: phi_; Eq7: theta_; Eq8: psi_; Eq9: theta_1; Eq1: r_; (%o53 ' (%o54 (%o55 (%o56 theta_1 (%o57 r (%i93 s: rk([eq1, Eq, Eq3, Eq6, Eq7, Eq8], [phi_, theta_, psi_, phi, theta, psi], [.1,.1,.1,., %pi/, ], [t,,.1,.1]$
9 369(9-rot.wxm 9 / 1 (%i15 Eq1;Eq;Eq3; (%o13 ( sin( + ( sin( (%o14 ( + +( ' ' 3 cos( + ' ' ' 3 ' +' cos( + ' sin( /( ( (%o15 ((( sin( +( ' ' ' +( + cos( + ' sin( + ' 3 ' ' ' /(( (%i116 trigsimp(enom(eq1;trigsimp(num(eq1; (%o115 (%o Graphics (%i99 c1: makelist([p[1],p[]],p,s$ c: makelist([p[1],p[3]],p,s$ c3: makelist([p[1],p[4]],p,s$ c4: makelist([p[1],p[5]],p,s$ c5: makelist([p[1],p[6]],p,s$ c6: makelist([p[1],p[7]],p,s$ 1. Plot theta, phi, psi (%i11 c1;c; (%o1 [[.,.1]] (%o11 [[.,.1]]
10 369(9-rot.wxm 1 / 1 (%i78 wxplot([[iscrete, c1], [iscrete, c]], [xlabel, "time"], [legen, "phi_{ot}", "theta_{ot}"]$ Message from maxima's sterr stream: Warning: empty x range [:], ajusting (%t78
11 369(9-rot.wxm 11 / 1 (%i79 wxplot([[iscrete, c3], [iscrete, c4]], [xlabel, "time"], [legen, "psi_{ot}", "theta_1_{ot}"]$ Message from maxima's sterr stream: Warning: empty x range [:], ajusting Warning: empty y range [:], ajusting to [-1:1] (%t79
12 369(9-rot.wxm 1 / 1 (%i8 wxplot([[iscrete, c5], [iscrete, c6]], [xlabel, "time"], [legen, "r_{ot}", "phi"]$ Message from maxima's sterr stream: Warning: empty x range [:], ajusting (%t8 ; ;
Gyroscope free falling with R
369(1b.wxm 1 / 13 Gyroscope free falling with R (%i1 kill(all; (%o0 one 1 Coorinates (%i1 epens([x,y,theta,phi,psi,r,v,t,u,l,phi,omega],t; (%o1 [x( t,y( t,( t,'( t,( t,r( t,v( t,t( t,u( t,l( t,phi( t,
More informationGyroscope with variable R
369(1a).wxm 1 / 11 Gyroscope with variable R (%i1) (%o0) kill(all); one 1 Coorinates (%i1) (%o1) t],!() epens([x,y,theta,phi,psi,r,v,t,u,l,phi,omega],t); t,y() t,ò() t,'() t,é() t,r() t,v() t,t() t,u()
More informationSpinning top - represented by dumbbell
Lagange-umbbell.wxm 1 / 10 Spinning top - epesente by umbbell (%i1) kill(all); (%o0) one 1 Cooinates (%i1) epens([x,y,z,x,y,z,theta,phi,theta_1,phi_1,,t,u,l],t); (%o1) [x( t ),y( t ),z( t ),X( t ),Y( t
More informationWarning, the protected names norm and trace have been redefined and unprotected
22M:28 Spring 5 J. Simon An example of computing a triple integral with several coordinate systems > > with(plots):with(linalg): Warning, the name changecoords has been redefined Warning, the protected
More informationPhysics 339 Euler Angles & Free Precession November 2017 Hamilton s Revenge
Physics 339 Euler Angles & Free Precession November 217 Hamilton s Revenge As described in the textbook, Euler Angles are a way to specify the configuration of a 3d object. Starting from a fixed configuration
More informationFiltering Color Mapped Textures and Surfaces
Filtering Color Mapped Textures and Surfaces Eric Heitz, Derek Nowrouzezahrai Pierre Poulin Fabrice Neyret Universite de Grenoble and CNRS Universite de Montre al Abstract This document provides additional
More informationME352 Project 4. Internal Combustion Engine Analysis. Jeremy Bourque Lab Division 9 April 24 th, 2018
ME352 Project 4 Internal Combustion Engine Analysis Jeremy Bourque Lab Division 9 April 24 th, 2018 Summary of Report This project aims to model an internal combustion engine from a lawnmower as a slidercrank
More informationR e p u b lic o f th e P h ilip p in e s. R e g io n V II, C e n tra l V isa y a s. C ity o f T a g b ila ran
R e p u b l f th e P h lp p e D e p rt e t f E d u t R e V, e tr l V y D V N F B H L ty f T b l r Ju ly, D V N M E M R A N D U M N. 0,. L T F E N R H G H H L F F E R N G F R 6 M P L E M E N T A T N T :,
More information2- Scalars and Vectors
2- Scalars and Vectors Scalars : have magnitude only : Length, time, mass, speed and volume is example of scalar. v Vectors : have magnitude and direction. v The magnitude of is written v v Position, displacement,
More informationJackson, Classical Electrodynamics, Section 14.8 Thomson Scattering of Radiation
High Energy Cross Sections by Monte Carlo Quadrature Thomson Scattering in Electrodynamics Jackson, Classical Electrodynamics, Section 14.8 Thomson Scattering of Radiation Jackson Figures 14.17 and 14.18:
More informationforce reduces appropriately to the force exerted by one point charge on another. (20 points)
Phsics III: Theor an Simulation Examination 3 December 4, 29 Answer All Questions Analtical Part: Due 5: p.m., M, 12/7/9 Name SOUTIONS 1. Two line charges A an B of the same length are parallel to each
More information. ) f. i J f. q!!;."" ~ ~ i. ~ Jl j. v ~ "' y; f:j. f' ~ ~ lj ~ e~ 1 t;i' c! f{ 1 tl. ~), jl ( 1 3 ~ ~ /' h 'I. ~ - ~; ~ / jcd-j 1 ~ ~ 1~ ~ ...
{I I _j'..j 1...J r; tl,, tl r t 0. 5 t f j J J _.. /; 4 1 j 'J J r v, 1 1 1 j tl 1 1 I :1 I),.1 (1 1 r- ' 1f.t j q!!;."" J,...Y' 3 1.:f,., v "' C> ' ev J\ j A D d l d. ) f q { j j < / ( 1 i i J f rjl
More informationAN ERROR ANALYSIS OF 3-D ELECTROMAGNETIC TRACKERS. A Dissertation Presented. Katherine Gallo. The Faculty of the Graduate College
AN ERROR ANALYSIS OF 3-D ELECTROMAGNETIC TRACKERS A Dissertation Presented by Katherine Gallo to The Faculty of the Graduate College of The University of Vermont In Partial Fullfillment of the Requirements
More informationAppendix A Solution of Waveform Design Optimization Problem
Appix A Solution of Waveform Design Optimization Problem A.1 Preliminaries This section presents some preliminary results used in the proofs throughout the paper. For proofs of the following theorems,
More informationWorkbook on STATICS Department of Mechanical Engineering LINGAYA S UNIVERSITY Nachauli, Old Faridabad Jasana Road Faridabad , Haryana, India
Workbook on STATICS Department of Mechanical Engineering LINGAYA S UNIVERSITY (Deemed to be University under Section 3 of UGC Act, 1956) Nachauli, Old Faridabad Jasana Road Faridabad 121 2, Haryana, India
More informationt t t ér t rs r t ét q s
rés té t rs té s é té r t q r r ss r t t t ér t rs r t ét q s s t t t r2 sé t Pr ss r rs té P r s 2 t Pr ss r rs té r t r r ss s Pr ss r rs té P r q r Pr ss r t r t r r t r r Prés t r2 r t 2s Pr ss r rs
More informationI I. R E L A T E D W O R K
A c c e l e r a t i n g L a r g e S c a l e C e n t r o i d - B a s e d C l u s t e r i n g w i t h L o c a l i t y S e n s i t i v e H a s h i n g R y a n M c C o n v i l l e, X i n C a o, W e i r u L
More informationProblem (a) θ 0 = 0
Problem 4.1.11 (a) θ 0 = 0 Beam Pattern: 3 For circular symmetric beam pattern, the weights a nm are given by ( ψ R J 1 ψ R n ) + m a(n, m) = π n + m In order to be consistent with Problem 4.1.5, we take
More information? D. 3 x 2 2 y. D Pi r ^ 2 h, r. 4 y. D 3 x ^ 3 2 y ^ 2, y, y. D 4 x 3 y 2 z ^5, z, 2, y, x. This means take partial z first then partial x
PartialsandVectorCalclulus.nb? D D f, x gives the partial derivative f x. D f, x, n gives the multiple derivative n f x n. D f, x, y, differentiates f successively with respect to x, y,. D f, x, x 2, for
More information5.2 - Euler s Method
5. - Euler s Method Consider solving the initial-value problem for ordinary differential equation: (*) y t f t, y, a t b, y a. Let y t be the unique solution of the initial-value problem. In the previous
More informationA L A BA M A L A W R E V IE W
A L A BA M A L A W R E V IE W Volume 52 Fall 2000 Number 1 B E F O R E D I S A B I L I T Y C I V I L R I G HT S : C I V I L W A R P E N S I O N S A N D TH E P O L I T I C S O F D I S A B I L I T Y I N
More informationMassachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.061/6.690 Introduction to Power Systems
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.061/6.690 Introduction to Power Systems Problem Set 9 Solutions April 18, 011 Problem 1: Chapter 9, Problem
More informationKinetic energy. Analytical Mechanics: Link Mechanisms. Kinetic energy. Agenda. Kinematics of two link open mechanism.
Analtical Mechanics: Link Mechanisms Shinichi Hirai Dept. Robotics, Ritsumeikan Univ. Kinetic energ velocit of the center of mass of link : angular velocit of link : kinetic energ of link : ẋ c l c θ S
More informationSystems of Linear Equations
Systems of Linear Equations As stated in Section G, Definition., a linear equation in two variables is an equation of the form AAAA + BBBB = CC, where AA and BB are not both zero. Such an equation has
More information2. MATLAB Basics. (d) array1([1 3],end) consists of the elements in the first and third rows on the last column of array1:
2. MATLAB Basics 2.1 (a) The size of array1 is 4 5. (b) The value of array1(1,4) is -3.5. (c) array1(:,1:2:5) is a 4 3 array consisting of the first, third, and fifth columns of array1: >> array1(:,1:2:5)
More information1.1.3 The narrowband Uniform Linear Array (ULA) with d = λ/2:
Seminar 1: Signal Processing Antennas 4ED024, Sven Nordebo 1.1.3 The narrowband Uniform Linear Array (ULA) with d = λ/2: d Array response vector: a() = e e 1 jπ sin. j(π sin )(M 1) = 1 e jω. e jω(m 1)
More informationAlso interested only in internal energies Uel (R) only internal forces, has symmetry of molecule--that is source of potential.
IV. Molecular Vibrations IV-1 As discussed solutions, ψ, of the amiltonian, (Schrödinger Equation) must be representations of the group of the molecule i.e. energy cannot change due to a symmetry operation,
More informationApplication 10.5D Spherical Harmonic Waves
Applicatio 10.5D Spherical Haroic Waves I probles ivolvig regios that ejoy spherical syetry about the origi i space, it is appropriate to use spherical coordiates. The 3-diesioal Laplacia for a fuctio
More informationy + 3y = 0, y(0) = 2, y (0) = 3
MATH 3 HOMEWORK #3 PART A SOLUTIONS Problem 311 Find the solution of the given initial value problem Sketch the graph of the solution and describe its behavior as t increases y + 3y 0, y(0), y (0) 3 Solution
More informationTransformations. Lars Vidar Magnusson. August 24,
Transformations Lars Vidar Magnusson August 24, 2012 http://www.it.hiof.no/~larsvmag/iti43309/index.html 2D Translation To translate an object is to move it in two-dimensinal space. If we have a point
More informationMassachusetts Institute of Technology Department of Electrical Engineering and Computer Science Electric Machines
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.685 Electric Machines Problem Set Solutions September 17, 5 Problem 1: First, note that λ 1 = L 1 i 1 +
More informationAnalytical Mechanics: Variational Principles
Analytical Mechanics: Variational Principles Shinichi Hirai Dept. Robotics, Ritsumeikan Univ. Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: Variational Principles 1 / 71 Agenda
More informationSTEEL PIPE NIPPLE BLACK AND GALVANIZED
Price Sheet Effective August 09, 2018 Supersedes CWN-218 A Member of The Phoenix Forge Group CapProducts LTD. Phone: 519-482-5000 Fax: 519-482-7728 Toll Free: 800-265-5586 www.capproducts.com www.capitolcamco.com
More informationUnit 4 Example Review. 1. Question: a. [A] b. [B] c. [C] d. [D]
Unit 4 Example Review 1. Question: Answer: When the line is wrapped around the circle, the length of the line will be equal to the length of the arc that is formed by the wrap starting from the point (1,0)
More informationComputation of Multibody Kinematics And Dynamics 1.1 Multibody Mechanical Systems. 2D Planar Mechanism. 3D Spatial Mechanism
Computation of Multibody Kinematics And Dynamics 1.1 Multibody Mechanical Systems D Planar Mechanism D Spatial Mechanism 1 Computer-Aided Design (CAD) Mechanics: Statics and dynamics. Dynamics : kinematics
More informationVICTORIA POWERS AND THORSTEN W ORMANN Suppose there exists a real, symmetric, psd matrix B such that f = x B x T and rank t. Since B is real symmetric
AN ALGORITHM FOR SUMS OF SQUARES OF REAL POLYNOMIALS Victoria Powers and Thorsten Wormann Introduction We present an algorithm to determine if a real polynomial is a sum of squares (of polynomials), and
More information1. State the amplitude and period for the function y = -3 sin 3. 3 * theta = 2* pi (Because, 2*pi is period for sine function
Unit 4 Review 1. State the amplitude and period for the function y = -3 sin 3. Amplitude: 3 (Always a positive quanitity) 3 * theta = 2* pi (Because, 2*pi is period for sine function theta = pi * (2/3)
More informationQuantitative Understanding in Biology Module III: Linear Difference Equations Lecture II: Complex Eigenvalues
Quantitative Understanding in Biology Module III: Linear Difference Equations Lecture II: Complex Eigenvalues Introduction In the previous section, we considered the generalized twovariable system of linear
More informationStatic Unbalance. Both bearing reactions occur in the same plane and in the. After balancing no bearing reaction remains theoretically.
BALANCING Static Unbalance A disk-shaft system on rigid rails An unbalanced disk-shaft system mounted on bearings and rotated For the unbalanced rotating system: Dynamic Bearing Reactions + Inertia Forces
More informationMarine Hydrodynamics Prof.TrilochanSahoo Department of Ocean Engineering and Naval Architecture Indian Institute of Technology, Kharagpur
Marine Hydrodynamics Prof.TrilochanSahoo Department of Ocean Engineering and Naval Architecture Indian Institute of Technology, Kharagpur Lecture - 10 Source, Sink and Doublet Today is the tenth lecture
More informationD Alembert s principle of virtual work
PH101 Lecture 9 Review of Lagrange s equations from D Alembert s Principle, Examples of Generalized Forces a way to deal with friction, and other non-conservative forces D Alembert s principle of virtual
More informationLecture 27: Generalized Coordinates and Lagrange s Equations of Motion
Lecture 27: Generalize Coorinates an Lagrange s Equations of Motion Calculating T an V in terms of generalize coorinates. Example: Penulum attache to a movable support 6 Cartesian Coorinates: (X, Y, Z)
More informationLinear and Nonlinear Properties of Pancharatnam s Phase in Polarization Optics. Hyrum Richardson
Linear and Nonlinear Properties of Pancharatnam s Phase in Polarization Optics Hyrum Richardson A senior thesis submitted to the faculty of Brigham Young University in partial fulfillment of the requirements
More informationGiven P(1,-4,-3), convert to cylindrical and spherical values;
CHAPTER 1 Poblems Pob. 1.1 Pob. 1.2 () Given P(1,-4,-3), convet to cylindicl nd spheicl vlues; 4 x y = + = + = = 1 ( 4) 17 4.123. 1 y 1 4 = tn = tn = 284.04. x 1 P(,, ) = (4.123, 284.04, 3). Spheicl :
More informationQuantum Mechanics for Scientists and Engineers. David Miller
Quantum Mechanics for Scientists and Engineers David Miller Background mathematics Basic mathematical symbols Elementary arithmetic symbols Equals Addition or plus 3 5 Subtraction, minus or less 3 or Multiplication
More informationHeat, Work, and the First Law of Thermodynamics. Chapter 18 of Essential University Physics, Richard Wolfson, 3 rd Edition
Heat, Work, and the First Law of Thermodynamics Chapter 18 of Essential University Physics, Richard Wolfson, 3 rd Edition 1 Different ways to increase the internal energy of system: 2 Joule s apparatus
More informationCMPSCI 250: Introduction to Computation. Lecture #11: Equivalence Relations David Mix Barrington 27 September 2013
CMPSCI 250: Introduction to Computation Lecture #11: Equivalence Relations David Mix Barrington 27 September 2013 Equivalence Relations Definition of Equivalence Relations Two More Examples: Universal
More informationUnit 4 Example Review. 1. Question: a. [A] b. [B] c. [C] d. [D]
Unit 4 Example Review 1. Question: Answer: When the line is wrapped around the circle, the length of the line will be equal to the length of the arc that is formed by the wrap starting from the point (1,0)
More informationPhysics 828 Problem Set 7 Due Wednesday 02/24/2010
Physics 88 Problem Set 7 Due Wednesday /4/ 7)a)Consider the proton to be a uniformly charged sphere of radius f m Determine the correction to the s ground state energy 4 points) This is a standard problem
More informationExamples: Solving nth Order Equations
Atoms L. Euler s Theorem The Atom List First Order. Solve 2y + 5y = 0. Examples: Solving nth Order Equations Second Order. Solve y + 2y + y = 0, y + 3y + 2y = 0 and y + 2y + 5y = 0. Third Order. Solve
More informationProject 10.5C Circular Membrane Vibrations
Projet 1.5C Cirular Membrane Vibrations In problems involving regions that enjoy irular symmetry about the origin in the plane (or the vertial z-axis in spae), the use of polar (or ylindrial) oordinates
More information4 Finite Element Method for Trusses
4 Finite Element Method for Trusses To solve the system of linear equations that arises in IPM, it is necessary to assemble the geometric matrix B a. For the sake of simplicity, the applied force vector
More informationChapter DEs with Discontinuous Force Functions
Chapter 6 6.4 DEs with Discontinuous Force Functions Discontinuous Force Functions Using Laplace Transform, as in 6.2, we solve nonhomogeneous linear second order DEs with constant coefficients. The only
More informationMathematics 22: Lecture 7
Mathematics 22: Lecture 7 Separation of Variables Dan Sloughter Furman University January 15, 2008 Dan Sloughter (Furman University) Mathematics 22: Lecture 7 January 15, 2008 1 / 8 Separable equations
More informationEngineering Mechanics Statics
Mechanical Systems Engineering_2016 Engineering Mechanics Statics 6. Moment of a Couple Dr. Rami Zakaria Moment of a Couple We need a moment (or torque) of (12 N m) to rotate the wheel. Notice that one
More informationPractice Exam 1 (with solutions)
Practice Exam 1 (with solutions) Multiple Choice Identify the choice that best completes the statement or answers the question. 1. The quantity with the same units as force times time, Ft, with dimensions
More informationMassachusetts Institute of Technology Department of Electrical Engineering and Computer Science Electric Machines
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.685 Electric Machines Problem Set solutions September 8, 3 e dλ dt Problem : Torque is, for this system,
More informationGAUTENG DEPARTMENT OF EDUCATION SENIOR SECONDARY INTERVENTION PROGRAMME MATHEMATICS GRADE 12 SESSION 20 (LEARNER NOTES)
MATHEMATICS GRADE SESSION 0 (LEARNER NOTES) TRIGONOMETRY () Learner Note: Trigonometry is an extremely important and large part of Paper. You must ensure that you master all the basic rules and definitions
More information10.4 The Cross Product
Math 172 Chapter 10B notes Page 1 of 9 10.4 The Cross Product The cross product, or vector product, is defined in 3 dimensions only. Let aa = aa 1, aa 2, aa 3 bb = bb 1, bb 2, bb 3 then aa bb = aa 2 bb
More informationEE364 Review Session 4
EE364 Review Session 4 EE364 Review Outline: convex optimization examples solving quasiconvex problems by bisection exercise 4.47 1 Convex optimization problems we have seen linear programming (LP) quadratic
More informationElectromagnetic Properties of Materials Part 2
ECE 5322 21 st Century Electromagnetics Instructor: Office: Phone: E Mail: Dr. Raymond C. Rumpf A 337 (915) 747 6958 rcrumpf@utep.edu Lecture #3 Electromagnetic Properties of Materials Part 2 Nonlinear
More informationThe AGM Function. Kurt Nalty. July 1, 2014
The AGM Function Kurt Nalty July 1, 014 Abstract These are some informal notes on the Arithmetic Geometic Mean function, and the connection to elliptic integrals of the first and second kind. The primary
More informationφ φ sin φ θ sin sin u = φθu
Project 10.5D Spherical Harmoic Waves I problems ivolvig regios that ejoy spherical symmetry about the origi i space, it is appropriate to use spherical coordiates. The 3-dimesioal Laplacia for a fuctio
More informationDevelopment Economics (PhD) Intertemporal Utility Maximiza
Development Economics (PhD) Intertemporal Utility Maximization Department of Economics University of Gothenburg October 7, 2015 1/14 Two Period Utility Maximization Lagrange Multiplier Method Consider
More informationEngineering Mechanics Statics
Mechanical Systems Engineering- 2016 Engineering Mechanics Statics 2. Force Vectors; Operations on Vectors Dr. Rami Zakaria MECHANICS, UNITS, NUMERICAL CALCULATIONS & GENERAL PROCEDURE FOR ANALYSIS Today
More informationLinear least squares problem: Example
Linear least squares problem: Example We want to determine n unknown parameters c,,c n using m measurements where m n Here always denotes the -norm v = (v + + v n) / Example problem Fit the experimental
More informationVolume 37, Issue 3. On the Likelihood of the Coincidence of Condorcet Committees
Volume 37, Issue 3 On the Likelihood of the Coincidence of Condorcet Committees Eric Kamwa Vincent Merlin Université des Antilles, LC2S UMR CNRS 8053, F-97275 CNRS and Normandie Univ, France ; UNAEN, CREM,
More informationFH6 GENERAL CW CW CHEMISTRY # BMC10-4 1BMC BMC10-5 1BMC10-22 FUTURE. 48 x 22 CART 1BMC10-6 2NHT-14,16,18 36 X x 22.
0 //0 :0: TR -00 A 0-W -00 00-W -00 A 00-W R0- W- & W- U RTAY TA- RAT WT AA TRATR R AT ' R0- B U 0-00 R0- B- B0- B0- R0- R0- R0- R0- B0- R0-0 B0- B0- A. B0- B0- UTR TAT 0-00 R0-,, T R-,,0 0-0 A R- R0-,
More informationThe discrete Fourier transform
3 The discrete Fourier transform Page 1 3.1 Introduction Let f(t) be a function describing a physical process in the time domain. An equivalent representation of the process, F (ν), is possible in the
More informationEquilibrium. Forward and Backward Reactions. Hydrogen reacts with iodine to make hydrogen iodide: H 2 (g) + I 2 (g) 2HI(g)
Equilibrium Forward and Backward Reactions Hydrogen reacts with iodine to make hydrogen iodide: H 2 (g) + I 2 (g) 2HI(g) forward rate = k f [H 2 ][I 2 ] 2HI(g) H 2 (g) + I 2 (g) backward rate = k b [HI]
More informationMATH 304 Linear Algebra Lecture 20: Review for Test 1.
MATH 304 Linear Algebra Lecture 20: Review for Test 1. Topics for Test 1 Part I: Elementary linear algebra (Leon 1.1 1.4, 2.1 2.2) Systems of linear equations: elementary operations, Gaussian elimination,
More informationLinear Algebra. Jim Hefferon. x x x 1
Linear Algebra Jim Hefferon ( 1 3) ( 2 1) 1 2 3 1 x 1 ( 1 3 ) ( 2 1) x 1 1 2 x 1 3 1 ( 6 8) ( 2 1) 6 2 8 1 Notation R real numbers N natural numbers: {, 1, 2,...} C complex numbers {......} set of... such
More informationLogic and Modelling. Introduction to Predicate Logic. Jörg Endrullis. VU University Amsterdam
Logic and Modelling Introduction to Predicate Logic Jörg Endrullis VU University Amsterdam Predicate Logic In propositional logic there are: propositional variables p, q, r,... that can be T or F In predicate
More information1 The concept of limits (p.217 p.229, p.242 p.249, p.255 p.256) 1.1 Limits Consider the function determined by the formula 3. x since at this point
MA00 Capter 6 Calculus and Basic Linear Algebra I Limits, Continuity and Differentiability Te concept of its (p.7 p.9, p.4 p.49, p.55 p.56). Limits Consider te function determined by te formula f Note
More informationMath156 Review for Exam 4
Math56 Review for Eam 4. What will be covered in this eam: Representing functions as power series, Taylor and Maclaurin series, calculus with parametric curves, calculus with polar coordinates.. Eam Rules:
More informationFind the equation of a plane perpendicular to the line x = 2t + 1, y = 3t + 4, z = t 1 and passing through the point (2, 1, 3).
CME 100 Midterm Solutions - Fall 004 1 CME 100 - Midterm Solutions - Fall 004 Problem 1 Find the equation of a lane erendicular to the line x = t + 1, y = 3t + 4, z = t 1 and assing through the oint (,
More informationIMPROVEMENT OF AN APPROXIMATE SET OF LATENT ROOTS AND MODAL COLUMNS OF A MATRIX BY METHODS AKIN TO THOSE OF CLASSICAL PERTURBATION THEORY
IMPROVEMENT OF AN APPROXIMATE SET OF LATENT ROOTS AND MODAL COLUMNS OF A MATRIX BY METHODS AKIN TO THOSE OF CLASSICAL PERTURBATION THEORY By H. A. JAHN {University of Birmingham) [Received 7 October 947]
More information144 Chapter 3. Second Order Linear Equations
144 Chapter 3. Second Order Linear Equations PROBLEMS In each of Problems 1 through 8 find the general solution of the given differential equation. 1. y + 2y 3y = 0 2. y + 3y + 2y = 0 3. 6y y y = 0 4.
More informationLesson 14: Solving Inequalities
Hart Interactive Algebra 1 Lesson 14 Classwork 1. Consider the inequality xx 2 + 4xx 5. a. Think about some possible values to assign to xx that make this inequality a true statement. Find at least two
More informationDerivation Of Lagrange's Equation Of Motion For Nonholonomic Constraints Using Lagrange s Multiplier
Derivation Of Lagrange's Equation Of Motion For Nonholonomic Constraints Using Lagrange s Multiplier We consier a path ABC in x, y, x., y.,t space which a system traverse from time to time. Now we consier
More informationECE251DN: Homework #4 Solutions
ECE51DN: Homework #4 Solutions Problem 4.1.5 (a) θ 0 = 0 Beam Pattern: The Kaiser window can be computed by ( ) n W kaiser [n] = I 0 β 1, N 1 n N 1 N Where β = 6, N = 11. For a broadside array, the separable
More informationElizabeth Gire Sanjay Rebello
Elizabeth Gire Sanjay Rebello What factors affect the difficulty of a physics problem? Structure of a Problem/ Type of Problem Textbook problem Ranking task Jeopardy yquestion Conceptual question Multiple
More informationSIMPLIFICATION OF FORCE AND COUPLE SYSTEMS & THEIR FURTHER SIMPLIFICATION
SIMPLIFICATION OF FORCE AND COUPLE SYSTEMS & THEIR FURTHER SIMPLIFICATION Today s Objectives: Students will be able to: a) Determine the effect of moving a force. b) Find an equivalent force-couple system
More informationNARAYANA I I T / N E E T A C A D E M Y. C o m m o n P r a c t i c e T e s t 1 3 XI-IC SPARK Date: PHYSICS CHEMISTRY MATHEMATICS
. (C). (C) 3. (A) 4. (B) 5. (C) 6. (A) 7. (B) 8. (A) 9. (B) 0. (B). (D). (D) 3. (A) 4. (C) 5. (D) NARAYANA I I T / N E E T A C A D E M Y XIS-IC-IIT-SPARK (8..7) C o m m o n P r a c t i c e T e s t 3 XI-IC
More informationMonte Carlo Model of Comet Dust
Monte Carlo Model of Comet Dust Computational Physics Monte Carlo Model of Comet Dust Project 1 About Comets... Coma (Head) Dust Tail Ion Tail Nucleus 10^4-10^5 km 10^5-10^6 km 10^5-10^6 km 1 km All the
More informationMOMENT OF A FORCE SCALAR FORMULATION, CROSS PRODUCT, MOMENT OF A FORCE VECTOR FORMULATION, & PRINCIPLE OF MOMENTS
MOMENT OF A FORCE SCALAR FORMULATION, CROSS PRODUCT, MOMENT OF A FORCE VECTOR FORMULATION, & PRINCIPLE OF MOMENTS Today s Objectives : Students will be able to: a) understand and define moment, and, b)
More informationChapter 3 + some notes on counting the number of degrees of freedom
Chapter 3 + some notes on counting the number of degrees of freedom Minimum number of independent parameters = Some number of dependent parameters minus the number of relationships (equations) you can
More informationEE4440 HW#8 Solution
EE HW#8 Solution April 3,. What frequency must be used for basis symbols φ = sin(ω c t) and φ = cos(ω c t) with a symbol period of /(96Hz) to contain 6 cycles of the sin and cos? Solution: The frequency
More informationDifferential Lengths, Surfaces and Volumes
Differential Lengths, Surfaces and Volumes When integrating along lines, over surfaces, or throughout volumes, the ranges of the respective variables define the limits of the respective integrations. In
More informationModeling the Yildiz Motor revisited
Modeling the Yildiz Motor revisited Duarte, J.L. DOI: 10.6100/6b67487a-bd78-4a24-bdd1-d7b9b3dbd5b7 Published: 22/07/2018 Document Version Other version Please check the document version of this publication:
More informationLie symmetry and Mei conservation law of continuum system
Chin. Phys. B Vol. 20, No. 2 20 020 Lie symmetry an Mei conservation law of continuum system Shi Shen-Yang an Fu Jing-Li Department of Physics, Zhejiang Sci-Tech University, Hangzhou 3008, China Receive
More informationODE Homework Solutions of Linear Homogeneous Equations; the Wronskian
ODE Homework 3 3.. Solutions of Linear Homogeneous Equations; the Wronskian 1. Verify that the functions y 1 (t = e t and y (t = te t are solutions of the differential equation y y + y = 0 Do they constitute
More informationDA 4: Maximum and Minimum Values
Differentiation Applications : Maximum and Minimum Values 7 DA : Maximum and Minimum Values Model : Stage of 0 Tour de France Summary Box DA. For a function f on a domain D, there is a global maximum at
More informationShan Feng
Shan Feng fengs2@mcmaster.caca Shan Feng: 2pm~5pm, Monday, ITB A302; fengs2@univmail.cis.mcmaster.ca Yong Zhao: 1pm~4pm, Wednesday, ITB A103; zhaoy49@grads.ece.mcmaster.ca Heng Wang:2pm~5pm, Thursday,
More information4 4 N v b r t, 20 xpr n f th ll f th p p l t n p pr d. H ndr d nd th nd f t v L th n n f th pr v n f V ln, r dn nd l r thr n nt pr n, h r th ff r d nd
n r t d n 20 20 0 : 0 T P bl D n, l d t z d http:.h th tr t. r pd l 4 4 N v b r t, 20 xpr n f th ll f th p p l t n p pr d. H ndr d nd th nd f t v L th n n f th pr v n f V ln, r dn nd l r thr n nt pr n,
More informationEquations with regular-singular points (Sect. 5.5).
Equations with regular-singular points (Sect. 5.5). Equations with regular-singular points. s: Equations with regular-singular points. Method to find solutions. : Method to find solutions. Recall: The
More information1 Equation set of 7 equations
380(4).wxm 1 / 8 (%i1) kill(all); (%o0) one (%i1) cross(a,b) := [a[2]*b[3] - a[3]*b[2], a[3]*b[1] - a[1]*b[3], a[1]*b[2] - a[2]*b[1]]; (%o1) cross( a,b ):=[a 2 b 3 a 3 b 2,a 3 b 1 a 1 b 3,a 1 b 2 a 2 b
More informationSOUTHWESTERN ELECTRIC POWER COMPANY SCHEDULE H-6.1b NUCLEAR UNIT OUTAGE DATA. For the Test Year Ended March 31, 2009
Schedule H-6.lb SOUTHWSTRN LCTRIC POWR COMPANY SCHDUL H-6.1b NUCLAR UNIT OUTAG DATA For the Test Year nded March 31, 29 This schedule is not applicable to SVvPCO. 5 Schedule H-6.1 c SOUTHWSTRN LCTRIC POWR
More informationA evaluation proposal for the Wuerth centrifugal pumper by W.D. Bauer,
A evaluation proposal for the Wuerth centrifugal pumper by W.D. Bauer, 2.7.02 1. Introduction In a previous article the author described Wuerth s rotor experiment [1]. The result of the experiment was
More informationLecture 7 Open-loop & closedloop experiments. The Bias Formula goes closed-loop and the returns
Lecture 7 Open-loop & closedloop experiments The Bias Formula goes closed-loop and the returns Closed-loop vs open-loop Open-loop system disturbance Feedback-free system input output Control input independent
More information