Scribd
Upload
16 views12 pages

01 Introduction To ENGG 501

Uploaded by

smnyqstns
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
16 views12 pages

01 Introduction To ENGG 501

Uploaded by

smnyqstns
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 12

01.

Introduction to
Computational Mathematics 1
ENGG 501
Class Professor: Dr. Oliver Lexter July A. Jose
Schedule: Sunday, 7:00 – 10:00 a.m.
Credit: 3 Units
Course Outline
• Introduction Modeling, Computer and Error
Analysis
• Roots of Equation
• Linear Algebraic Equation
• Optimization
• Curve Fitting
• Numerical Differentiation/Integration
• Ordinary DE
Strategy and Final requirement
• In lectures, the instructor will discuss in English only.
• Students (by pair) need to complete one final report.
• The final project is the application of numerical methods to some
existing real-world problem/s (preferably relevant to your thesis or
dissertation).
• The final report should be in:
– IEEE conference format (4 to 5 pages)
– Similarity report, less than 25 %
– Grammarly, 90% above
– Codes using Phyton, Matlab, or any similar software (separate from
the conference paper)
• Grading
– Not submitted 1.25 – 2.0
– Submitted 1.25 – 1.75
– Submitted and accepted 1.0 – 1.5
Simple Mathematical Model
• equation that expresses the essential
features of a physical system in
mathematical terms
𝑑𝑒𝑝𝑣𝑎𝑟 = 𝑓(𝑖𝑛𝑑𝑒𝑝𝑣𝑎𝑟, 𝑝𝑎𝑟𝑎𝑚, 𝑓𝑜𝑟𝑓𝑢𝑛𝑐)
where
𝑑𝑒𝑝𝑒𝑛𝑑𝑒𝑛𝑡 𝑣𝑎𝑟𝑖𝑎𝑏𝑙𝑒 − 𝑏𝑒ℎ𝑎𝑣𝑖𝑜𝑟 𝑜𝑓 𝑡ℎ𝑒 𝑠𝑦𝑠𝑡𝑒𝑚
𝑖𝑛𝑑𝑒𝑝𝑒𝑛𝑑𝑒𝑛𝑡 𝑣𝑎𝑟𝑖𝑎𝑏𝑙𝑒 − 𝑑𝑖𝑚𝑒𝑛𝑠𝑖𝑜𝑛𝑠
𝑝𝑎𝑟𝑎𝑚𝑒𝑡𝑒𝑟𝑠 − 𝑝𝑟𝑜𝑝𝑒𝑟𝑡𝑖𝑒𝑠
𝑓𝑜𝑟𝑐𝑖𝑛𝑔 𝑓𝑢𝑛𝑐𝑡𝑖𝑜𝑛𝑠 − 𝑒𝑥𝑡𝑒𝑟𝑛𝑎𝑙
• From Newton 2nd law
𝐹
𝑎=
𝑚
• Expressing the 𝑎 as the time rate of change
𝑑𝑣 𝐹
=
𝑑𝑡 𝑚
𝐹 = 𝐹𝐷 𝑑𝑜𝑤𝑛𝑤𝑎𝑟𝑑 + 𝐹𝑈 𝑢𝑝𝑤𝑎𝑟𝑑
𝑑𝑣 𝑚𝑔 − 𝑐𝑑 𝑣 2 𝑐𝑑 𝑣 2
= =𝑔−
𝑑𝑡 𝑚 𝑚
• For analytical solution, initially at rest (𝑣 =
0 @ 𝑡 = 0)
𝑔𝑚 𝑔𝑚
𝑣 𝑡 = tanh 𝑡
𝑐𝑑 𝑐𝑑

Where
𝑒 𝑥 − 𝑒 −𝑥
𝑡𝑎𝑛ℎ𝑥 = 𝑥
𝑒 + 𝑒 −𝑥
Example 1
• A bungee jumper with a mass of 68.1 kg
leaps from a stationary hot air balloon (use
drag coefficient of 0.25 kg/m).
– compute velocity for the first 12 (increment of
2) s of free fall.
– determine the terminal velocity that will be
attained for an infinitely long cord
Numerical Methods
• Analytical method satisfies the original DE
• Numerical methods solved by arithmetic
operations by reformulating mathematical
problems
𝑑𝑣 Δ𝑣 𝑣 𝑡𝑖+1 − 𝑣 𝑡𝑖
≅ =
𝑑𝑡 Δ𝑡 𝑡𝑖+1 − 𝑡𝑖
𝑑𝑣 Δ𝑣
≅ lim
𝑑𝑡 Δ𝑡→0 Δ𝑡
𝑣 𝑡𝑖+1 − 𝑣 𝑡𝑖 𝑐𝑑 2
=𝑔− 𝑣
𝑡𝑖+1 − 𝑡𝑖 𝑚
Manipulating the eqns
𝑐𝑑
𝑣 𝑡𝑖+1 = 𝑣 𝑡𝑖 + 𝑔 − 𝑣 𝑡𝑖 2 (𝑡𝑖+1 − 𝑡𝑖 )
𝑚
𝑑𝑣𝑖
𝑣𝑖+1 = 𝑣𝑖 +
𝑑𝑡
Euler’s Method
𝑉𝑎𝑙𝑢𝑒𝑛𝑒𝑤 = 𝑉𝑎𝑙𝑢𝑒𝑜𝑙𝑑 + (𝑠𝑙𝑜𝑝𝑒 × 𝑠𝑡𝑒𝑝 𝑠𝑖𝑧𝑒)
Example 2
• Perform the same computation as in
Example 1 using Euler’s method
Step size = 2, 1, and 0.5 VS actual
Thank you for your attention

You might also like

pFad - Phonifier reborn

Pfad - The Proxy pFad of © 2024 Garber Painting. All rights reserved.

Note: This service is not intended for secure transactions such as banking, social media, email, or purchasing. Use at your own risk. We assume no liability whatsoever for broken pages.


Alternative Proxies:

Alternative Proxy

pFad Proxy

pFad v3 Proxy

pFad v4 Proxy