м |
|
Строка 1: |
Строка 1: |
- | [[Численные Методы, 14 лекция (от 02 апреля)|Предыдущая лекция]] | [[Численные Методы, 16 лекция (от 09 апреля)|Следующая лекция]]
| + | == From Ebaums Inc to MurkLoar. == |
- | | + | We at EbaumsWorld consider you as disgrace of human race. |
- | = Глава 4. Разностные методы решения задач математической физики = | + | Your faggotry level exceeded any imaginable levels, and therefore we have to inform you that your pitiful resourse should be annihilated. |
- | == Параграф 1. Разностные схемы для первой краевой задачи уроавнения теплопроводности ==
| + | Dig yourself a grave - you will need it. |
- | === Пункт 3. Симметричная разностная схема (схема Кранка-Никольсена) ===
| + | |
- | | + | |
- | δ<sup>2</sup>u/δx<sup>2</sup> + f(x, t) (1)
| + | |
- | * u(0, t) = μ<sub>1</sub>(t); u(1, t) = μ<sub>2</sub>(t), 0 ≤ t ≤ T (2)
| + | |
- | * u(x, 0) = u<sub>0</sub>(x), 0 ≤ x ≤ 1 (3)
| + | |
- | * y<sub>x_x, i</sub><sup>m</sup> = (y<sub>i + 1</sub><sup>m</sup> − 2y<sub>i</sub><sup>m</sup> + y<sub>i − 1</sub><sup>m</sup>)
| + | |
- | * (y<sub>i</sub><sup>n + 1</sup> − y<sub>i</sub><sup>n</sup>)/τ = 0,5(y<sub>x_x, i</sub><sup>n + 1</sup> + y<sub>x_x, i</sub><sup>n</sup>) + f(x<sub>i</sub>, t<sub>т + 1/2</sub>), (x<sub>i</sub>, t<sub>j</sub>) ∈ ω<sub>τh</sub> (4)
| + | |
- | * y<sub>0</sub><sup>n + 1</sup> = μ<sub>1</sub>(t<sub>n + 1</sub>); y<sub>N</sub><sup>n + 1</sup> = μ<sub>2</sub>(t<sub>n + 1</sub>), t<sub>n + 1</sub> ∈ ω_<sub>τ</sub> (5)
| + | |
- | * y<sub>i</sub><sup>0</sup> = u<sub>0</sub>(x<sub>i</sub>), x<sub>i</sub> ∈ ω_<sub>h</sub> (6)
| + | |
- | | + | |
- | Шаблон (типа ящик)
| + | |
- | | + | |
- | Эта схема неявная. Решается она так же, как и чисто неявная схема — методом прогонки.
| + | |
- | | + | |
- | Сходимость будем рассматривать среднеквадратичную, а не в С, ибо в С очень сложно. Будем использовать новый аппарат. Ради чего напрягаться? Оказывается, эта схема имеет второй порядок и по тау, и по h, то есть лучше. Сегодня освоим аппарат Фурье метода оценок, и его тоже можно здесь применить. Можно строить функцию Грина, можно энергетическими неравенствами. Способов очень много методов оценок.
| + | |
- | | + | |
- | ==== Погрешность ====
| + | |
- | z<sub>i</sub><sup>n</sup> = y<sub>i</sub><sup>n</sup> − u<sub>i</sub><sup>n</sup>
| + | |
- | | + | |
- | * (z<sub>i</sub><sup>n + 1</sup> − z<sub>i</sub><sup>n</sup>)/τ = 0,5(z<sub>x_x, i</sub><sup>n + 1</sup> + z<sub>x_x, i</sub><sup>n</sup>) + ψ<sub>i</sub><sup>n</sup>
| + | |
- | | + | |
- | Здесь принципиально, как брать, и это принципиально, иначе не получим второй порядок.
| + | |
- | | + | |
- | * z<sub>0</sub><sup>n + 1</sup> = 0, z<sub>N</sub><sup>n + 1</sup> = 0, z<sub>i</sub><sup>0</sup> = 0
| + | |
- | * ψ<sub>i</sub><sup>n</sup> = (u<sub>x_x</sub><sup>n + 1</sup> + u<sub>x_x</sub><sup>n</sup>) + f(x<sub>i</sub>, t<sub>n + 1/2</sub>) − (u<sub>i</sub><sup>n + 1</sup> − u<sub>i</sub><sup>n</sup>)/τ (8)
| + | |
- | | + | |
- | '''Задача'''. Доказать, что ψ<sub>i</sub><sup>n</sup> = O(τ<sup>2</sup> + h<sup>2</sup>).
| + | |
- | | + | |
- | ==== Задача Штурма-Лиувилля ====
| + | |
- | | + | |
- | * δ<sup>2</sup>u/δx<sup>2</sup> + λu(x) = 0, 0 < x < 1 (9)
| + | |
- | * u(0) = u(1) = 0, u(x) !== 0
| + | |
- | * λ — собственное значение, u(x) — собственная функция
| + | |
- | * &lambda<sub>k</sub> = (πk)<sup>2</sup>, k = 1, 2, …
| + | |
- | * 0 < &lambda<sub>1</sub> < λ<sub>2</sub> < λ<sub>k</sub> < …
| + | |
- | * u<sub>k</sub>(x) = sqrt(2)sin πkx
| + | |
- | * {u<sub>k</sub>(x)}<sub>1</sub><sup>∞</sup> — ортонормированный базис
| + | |
- | * f(x) ∈ L<sub>2</sub>(0, 1)
| + | |
- | * интеграл от 0 до 1 f<sup>2</sup>(x) dx < infin;,
| + | |
- | ... | + | |
- | | + | |
- | ||f||<sub>L<sub>2</sub></sub><sup>2</sup> = ∑<sub>k = 1</sub><sup>∞</sup> c<sub>k</sub><sup>2</sup> — уравнение Парсеваля
| + | |
- | Разностная задача Штурма-Лиувилля
| + | |
- | {|
| + | |
- | |rowspan="2"|{
| + | |
- | |y<sub>x_x, i</sub> + λy<sub>i</sub> = 0, x<sub>i</sub> ∈ ω_<sub>h</sub>
| + | |
- | |rowspan="2"|(10)
| + | |
- | |--
| + | |
- | |y<sub>0</sub> = y<sub>N</sub> = 0
| + | |
- | |}
| + | |
- | | + | |
- | * y<sub>i + 1</sub> − 2y<sub>i</sub> + y<sub>0i − 1</sub> + λh<sup>2</sup>y<sub>i</sub> = 0
| + | |
- | * y<sub>i + 1</sub> + y<sub>i − 1</sub> − (2 − λh<sup>2</sup>)y<sub>i</sub> = 0
| + | |
- | * y<sub>i</sub> = sin αx<sub>i</sub> ≠ 0
| + | |
- | * y<sub>i + 1</sub> = sin α(x<sub>i</sub> + h) ≠ 0
| + | |
- | * y<sub>i − 1</sub> = sin α(x<sub>i</sub> − h) ≠ 0
| + | |
- | * y<sub>i + 1</sub> + y<sub>i − 1</sub> = sin (αx<sub>i</sub> + αh) + sin (αx<sub>i</sub> − αh) = 2sin αx<sub>i</sub> cos αx
| + | |
- | * 2cos αh sin αx<sub>i</sub> − (2 − λh<sup>2</sup>)y<sub>i</sub> = 0
| + | |
- | * 2cos αh − 2 + λh<sup>2</sup> = 0
| + | |
- | * λ = 2(1 − cos αh)/h<sup>2</sup> = 4/h<sup>2</sup>sin(αh/2)
| + | |
- | * λ<sub>k</sub> = 4/h<sup>2</sup> sin<sup>2</sup> πkh/2, k = 1…N − 1
| + | |
- | * y<sub>k</sub>(x<sub>i</sub>) = sqrt(2) πkx<sub>i</sub>, x<sub>i</sub> ∈ ω<sub>h</sub>, k = 1…N − 1
| + | |
- | | + | |
- | равенство Парсеваля:
| + | |
- | dim H<sub>N − 1</sub> = N − 1, y<sub>0</sub> = y<sub>n</sub> = 0
| + | |
- | | + | |
- | f, g ∈ H<sub>N − 1</sub>, (f, g) = &sum ;<sub>i = 1</sub><sup>N − 1</sup> f<sub>i</sub>g<sub>i</sub>h
| + | |
- | * ||f||<sub>L<sub>2</sub>(ω<sub>n</sub>)</sub> = ||f||<sub>L<sub>2</sub></sub> = (f, f)<sup>1/2</sup> = ∑<sub>i = 1</sub><sup>N − 1</sup>f<sub>i</sub><sup>2</sup>h
| + | |
- | | + | |
- | * (y<sub>k</sub>(x<sub>i</sub>) y<sub>l</sub>(x<sub>i</sub>)) = δ<sub>k, l</sub>
| + | |
- | * f(x<sub>i</sub>) = ∑<sub>k = 1</sub><sup>N − 1</sup> c<sub>k</sub>y<sub>k</sub>(x<sub>i</sub>)
| + | |
- | | + | |
- | Тождественно обозначим: y<sub>k</sub>(x<sub>i</sub>0 = μ<sub>k</sub>(x<sub>i</sub>)
| + | |
- | | + | |
- | * ||z<sup>n</sup>||<sub>L<sub>2</sub>(ω<sub>n</sub>)</sub> → 0, &tau, h → r
| + | |
- | * z<sub>i</sub><sup>n</sup> = ∑<sub>k = 1</sub><sup>N − 1</sup> c<sub>k</sub>(t<sub>n</sub>)μ<sup>k</sup>(x<sub>i</sub>)
| + | |
- | * ψ<sub>i</sub><sup>n</sup> = ∑<sub>k = 1</sub><sup>N − 1</sup> ψ<sup>k</sup>(t<sub>n</sub>)μ<sup>k</sup>(x<sub>i</sub>)
| + | |
- | * ∑<sub>k = 1</sub><sup>N − 1</sup> (c<sub>k</sub>(t<sub>n + 1</sub>) − c<sub>k</sub>(t<sub>n</sub>))/τ × μ<sup>k</sup>(x<sub>i</sub>) = 0,5 ∑<sub>k = 1</sub><sup>N − 1</sup>(c<sub>k</sub>(t<sub>n + 1</sub>) + c<sub>k</sub>(t<sub>n</sub>))μ<sub>x_x, i</sub><sup>k</sup> + ∑<sub>k = 1</sub><sup>N − 1</sup> ψ<sup>k</sup>(t<sub>n</sub>)μ<sup>k</sup>(x<sub>i</sub>)
| + | |
- | * ∑<sub>k = 1</sub><sup>N − 1</sup> ((c<sub>k</sub>(t<sub>n + 1</sub>) − c<sub>k</sub>(t<sub>n</sub>))/τ + 0,5λ<sub>k</sub>(c<sub>k</sub>(t<sub>n + 1</sub>) + c<sub>k</sub>(t<sub>n</sub>)))μ<sup>k</sup>(x<sub>i</sub>) = ∑<sub>k = 1</sub><sup>N − 1</sup> ψ<sup>k</sup>(t<sub>n</sub>)μ<sup>k</sup>(x<sub>i</sub>)
| + | |
- | * (c<sub>k</sub>(t<sub>n + 1</sub>) − c<sub>k</sub>(t<sub>n</sub>))/τ + 0,5λ<sub>k</sub>(c<sub>k</sub>(t<sub>n + 1</sub>) + c<sub>k</sub>(t<sub>n</sub>)) = ψ<sup>k</sup>(t<sub>n</sub>)
| + | |
- | * c<sub>k</sub>(t<sub>n + 1</sub>) − c<sub>k</sub>(t<sub>n</sub>) + 0,5τλ<sub>k</sub>(c<sub>k</sub>(t<sub>n + 1</sub>) + c<sub>k</sub>(t<sub>n</sub>)) = τψ<sup>k</sup>(t<sub>n</sub>)
| + | |
- | * (1 + 0,5τλ)c<sub>k</sub>(t<sub>n + 1</sub>) &minus (1 − 0,5τλ)c<sub>k</sub>(t<sub>n</sub>) = τψ<sup>k</sup>(t<sub>n</sub>)
| + | |
- | * c<sub>k</sub>(t<sub>n + 1</sub>) = (1 − 0,5τλ)/(1 + 0,5τλ) × c<sub>k</sub>(t<sub>n</sub>) + τ/(1 + 0,5τλ)ψ<sup>k</sup>(t<sub>n</sub>)
| + | |
- | * q<sub>k</sub> = (1 − 0,5τλ)/(1 + 0,5τλ)
| + | |
- | * |q<sub>k</sub>| < 1
| + | |
- | * z<sub>i</sub><sup>n + 1</sup> = ∑<sub>k = 1</sub><sup>N − 1</sup>(q<sub>k</sub>c<sub>k</sub>(t<sub>n</sub>) + τ/(1 + 0,5τλ)ψ<sup>k</sup>(t<sub>n</sub>))μ<sup>k</sup>(x<sub>i</sub>) = ∑<sub>k = 1</sub><sup>N − 1</sup>q<sub>k</sub>c<sub>k</sub>(t<sub>n</sub>)μ<sup>k</sup>(x<sub>i</sub>) + ∑<sub>k = 1</sub><sup>N − 1</sup>τ/(1 + 0,5τλ)ψ<sup>k</sup>(t<sub>n</sub>)μ<sup>k</sup>(x<sub>i</sub>)
| + | |
- | * v<sub>i</sub><sup>n</sup> = ∑<sub>k = 1</sub><sup>N − 1</sup>q<sub>k</sub>c<sub>k</sub>(t<sub>n</sub>)μ<sup>k</sup>(x<sub>i</sub>)
| + | |
- | * ω<sub>i</sub><sup>n</sup> = ∑<sub>k = 1</sub><sup>N − 1</sup>τ/(1 + 0,5τλ)ψ<sup>k</sup>(t<sub>n</sub>)μ<sup>k</sup>(x<sub>i</sub>)
| + | |
- | * z<sub>i</sub><sup>n + 1</sup> = v<sub>i</sub><sup>n</sup> + ω<sub>i</sub><sup>n</sup>
| + | |
- | * ||z<sup>n + 1</sup>||<sub>L<sub>2</sub>(ω<sub>n</sub>)</sub> = ||v<sup>n</sup>||<sub>L<sub>2</sub>(ω<sub>n</sub>)</sub> + ||ω<sup>n</sup>||<sub>L<sub>2</sub>(ω<sub>n</sub>)</sub>
| + | |
- | * ||v<sup>n</sup>||<sub>L<sub>2</sub>(ω<sub>n</sub>)</sub> = ∑<sub>k = 1</sub><sup>N − 1</sup>q<sub>k</sub><sup>2</sup>c<sub>k</sub><sup>2</sup>(t<sub>n</sub>) ≤ ∑<sub>k = 1</sub><sup>N − 1</sup>c<sub>k</sub><sup>2</sup>(t<sub>n</sub>) = ||z<sup>n</sup>||<sub>L<sub>2</sub>(ω<sub>n</sub>)</sub>
| + | |
- | * ||ω<sup>n</sup>||<sub>L<sub>2</sub>(ω<sub>n</sub>)</sub> ≤ τ<sup>2</sup>||ψ<sup>n</sup>||<sub>L<sub>2</sub>(ω<sub>n</sub>)</sub>
| + | |
- | * ||z<sup>n + 1</sup>||<sub>L<sub>2</sub>(ω<sub>n</sub>)</sub> ≤ ||z<sup>n</sup>||<sub>L<sub>2</sub>(ω<sub>n</sub>)</sub> + τ||ψ<sup>n</sup>||<sub>L<sub>2</sub>(ω<sub>n</sub>)</sub>
| + | |
- | | + | |
- | Получили ту же самую формулу, что и в прошлый раз, следовательно, её можно сделать рекуррентной.
| + | |
- | | + | |
- | Мы не использовали связь шагов, поэтому сходимость абсолютная.
| + | |
- | | + | |
- | * ||z<sup>n + 1</sup>||<sub>L<sub>2</sub>(ω<sub>n</sub>)</sub> ≤ ||z<sup>0</sup>|| { = 0 } + ∑<sub>j = 0</sub><sup>n</sup> τ||ψ<sup>n</sup>||<sub>L<sub>2</sub>(ω<sub>n</sub>)</sub>
| + | |
- | * ψ<sub>i</sub><sup>n</sup> = O(τ<sup>2</sup> + h<sup>2</sup>)
| + | |
- | * ∃ M ≥ 0 не зависит от τ и h
| + | |
- | * ||ψ<sup>n</sup>||<sub>L<sub>2</sub>(ω<sub>n</sub>)</sub> ≤ M(τ<sup>2</sup> + h<sup>2</sup>)
| + | |
- | * ||z<sup>n + 1</sup>||<sub>L<sub>2</sub>(ω<sub>n</sub>)</sub> ≤ M(τ<sup>2</sup> + h<sup>2</sup>)∑<sub>j = 0</sub><sup>n</sup> τ = Mt<sub>n</sub>(τ<sup>2</sup> + h<sup>2</sup>) ≤ MT(τ<sup>2</sup> + h<sup>2</sup>)
| + | |
- | * ||z<sup>n + 1</sup>||<sub>L<sub>2</sub>(ω<sub>n</sub>)</sub> ≤ M<sub>1</sub>(τ<sup>2</sup> + h<sup>2</sup>), M<sub>1</sub> не зависит от τ и h
| + | |
- | | + | |
- | ==== Оценка устойчивости ====
| + | |
- | Если бы мы рассматривал нашу оценку и при нулевых краевых условиях, то мы бы получили совершенно аналогично:
| + | |
- | * ||y<sup>n + 1</sup>||<sub>L<sub>2</sub>(ω<sub>n</sub>)</sub> ≤ ||u<sub>0</sub>||<sub>L<sub>2</sub>(ω<sub>n</sub>)</sub> + ∑<sub>j = 0</sub><sup>n</sup> τ||f(t)||<sub>L<sub>2</sub>(ω<sub>n</sub>)</sub> — априорная оценка, означающая устойчивость
| + | |
- | | + | |
- | Симметричная схема обладает более высоким порядком точности, но мы её доказали пока только для среднеквадратичной схемы.
| + | |
- | | + | |
- | {{Численные Методы}}
| + | |
- | {{Lection-stub}}
| + | |
We at EbaumsWorld consider you as disgrace of human race.
Your faggotry level exceeded any imaginable levels, and therefore we have to inform you that your pitiful resourse should be annihilated.
Dig yourself a grave - you will need it.