Scribd
Upload
7 views

01 Introduction To ENGG 501

Uploaded by

smnyqstns
Copyright
© © All Rights Reserved
Available Formats
Download as PDF, TXT or read online on Scribd
7 views

01 Introduction To ENGG 501

Uploaded by

smnyqstns
Copyright
© © All Rights Reserved
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