Metody Numeryczne (III INF)
Egzamin, 3 luty 2015 Uwaga.
Rozwiazanie ka˙zdego z dziesi, eciu zada´, n nale˙zy pisa´c na osobnej kartce podpisanej imie- miem, nazwiskiem i numerem indeksu. Ka˙zde zadanie warte jest 5 punkt´ow.
1. Niech
f (x) =√
2x + 30 −√
2x + 10, x ≥ −5.
Rozpatrzmy dwa algorytmy, A1 i A2, obliczania f (x).
A1: zgodnie z powy˙zszym wzorem, A2: korzystajac z r´, ownowa˙znego wzoru
f (x) = 20
√2x + 30 +√
2x + 10.
Kt´ory z tych algorytm´ow nale˙zy zastosowa´c do obliczenia f (x) w arytmetyce fl po- jedynczej precyzji dla x = 108? Podaj kr´otkie uzasadnienie.
2. Niech
A =
3 2 1
30 22 11 0 20 110
∈ R3,3.
Znajd´z metoda eliminacji Gaussa z wyborem elementu g l´, ownego w kolumnie macierze L, R i P takie, ˙ze P A = LR, gdzie L ∈ R3,3 jest macierza tr´, ojkatn, a doln, a z jedyn-, kami na g l´ownej przekatnej i wszystkimi elementami co do modu lu nie wi, ekszymi od, jedno´sci, R ∈ R3,3 jest tr´ojkatna g´, orna, a P ∈ R3,3 jest macierza permutacji.,
3. Dla danych w tabeli punkt´ow
xi −1 −1 −2 0
yi 1 0 2 0
poszukujemy wielomianu postaci w(x) = a0+ a1x minimalizujacego wielko´s´, c
4
X
i=1
(w(xi) − yi)2.
Sformu luj odpowiednie Liniowe Zadanie Najmniejszych Kwadrat´ow z macierza A i, wektorem ~b, a nastepnie rozwi, a˙z je dowoln, a metod, a.,
4. Niech
p0(x) = 1, p1(x) = x − 1,
p2(x) = (x − 1)(x − 2),
p3(x) = (x − 1)(x − 2)(x − 3),
oraz w(x) = p0(x) + p1(x) − 2p2(x) + 3p3(x). Oblicz r´o˙znice dzielon, a w[2, 3, 4]., Dalszy ciag na nast, epnej stronie!,
5. Niech
f (x) = 32x3− 48x2+ 18x − 1.
Wyznacz wielomian stopnia ≤ 2 najlepiej przybli˙zajacy f na przedziale [0, 1] w nor-, mie jednostajnej, tzn. w normie
kgk∞= sup
0≤t≤1
|g(t)|.
6. Niech {pk}k≥0 bedzie ci, agiem wielomian´, ow ortogonalnych, deg pk = k, wzgledem, iloczynu skalarnego hf, gi = R1
−1f (x)g(x)%(x) dx, gdzie waga % jest parzysta, tzn.
%(x) = %(−x) ∀x. Wyka˙z, ˙ze pk(0) = 0 wtedy i tylko wtedy gdy k jest nieparzyste.
7. Znajd´z maksymalny rzad kwadratury postaci,
Q(f ) = a0f (1/3) + a1f0(1/3) + a2f (c), gdzie a0, a1, a2 ∈ R oraz c ∈ [0, 1], przybli˙zajacej ca lk, e, R1
0 f (x) dx.
8. Do macierzy
A = −11 45 45 −11
zastosowano odwrotna metod, e pot, egow, a,,
~
yk = (A − σI)−1~xk−1,
~
xk = ~yk/k~ykk2, k = 1, 2, 3, . . .
z parametrem σ = 100 oraz wektorem poczatkowym ~, x0 = [10, 1]T. Czy otrzymany ciag jest zbie˙zny, a je´sli tak to jaka jest jego granica? Jak zmieni si, e odpowied´, z gdy σ = −100?
9. Niech
f (x) = (x − 1)2(x − 2)(x − 3)3.
Rozwa˙zmy r´ownanie f (x) = 0. Dla jakich warto´sci parametru τ metoda iteracji prostej dana wzorem xk+1 = ϕ(xk), gdzie
ϕ(x) = x − τ f (x),
jest lokalnie zbie˙zna do rozwiazania α = 2? Odpowied´, z uzasadnij.
10. Stosujemy metode Newtona dla znalezienia zera funkcji, f (x) = (x − 2)2+ ex−2− 2.
Poka˙z, ˙ze r´ownanie f (x) = 0 ma dok ladnie jedno rozwiazanie x, ∗ w przedziale [2, 3].
Czy metoda jest zbie˙zna lokalnie do x∗, a je´sli tak to z jakim rzedem? Czy metoda, jest zbie˙zna dla wszystkich punkt´ow poczatkowych x, 0 ≥ x∗?
KONIEC ZADA ´N