2008

1.Метод простой итерации.

Решим систему Ах=В (1), А= х= В= пусть аij ≠ 0 тогда Х₁=β₁+£₁₂Х₂+£₁₃Х₃+…+£_1nХ_n. Хn= βn+£_n1Х₁+£_n2Х₂+ …+£nn-1An-1, где βi=вi/Aii £ij=Aij/Aii Х=β+£х – эту систему будем наз. методом последовательных приближений. За 0-е приближение можно взять столбец свобод. членов Х⁰=β, тогда Х¹=β+£, Х¹ – 1-ое прибл. Если посл-ть х0,х1… имеет предел, то этим пределом будет решение сис. т.е. Х^k=∑(i=1 до n) £*Xi^(K-1) +βi. Дост. услов. сходим. 1) ∑(от i=1 до n)|Aij|<1 для люб.Y 2) ∑(от j=1 до n) |Aij|<1 для любого i Следствие: для сист. метод итер. сходит. если все модули диагонал. коэффициентов для каждого Ур-ия >, суммы модулей всех остал. коэффиц.|Aii|<∑(i ≠ j)|Aij|.

2.Метод Зейделя.

Основная идея: при вычислении (К+1)-го приближения неизвестной Хi учитываются уже вычисленные ранее (К+1)-е приближение х1,х2,… Пусть дана приведенная лин. система Хi=βi+(∑ от j=1 до n)£ij*xi^K Предпол. что К-е приближ. известны тогда (К+1) имеет вид Xn^(K+1)=βn+(∑ отj=1 до n-1) £niXj^(K+1) + £nnXn^(K). Достаточное условие: 1) ∑(от i=1 до n)|Aij|<1 для любого y 2) ∑(от j=1 до n)|Aij|<1 для любого i Следствие: для сист. метод итер. сходит. если все модули диагонал. коэффициентов для каждого Ур-ия >, суммы модулей всех остал. коэффиц. |Aii| <(i≠j)|Aij|.

3.Достаточные условия сходимости процесса итерации (доказать)

Теорема: процесс итерации для системы Х=£x+β сх. к единств. решению если какая либо норма матр. £<1 или ||£||<1, x^K=£x^(K-1)+β,

Док-во: постр. последов. приближен., таких, что X^K=(£^(K-1) +£^(K-2)+…+ £+E)β+£^KX⁰ Т.к. норма матр.<1, то норма £к→0, при к→∞ E+£+£²+…+£^K-1=(∑от n=0 до K-1) £n| ||£||<1 |=(E-£)^-1 X=lim(k→∞)X^к=lim[ (£^(K-1) +£^(K-2)+…+ £+E)β+£^KX⁰ ]=(E-£ )^-1β Сходимость доказана.

4. Отделение корней уравнения.

Отделение корней может происходить аналитически или графически.

F(x)=0 X^*- решение( F(X^*) =0). Т. Если непрерывная ф-ция F(x) принимает значения различных знаков на [a,b], т.е F(a)*F(b)<0,то внутри этого отрезка содержится по крайней мере один корень уравнения F(x)=0

Метод половинного деления.

Будем решать ур-ние f(x)=0 на [a,b], известно, что f(a)*f(b)<0. Найти X такое [x-x*] < ε. Метод состоит в последовательном построении интервалов [an, bn], вложенных друг в друга и содерж. реш. X*. Пусть a₀=a, b₀=b предпол. что [ai,bi]промежуток построения, причем f(ai)*f(bi)<0 т.е. X* Є[ai,bi] найдем (.)Ci=(ai+bi)/2. Если f(ci)=0 то (ci) точное решение, если f(ci) ≠ 0 то: либо 1) f(ai)*f(сi) < 0, тогда a(i+1)=ai*b(i+1)=ci, либо 2) f(ci)*f(bi)<0, тогда то a(i+1)=ci, b(i+1)=bi В любом случае получ. интервал вдвое меньше исходного, причем f(ai+1)*f(bi+1) < 0 Процесс заканчив. либо 1)|ai-bi|<E,тогда (F(ci)=0),ci- –выбир. любую (.) и приним. ее за решен. ОЦЕНКА точности. f(ai)*f(bi)<0; bn-an=1/2ⁿ X(с чертой)-an<=1/2ⁿ (b-a) k-кол-во шагов k=log₂((b-a)/E)+1. Решение нелин. ур-ий . Решен. осущ. в 2 этапа 1) локализация корней , т.е. нахожд. промежутков [a,b] котор. принадл. корень ур-ия 2)уточнение корней ,т.е. решение с задан. точностью.

Метод хорд.

Пусть дана ф-ция F(x) = 0 на [a, b], причем F(a)*F(b) < 0. Для определенности предпол., что F(a)<0 а F(b)>0, тогда поделим отрезок [a,b] на F(a)/F(b). (x-a)/(b-a)=(y-f(a))/(f(b)-f(a)); Допустим X=X₁, y=0; X₁=a – (f(a) / (f(b)-f(a)))*(b-a) по этой ф-ле можно записать итерац. процесс X₁=a+h, где h₁=-f(a)/(f(b)-f(a)) Докаж. сходим. итер. процесса будем предполагать, что корень f(x) определен и сохр. знак. на [a,b]. Имеем 2 сл. f(a)>0 и f(a)<0

1)a>0 ;

2) a<0

Т.е.1) неподвиж. тот конец , для кот. знак ф-ции совпад. со знаком втор. производ. 2)Послед. приближен. Xn лежат по ту сторон. X*, где f(x) имеет знак противоп знаку ее 2-ой произв. Критерий остановки |X_n+1 - X_n|< ε

5. Метод касательных (Ньютона).

y-y0=f ^‘ (X)(X-X0) пусть f ^‘ (X), f ^{‘ ‘} (X) –непрер. и сохр. опред. знаки на отрезке [a,b] пусть f ^{‘ ‘} (X)>0 и f(b)>0 (рис.) Ур-ие касат. y-f(Xn)= f ^‘ (Xn)(X-Xn) y=0, X=Xn+1; X(n+1)=Xn-(f(Xn)/ f ^‘ (Xn)) Исходя из геом. соображен. , вывед. тоже самое аналитич. X*=Xn+h_n Прим. ф-лу Тэйлора 0=f(Xn+ h_n)=f(Xn)+ h_nf ^‘ (Xn); h_n =-f(Xn)/f ^‘ (Xn) Выбор начал. точки f(X0)*f ^{‘ ‘} (X0)>0. Теорема: если f(a)*f(b)<0 причем втор. произв. сохр. знак на [a,b] , то исходя из любого начал. приближен. X0c[a,b] удовл. условию f(X0)* f ^{‘ ‘} (X0)>0 можно вычисл. корень ур-ия с любой степенью точности. Д-во: пусть f(a)<0 , f(b)>0 f ^‘ (x)>0; f ^{‘ ‘} (x)>0(вогнута) т.е. f(b)* f ^{‘ ‘} (b)>0 ; X0=b-max точка.

X(n+1)=Xn-(f(Xn)/ f ^‘ (Xn))>X*.

Геометрически прямая параллельная первой касательной . т.е. здесь касательная заменяется на прямую, паралл. первой касательной

Теорема о методе касательных.

Если ф-ия f(a)*f(b)>0, причем вторая призводная сохраняет знак на отрезке [a,b], то исходя из любого начального приближения х₀Є[a,b], удовл. условию f(х₀) f”( х₀)>0 то можно вычислить корень ур-я с любой степенью точности. Док-во: f(a)<0; f(b)>0; f”(x)>0; f’(x)>0, т.е. f(b) f”(b)>0 x₀=b. Методом индукции докажем что, все x_n>x*. Т.е. f(x_n)>0, очевидно х₀>x*. Путь х_n>х*. Представим х*=х_n+(х*-х_n) по форм. Тейлора f(x*)=f(x_n)+f’(x_n) (х*-х_n)+1/2 f”(x_n) (х*-х_n)², тогда f(x_n)-f’(x_n) (х*-х_n)<0 х_n+1=x_n(+)-f(x_n)/f’(x_n)>x*. Ч.т.д.

Видоизмененный метод касательных.

X(n+1)=Xn-f(Xn)/ (Xn). Если производная исходной ф-ции меняется мало на [a,b] то можно предположить f(Xn)≈f(X0) X(n+1)=Xn-f(Xn)/ (X0) Это дает воз-ть не вычисл. зн-е произв. в каждой (.).

Комбинированный метод.

f(x)=0 : (a,b) : f(a)*f(b)<0 имеем 4 случая:

Расм. только 1-ой случай т.к. оставшиеся случаи либо сводяться к нему либо аналогичны. _

X* Є (a,b) x₀=a x₀=b

X* Є

Если допустимая погрешность задана u=E, то как только , то процесс заканчивается. За решение ур-ия удобно взять среднее

6. Метод итерации.

f(x)=0, заменим это ур-е равносильным x=f(x) и построим приближение x₂=f(х₁),…x_n+1=f(x). Если эта последовательность существует и сходиться, т.е. сущ. Предел x*=f(x*) с любой степенью точности.

Если φ(х) возр. то получается лестницей, если убыв. то строится по спирали.

Теорема о сходим. метода итерации.

Пусть φ(х) опред. и диф. на [a,b] тогда, если сущ. дробь q такая |φ ^‘ (x)|<=q<1 Тогда 1) процесс итерации А(n+1)=φ(Xn) сход. независ. от Х0.2) Предельн. знач-е явл. единств. корнем ур-ия LimXn=X*(n→бескон).Д-во: рассмотр. Хn=φ(X(n-1)); X(n+1)=φ(Xn); вычтем, по теорем. Лагранжа будем иметь a(n+1)-Xn=(Xn-X(n-1))*φ’(Xn( с чер.));|X(n+1)-Xn|<=

q |Xn-X(n-1)|. Пусть n=1,2… Тогда.|X(n+1)-Xn|<=qⁿ|Xn-X0|. Рассмотр. ряд Х0+(Х1-Х0)+(Х2-Х1)+… для этого ряда посл-ти приближен. явл. (n+1)-ми частными суммами Т.е. Xn=S(n+1) В силу неравенств все члены ряда < соотв. членов геом. прогрессии .У котор. знаменат. <1 поэт. ряд сход. абсолютно. LimSn+1=LimXn=X*(n-бескон.)-единств. решение.X(n+1)=φ(Xn), X*. Итер. процесс сход. тем быстрее, чем |φ ^‘ (x)|<1.

Оценка итер. процесса:

|X*-Xn|<=q ⁿ /(1-q)*|X1-X0|. Рассм. как свести ур-ие к виду, удобн. для итерации. f(x)=0; x=x-λf(x)-ур-ие равносил. исходному. x-λf(x)=φ(x); [x=φ(x)];Пусть f ’(x) удовл. неравенству 0<m<f ‘(x)<M. Нужно подобр. такое λ чтоб |φ’(x)|=|1-λf ‘’(x)|<1;Т.к.рассматр. положит. знач-я, то модули можно снять. φ’(x)=1-λf ‘(x); 1-λM₁<1-λm₁<=q<1;Т.о. можно положить λ=1/M₁ . Исходя из этого вычисл. q; q=1-m/M; Этот спос. удобен если извест. интервал, в котор. нах. точное знач-е корня. Если задан. Х0, то λ наход. из ур-ия 1-λf ^‘ (X0)=0.

7. Выбор коэффициента в методе итераций.

Нужно учитывать что в у-ии х= φ(х), φ(х) должно удовлет усл φ(х)<1. Этого можно добиться таким способом Рассмотр ур-е F(x)=0 это ур-е равносильно х=х-ЛF(x)=> φ(х)=ч-ЛF(x), где парам Л выбирают таким образом что бы 0=<1-ЛF’(x)=<q<1 По условию задачи m1=<F’(x)=<M1(наиб и наим значение) на [a;b] => Л=1/M1

(Примечание (Л – лямбда))

Метод итераций решения нелинейных систем второго порядка (x,y)-реш. сист

При локализации корня графич. способом, удобно эти значения применять за нулевое приблежение. Из 1-ого ур-я явно выразим х, а из 2-ого у.

Lim(n  )Xn=x* Lim(n  )Yn=y*

8.Метод Ньютона.

f(c чер.)(x)=0 f1(x1..xn)=0, fn(x1..xn)=0 предпол. что Xp(c чер.)=(X1p,…Xnp) найдено 1-ого из корней, тогда этот корень X=Xp(c чер.)+E(c чер.)p, где E(c чер.)p=(E1p…Enp). E-погрешн. f(x)= f(a)+ f ^‘ (a)(x-a). Выпиш. это ур-ие в разв. виде fn(…Xnp+Enp)=fn(…Xnp)+f ^‘ n(X1)*(…Xnp)*E1p+…

+f ^‘n(Xn)*(X1p…Xnp)*Enp. т.о f ^‘(х)=W(X)=(матрица)( 1) f ^‘1(Х1)… f ^‘1Хn 2) f ^‘nX1… f ^‘nXn)- Якобиан

Исход. сист. можно зап. в виде f(Xp)+W(Xp)Ep=0; Ep=-W(Xp)f(Xp). Итерац. процесс X(p+1)=Xp-W^–1 (Xp)f(Xp). Теорема: пусть Х0-нач. приближ. причем вып. условия 1) W(X)-якобиан имеет обр. матр. Г0= W^-1(x0), ||Г0||<=A 2)норма ||Г0f(x0)||<=B 3) сумма|δ² fi/δXiδXj|<=C 4)множитель μ0=2nABC<1 . При выпол. этих услов. процесс сходится. X(p+1)=Xp-W^-1(Xp)f(Xp).

Модиф. метод Ньютона.

На каждом шаге процесса Ньют. приход. обращать матрицу. Предпол. что W^-1(Xp)(примерно =) W^-1(x0). Тогда X(p+1)=Xp- W^-1(x0)f(Xp)-эта ф-ла в точности совпад. с фор-лой метода прост. итерации.Теорема: при выполн. услов. 1) W(X)-якобиан имеет обр. матр. Г0= W^-1(x0), ||Г0||<=A 2)норма ||Г0f(x0)||<=B 3) сумма|δ² fi/δXiδXj|<=C 4)множитель μ0=2nABC<1 , модиф. метод Ньютона сход. к решению системы.

9. Метод скорейшего спуска решения линейных систем

Есть система

Предположим, что все f_i непр. дифф., рассм. Ф-цию U(x) = (,)=

. Каждое решение сис. обращает в 0 ф-цию U и наоборот

U()  0 => будем искать минимум ф-ции U, кот. и будет решением

Пусть - 0-ое приближение. Проведем пов. уровня U() и из (.)будем двигаться по нормали к пов-ти уровня (в напр., противоп. grad U = U = (), пока не “коснемся” следующей пов-ти уровня. Точка касания –. И т д. Очевидно, что , таким образом придем в (.) min U, а этот min – решение сист. Можно записать , где _p – это нек-й коффициет, вычисляемый на каждом шаге. Задача в том, чтобы его найти. Рассм ф-ю (_p) = U(). Она задает новый уровень ф-ции U который получен продвижением вдоль соответствующей нормали к пов-ти уровня в (.) x^p и зависит от _p. Множители _p нужно выбирать т.о. чтобы (_p) имела минимум при этом знач. _p, т.е. (_p)/_p = 0. Наименьший положит корень это и есть _p. (_p) = Разложим функции f_i в ряд Тейлора до линейных членов. (_p) = =>

выразим отсюда _p и подставим значение U = 2W^Tf(x). В итоге получим:

x^p+1=x^p-_pW^T(x^p)f(x^p), где _p = 2_p = (f_p,W_pW_p^Tf_p)/( W_pW_p^Tf_p, W_pW_p^Tf_p)

Для линейных систем видаW=A и f_p = r_k =- называется вектором невязки линейной системы.

Метод скорейшего спуска решения систем.(grad)

Grad-задаёт направление. Предположим что все fi непрерыв диф-мы Рассмотрим ф-ю u(x)=(_f(_x);_f(_x))= ф-я сама на себя скалярная = (i от 1 до n)(fi(_x))². Очевидно что каждое решение исход сист обращ в 0 ф-ю u(x) и наоборот те числа обращают в 0 ф-ю u(x) явл решеним исх сист. Предположим что исх сист имеет лишь решения которые явл (.) мин ф-ии u таким образом задача сводится к нахождения (.) мин ф-ии u. Пусть _х – решение сист _х(0) – нулевое приблежение. В (.) _х(0) проведём поверхность уровня u(x)=u(xi0), если нач (.) близка к решению, то поверхность похожа на эллипсоид . Из(.) х(0) будем двигаться по нормали поверхности u(x)=u(x0) до тех пор пока нормаль не каснётся другой поверхности уровня u(x)=u(x(1)) Затем отправляемся из (.) х(1) по нормали поверхности уровня u=u(x(1)) до тех пор пока эта нормаль не каснётся поверхности уровня в (.) х(2) и т.д.

Очевидно что u(x(0))>(u(x(1)))>(u(x(n))), т.е. прийдём к (.) мин ф-ии, а этот минимум явл реш исх сист (рис1), grad U. X(p+1)=x(p) - Лp U(x(p)), Лр – нек коэф на каждом шаге задача состоит в том что бы его найти. Ф(Л)=U(x(p) - ЛU(x(p)) эта ф-я задет изменения ф-ии уровня и вдоль соотв нормали к поверхности уровня в (.) х(р) Множ Л=Лр нужно выбирать таким образом, что бы ф-я φ(Л) имела мин т.е. U(x(p) – ЛNU(x(p)) наим положит корень это и есть Л. Получ сист ур-й будем решать численно, поскольку предпол., что Л мало, то будем брать только лин. члены разлож ф-ии в ряд Тейлора x(p)=x(p) - μpwT (xP)f(x(p)) где μр=2Лр=

(Примечание (Л – лямбда))

10. Формула ТРАПЕЦИИ.

Разобьем отрезок [a; b]

на n равных отрезков

[x₀,x₁], [x₀,x₂], …, [x₀,x_n]

шаг h = (b-a)/n

y_i=f(x_i)

Ф-ла трапеции: = (y₀+y₁)*h/2+(y₁+y₂)*h/2+…+(y_n-1+y_n)*h/2= (y₀+y_n)*h/2+h*(y₁+y₂+…+y_n-1)

Нахождение остаточного члена формулы трапеций.

Рассмотрим остаточный член для 1-го звена, а остальные просуммируем

Предположим, что ф-ция f дважды дифференцируема. Получим ф-цию

R(h) =- h/2(y(x₀+h)+ y(x₀)). Продифференцируем это рав-во по h

; (y(x₀))’_h = 0; R’(h) = y(x₀+h) – (1/2)*(y₀(x₀+h) + y(x₀)) -

- (h/2)*y’(x₀+h) = (1/2)*(y(x₀+h) - y(x₀)) - (h/2)*y’(x₀+h)

R”(h) = (1/2)*y’(x₀ + h) - (1/2)*y’(x₀ + h) - (h/2)*y”(x₀ + h) = - (h/2)*y”(x₀ + h)

R(0)=R’(0)=R”(0)=0. Проинтегрируем по h и воспользуемся теоремой о среднем

, c₂(x₀, x₀+h)

, c₁(x₀, x₀+h)

Теперь просуммируем остаточные члены для одного звена и получим полный

R(h) = - (h³/12)*y”(c_i), c_i(x_i-1, x_i)

Далее,  (.) c[a, b], такая что (y”(c_i))/n = y”(c) => R(h) = - (h³/12)*n*y”(c) =>

R(h) = - (b-a)*h²*y^”(c)/12, c[a, b]

На практике: I_h = _h + Mh², I_2h = _2h + 4Mh² => R = (_h - _2h)/3

11.Формула СИМПСОНА (парабола).

∫f(x)dx (рис.).

Число разбиен. должно делиться на 4, чтобы можно было вычислить суммы с шагом 2h. Пары интервалов, следующие друг за другом, начиная с первой, будем “накрывать” параболой, проход. через 3 (.)-ки y=A₀+A₁x+A₂x²; (x₀,y₀), (x₁,y₁), (x₂,y₂). Предпол что x₀ = 0 => x₁ = h, x₂ = 2h. Для нахожд. A₀,A₁,A₂ подстав. эти (.)-ки в ур-ие параб. =>

y₀ = A₀; y₁ = A₀+A₁h+A₂h²; y₂=A₀+2A₁h+4A₂h²;

Получим, A₀ = y₀; A₁ = (-3y₀ + 4y₁ - y₂)/2h; A₂ = (y₀ - 2y₁ + y₂)/2h²

Площадь одного эл-та (зависит от h и значений ф-ции в выбранных точках) = = A₀2h+A₁(4h²)/2+A₂h³8/3;

Подставл знач-я А₀, А₁, А₂ получ. I_n=(h/3)*(y_n-2 + 4y_n-1 + y_n);

 I_n = = (h/3)*[y₀ + y_n + 2*y_2k + 4* y_2k-1] – ф-ла Симпсона

Оценка погрешности в методе Симпсона.

Предположим, что ф-я y непрерывна и трижды дифференцируема.

Запишем ошибку для первого интервала в виде:

- (h/3)*[y(x₁-h) + 4y(x₁) + y(x₁+h)]

Продифференцируем эту ф-ю 3 раза.

R’(h) = y(x₁+h) – y(x₁-h) – (1/3)[y(x₁-h) + y(x₁+h) + 4y(x₁)] – (h/3)*[-y’(x₁ – h) +

+y’(x₁+h)]

R”(h) = - (1/3)y’(x₁-h) + (1/3)y’(x₁+h) – (h/3)[y”(x₁-h) + y”(x₁+h)]

R”’(h) = -(h/3)[y”’(x₁-h) - y”’(x₁+h)] = -(2h²/3)*y^IV(c₄), c₄(x₁-h, x₁+h)

Теперь будем интегрировать:

R”(h) = , c₃(x₁-h, x₁+h)

R’(h) = = -(1/18)*h⁴*y⁴(c₂), c₂(x₁-h, x₁+h)

R(h) = = -(1/90)*h⁵*y^IV(c₁), c₁(x₁-h, x₁+h) – ошибка для 1-ой пары интервалов; Для всего отрезка ошибка R = -((b-a)/180)*h⁴*y^IV(c), c(a, b);

На практике: R_h = Mh⁴, R_2h = 16Mh⁴; R = (_h - _2h)/15

12. Метод Эйлера решения дифференциальных уравнений

y’=f(x,y); y(x₀) = y₀ - задача Коши. Разобьем ось х на промежут. x_i = x₀+i*h (рис.) В (.)M₀-пров. касател. до пересеч. с прямой х = х₁, получим (.)M₁. y₁ = х₁ M₁

Через (.) (x₁, f(x₁, y₁)) пров касат до пересечения с прямой x = x₂, получим M₂ и тд Получ. Ломаную M₀M₁M₂… считаем решением.

Аналитически: f(x_i,y_i) = y’ (y_i+1 - y_i)/h => y_i+1 = y_i + h*f(x_i,y_i)-ф-ла облад. малой точностью.(порядка h²); y(xi+h)=y(xi)+y’(xi)h остал. члены разлож. в ряд Тэйлора (y(x)=y(x₀)+y’(x₀)(x-x₀)+y”(x₀)(x-x₀)²/2!-поэт. точность порядка h².

Модифицированный метод Эйлера.

В отличии от метода Эйлера, когда для вычисления следующей точки (X_i+1, Y_i+1) требуется информация только Y предыдущей точки. Мод. метод предполагает знание о нек-ой промежуточной точки x_i+1/2=x_i+h/2 и x_i+1/2=x_i+(h/2)*f_i. Метод заключается: 1) через точку (x_i,y_i) проводится кас-ая с tg наклона tgα=f(x_i,y_i) до х с прямой x=x_i+1/2. В полученной точке х по методу Эйлера выч-ся знач-я ф-ии y=y_i+1/2 и выч-ся знач-е производной f_i+1/2=f(x_i+1/2, y=y_i+1/2). Знач-е этой производной определяет tg угла наклона 2-ой производной, кот. проводится из полученной точки. 2) Затем возвр-ся в исходную точку и через неё проводим прямую, парал 2-й кас-ой.

Усовершенствованный метод Эйлера.

1. через точку (x_i,y_i) проводится касс-ая до х с прямой x_i+1/2=x_i+h. Угл. коэф-т этой кривой tgα₁=f(x_i,y_i). 2) В полученной точке по методу Эйлера выч-ся знач-е ф-ии y_i+1=y_i+h*( f(x_i,y_i)) и выч-ся новая производная tg α₂= f(x_i,y_i+1). 3) Далее происходит возврат в точку (x_i,y_i) и через нее проводится новая касательная, где tg угла наклона есть среднее арифметическое 2-ух предыдущих tg-енсов, т.е. tgα = (tgα₁+ tgα₂)/2 = (f(x_i,y_i)+F(x_i+1,y_i+1))/2. Новое значение y_i+1 = y_i+h*((f(x_i,y_i)+F(x_i+1,y_i+1))/2).

Метод Рунге-Кутта.

Обычно говорят, что метод Р-К явл. Методом 4 и 5 порядка точности, т.е. (h⁴,h⁵). Точно вычислить погрешность этого метода затруднительно, т.к. исходя из текущего значения y(x_i) вычисляют величину y(x_i+2h) 2-мя способами. 1 раз с шагом h, другой с шагом 2h. Если расхождение полученных значений не превышает допустимой погрешности, то шаг выбран правильно и полученное значение можно принять за верное, в противном случае шаг уменьшают в 2 раза. Мы решаем ур-е

Выбирается шаг h, наносится сетка x_i+1=x₀+hi. Рассматривают числа k₁⁽ⁱ⁾=h*f(x_i,y_i); k₂⁽ⁱ⁾=h*f(x_ih/2,y_i+k₁⁽ⁱ⁾/2); k₃⁽ⁱ⁾=h*f(x_ih/2,y_i+k₂⁽ⁱ⁾/2); k₄⁽ⁱ⁾=h*f(x_ih,y_i+k₃⁽ⁱ⁾). Вычисляются знач-я k, соотв-но y_i+1=y_i+Δy_i. Шаг Δ выч-ся т.о., что Δy_i=1/6 (k₁⁽ⁱ⁾+ 2k₂⁽ⁱ⁾+ 2k₃⁽ⁱ⁾+ k₄⁽ⁱ⁾). Удобно на каждом шаге записывать таблицу

13 .Метод прогонки

Будем рассматривать д.у 2 порядка: y” +p(x)y’+q(x)y =f(x) при заданных условиях α0y(a)+ α1y’(a)=A β0y(b)+ β1y’(b)=B │α0│+│α1│≠0 │ β0│+│ β1│≠0 Исходная система может быть представлена в виде разностной схемы (yi+1–2yi+yi-1)/h2+pi(yi+1-yi-1)/2h +qiyi=fi (1), где pi, qi, fi –значения коэффициентов в точках xi. α0y0+ (α1(y1-y0))/h=A β0yn+ (β1(yn-1-yn))/-h=B Будем использовать метод прогонки Из ур (1) можно получить yi+1+miyi+niyi-1=f^ih2 (4) где коэффициенты имеют вид: mi=(2-qih2)/(1+pih/2) ni= (1-pih/2)/(1+pih/2) f^i=fi/(1+pih/2) Из уравнения (4) выразим yi yi=f^ih2/mi – (1/mi)(yi+1)-(ni/mi)(yi-1) (5) Предположим, что с помощью системы (2), (3), (4) из ур-ния (5) исключено yi-1,тогда yi= ci(di-yi+1) (6) Из ур (6) можно записать yi-1= ci-1(di-1-yi) Подставляя это уравнение в ур-ние (4) получим yi+1+miyi+nici-1di-1-yi=f^ih2

Отсюда выразим yi

yi=( f^ih2- nici-1di-1-yi+1)/(mi-nici-1) Исходя из сравнения полученных формул получим ci=1/mi-nici-1 di=f^ih2-nici-1di-1 c0= α1/ (α0h - α1) d0=Ah/ α1 На основании этих формул последовательно определяются сi и di при i=1 до n-1 включительно Обратный ход начинается с вычисления yn Вычислив yn по ф-ле (6) определяются все остальные значения y

14.Уравнение Лапласа в конечных разностях

Будем решать δ²u/δx²+ δ²u/δy²=0 Будем решать задачу Дирихле, т.е найдём решение уравнения Лапласа в заданной области, ограниченной контуром гамма причём u(p)=φ(p) PГ для всех точек p принадлежащих границе.Т.е на границе области задана непрерывная функция. Перейдём от дифура к его разностному аналогу δ²u/δx²= =(u(x+h,y)-2u(x,y)+u(x-h,y))/h² δ²u/δy²=(u(x,y+h)-2u(x,y)+u(x,y-h))/h² Подставим u(x,y)=1/4(u(x+h,y)+u(x-h,y)+ u(x,y+h) +u(x,y-h)) Это уравнение соответствует первой разностной схеме. 1 осн схема изображена на левом рисунке. 2 осн сх: u(x,y)=1/4(u(x+h,y+h)+u(x+h,y-h)+ u(x-h,y-h) +u(x-h,y+h)) И 1 и 2 схема являются точными до порядка h²

Метод сеток. Построение шаблонов

Пусть в плоскости XOY имеется область G с границей Г, построим на плоскости два семейства параллельных прямых x= x₀+ih y=y₀+kh

В каждой внутренней точке заменим исходное уравнение конечно-разностным по 1 или 2 схеме. В граничных точках u(Bh)=u(B)=φ(B),где B-ближайшая к границе точка. Для внутренних узлов:u^k_ij =1/4(u^k-1_i-1,j+ u^k-1_i+1,j+ u^k-1_i,j-1+ u^k-1_i,j+1) В граничных точках значения могут быть поправлены по формулам линейной интерполяции u⁰(Ah)= u(a)=φ(A) u^k(Ah)= u(a)+(( u^k-1(B)-u(A))/(h+δ))δ точка A ближайшая к узловой точке Ah границы Г B-ближайший к Ah внутренний узел сетки δ- удаление узла Ah от точки A. Значение δ может быть как + так и – Если узел расположен на границе u^k(Ah)= u(a)=φ(A) При решении задачи для выбора u⁰_ij используют линейную интерполяцию. Иногда для решения задачи используют крупную сетку, а затем решают её же но на более мелкой сетке. Построение шаблонов. Для построения шаблонов строим вторую сетку, линии которой проходят посередине между линиями первой, причём так, что узлы первой сетки попадают в центры клеток 2 сетки

16. Интерполирование функций. Постановка задачи.

Пусть на [a,b] заданы (n+1) точек x₀,x₁,…,x_n (эти точки наз-ся узлами интерполяции). Известны значения нек-ой ф-ии в этих точках: f(x₀)=y₀, f(x₁)=y_1, f(x_n)=y_n. По этим точкам требуется построить ф-ю F(x), принадлежащую известному классу, принимающую в узлах интерполяции известные знач-я. В такой интерполяции ф-ия F(x) либо не однозначна, либо не сущ-ет вовсе. Если же ф-ия F(x)=P_n(x) имеет вид полинома степени не выше степени n, то такая задача имеет единственное решение. Полученную интерполяционную ф-ию F(x) используют для нахождения значений ф-ии в точках, отличных от узлов интеропляции. Если эти точки принадлежат промежутку [a,b], то такая задача наз-ся интерполированием в узком смысле. Если же точки выходят за границы [a,b], то такая задача наз-ся экстраполированием.

Метод наименьших квадратов.

(матем. обработка данных).Пусть в рез-те исслед. некоторой величины Х знач-я Х1,Х2…Хn поставлены в соответствие др. знач-я …Xn→y1,y2…yn; Треб. подобрать вид аналитич. записи y=y(x) связыв. перемен. х и у; В задаче треб. найти некую прямую y=bx; Задача закл. в определен. пост. b; В методе наим. квадр. параметр b опред. из условий min суммы квадратов отклонений эксперим. данных yi от расчетных величин yi*; F=∑(от 1 до бесконеч.)(yi-bxi)²→min; δF/δb =2∑((yi-bxi)*(-xi))=0; b=(∑xiyi)/(∑xi²); В нек. случаях иском. прям. имеет вид y=a+bx; F=∑(yi-a-bxi)²→min; δF/δb=2∑[(yi-a-bxi)(-xi)]; ∑(yixi-axi-bxi²)=0; ∑yixi=∑(axi)+∑bxi²; ∑xiyi=a∑xi+b∑xi²; a= (∑yi-b∑xi)/n; n∑xiyi=∑xiyi-b(∑xi)²+bn∑xi²; b=(n∑xiyi-∑xi∑yi)/(n∑xi²-(∑xi)²).

17. Формула Лагранжа.

Будем считать данную ф-ию f(x) и полином Q_m(x)=a₀+a₁x+a₂x²…+a_mx^m близкими, если они совпадают на заданной системе точек x₀,x₁,x₂…x_n. Задача состоит в том, чтобы построить многочлен возможно низшей степени m, принимающий в данных точках известные значения. По основной теореме алгебры можно предположить, что m=n. Исходя из условия задачи, можно записать систему линейных уравнений.. Эта система линейная, её можно решать по формуле Крамера (определитель Вандермода): Δ==П(x_q-x_p) – произведение, где <q. Исходя из этого, можно выписать интерполяционный многочлен Лагранжа:

18.Формула Ньютона для разностных узлов.

∆y₀ = y₁-y₀ ∆y₁= y₂-y₁ ∆y₂= y₃-y₂ ∆y_n-1 = y_n - y_n-1

∆y₀ = y₁-y₀ ∆²y₀ = ∆y₁- ∆y₀ = y₂ – y₁ – y₁ + y₀ = y₂ – 2y₁ + y₀ ∆³y₀= ∆²y₁- ∆²y₀ = y₃ – 3y₂ +3y₁ – y₀

∆^ky₀ = y_k – ky_k-1 + ((k(k-1))/2!)y_k-2 + ((k(k-1)(k-2))/3!)y_k-3 + ...+ (-1)^Ky₁

Запишем эту формулу для значения разности в узле x_i: ∆^ky_i = y_k+i – y_k+i-1 + ((k(k-1))/2!)y_k+i-2 y_k=y₀ +k∆y₀ +((k(k-1)/2!)∆²y₀ + … ∆^ky₀

Построим интерполяционный многочлен Ньютона: N(x)= a₀ +a₁(x-x₀) +a₂(x-x₀)(x-x₁) +a₃(x-x₀)(x-x₁)(x-x₂) + … a_n(x-x₀)(x-x₁)…(x-x_n-1) Этот многочлен должен проходить через заданные узлы, поэтому: N(x₀)= a₀=y₀ N(x₁)= a₀ +a₁(x-x₀) = a₀ + a₁h=y₁ N(x₂)= a₀ +2a₁h+ 2a₂h² = y₂ Отсюда найдём коэффициенты: a₀=y₀ a₁= (y₁-y₀)/h=∆y₁/h a₂=∆²y₀/2h² Следовательно, любой коэффициент имеет вид: a_k = ∆^ky₀/k!h^K и полином Ньютона будет иметь вид: N(x)=y₀+ (∆y₀(x-x₀))/h + (∆²y₀(x-x₀)(x-x₁))/2!h² + … + (∆ⁿy₀(x-x₀)(x-x₁)…(x – x_n-1))/n!hⁿ

Запишем его иначе: пусть (x-x₀)/h =q тогда N(x) =(y₀ +q∆y₀∆²y₀q(q-1))/2! + … + ∆ⁿy₀(1/n!)q(q-1)(q-2)…(q-n+1) Многочленами Ньютона пользуются в случае равноотстоящих узлов.

Приближенное дифференцирование. Постановка краевых задач

Пусть имеем функцию y(x) заданную в равноотстоящих точках. Запишем многочлен Ньютона для этой функции y(x)= y₀+q∆y₀+(q(q-1) ∆²y₀)/2 + (q(q-1)(q-2)∆³y₀)/3! + (q(q-1)(q-2)(q-3)∆⁴y₀)/4! Для нахождения производной функции будем искать производную этого многочлена h=x_i+1-x_i; q=(x-x₀)/h y(x)= y₀+q∆y₀+ ((q²-q)∆²y₀)/2+((q³-3q²+2q)∆³y₀)/3! + ((q⁴-6q³+4q²-6) ∆⁴y₀)/4!+… dy/dx=(dy/dq)(dq/dx)=(1/h)(dy/dq) y’(x)=1/h(∆y₀+((2q-1)∆²y₀)/2+((3q²-6q+2)∆³y₀)/3! + ((4q³-18q²+22q-6)∆⁴y₀)/4!+…) Аналогично вычисляются производные любых порядков. Иногда требуеся находить производные функции в основных табличных точках x_i, тогда производная выглядит следующим образом: y’(x₀)= (1/h)(∆y₀ - ∆²y₀/2 + ∆³y₀/3 -∆⁴y₀/4+…) Постановка краевых задач: пусть дано: y”=f(x,y,y’) y(a)=A y(b)=B Геометрически это означает, что требуется найти интегральную кривую, проходящую через заданные точки. В некоторых случаях краевые условия сводятся к заданиям значений производных y’(a)=A y’(b)=B, т.е нужно найти такую интегральную кривую, которая в заданных точках проходит под известным углом.

19. Абсолютная и относительная погрешность.

Абсолютной погрешностью приближенного значения а* называется величина . Относительной погрешностью значения а* (при a не равно 0) называется величина   .