билеты по матмоду (1239874), страница 2
Текст из файла (страница 2)
Прямой ход:предп., что в СЛАУ а11≠0; m - 1 шагов. k – ый шаг: мастаб коэфф.: mik = aik(k-1)/akk(k-1), i = k + 1, … m. aij(k) = aij(k-1) - mikakj(k-1), i = k + 1, … m; j = k, … m. bi(k) = bi(k-1) - mikbk(k-1). Обратный ход: xm = bm(m-1)/amm(m-1). xk = (1/akk)[bk(k-1) - Sj=k+1makjxi(k-1)]. Труд-ть м: 1 шаг: m – 1 дел-е, (m – 1)2 умн-е, (m – 1)2 выч-е. 2 шаг: m – 2 дел-я, (m – 2)2 умн-е, (m – 2)2 выч-е. 3. m – 1 шаг: все по 1 разу. Итого Sk=1m-1k + 2Sk=1m-1k2 » 2m3/3 арифм. действий.
15. LU-разл-е м-цы. Дост-е усл-е разл-ти м-цы.
Это тот же м.Г., но для самой м-цы А, не трогая b. A = LU, L – НТМ, сост. из масштаб. коэфф. постр: L = (1, 0, … , 0;m21, 1, 0, …, 0;… ; mm1, mm2, … , mm,m-1, 1). U – ВТМ, то, что ост. от А в проц. преобраз. Постр.: U = (a11, a12, …, a1m; 0, a21(1), …, a2m(1); 0, 0, …, amm(m-1)). Р-е: Ly = b ® y, Ux = y ® x – ответ. Дост.усл.разл-ти: Если все гл. мин. матр. А ¹ 0, то $ ! LU-разл. м
16. Зад., реш-е на осн. LU: A-1, det(A).
Примен. для реш. систем с b, т. к. LU – разл. требует (2/3)m3 арифм. действ., а ост-е д-я лишь 2m2. detA = (-1)sa11(0)a22(1)…amm(m-1), где s – перест.строк. A A-1=E, реш. как AX=E
17. МГ(сх. част.выб.): опис-е м., выч.уст-ть.
Усл-е:|mik|£ 1(выч.уст-ть). на каж.(k-ом) ш. м. выб. макс. по мод.коэфф. из стб, и мен.стр местами, макс. кофф. аij попала на 1место. И.т.д. Это чтобы уменьш. mik, погр-ть нар. <
18. МХ(СЛАУ): опис-е м., его пр-ва.
Для симм. пол-но опр-х м-ц. Положит. опр. – значит все собств. числа > 0, в нек-ром. случ. достат. проверять, что |aii|£ Si¹j|aij|. А = LLT. Матрица L: L = (l11, 0, …, 0; l21, l22, 0, …, 0; lm1, lm2, …, lmm). Эл-ты матр. LLT =матр. А. Ly = b ® y, LTx = y ® x. Преимущ.: труд-ть » m3/3 + 2m2 , в 2р.<МГ. Симм-я матр. позв. ком-но хр-ть и зам-ть эл-ты аij эл-ми lij, м. гар-но уст-в.
19. М.прог-ки с 3-диаг.м-й: оп-е м., усл-я его прим-ти и дост-ва.
построчно: (b1, c1, 0, …, 0; a2, b2, c2, 0, …, 0; 0, a3, b3, c3, 0, …, 0; 0, …, 0, am, bm). Пр.ч.: d1, d2, …, dm. Пр. ход – выч-е. прог-х. коэфф.: При i = 1:a1 = - с1/b1;b1 = d1/b1. При i = 2, …, m-1; gi = (aiai-1 + bi); ai = - ci/gi; bi = (di - aibi-1)/gi. Обратный ход – вычисл. х: xm = bm = (dm – ambm-1)/(amam-1 + bm); далее: xi = aixi+1 + bi, при i = m-1, m-2, …, 1. Труд-ть » 8m арифм. д-й. Усл-я прим-ти: Если bk|³|ak|+|ck|;|bk|>|ak|. Тогда, gk ¹ 0, т. е. пр-ка м.б. довед. до конца и |ak|£1 – т.е. погр-ть не нак-ся.
20. МЯ (СЛАУ): оп-ие м., ус-е сх-ти, оц.погр.
Ах = b ® x = Bх + c, причем, bij = - aij/aii; ci = bi/aii. Т. Усл-е сх-ти: Пусть Ax = b ® x = Bx + c (1); x(n+1) = Bx(n) + c (2). Если ||B|| < 1, то итер. процесс сход. к корню ур-я и справедл. след. оценки: ||x(n) - || £ ||B||n * ||x(0) -
|| (3);||x(n) -
|| £ (||B||/(1 - ||B||))*||x(n) – x(n-1)|| (4). Д.(3). Имеем системы: x(n+1) = Bx(n) + c и
= B
+ c, вычтем их: x(n+1) -
= Bx(n) - B
= В(x(n) -
). ||x(n+1) -
|| = ||В(-
)|| £ ||В|||| (x(n) -
)|| £ ||В2|||| (x(n-1) -
)|| £||В||n+1 * ||x(n-1) -
||. Д.(4). x(n+1) -
= В(x(n-1) -
) = В(x(n-1) - x(n) + x(n) -
) ® ||x(n+1) -
||=||В(x(n) -
)+ В(x(n-1) - x(n))|| £ ||В||||(x(n) -
)|| + ||В||*||(x(n-1) - x(n))||. (1 - ||В||)||(x(n) -
)|| £ ||В||||(x(n-1) - x(n))|| ч.т.д. ->Крит-й ок-я ит-ий: ||(x(n-1) - x(n))|| < e(1 - ||В||)/||B||. Апостер.оц.погр-ти:
21. МЗ(СЛАУ): оп-е м., усл-е сх-ти, оц. погр.
Найд-е на n-ом ш. зн-я неиз-х x1(n), x2(n) … сразу же вкл. в расчет посл-х неизв-х. x(n+1) = B1x(n+1) + B2x(n) + c. В1 – НТМ, на гл. диаг – 0. В2 – ВТМ, на гл. диаг – 0. В = В1 + В2. Коэфф. по тем же ф-лам, что и в МПИ: bij = - aij/aii; ci = bi/aii. Расч-я ф-ла: xi(n+1) = Si=1i-1bijxj(n+1) + Sj=i+1mbijxj(n) + c1, i = 1, …, m.
22. М.рел-ции: оп-е м., усл-е сх-ти.
Модиф. МПИ. как и в МЗ, наход-е Х исп-ся на этом же шаге,:xn=(1-ω)xn-1+ω(B1xn+B2xn-1+c). ω=0..2–коэфф.рел-ции. После выч-я i-коорд.: (ω-1)(xin-xin-1) т. Дост.ул.сх-ти.: A=Aт>0 при ω[0,2] м.р. сход. Ɐx0.
23. Аппр-я ф-й по МНК. Пост. Зад.. Вывод норм. сист. м. Выбор ст.ап-го мн-а.
В случае, если: 1. f(x) зад.табл. 2. f(x) – вычисл. сложно. Задача: прибл-ть ф-ю f(x) ф-ей F(x), т. е. f(x) F(x). - средеквадр. отклон. должно быть миним. = ((1/(n + 1))i=0n(yi – F(xi))2) min. Прибл-ем мн-ом Рn(x) = a0 + a1x + a2x2 + … + amxm = 0. = ((1/(n + 1))i=0n(a0 + a1xi + a2xi2 + … + amxim - yi)2) min. i=0n(a0 + a1xi + a2xi2 + … + amxim - yi)2 min по a0, a1, a2, …, am, т. е. g/a0 = 2i=0n(a0 + a1xi + a2xi2 + … + amxim - yi)*1) = 0; … g/am = 2i=0n(a0 + a1xi + a2xi2 + … + amxim - yi)*xim) = 0; Получаем нормальную систему ур-ий МНК: (i=0nxi0)a0 + (i=0nxi1)a1 + … + (i=0nxim)am = i=0nyi; … ; (i=0nxim)a0 + (i=0nxim+1)a1 + … + (i=0nxi2m)am = i=0nyixim. Система симм. и плохообусл. при m > 5 велика погр-ть. Из нее нах. a0, a1, a2, …, am. Степень аппр. мн. m выб. такой, когда ф-ия (m) имеет 1й лок. Мин.
24. Пост. зад. интерп-ии. Т. о ⱻ! интерп-го мн-на.Дана n + 1 точка, постр. интерп. мн. ст. n, так, чтобы Pn(xi) = yi. Т. интерп. мн-н и !. Д. 1.ⱻ:. Рассм. м-н Лагранжа: Ln(x) = i=0nyi((x – x0)(x-x1)…(x-xi-1)(x-xi+1)…(x-xn)/(xi-x0)(xi-x1)…(xi-xi-1)(xi-xi+1)…(xi-xn)). Ln(x) - мн-ов n-ой степени, сам он тоже n-ой степени. Ln(xk) = yk, k = 0, 1, …, n. 2. !. Рассм. Pn1, Pn2 – 2 интерп. мн-на n – ой ст. Pn(x) = Pn1 - Pn2. Pn(xk) = 0; k = 0, …, n
25. Мн-н Лагранжа. Оц. погр-ти интерп-ии.
Ln(x) = i=0nyi((x – x0)(x-x1)…(x-xi-1)(x-xi+1)…(x-xn)/(xi-x0)(xi-x1)…(xi-xi-1)(xi-xi+1)…(xi-xn)). Ln(x) - мн-ов n-ой степени, сам он тоже n-ой степени.
Погрешность интерп.: Rn(x)=y(x)–Pn(x). T. y(x)cn+1[x0,xn]. Rn(x)(Mn+1/(n+1))!n+1(x), где: Mn+1=max[x0, xn]y(n + 1)(x); n+1(x) = (x – x0)(x - x1)…(x – xn). без. док. maxn+1(x) невелико внутри [x0, xn], а вне быстро растет.
26. Мн-н Ньютона с к.р. Оц. погр-ти интерп.
Прим., когда сетка равном. Сост. диаг. Табл. К.р. строка заголов.: xi; yi; yi; 2yi;… Возьмем интерп. мн. в след. виде: n(x) = a0 + a1(x – x0) + a2(x – x0)(x – x1) + … + an(x – x0)(x – x1)…(x – xn). Условие: n(xi) = yi. (xk – xk-1) = h. Рассм. n(x0) = a0 = y0; n(x1) = … и выв. ф-лу: ak = ky0/(k!hk). Интерп. м-н Н. при интерп. вперед от точки x0: n(x) = y0 + (y0/1!h)(x – x0) + … + (ny0/n!hn)(x – x0)(x – x1)…(x – xn). Интер. МН при инт. назад от точки xn: n(x) = yn + (yn-1/1!h)(x – x0) + … + (ny0/n!hn)(x – x0)(x – x1)…(x – xn). Оц. погр.: как для Лагранжа Rn(x)(Mn+1/(n+1)!)hn+1. Преимущ.: при доб. новых узлов надо только доб.новое слаг-е, можно строить от точки табл.
27. Интерп. с кр. узлами. Мн-н Эрмита.
Изв. в хi не только зн-я f(xi), но и f '(xi), … f(k-1)(xi). (xi – узел кр-ти k), Ɐ(i = 0, …, m). k1 – кр-ть 1-ого узла, …, km – кр-ть m-ого узла. n = k0 + … + km. ! Мн-н ст. n, такой, что Pn(xi) = yi; P'(xi) = y'(x), … до (k – 1) для всех i = 0, …, m. 1. Задан 1уз. Кр-ти k. Pn(x) = i=0k-1y0(i)(x – x0)i/(i!). 2. Зад. неск. точек кр-ти 2, т. е. мы знаем y(xi) и y'(xi). По кажд. 2 сос. Т-м в этом случае стр мн-н Эрмита: P3(x) = yi-1((x – xi)2(2(x – xi-1) + hi)/hi3) + yi((x – xi-1)2(2(xi – x) + hi)/hi3) + y'i-1((x – xi-1)(x – xi)2/hi2) + y'i((x – xi-1)2(x – xi)/hi2). (лок. сплайн). Погр-ть max[x0, x1]f(x) – P3(x) (M4/384)h4, M4 = max[x0, x1]f(4)(x). h = (x1 – x0).
28. Понятие сплайна. Построение линейного сплайна.
Сплайн – ф-я, удовл.след.усл.: 1S(x)предст.на отр. [a,b]нек-рый мн-н ст.n; 2) S(xi)=yi, i=0…n; 3)S`(x)непр. до р-го пор-вкл. m-p – дефект с. Лин. сплайн: yi+((yi+1–yi)/(xi+1–xi))(x–xi), x [xi, xi+1]. Дефект =1.
29 Постр-е куб. сплайна. Разл-е виды гр-х усл-й. Оц. погр-ти прибл-я ф-ции куб.с.
Куб. сплайн.: S3(x)=P3,i(x)=yi-1((x–xi)2(2(x–xi-1)+hi)/hi3)+yi((x–xi-1)2(2(xi–x)+ hi)/hi3) + si-1((x – xi-1)(x - xi)2/hi2) + s'i((x – xi-1)2(x – xi)/hi2). Необходимо, чтобы P''3,i(xi) = P''3,i+1(xi), i = 1, 2, …, n – 1. P''3,i(xi) = 2si-1/hi + 4si/hi – 6(yi – yi-1)/hi2; P''3,i+1(xi) = - 4si/hi+1 – 2si+1/hi+1 + 6(yi+1 – yi)/hi+12. Эти рав-ва приводят к системе ур-ий отн. si: hi-1si-1+2(hi-1+hi+1-1)si+hi+1-1si+1 = 3[hi-2(yi – yi-1) + hi+1-2(yi+1 – yi)]. Доп-яя усл-ми s0 = y'(x0); sn = y'(xn) получ. Фунд-ый куб. с. Погр-ть куб. с.: y(x) – S3(x) cM4hmax4, M4 = max[x0, xn]y(4)(x). Гр-е усл-я: 1. Изв. y'(x0), y'(xn), то пол-т, что s0 = y'(x0); sn=y'(xn), где si – наклон с., si = S'm(xi). 2. Изв. y''(x0), y''(xn) - 4s0/h1 – 2s1/h1 + 6(y1 – y0)/h12 = y''(x0); 2sn-1/hn + 4sn/hn – 6(yn – yn-1)/hn2 = y''(xn). 3. Естеств. С.: y'(x0) = y'(xn) = 0. 4. Отс-е узла. выбор наклонов si произв-т так, чтобы P3,1(x)=P3,2(x); P3,n-1(x)=P3,n(x). для этого потреб. совпад. соотв. 3-их произв.: P'''3,1(x1) = P'''3,2(x1); P'''3,n-1(xn-1) = P'''3,n(xn-1). 2h1-3(y0 – y1) + h1-2(s0 – s1) = 2h2-3(y1 – y2) + h2-2(s1 + s2), 2hn-1-3(yn-2 – yn-1) + hn-1-2(sn-2 + sn-1) = 2hn-3(yn-1 – yn) + hn-2(sn-1 + sn).
30. Реш-е сист. нел-х ур-й м. Ньютона.