Scribd Logo
Upload

Welcome to Scribd!

  • 26 views

    SI 211: Numerical Analysis HW1 2018231036 顾怡恬

    Uploaded by

    gabrielsi
    This document contains the solutions to homework problems from a numerical analysis class. 1. The first problem examines numerical approximation error when calculating sin(10000π) due to limitations in floating point numbers. It also derives expressions for error propagation in a series of calculations. 2. The second problem provides code to numerically calculate the derivative of a function using a forward difference formula. 3. The third problem derives Taylor series expressions for a function f(x) and its derivative f'(x), then calculates the numerical errors associated with evaluating f'(x) at x=0, finding the error is approximately 0.

    Copyright:

    © All Rights Reserved
    Download as DOCX, PDF, TXT or read online from Scribd

    SI 211: Numerical Analysis HW1 2018231036 顾怡恬

    Uploaded by

    gabrielsi
    26 views2 pages
    This document contains the solutions to homework problems from a numerical analysis class. 1. The first problem examines numerical approximation error when calculating sin(10000π) due to limitations in floating point numbers. It also derives expressions for error propagation in a series of calculations. 2. The second problem provides code to numerically calculate the derivative of a function using a forward difference formula. 3. The third problem derives Taylor series expressions for a function f(x) and its derivative f'(x), then calculates the numerical errors associated with evaluating f'(x) at x=0, finding the error is approximately 0.

    Original Description:

    shanghaitech SI211 HW1

    Original Title

    SI211HW1

    Copyright

    © © All Rights Reserved

    Available Formats

    DOCX, PDF, TXT or read online from Scribd

    Did you find this document useful?

    Is this content inappropriate?

    This document contains the solutions to homework problems from a numerical analysis class. 1. The first problem examines numerical approximation error when calculating sin(10000π) due to limitations in floating point numbers. It also derives expressions for error propagation in a series of calculations. 2. The second problem provides code to numerically calculate the derivative of a function using a forward difference formula. 3. The third problem derives Taylor series expressions for a function f(x) and its derivative f'(x), then calculates the numerical errors associated with evaluating f'(x) at x=0, finding the error is approximately 0.

    Copyright:

    © All Rights Reserved
    Download as DOCX, PDF, TXT or read online from Scribd
    Download as docx, pdf, or txt
    26 views2 pages

    SI 211: Numerical Analysis HW1 2018231036 顾怡恬

    Uploaded by

    gabrielsi
    This document contains the solutions to homework problems from a numerical analysis class. 1. The first problem examines numerical approximation error when calculating sin(10000π) due to limitations in floating point numbers. It also derives expressions for error propagation in a series of calculations. 2. The second problem provides code to numerically calculate the derivative of a function using a forward difference formula. 3. The third problem derives Taylor series expressions for a function f(x) and its derivative f'(x), then calculates the numerical errors associated with evaluating f'(x) at x=0, finding the error is approximately 0.

    Copyright:

    © All Rights Reserved
    Download as DOCX, PDF, TXT or read online from Scribd
    Download as docx, pdf, or txt
    You are on page 1of 2

    SI 211: Numerical Analysis HW1 2018231036 顾怡恬

    1. (a) Numerical approximation error is -1.5460e-13. This error comes from the
    limitation of the floating numbers. sin(10000π) cannot be exactly 0 because
    π is an infinite non-repeating decimal.
    (b)
    a0 = 𝑥 ∆a0 = 𝑒𝑝𝑠 ∆a0 |𝑥=10−10 ≈ 𝑒𝑝𝑠
    4 4
    a1 = 10 ∗ 𝑎0 ∆a1 = 10 ∗ ∆a0 + 𝑒𝑝𝑠 ∆a1 |𝑥=10−10 ≈ 1.0001 × 104 𝑒𝑝𝑠
    a2 = sin⁡(a1 ) ∆a2 = |cos(𝑎1 )| ∗ ∆a1 + 𝑒𝑝𝑠 ∆a2 |𝑥=10−10 ≈ 1.0002 × 104 𝑒𝑝𝑠
    𝑎2 𝑎2 ∆a2
    a3 = ∆a3 = |− ∗ ∆a0 | + | | + 𝑒𝑝𝑠 ∆a3 |𝑥=10−10 ≈ 2.0002 × 1014 𝑒𝑝𝑠
    𝑎0 𝑎02 𝑎0

    2. (a) The code is followed.


    function⁡df=diff(f1,x1,h)
    syms⁡x;
    f(x)=f1;
    f1=double(f(x1+h));
    f2=double(f(x1-h));
    df=(f1-f2)/(2*h);
    end
    (b)

    3. (a)
    a0 = 𝑥 b0 = 1
    a1 = cos⁡(𝑥) b1 = −sin⁡(𝑥)
    a2 = 𝑎1 ∗ 𝑎1 b2 = 𝑎1 ∗ 𝑏1 + 𝑏1 ∗ 𝑎1
    a3 = sin⁡(𝑎1 ) b3 = cos(𝑎1 ) ∗ 𝑏1
    a4 = 𝑎2 ∗ 𝑎3 b4 = 𝑏2 ∗ 𝑎3 + 𝑎2 ∗ 𝑏3
    𝑓(𝑥) = a4 f ′ (𝑥) = 𝑏4
    (b)
    a0 = 0 ∆a0 = 𝑒𝑝𝑠(0) = 0
    a1 = 1 ∆a1 = − sin(0) + 𝑒𝑝𝑠 ≈ 2.0 ∗ 10−16
    a2 = 1 ∆a2 = 𝑎1 ∗ ∆a1 + 𝑏1 ∗ ∆a1 + 𝑒𝑝𝑠 ≈ 6.0 ∗ 10−16
    SI 211: Numerical Analysis HW1 2018231036 顾怡恬

    a3 = 0.841470985 ∆a3 = cos(𝑎1 ) ∗ ∆a1 + 𝑒𝑝𝑠 ≈ 3.08060461 ∗ 10−16


    a4 = 0.841470985 ∆a4 = ∆a2 ∗ 𝑎3 + 𝑎2 ∗ ∆a3 ≈ 8.12943052 ∗ 10−16

    b0 = 1 ∆b0 = 𝑒𝑝𝑠(0) = 0
    b1 = 0 ∆b1 = |−cos⁡(𝑎0 )| ∗ ∆a0 + 𝑒𝑝𝑠(0) = 0
    b2 = 2 ∗ 𝑎1 ∗ 𝑏1 = 0 ∆b2 = 2 ∗ |𝑏1 | ∗ ∆a1 + 2 ∗ |𝑎1 | ∗ ∆b1 + 𝑒𝑝𝑠(0) = 0
    b3 = cos(𝑎1 ) ∗ 𝑏1 ∆b3 = |cos⁡(𝑎1 )| ∗ ∆b1 + |𝑏1 | ∗ |sin(𝑎1 )| ∗ ∆a1 + 𝑒𝑝𝑠(0) = 0
    b4 = 𝑏2 ∗ 𝑎3 + 𝑎2 ∗ 𝑏3 ∆b4 = |𝑏3 | ∗ ∆a2 + |𝑎2 | ∗ ∆b3 + |𝑏2 | ∗ ∆a3 + |𝑎3 | ∗ ∆b2 + 𝑒𝑝𝑠(0) =
    0
    So the numerical error that is associated with evaluating the derivative is 0.

    You might also like