API tools faq
paste
Advertisement
Guest User

Untitled

a guest
May 23rd, 2018
84
0
Never
Add comment
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 47.25 KB | None | 0 0
raw download clone embed print report
  1. \documentclass[11pt]{article}
  2. %\usepackage[a5paper]{geometry}
  3. \usepackage[left=20mm, right=20mm, top=20mm, bottom=20mm, includefoot]{geometry}
  4. %\geometry{left=20mm, right=20mm, top=20mm, bottom=20mm, includefoot}
  5.  
  6. \usepackage[cp1251]{inputenc}
  7.  
  8. \usepackage[russian]{babel}
  9. \usepackage[intlimits]{amsmath}
  10. \interdisplaylinepenalty=2500
  11. \usepackage{amsfonts}
  12. \usepackage{amsthm}
  13. \usepackage{amssymb}
  14. \usepackage{euscript}
  15. \usepackage{amsmath}
  16. \usepackage{latexsym}
  17. \usepackage{graphicx}
  18.  
  19. \DeclareMathOperator{\sech}{sech}
  20. %нужен, чтобы вставлять рисунки
  21. %\usepackage{caption} %нужен, чтобы убрать автоматическую нумерацию рисунков
  22. \theoremstyle{plain}
  23. \newtheorem{theorem}{Теорема}[section]
  24. \newtheorem{lemma}{Лемма}[section]
  25. \newtheorem{propos}{Предложение}[section]
  26. \newtheorem{coll}{Следствие}[section]
  27. \newtheorem{statement}{Утверждение}[section]
  28.  
  29.  
  30.  
  31.  
  32. \begin{document}
  33.  
  34. \begin{titlepage}
  35. \begin{center}
  36.  
  37. Московский государственный университет им. М.В. Ломоносова
  38. \vspace{0.25cm}
  39.  
  40. Механико-математический факультет
  41.  
  42. Кафедра теоретической механики и мехатроники
  43. \vfill
  44.  
  45.  
  46.  
  47. \textsc{Отчет}\\[5mm]
  48.  
  49. {\LARGE По задаче 4 вычислительного практикума №2}
  50. \bigskip
  51.  
  52. Каргинова Екатерина Евгеньевна
  53. \vfill
  54. 4 курс, группа 422
  55. \end{center}
  56. \vfill
  57.  
  58.  
  59.  
  60. \begin{center}
  61. Москва, 2017 г.
  62. \end{center}
  63. \end{titlepage}
  64. \begin{center}
  65.  
  66. \end{center}
  67.  
  68.  
  69. \section{Постановка задачи.} \subsection{Линейный случай.}
  70. Рассмотрим уравнение в частных производных
  71. \begin{equation*}\frac{\partial u}{\partial t}+\frac{\partial F (u)}{\partial x }=0; \qquad t\in[0,1]\end{equation*}
  72.  
  73. Исследуем два случая: $F(u)$, называемая функцией состояния, линейна $F(u) = -\frac{1}{2}u$, и квадратична $F(u) = -\frac{1}{2}u^2$.
  74.  
  75. Таким образом, решаем уравнения:
  76. \begin{equation}\label{1}\frac{\partial u}{\partial t}-\frac{1}{2}\frac{\partial u}{\partial x }=0; \qquad t\in[0,1]\end{equation}
  77.  
  78.  
  79. \begin{equation}\label{nel}\frac{\partial u}{\partial t}-u\frac{\partial u}{\partial x }=0; \qquad t\in[0,1]\end{equation}
  80.  
  81.  
  82.  
  83. Для численного решения предлагаются следующие схемы:
  84.  
  85. $$ v_t + 0,5(1-sign(F_v'(v)))(F(v))_x+0,5(1+sign(F_v'(v)))(F(v))_{\bar x} = 0$$
  86. $$v_t +(F(\hat v))_{\overset{\circ}{x}}=0$$
  87. $$v_t +(F(\hat v))_{\overset{\circ}{x}}=\frac{1}{8}(v^{n}_{m+2}-4v^{n}_{m+1}+6v^{n}_{m}-4v^{n}_{m-1}+v^{n}_{m-2})$$
  88. Здесь используются обозначения разностных операторов:
  89. $$g_t = \frac{g^{n+1}_m-g^n_m}{\tau} \qquad g_x = \frac{v^{n}_{m+1}-v^n_m}{h} \qquad g_{\overset{\circ}{x}} =
  90. \frac{v^{n}_{m+1}-v^n_{m-1}}{2h}$$ и обозначение $$g^{n+1}_m = \hat g$$
  91. Тогда расчетные схемы в линейном случае имеют вид:
  92. \begin{equation}\label{2}\frac{v^{n+1}_m-v^n_m}{\tau} -\frac{v^{n}_{m+1}-v^{n}_m}{2h}=0\end{equation}
  93. \begin{equation}\label{3}\frac{v^{n+1}_m-v^n_m}{\tau} -\frac{v^{n+1}_{m+1}-v^{n+1}_{m-1}}{4h}=0\end{equation}
  94. \begin{equation}\label{4}\frac{v^{n+1}_m-v^n_m}{\tau}
  95. -\frac{v^{n+1}_{m+1}-v^{n+1}_{m-1}}{4h}=\frac{1}{8}(v^{n}_{m+2}-4v^{n}_{m+1}+6v^{n}_{m}-4v^{n}_{m-1}+v^{n}_{m-2})\end{equation}
  96.  
  97. В нелинейном случае:
  98.  
  99. \begin{equation}\label{nel1}\frac{v^{n+1}_m-v^n_m}{\tau} -\frac{(v^{n}_{m+1})^2-(v^{n}_m)^2}{2h}=0\end{equation}
  100. \begin{equation}\label{nel2}\frac{v^{n+1}_m-v^n_m}{\tau} -\frac{(v^{n+1}_{m+1})^2-(v^{n+1}_{m-1})^2}{4h}=0\end{equation}
  101. \begin{equation}\label{nel3}\frac{v^{n+1}_m-v^n_m}{\tau}
  102. -\frac{(v^{n+1}_{m+1})^2-(v^{n+1}_{m-1})^2}{4h}=\frac{1}{8}(v^{n}_{m+2}-4v^{n}_{m+1}+6v^{n}_{m}-4v^{n}_{m-1}+v^{n}_{m-2})\end{equation}
  103.  
  104.  
  105. Начальные условия задаются функцией
  106. \begin{equation*} u_0(x) = \left\{
  107. \begin{aligned}
  108. & 0 , \qquad x\le0\\
  109. &4x,\qquad 0<x\le0,25\\
  110. &1, \qquad x>0,25\\
  111. \end{aligned}\right. \end{equation*}
  112. %\begin{equation*} u_0(x) = \left\{
  113. %\begin{aligned}
  114. %& 1 , \qquad x\le0\\
  115. %&0, \qquad x>0,25\\
  116. %\end{aligned}\right. \end{equation*}
  117.  
  118.  
  119. Характеристическая система:
  120. \begin{equation}\label{char}\frac{d X}{d t} = -\frac{1}{2} \qquad \frac{d U}{d t} = 0\end{equation}. Обозначим ее решение
  121. (уравнение характеристик) за $x = X(t), u = U(t)$. Тогда характеристики имеют вид $x = -\frac{1}{2}t + C$, $u = u_0$, где $C = const$. Отсюда в силу постоянства решения на характеристиках заключаем, что достаточно рассматривать задачу в области $Q_T = \{(t,x) | 0 \le t \le 1; -1\le x \le 1\}$.
  122.  
  123. Граничные условия следующие: $$u(0,x) = u_0 \qquad u(t, -1) = 0 \qquad u(t, 1) = 1$$
  124.  
  125.  
  126.  
  127.  
  128.  
  129. \section{Точное решение.Линейный случай.}
  130.  
  131. Итак, решаем уравнение $\frac{\partial u}{\partial t}-\frac{1}{2}\frac{\partial u}{\partial x }=0 $ с начальными условиями, задаваемыми функцией $u(0,x) = u_0(x)$. \begin{equation*} u_0(x) = \left\{
  132. \begin{aligned}
  133. & 0 , \qquad x\le0\\
  134. &4x,\qquad 0<x\le0,25\\
  135. &1, \qquad x>0,25\\
  136. \end{aligned}\right. \end{equation*}
  137. Характеристическая система (\ref{char}) имеет решение $x = X(t) = -0,5t + x_0$, $u = U(t) = u_0$. Через каждую точку $(x,t)$ полуплоскости $t\ge 0$ проходит единственная характеристика $x = X(t, x_0, u_0(x_0))$, то есть уравнение характеристики однозначно разрешается относительно $x_0$.
  138.  
  139. Функция $u(t,x)$ постоянна вдоль характеристик, то есть вдоль каждой характеристики имеем $u(t,x)=u_0(x_0(t,x))$.
  140.  
  141. Таким образом, при $x\le -0,5t$ имеем $u(t, x) = 0$, а при $x>0,25-0,5t$ имеем $u(t,x) = 1$.
  142. В случае, когда $ -0,5t<x\le 0,25-0,5t$, имеем $x_0 = x+0,5t$. Таким образом $u(t,x) = u_0(x_0(t,x)) = 4x+2t.$
  143. Итак, решение имеет вид:
  144.  
  145. \begin{equation*} u(t,x) = \left\{
  146. \begin{aligned}
  147. & 0 , \qquad x\le -0,5t\\
  148. &4x+2t,\qquad -0,5t<x\le 0,25-0,5t\\
  149. &1, \qquad x>0,25-0,5t\\
  150. \end{aligned}\right. \end{equation*}
  151.  
  152. \section{Точное решение. Нелинейный случай.}
  153. решаем уравнение $\frac{\partial u}{\partial t}-u\frac{\partial u}{\partial x }=0 $ с начальными условиями, задаваемыми функцией $u(0,x) = u_0(x)$. \begin{equation*} u_0(x) = \left\{
  154. \begin{aligned}
  155. & 0 , \qquad x\le0\\
  156. &4x,\qquad 0<x\le0,25\\
  157. &1, \qquad x>0,25\\
  158. \end{aligned}\right. \end{equation*}
  159. Характеристическая система (\ref{char}) имеет решение $x = X(t) = -u(x_0(t)) t + x_0$, $u = U(t) = u_0$. Через каждую точку $(x,t)$ полуплоскости $t\ge 0$ проходит единственная характеристика $x = X(t, x_0, u_0(x_0))$, то есть уравнение характеристики однозначно разрешается относительно $x_0$.
  160.  
  161. Функция $u(t,x)$ постоянна вдоль характеристик, то есть вдоль каждой характеристики имеем $u(t,x)=u_0(x_0(t,x))$.
  162.  
  163. Таким образом, при $x\le 0$ имеем $u(t, x) = 0$, а при $x>-t+0,25$ имеем $u(t,x) = 1$.
  164. В случае, когда $ 0<x\le -t+0.25$, имеем $x_0 = \frac{4}{1-4t}$. Таким образом $u(t,x) = u_0(x_0(t,x)) = \frac{4x}{1-4t}.$
  165. Итак, решение имеет вид:
  166.  
  167. \begin{equation*} u(t,x) = \left\{
  168. \begin{aligned}
  169. & 0 , \qquad x\le 0\\
  170. &\frac{4x}{1-4t},\qquad 0<x\le -t+0.25\\
  171. &1, \qquad x>-t+0.25\\
  172. \end{aligned}\right. \end{equation*}
  173.  
  174.  
  175. \section{Теоретическое исследование разностных схем. Линейный случай.}
  176. \subsection{Аппроксимация.} \label{app}
  177. Найдем порядок аппроксимации уравнения (\ref{1}) каждой из разностных схем (\ref{2}), (\ref{3}), (\ref{4}) на гладких его решениях.
  178.  
  179. Пусть функция $u(t,x)$ является решением уравнения и дифференцируема по каждой из переменных при $t \ge 0$. Подставим проекцию на сетку
  180. этой функции в разностное уравнение (\ref{2}). Получим некоторую сеточную функцию $\phi$:
  181.  
  182. $$\phi^n_m = \frac{u^{n+1}_m-u^n_m}{\tau} -\frac{u^{n}_{m+1}-u^{n}_m}{2h}$$
  183. Эта функция может быть отлична от нуля. Решение задачи сводится к установлению порядка значений функции $\phi^n_m$ в зависимости от шагов сетки $\tau, h$ для гладких функций $u$.
  184.  
  185. Представим значение функции $u$ в узлах $(n+1,m)$ и $(n, m+1)$ в виде разложения в ряд Тейлора в окрестности точки $(n, m)$. Далее для сокращения записи используем обозначение $v^n_m = v$.
  186.  
  187. $$ u^{n+1}_m = u+\tau \dot u+\frac{\tau^2}{2} \ddot u + \overset{}{\underset{=}{O}}(\tau^3)$$
  188. $$ u^n_{m+1} = u+ h u' + \frac{h^2}{2} u'' + \overset{}{\underset{=}{O}}(h^3)$$
  189.  
  190.  
  191.  
  192. Подставим эти разложения в функцию $\phi$ и учтем, что $u$ - точное решение уравнения). Тогда получим: $$ \phi^n_m = \frac{u+\tau \dot u+\frac{\tau^2}{2} \ddot u + \overset{}{\underset{=}{O}}(\tau^3) - u}{\tau} -\frac{u+ h u' + \frac{h^2}{2} u'' + \overset{}{\underset{=}{O}}(h^3) - u}{2h}
  193. = \dot u -\frac{1}{2} u' + \frac{\tau}{2} \ddot u -\frac{h}{4} u'' +\overset{}{\underset{=}{O}}(\tau^2+h^2) = \overset{}{\underset{=}{O}}(\tau + h)$$.
  194.  
  195. Итак, разностная схема (\ref{2}) аппроксимирует с линейным порядком аппроксимации.
  196.  
  197. Проделаем аналогичную операцию со схемами (\ref{3}),(\ref{4}).
  198.  
  199. Для схемы \ref{3} сначала сдвинем индексы. Заменим индекс $n$ на индекс $p-1$. Тогда $u^{n+1}_m$ перепишется в виде $u^p_m$, а
  200. $u^{n}_m$ в виде $u^{p-1}_m$
  201.  
  202. Тогда:
  203. %\begin{multline}
  204. \begin{multline*}
  205. \phi^n_m = \frac{u^{p}_m-u^{p-1}_m}{\tau} -\frac{u^{p}_{m+1}-u^{p}_{m-1}}{4h}=\\
  206. = \frac{u - (u-\tau \dot u+\frac{\tau^2}{2} \ddot u -\frac{\tau^3}{6}\dddot u + \overset{}{\underset{=}{O}}(\tau^4))}{\tau} - \frac{ u+ h u' + \frac{h^2}{2} u'' +\frac{h^3}{6} u'''+ \overset{}{\underset{=}{O}}(h^4)-( u- h u' + \frac{h^2}{2} u'' -\frac{h^3}{6}u'''+ \overset{}{\underset{=}{O}}(h^4))}{4h} = \\
  207. = \dot u - \frac{1}{2}u' -\frac{\tau}{2}\ddot u -\frac{h^2}{12}u''' + \overset{}{\underset{=}{O}}(\tau^2+h^4) = = \dot u - \frac{1}{2}u'+
  208. \overset{}{\underset{=}{O}}(\tau+h^2)
  209. \end{multline*}
  210. Учитывая, что $u$ - точное решение уравнения (\ref{1}), получим, что схема (\ref{3}) аппроксимирует с порядком аппроксимации линейным по $\tau$ и квадратичным по $h$
  211.  
  212. Рассмотрим схему (\ref{4}). Нетрудно заметить, что левая часть имеет порядок аппроксимации такой же, как и в случае схемы (\ref{3}).
  213. Посмотрим, как ведет себя правая часть при разложении входящих в нее сеточных функций в ряд Тейлора в окрестности точки $(n,m)$
  214.  
  215. \begin{multline*}
  216. \frac{1}{8}\left(u^{n}_{m+2}-4u^{n}_{m+1}+6u^{n}_{m}-4u^{n}_{m-1}+u^{n}_{m-2}\right) =\\
  217. \frac{1}{8}\left( u+ 2h u' + \frac{(2h)^2}{2} u'' +\frac{(2h)^3}{6}u'''+\frac{(2h)^4}{24}u''''+ \overset{}{\underset{=}{O}}(h^5)\right)
  218. -4\left( u+ h u' + \frac{h^2}{2} u'' +\frac{h^3}{6}u'''+\frac{(h)^4}{24}u''''+ \overset{}{\underset{=}{O}}(h^5)\right)+\\+6u
  219. -4\left(u- h u' + \frac{h^2}{2} u'' -\frac{h^3}{6}u'''+ \frac{(h)^4}{24}u''''+\overset{}{\underset{=}{O}}(h^5)\right)+
  220. \left(u- 2h u' + \frac{(2h)^2}{2} u'' -\frac{(2h)^3}{6}u'''+ \frac{(2h)^4}{24}u''''+ \overset{}{\underset{=}{O}}(h^5)\right) = \\
  221. = \frac{1}{8}(h^4+\overset{}{\underset{=}{O}}(h^5)) = \overset{}{\underset{=}{O}}(h^4)
  222. \end{multline*}
  223. Нетрудно заметить, что при осуществлении аналогичной операции поиска порядка аппроксимации схемы (\ref{4}), получим, что эта схема, как и схема (\ref{3}), аппроксимирует уравнение (\ref{1}) с порядком аппроксимации линейным по $\tau$ и квадратичным по $h$.
  224.  
  225.  
  226.  
  227.  
  228. \subsection{Дифференциальное приближение.}
  229.  
  230. Выпишем дифференциальное приближение каждой из трех разностных схем (\ref{2}), (\ref{3}), (\ref{4}) с точностью до членов порядка $\overset{}{\underset{=}{O}}(\tau^3+h^3)$.
  231.  
  232. Пусть функция $v$ - гладкое решение уравнения (\ref{1}). Тогда в схему подставим вместо $v^{n+1}_m$ и $v^n+{m+1}$ их представления в виде рядов Тейлора в окрестности точки $(n, m)$ с точностью до членов порядка $\overset{}{\underset{=}{O}}(\tau^4)$ и
  233. $\overset{}{\underset{=}{O}}(h^4)$. Далее для сокращения записи используем обозначение $v^n_m = v$.
  234.  
  235. $$ \frac{v+\tau \dot v+\frac{\tau^2}{2} \ddot v +\frac{\tau^3}{6}\dddot v + \overset{}{\underset{=}{O}}(\tau^4) - v}{\tau} -\frac{v+ h v' + \frac{h^2}{2} v'' + \frac{h^3}{6}v'''+
  236. \overset{}{\underset{=}{O}}(h^4) - v }{2h} =0
  237. $$
  238.  
  239. После несложных преобразований, уравнение приводится к виду
  240. \begin{equation} \label{wuf}
  241. \dot v - \frac{1}{2}v' = -\frac{\tau}{2}\ddot v -\frac{\tau^2}{6}\dddot v +\frac{h}{4}v''+\frac{h^2}{12}v'''+ \overset{}{\underset{=}{O}}(\tau^3+h^3) \end{equation}
  242. В левой части равенства - иcходное волновое уравнение, а в правой погрешность аппроксимации, которая в общем случае отлична от нуля.
  243. Для более удобного анализа характера возмущений производные по времени заменим на пространственные производные с требуемой точностью.
  244. Для этого дифференцируем по времени (\ref{wuf})
  245.  
  246. \begin{equation}
  247. \ddot v - \frac{1}{2}\dot v' = -\frac{\tau}{2}\dddot v -\frac{\tau^2}{6}\ddddot v +\frac{h}{4}\dot v''+\frac{h^2}{12}\dot v'''+ \overset{}{\underset{=}{O}}(\tau^3+h^3)
  248. \end{equation}
  249.  
  250. Дифференцируем (\ref{wuf}) по пространственной переменной и разделим на 2
  251.  
  252. \begin{equation}
  253. \frac{\dot v'}{2} - \frac{1}{4} v'' = -\frac{\tau}{4}\ddot v' -\frac{\tau^2}{12}\dddot v' +\frac{h}{8} v'''+\frac{h^2}{24}v''''+
  254. \overset{}{\underset{=}{O}}(\tau^3+h^3) \end{equation}
  255.  
  256. Сложим два полученных уравнения и выразим $\ddot v$:
  257. \begin{equation} \ddot v = \frac{1}{4}v''-\frac{\tau}{2}\dddot v +\frac{h}{4}\dot v''-\frac{\tau}{4}\ddot v'+\frac{h}{8}v'''+ \overset{}{\underset{=}{O}}(\tau^2+h^2)
  258. \end{equation}
  259.  
  260. Аналогичными действиями получаем:
  261.  
  262. $$ \dot v''=\frac{1}{2}v'''+ \overset{}{\underset{=}{O}}(\tau^2+h^2) $$
  263. $$ \ddot v' = \frac{1}{4}v'''+ \overset{}{\underset{=}{O}}(\tau^2+h^2) $$
  264. $$\dddot v = \frac{1}{8}v'''+\overset{}{\underset{=}{O}}(\tau^2+h^2) $$
  265.  
  266. Тогда $$\ddot v = \frac{1}{4}v''-\frac{\tau}{16}v'''+\frac{h}{8}v'''-\frac{\tau}{16}v'''+\frac{h}{8}v''' = \frac{1}{4}v''+v'''(\frac{h}{4}-\frac{\tau}{8})$$
  267.  
  268. Подставляя полученные выражения в искомое, получим: $$\dot v - \frac{1}{2}v' = \frac{v''}{8}(2h - \tau) +\frac{v'''}{24}(2h^2 - 3\tau h + \tau^2)$$
  269.  
  270. В схеме (\ref{3}) сделаем замену индекса $n \longmapsto p-1$
  271.  
  272. $$ \frac{v^{p}_m-v^{p-1}_m}{\tau} -\frac{v^{p}_{m+1}-v^{p}_{m-1}}{4h}=0$$
  273.  
  274. Аналогично, считаем $v$ гладкой функцией, представим ее значения в узлах сетки в виде ряда Тейлора с центом в точке $(p, m)$.
  275. $$\frac{v - (v-\tau \dot v+\frac{\tau^2}{2} \ddot v -\frac{\tau^3}{6}\dddot v + \overset{}{\underset{=}{O}}(\tau^4))}{\tau}- \frac{ v+ h v' + \frac{h^2}{2} v'' +\frac{h^3}{6} v'''+
  276. \overset{}{\underset{=}{O}}(h^4)-( v- h v' + \frac{h^2}{2} v'' -\frac{h^3}{6}v'''+ \overset{}{\underset{=}{O}}(h^4))}{4h}
  277. $$
  278.  
  279. $$\dot v -\frac{1}{2}v' =\frac{\tau}{2}\ddot v - \frac{\tau^2}{6}\dddot v + \frac{h^2}{12}v'''+ \overset{}{\underset{=}{O}}(\tau^3+h^3)$$
  280.  
  281. Заменим производные по времени производными по координате
  282.  
  283. $$\dot v -\frac{1}{2}v' =\frac{\tau}{8}v''+\frac{v'''}{24}(2\tau^2+h^2) + \overset{}{\underset{=}{O}}(\tau^3+h^3)$$
  284.  
  285. Рассмотрим схему (\ref{4}). Заметим, что в правой части все линейные и квадратичные члены взаимоуничтожаются (это было подробно показано в пункте \ref{app}),
  286. поэтому дифференциальное приближение схемы (\ref{4}) такое же, как и для схемы (\ref{3}).
  287.  
  288. \subsection{Устойчивость.}
  289. Применим для исследования спектральный признак устойчивости.
  290.  
  291. \begin{statement}
  292. Пусть на сетке с узлами $(n,m)$ построена некоторая аппроксимация вида $$L_{\tau,h}v|_{n,m} = 0$$
  293. Выпишем все частные решения уравнения $L_{\tau,h}v_h = 0$, имеющие вид $v^n_m = \lambda(\varphi)^n e^{i m \varphi }$
  294. Если при заданном законе стремления шагов $\tau$ и $h$ к нулю существует $C<\infty$ такое, что $$|\lambda(\varphi)|\le e^{C \tau} \qquad \forall \varphi,$$
  295. то разностная схема может быть применена для численного решения соответствующей задачи Коши. В противном случае от применения разностной схемы следует воздержаться. \end{statement}
  296.  
  297. В каждую из схем (\ref{2}), (\ref{3}), (\ref{4}) подставим вместо $v^n_m$ выражение $\lambda(\varphi)^n e^{i m \varphi }$.
  298.  
  299. Тогда первая схема перепишется в виде:
  300.  
  301. $$ \frac{\lambda(\varphi)^{n+1} e^{i m \varphi } - \lambda(\varphi)^n e^{i m \varphi }} {\tau} - \frac{\lambda(\varphi)^n e^{i (m+1) \varphi }- \lambda(\varphi)^n e^{i m \varphi }}{2h} = 0$$
  302.  
  303. Разделим уравнение на $\lambda(\varphi)^n e^{i m \varphi }$. Тогда получим:
  304. $$ \frac{\lambda(\varphi) - 1} {\tau} - \frac{ e^{i m \varphi }- 1}{2h} = 0$$
  305.  
  306. Откуда получаем, что $\lambda(\varphi) = 1+ \frac{\tau}{2h}(e^{i \varphi}- 1)$.
  307.  
  308. Это график окружности с центром в точке $(1-\frac{\tau}{2h}, 0)$ и радиусом $\frac{\tau}{2h}$. Необходимое условие устойчивости будет выполнено, если окружность лежит целиком в единичной окружности с центром в начале координат.
  309.  
  310. Итак, выбор шага сетки (\ref{2}) стоит осуществлять так, чтобы выполнялось неравенство $0<\frac{\tau}{2 h}<1$.
  311.  
  312. Вторая схема перепишется в виде:
  313. $$ \frac{\lambda(\varphi)^{n+1} e^{i m \varphi } - \lambda(\varphi)^n e^{i m \varphi }} {\tau} - \frac{\lambda(\varphi)^{n+1} e^{i (m+1) \varphi }- \lambda(\varphi)^{n+1} e^{i (m-1) \varphi }}{4h} = 0$$
  314.  
  315. Разделим на $\lambda(\varphi)^{n+1}e^{i m \varphi}$, получим:
  316. $$ \frac{1 - \frac{1}{\lambda(\varphi)}} {\tau} - \frac{e^{i
  317. \varphi }- e^{-i \varphi }}{4h} = 0$$
  318.  
  319. Откуда получаем следующее выражение: $\lambda(\varphi) = \cfrac{1}{1 - i \cfrac{\tau}{2h}\sin{\varphi}} = \cfrac{1}{1 - i a\sin{\varphi}}$
  320. На комплексной плоскости число $1 - i \cfrac{\tau}{2h}\sin{\varphi}$ при любом значении $\varphi$ лежит на вертикальной прямой, проходящей через точку $(1,0)$. Тогда модуль этого числа больше единицы при любом $a$, следовательно, при любом соотношении между $\tau$ и $h$ число $\cfrac{1}{1 - i \cfrac{\tau}{2h}\sin{\varphi}}$ лежит внутри единичного круга с центром в начале координат, таким образом, необходимое условие устойчивости выполнено при любых $\tau$ и $h$.
  321.  
  322. Теперь рассмотрим схему (\ref{4}). При подстановке $v^n_m = \lambda(\varphi)^{n}e^{i m \varphi}$ получим следующее выражение:
  323. \begin{multline*}
  324. \frac{\lambda(\varphi)^{n+1} e^{i m \varphi } - \lambda(\varphi)^n e^{i m \varphi }} {\tau} - \frac{\lambda(\varphi)^{n+1} e^{i (m+1)
  325. \varphi }- \lambda(\varphi)^{n+1} e^{i (m-1) \varphi }}{4h} =\\
  326. = \frac{1}{8}(\lambda(\varphi)^{n}e^{i (m+2) \varphi}-4\lambda(\varphi)^{n}e^{i (m+1) \varphi}+6\lambda(\varphi)^{n}e^{i m \varphi}-4\lambda(\varphi)^{n}e^{i (m-1) \varphi}+\lambda(\varphi)^{n}e^{i (m-2) \varphi})
  327. \end{multline*}
  328.  
  329. Разделим на $\lambda(\varphi)^{n}e^{i m \varphi}$
  330.  
  331.  
  332. $$ \frac{\lambda(\varphi) -1} {\tau} - \frac{\lambda(\varphi)( e^{i
  333. \varphi }- e^{-i \varphi })}{4h}
  334. = \frac{1}{8}(e^{2i \varphi}-4e^{i \varphi}+6-4e^{-i \varphi}+e^{-2i \varphi})
  335. $$
  336.  
  337. Откуда получаем $$\lambda(\varphi) = \frac{\frac{\tau}{4}(\cos{2\varphi} - 4\cos{\varphi}+3)+1} {1-\frac{\tau}{2h}i \sin{\varphi}}$$
  338.  
  339. Заметим, что при любом выборе $\tau$ и $h$ для $\varphi = \pi$ получаем, что $\lambda(\pi) = 1+2\tau >1$. Значит, необходимое условие устойчивости не выполнено и схема неустойчива.
  340.  
  341. \section{Теоретическое исследование разностных схем. Нелинейный случай.}
  342. \subsection{Аппроксимация.}
  343.  
  344.  
  345. Установим порядок аппроксимации на гладких решениях уравнения (\ref{nel}) для схем (\ref{nel1}),(\ref{nel2}), (\ref{nel3}).
  346.  
  347. Пусть $u(t,x)$ - гладкое решение уравнения (\ref{nel}) и дважды дифференцируема по переменным $t$ и $x$. Тогда рассмотрим проекцию этой функции на ее сетку и подставим эту проекцию в каждую из схем.
  348.  
  349. Получим сеточную функцию $\phi^n_m$, имеющую вид:
  350. $$\phi^n_m = \frac{u^{n+1}_m-u^n_m}{\tau} -\frac{(u^{n}_{m+1})^2-(u^{n}_m)^2}{2h}$$
  351.  
  352. Эта функция в общем случае отлична от нуля. Найдем порядок этой функции в зависимости от шагов сетки $\tau$ и $h$.
  353.  
  354. Введем функцию $s(t,x) = u^2(t,x)$. Поскольку эта функция есть композиция двух гладких функций, дважды дифференцируемых по обоим своим переменным, то эту функцию можем разложить в ряд Тейлора в окрестности точки $(n,m)$.
  355.  
  356. Тогда, как и в линейном случае, подставим в $\phi^n_m$ представления входящих в нее значений функций в узлах в виде соответствующих рядов Тейлора с центром в точке $(n,m)$. Учитывая, что $u(t,x)$ - точное решение уравнения (\ref{nel}), получим:
  357.  
  358. \begin{multline*}\phi^n_m = \cfrac{u^{n}_m+\tau \dot n^n_m +\cfrac{\tau^2}{2}\ddot u^n_m + \overset{}{\underset{=}{O}}(\tau^3)-u^n_m}{\tau} -\cfrac{s^n_m+h s^n_m{}' +\cfrac{h^2}{2}s^n_m{}''+ \overset{}{\underset{=}{O}}(h^3)-u^n_m}{2h} =\\
  359. = \dot u^n_m + \frac{s^n_m{}'}{2}+ \overset{}{\underset{=}{O}}(\tau+h)
  360. = \frac{\partial u(n\tau, mh)}{\partial t} +\frac{\partial u^2(n\tau, mh)}{2 \partial x}+\overset{}{\underset{=}{O}}(\tau+h) = \overset{}{\underset{=}{O}}(\tau+h)
  361. \end{multline*}
  362.  
  363. Таким образом, схема (\ref{nel1}) имеет линейный порядок аппроксимации.
  364.  
  365. Со схемами (\ref{nel2}) и (\ref{nel3}) проделываем аналогичную операцию.
  366. В случае схемы (\ref{nel2}) заменим индекс $n \longmapsto p-1$, тогда схема перепишется в виде:
  367.  
  368. $$\frac{u^p_m-u^{p-1}_m}{\tau} -\frac{s^{p}_{m+1}-s^{p}_{m-1}}{4h}=0$$
  369.  
  370. После подстановки вместо значений функций в узлах сетки соответствующих рядов Тейлора, имеем:
  371.  
  372. \begin{multline*}
  373. \phi^n_m = \cfrac{u^p_m - (u^p_m-\tau \dot u^p_m+\cfrac{\tau^2}{2} \ddot u^p_m -\cfrac{\tau^3}{6}\dddot u^p_m + \overset{}{\underset{=}{O}}(\tau^4))}{\tau}- \\
  374. - \cfrac{ s^p_m+ h s^p_m{}' + \cfrac{h^2}{2} s^p_m{}'' +\cfrac{h^3}{6} s^p_m{}'''+ \overset{}{\underset{=}{O}}(h^4)-( s^p_m- h s^p_m{}' + \cfrac{h^2}{2} s^p_m{}''
  375. -\cfrac{h^3}{6}s^p_m{}'''+ \overset{}{\underset{=}{O}}(h^4))}{4h} = \\
  376. = \dot u^p_m - \frac{1}{2}s^p_m{}' -\frac{\tau}{2}\ddot u^p_m -\frac{h^2}{12}s^p_m{}''' + \overset{}{\underset{=}{O}}(\tau^2+h^4)
  377. =\\
  378. = \frac{\partial u(n\tau,mh)}{\partial t} - \frac{\partial u(n\tau,mh)^2}{2\partial x}+
  379. \overset{}{\underset{=}{O}}(\tau+h^2)
  380. \end{multline*}
  381.  
  382. С учетом того, что $u(t,x)$ - точное решение уравнения (\ref{nel}), получим, что схема (\ref{nel2}) имеет порядок аппроксимации линейный по $\tau$ и квадратичный по $h$.
  383.  
  384. Рассмотрим схему (\ref{nel3}). Нетрудно заметить, что левая ее часть имеет такой же порядок аппроксимации, как и схема (\ref{nel2}). Тогда подсчет порядка аппроксимации правой части этой схемы в точности повторяет подсчет порядка аппроксимации правой части схемы (\ref{4}), подробно изложенный в пункте \ref{app}. Как и в том случае, получим, что правая часть имеет порядок аппроксимации $\overset{}{\underset{=}{O}}(h^4)$, поэтому вся схема (\ref{nel3}) имеет порядок аппроксимации такой же, как и схема (\ref{nel2}): линейный по $\tau$ и квадратичный по $h$.
  385. \section{Численное решение. Линейный случай.}
  386.  
  387. Решаем уравнение
  388. $\frac{\partial u}{\partial t}-\frac{1}{2}\frac{\partial u}{\partial x }=0$ в области $Q_T = \{(t,x) | 0 \le t \le
  389. 1; -1\le x \le 1\}$
  390. при помощи явной схемы
  391.  
  392. $$\frac{v^{n+1}_m-v^n_m}{\tau} -\frac{v^{n}_{m+1}-v^{n}_m}{2h}=0$$
  393. и двух неявных схем
  394.  
  395. $$\frac{v^{n+1}_m-v^n_m}{\tau}
  396. -\frac{v^{n+1}_{m+1}-v^{n+1}_{m-1}}{4h}=\frac{\omega}{8}(v^{n}_{m+2}-4v^{n}_{m+1}+6v^{n}_{m}-4v^{n}_{m-1}+v^{n}_{m-2}), \qquad \omega = 0, 1$$.
  397.  
  398.  
  399.  
  400.  
  401. Начальные и граничные условия: $$\cfrac{\mbox{числитель}}{\mbox{знаменатель}}$$
  402.  
  403. $$u(0,x) = u_0 \qquad u(t, -1) = 0 \qquad u(t, 1) = 1$$
  404.  
  405. \begin{equation*} u_0(x) = \left\{
  406. \begin{aligned}
  407. & 0 , \qquad x\le0\\
  408. &4x,\qquad 0<x\le0,25\\
  409. &1, \qquad x>0,25\\
  410. \end{aligned}\right. \end{equation*}
  411.  
  412. Расчет по неявным схемам реализуется следующим образом. Точки на $n+1$ слое образуют вектор $\hat v = (v^{n+1}_1, v^{n+1}_2, \dots ,v^{n+1}_m)^T$, и он находится как решение следующей системы.
  413.  
  414.  
  415. \begin{equation*}\label{lag}\left\{
  416. \begin{aligned}
  417. &v^{n+1}_2-v^n_2-\frac{\tau}{4h}(v^{n+1}_{3}-v^{n+1}_{1}) = 0\\
  418. &v^{n+1}_3-v^n_3-\frac{\tau}{4h}(v^{n+1}_{4}-v^{n+1}_{2}) = 0\\
  419. &\vdots\\
  420. &v^{n+1}_{m-1}-v^n_{m-1}-\frac{\tau}{4h}(v^{n+1}_{m}-v^{n+1}_{m-2}) = 0\\
  421. \end{aligned}\right. \end{equation*}
  422.  
  423. Или, учитывая граничные условия $v^n_1 = 0$ и $v^n_m = 1$ $\forall n$,
  424.  
  425. \begin{equation*}\begin{pmatrix}
  426. 1 &-\frac{\tau}{4h}& 0 &0 &\ldots & 0\\
  427. \frac{\tau}{4h}& 1 &-\frac{\tau}{4h}&0 &\ldots & 0\\
  428. \vdots & \vdots &\ddots &-\ddots&\ddots &\vdots\\
  429. 0 & \ldots &0 &\frac{\tau}{4h} &1 &-\frac{\tau}{4h}\\
  430. 0 & \ldots &0 &0 &-\frac{\tau}{4h} &1
  431. \end{pmatrix}
  432. \begin{pmatrix}
  433. v^{n+1}_2\\
  434. v^{n+1}_3\\
  435. \vdots\\
  436. v^{n+1}_{m-2}\\
  437. v^{n+1}_{m-1}
  438. \end{pmatrix}
  439. =
  440. \begin{pmatrix}
  441. v^n_2\\
  442. v^n_3\\
  443. \vdots\\
  444. v^n_{m-2}\\
  445. v^n_{m-1} + \frac{\tau}{4h}
  446. \end{pmatrix}
  447. \end{equation*}
  448.  
  449. ~\
  450. Аналогично, выпишем систему для неявной схемы при значении параметра $\omega = 1$.
  451.  
  452. \begin{equation*}\begin{pmatrix}
  453. 1 &-\cfrac{\tau}{4h}& 0 &0 &\ldots & 0\\
  454. \cfrac{\tau}{4h}& 1 &-\cfrac{\tau}{4h}&0 &\ldots & 0\\
  455. 0& \cfrac{\tau}{4h}& 1 &-\cfrac{\tau}{4h}& \ldots & 0\\
  456. \vdots & \vdots &\ddots &\ddots&\ddots &\vdots\\
  457. 0 & \ldots &0 &\cfrac{\tau}{4h} &1 &-\cfrac{\tau}{4h}\\
  458. 0 & \ldots &0 &0 &\cfrac{\tau}{4h} &1
  459. \end{pmatrix}
  460. \begin{pmatrix}
  461. v^{n+1}_2\\
  462. ~\ \\
  463. v^{n+1}_3\\
  464. \\
  465. v^{n+1}_4\\
  466. \cfrac{\vdots}{}\\
  467.  
  468.  
  469. v^{n+1}_{m-2}\\
  470. \\
  471. v^{n+1}_{m-1}
  472. \end{pmatrix}
  473. =
  474. \begin{pmatrix}
  475. v^n_2 + \cfrac{\tau}{8}(v^{n}_{3}-4v^{n}_{2}+6v^{n}_{1})\\
  476. v^n_3+ \cfrac{\tau}{8}(v^{n}_{4}-4v^{n}_{3}+6v^{n}_{2}-4v^{n}_{1})\\
  477. v^n_4+\cfrac{\tau}{8}(v^{n}_{6}-4v^{n}_{5}+6v^{n}_{4}-4v^{n}_{3}+v^{n}_{2})\\
  478. \vdots\\
  479. v^n_{m-2}+\cfrac{\tau}{8}(1-4v^{n}_{m-1}+6v^{n}_{m-2}-4v^{n}_{m-3}+v^{n}_{m-4})\\
  480. v^n_{m-1}+ \cfrac{\tau}{4h}+ \cfrac{\tau}{8}(-3 + 6v^{n}_{m-1}-4v^{n}_{m-2}+v^{n}_{m-3})
  481. \end{pmatrix}
  482. \end{equation*}
  483.  
  484.  
  485. Далее используем обозначения $||v||_{C_h} = \max\limits_{x_i\in\omega_h} {|v_i|}$, $||v||_{L_{1,h}} = h \sum\limits_{x_i \in\omega_h} |v_i|$ \qquad (здесь $\omega_h$ - сетка),
  486.  
  487. $\Delta (v)_\alpha = ||v - u||_\alpha$, $\delta (v)_\alpha = \cfrac{||v-u||_\alpha}{||v||_\alpha}$ (здесь $u$ - точное решение)
  488.  
  489.  
  490.  
  491. \newpage
  492. \begin{center}
  493. \begin{tabular}[t]{|c|c|c|c|c|c|}
  494. \hline
  495. $\tau$ &$h$ & $\Delta (v)_{C_h}$ & $\Delta (v)_{L_{1,h}}$& $\delta (v)_{C_h}$ & $\delta (v)_{L_{1,h}}$ \\
  496. \hline
  497. 0.1 & 0.1 & 2.375000e-001 & 4.730469e-001 & 2.375000e-001 & 4.320063e-002\\
  498. \hline
  499. 0.01 & 0.1 & 3.068273e-001 & 6.972739e+000 & 3.068273e-001 & 6.502012e-002 \\
  500. \hline
  501. 0.001 & 0.1 & 3.121874e-001 & 7.151915e+001 & 3.121874e-001 & 6.683354e-002\\
  502. \hline
  503. 0.1 & 0.01 & 9.312502e+006 & 7.845286e+005 & 9.999999e-001 & 9.999895e-001\\
  504. \hline
  505. 0.01 & 0.01 & 7.958923e-002 & 4.996756e-001 & 7.958923e-002 & 4.451453e-003\\
  506. \hline
  507. 0.001 & 0.01 & 1.097955e-001 & 9.369601e+000 & 1.097955e-001 & 8.363848e-003\\
  508. \hline
  509. 0.1 & 0.001 & 2.471298e+016 & 1.827221e+014 & 1.000000e+000 & 1.000000e+000\\
  510. \hline
  511. 0.01 & 0.001 & 2.142381e+091 & 5.976315e+089 & 1.000000e+000 & 1.000000e+000\\
  512. \hline
  513. 0.001 & 0.001 & 2.522502e-002 & 5.000000e-001 & 2.522502e-002 & 4.445432e-004\\
  514. \hline
  515. \end{tabular}%\end{center}
  516. {\center Нормы погрешности расчетов: явная схема}
  517. %\end{table}
  518.  
  519. ~\
  520.  
  521. %\begin{table}[h!]\begin{center}
  522. \begin{tabular}[t]{|c|c|c|c|c|c|}
  523. \hline
  524. $\tau$ &$h$ & $\Delta (v)_{C_h}$ & $\Delta (v)_{L_{1,h}}$& $\delta (v)_{C_h}$ & $\delta (v)_{L_{1,h}}$ \\
  525. \hline
  526. 0.1 & 0.1 & 4.394995e-001 & 2.014531e-001 & 4.213793e-001 & 1.724507e-001\\
  527. \hline
  528. 0.01 & 0.1 & 4.996988e-001 & 1.838241e-001 & 4.352544e-001 & 1.513776e-001 \\
  529. \hline
  530. 0.001 & 0.1 & 5.102556e-001 & 1.876195e-001 & 4.359180e-001 & 1.539161e-001\\
  531. \hline
  532. 0.1 & 0.01 & 2.669414e-001 & 1.277566e-001 & 2.669190e-001 & 9.753321e-002\\
  533. \hline
  534. 0.01 & 0.01 & 1.029742e-001 & 2.142223e-002 & 1.028922e-001 & 1.580977e-002\\
  535. \hline
  536. 0.001 & 0.01 & 6.186953e-002 & 1.436181e-002 & 6.065315e-002 & 1.056404e-002\\
  537. \hline
  538. 0.1 & 0.001 & 2.479751e-001 & 1.258892e-001 & 2.479748e-001 & 9.512731e-002\\
  539. \hline
  540. 0.01 & 0.001 & 8.178958e-002 & 1.524753e-002 & 8.178893e-002 & 1.114178e-002\\
  541. \hline
  542. 0.001 & 0.001 & 2.733776e-002 & 2.311549e-003 & 2.733707e-002 & 1.683575e-003\\
  543. \hline
  544. \end{tabular}
  545. {\center Нормы погрешности расчетов: неявная схема, параметр $\omega = 0$}%\end{center}
  546. %\end{table}
  547.  
  548. ~\
  549.  
  550. %\begin{table}[h!]\begin{center}
  551. \begin{tabular}[t]{|c|c|c|c|c|c|}
  552. \hline
  553. $\tau$ &$h$ & $\Delta (v)_{C_h}$ & $\Delta (v)_{L_{1,h}}$& $\delta (v)_{C_h}$ & $\delta (v)_{L_{1,h}}$ \\
  554. \hline
  555. 0.1 & 0.1 & 5.971631e-001 & 3.870527e-001 & 4.286385e-001 & 3.284659e-001\\
  556. \hline
  557. 0.01 & 0.1 & 8.353191e-001 & 5.177420e-001 & 4.551356e-001 & 4.186652e-001\\
  558. \hline
  559. 0.001 & 0.1 & 8.990383e-001 & 5.451273e-001 & 4.734177e-001 & 4.386622e-001\\
  560. \hline
  561. 0.1 & 0.01 & 2.674390e-001 & 1.286974e-001 & 2.669338e-001 & 9.834490e-002\\
  562. \hline
  563. 0.01 & 0.01 & 1.102242e-001 & 2.667720e-002 & 1.069560e-001 & 1.969491e-002\\
  564. \hline
  565. 0.001 & 0.01 & 1.390244e-001 & 3.100314e-002 & 1.228700e-001 & 2.281165e-002\\
  566. \hline
  567. 0.1 & 0.001 & 2.479783e-001 & 1.259626e-001 & 2.479781e-001 & 9.518849e-002\\
  568. \hline
  569. 0.01 & 0.001 & 8.178958e-002 & 1.525727e-002 & 8.178484e-002 & 1.114899e-002\\
  570. \hline
  571. 0.001 & 0.001 & 2.733772e-002 & 2.316349e-003 & 2.733265e-002 & 1.687077e-003\\
  572. \hline
  573. \end{tabular}
  574. {\center Нормы погрешности расчетов: неявная схема, параметр $\omega = 1$}%\end{center}
  575. %\end{table}
  576.  
  577. \end{center}
  578.  
  579. \subsection{Измельчение сетки, линейный случай.}
  580.  
  581. Пусть $\tau_k = \frac{\tau}{2^k}$, $h_k = \frac{h}{2^k}$, где $ k = 1, 2, 3, 4$. Нормы погрешности становятся меньше при измельчении сетки, что согласуется с теорией.
  582. \begin{center}
  583. \begin{tabular}[t]{|c|c|c|c|c|c|}
  584. \hline
  585. $\tau_k$ &$h_k$ & $\Delta (v)_{C_h}$ & $\Delta (v)_{L_{1,h}}$& $\delta (v)_{C_h}$ & $\delta (v)_{L_{1,h}}$ \\
  586. \hline
  587. 0.05 & 0.05 & 1.747131e-001 & 4.507751e-002 & 1.747131e-001 & 3.339075e-002\\
  588. \hline
  589. 0.025 & 0.025 & 1.253250e-001 & 2.411326e-002 & 1.253250e-001 & 1.769780e-002 \\
  590. \hline
  591. 0.0125 & 0.0125 & 8.892778e-002 & 1.243342e-002 & 8.892778e-002 & 9.083777e-003\\
  592. \hline
  593. 0.00625 & 0.00625 & 6.297983e-002 & 6.248743e-003 & 6.297983e-002 & 4.554892e-003\\
  594. \hline
  595. \end{tabular}%\явный случай
  596. {\center Нормы погрешности при измельчении сетки $\tau = 0.1, h = 0.1$. Явная схема.}
  597.  
  598. ~\
  599.  
  600. \begin{tabular}[t]{|c|c|c|c|c|c|}
  601. \hline
  602. $\tau_k$ &$h_k$ & $\Delta (v)_{C_h}$ & $\Delta (v)_{L_{1,h}}$& $\delta (v)_{C_h}$ & $\delta (v)_{L_{1,h}}$ \\
  603. \hline
  604. 0.05 & 0.05 & 3.232995e-001 & 1.034413e-001 & 3.178564e-001 & 8.114104e-002\\
  605. \hline
  606. 0.025 & 0.025 & 1.869422e-001 & 5.235464e-002 & 1.856677e-001 & 3.951274e-002\\
  607. \hline
  608. 0.0125 & 0.0125 & 1.185948e-001 & 2.661411e-002 & 1.184491e-001 & 1.971416e-002\\
  609. \hline
  610. 0.00625 & 0.00625 & 7.718552e-002 & 1.359971e-002 & 7.715506e-002 & 9.981435e-003\\
  611. \hline
  612. \end{tabular}%\неявный случай, омега 0
  613. {\center Нормы погрешности при измельчении сетки $\tau = 0.1, h = 0.1$. Неявная схема, параметр $\omega = 0$.}
  614.  
  615. ~\
  616.  
  617. \begin{tabular}[t]{|c|c|c|c|c|c|}
  618. \hline
  619. $\tau_k$ &$h_k$ & $\Delta (v)_{C_h}$ & $\Delta (v)_{L_{1,h}}$& $\delta (v)_{C_h}$ & $\delta (v)_{L_{1,h}}$ \\
  620. \hline
  621. 0.05 & 0.05 & 6.953890e-001 & 3.088180e-001 & 4.521346e-001 & 2.451652e-001\\
  622. \hline
  623. 0.025 & 0.025 & 3.134525e-001 & 9.455170e-002 & 2.640357e-001 & 7.151918e-002\\
  624. \hline
  625. 0.0125 & 0.0125 & 1.410937e-001 & 3.576457e-002 & 1.341519e-001 & 2.650691e-002\\
  626. \hline
  627. 0.00625 & 0.00625 & 7.718227e-002 & 1.459929e-002 & 7.667829e-002 & 1.071654e-002\\
  628. \hline
  629. \end{tabular}
  630. {\center Нормы погрешности при измельчении сетки $\tau = 0.1, h = 0.1$. Неявная схема, параметр $\omega = 1$.}%\неявный случай, омега 1
  631. \end{center}
  632. \subsection{Графики, линейный случай.}
  633.  
  634. \begin{figure}[h!]
  635. \center{\includegraphics[height=60mm]{yav.png} \\График численного решения для $\tau = h = 0.001$. Явная схема.}
  636. \end{figure}
  637.  
  638. \begin{figure}[h!]
  639. \begin{minipage}[h]{0.49\linewidth}\center{\includegraphics[height=55mm]{del_yav_0.png} \\График модуля разности точного и численного решений для $\tau = h = 0.001$. Явная схема.}
  640. \end{minipage}
  641. \hfill
  642. \begin{minipage}[h]{0.49\linewidth}\center{\includegraphics[height=55mm]{del_yav_1.png} \\Деталь графика модуля разности точного и численного решений для $\tau = h = 0.001$. Явная схема.}
  643. \end{minipage}
  644. \end{figure}
  645. \newpage
  646. \begin{figure}[h!]
  647. \center{\includegraphics[height=60mm]{w0.png} \\График численного решения для $\tau = h = 0.001$. Неявная схема, параметр $\omega = 0$.}
  648. %\end{figure}
  649.  
  650. %\begin{figure}[h!]
  651. \begin{minipage}[h]{0.49\linewidth}\center{\includegraphics[height=55mm]{delw0_0.png} \\График модуля разности точного и численного решений для $\tau = h = 0.001$. Неявная схема, параметр $\omega = 0$.}
  652. \end{minipage}
  653. \hfill
  654. \begin{minipage}[h]{0.49\linewidth}\center{\includegraphics[height=55mm]{delw0_1.png} \\Деталь графика модуля разности точного и численного решений для $\tau = h = 0.001$. Неявная схема, параметр $\omega = 0$.}
  655. \end{minipage}
  656. \end{figure}
  657.  
  658. \newpage
  659.  
  660. \begin{figure}[h!]
  661. \center{\includegraphics[height=60mm]{w1.png} \\График численного решения для $\tau = h = 0.001$. Неявная схема, параметр $\omega = 1$.}
  662. %\end{figure}
  663. %\begin{figure}[h!]
  664. \begin{minipage}[h]{0.49\linewidth}\center{\includegraphics[height=55mm]{delw1_0.png} \\График модуля разности точного и численного решений для $\tau = h = 0.001$. Неявная схема, параметр $\omega = 1$.}
  665. \end{minipage}
  666. \hfill
  667. \begin{minipage}[h]{0.49\linewidth}\center{\includegraphics[height=55mm]{delw1_1.png} \\Деталь графика модуля разности точного и численного решений для $\tau = h = 0.001$. Неявная схема, параметр $\omega = 1$.}
  668. \end{minipage}
  669. \end{figure}
  670.  
  671. \newpage
  672.  
  673. \subsection{Нелинейный случай: численное решение.}
  674. %\begin{table}
  675. Решаем уравнение
  676. $\frac{\partial u}{\partial t}-u\frac{\partial u}{\partial x }=0$ в области $Q_T = \{(t,x) | 0 \le t \le
  677. 1; -1\le x \le 1\}$
  678. при помощи явной схемы
  679.  
  680. $$\frac{v^{n+1}_m-v^n_m}{\tau} -\frac{(v^{n}_{m+1})^2-(v^{n}_m)^2}{2h}=0$$
  681. и двух неявных схем
  682.  
  683. $$\frac{v^{n+1}_m-v^n_m}{\tau}
  684. -\frac{(v^{n+1}_{m+1})^2-(v^{n+1}_{m-1})^2}{4h}=\frac{\omega}{8}(v^{n}_{m+2}-4v^{n}_{m+1}+6v^{n}_{m}-4v^{n}_{m-1}+v^{n}_{m-2}), \qquad \omega = 0, 1$$.
  685.  
  686.  
  687.  
  688.  
  689. Начальные и граничные условия: $$u(0,x) = u_0 \qquad u(t, -1) = 0 \qquad u(t, 1) = 1$$
  690.  
  691. \begin{equation*} u_0(x) = \left\{
  692. \begin{aligned}
  693. & 1 , \qquad x\le0\\
  694. &0, \qquad x>0\\
  695. \end{aligned}\right. \end{equation*}
  696.  
  697. Расчет по явным схемам реализуется аналогично линейному случаю, а именно: каждая точка на следующем слое явно выражается через некоторую комбинацию точек на предыдущем слое.
  698.  
  699. Расчет по неявным схемам реализуется следующим образом. Точки на $n+1$ слое образуют вектор $\hat v = (v^{n+1}_0, v{n+1}_2, \dots , v{n+1}_m)^T$, и он находится как решение следующей системы.
  700.  
  701. \begin{equation*} \left\{
  702. \begin{aligned}
  703. &G_1(\hat v) = \hat v_1 - v_1 - a((\hat v_2)^2 - (\hat v_0)^2) = 0,\\
  704. &G_2(\hat v) = \hat v_2 - v_2 - a((\hat v_3)^2 - (\hat v_2)^2) = 0,\\
  705. &\vdots\\
  706. &G_{m-1}(\hat v) = \hat v_{m-1} - v_{m-1} - a((\hat v_{m})^2 - (\hat v_{m-2})^2) = 0,\\ \end{aligned}\right.
  707. \end{equation*}
  708. Здесь $v_i$ - точки на предыдущем слое, $a = \frac{\tau}{2h}$.
  709.  
  710. Используем метод Ньютона для решения этой системы.
  711. Обозначим как $J(\hat v)$ матрицу Якоби оператора $G$.
  712. Из разложения в ряд Тейлора находим $G(\hat v^k) + J(\hat v) (\hat v^{k+1} - \hat v_{k}) = 0$, за начальное приближение $\hat v^0$ берем $v$.
  713.  
  714.  
  715. \begin{multline*}\begin{pmatrix}
  716. \hat v_1 - v_1 - a((\hat v_2)^2 - (\hat v_0)^2)\\
  717. \hat v_2 - v_2 - a((\hat v_3)^2 - (\hat v_2)^2) = 0\\
  718. \vdots\\
  719. \hat v_{m-2} - v_{m-2} - a((\hat v_{m-1})^2 - (\hat v_{m-3})^2)\\
  720. \hat v_{m-1} - v_{m-1} - a((\hat v_{m})^2 - (\hat v_{m-2})^2)
  721. \end{pmatrix}+\\
  722. +\begin{pmatrix}
  723. 1 &-2 a \hat v^k_2& 0 &0 &\ldots & 0\\
  724. 2 a \hat v^k_1& 1 &-2 a \hat v^k_3&0 &\ldots & 0\\
  725. \vdots & \vdots &\ddots &\ddots&\ddots &\vdots\\
  726. 0 & \ldots &0 &2 a \hat v^k_{m-3} &1 &-2 a \hat v^k_{m-1}\\
  727. 0 & \ldots &0 &0 &2 a \hat v^k_{m-2} &1
  728. \end{pmatrix}
  729. \begin{pmatrix}
  730. \hat v^{k+1}_1 - \hat v^k_1\\
  731. \hat v^{k+1}_2 - \hat v^k_2\\
  732. \vdots\\
  733. \hat v^{k+1}_{m-2} - \hat v^k_{m-2}\\
  734. \hat v^{k+1}_{m-1} - \hat v^k_{m-1}
  735. \end{pmatrix}
  736. = 0
  737. \end{multline*}
  738.  
  739.  
  740. При расчетах учитываем граничные условия $v_0 = 0$ и $v_m = 1$.
  741.  
  742. ~\
  743. Аналогично, выпишем систему для неявной схемы при значении параметра $\omega = 1$.
  744.  
  745. \begin{equation*} \left\{
  746. \begin{aligned}
  747. &G_1(\hat v) = \hat v_1 - v_1 - a((\hat v_2)^2 - (\hat v_0)^2) - v_1 + \cfrac{\tau}{8}(v_{3}-4v_{2}+6v_{1}-4v_{0}+1) = 0,\\
  748. &G_2(\hat v) = \hat v_2 - v_2 - a((\hat v_3)^2 - (\hat v_2)^2) - v_2 + \cfrac{\tau}{8}(v_{4}-4v_{3}+6v_{2}-4v_{1}+v_0)= 0,\\
  749. &\vdots\\
  750. &G_{m-2}(\hat v) = \hat v_{m-2} - v_{m-2} - a((\hat v_{m-1})^2 - (\hat v_{m-3})^2) - v_{m-2} + \cfrac{\tau}{8}(v_{m}-4v_{m-1}+6v_{m-2}-4v_{m-3}+v_{m-4})= 0,\\
  751. &G_{m-1}(\hat v) = \hat v_{m-1} - v_{m-1} - a((\hat v_{m})^2 - (\hat v_{m-2})^2) - v_{m-1} + \cfrac{\tau}{8}(-4v_{m}+6v_{m-1}-4v_{m-2}+v_{m-3})= 0.\\ \end{aligned}\right.
  752. \end{equation*}
  753.  
  754. И решаем ее методом Ньютона.
  755.  
  756. Далее используем обозначения $||v||_{C_h} = \max\limits_{x_i\in\omega_h} {|v_i|}$, $||v||_{L_{1,h}} = h \sum\limits_{x_i \in\omega_h} |v_i|$ \qquad (здесь $\omega_h$ - сетка),
  757.  
  758. $\Delta (v)_\alpha = ||v - u||_\alpha$, $\delta (v)_\alpha = \cfrac{||v-u||_\alpha}{||v||_\alpha}$ (здесь $u$ - точное решение)
  759. \begin{center}
  760. \begin{tabular}[t]{|c|c|c|c|c|c|}
  761. \hline
  762. $\tau$ &$h$ & $\Delta (v)_{C_h}$ & $\Delta (v)_{L_{1,h}}$& $\delta (v)_{C_h}$ & $\delta (v)_{L_{1,h}}$ \\
  763. \hline
  764. 0.1 & 0.1 & 5.000000e-001 & 5.585476e-001 & 5.000000e-001 & 7.704105e-002\\
  765. \hline
  766. 0.01 & 0.1 & 3.427265e-001 & 5.951159e+000 & 3.427265e-001 & 7.864800e-002\\
  767. \hline
  768. 0.001 & 0.1 & 3.342380e-001 & 5.999320e+001 & 3.342380e-001 & 7.893810e-002\\
  769. \hline
  770. 0.0001 & 0.1 & 3.334235e-001 & 6.004013e+002 & 3.334235e-001 & 7.896525e-002\\
  771. \hline
  772. 0.01 & 0.01 & 5.000000e-001 & 9.028638e-001 & 5.000000e-001 & 1.207845e-002\\
  773. \hline
  774. 0.001 & 0.01 & 3.427265e-001 & 1.559717e+001 & 3.427265e-001 & 2.079356e-002\\
  775. \hline
  776. 0.0001 & 0.01 & 3.342380e-001 & 1.629826e+002 & 3.342380e-001 & 2.172037e-002\\
  777. \hline
  778. 0.001 & 0.001 & 5.000000e-001 & 1.409895e+000 & 5.000000e-001 & 1.880487e-003\\
  779. \hline
  780. 0.0001 & 0.001 & 3.427265e-001 & 2.897675e+001 & 3.427265e-001 & 3.863633e-003\\
  781. \hline
  782. \end{tabular}
  783. {\center Нормы погрешности расчетов: явная схема.}
  784.  
  785. ~\
  786.  
  787. \begin{tabular}[t]{|c|c|c|c|c|c|}
  788. \hline
  789. $\tau$ &$h$ & $\Delta (v)_{C_h}$ & $\Delta (v)_{L_{1,h}}$& $\delta (v)_{C_h}$ & $\delta (v)_{L_{1,h}}$ \\
  790. \hline
  791. 0.1 & 0.1 & 1.170820e+000 & 2.747209e+000 & 1.170834e+000 & 4.002677e-001\\
  792. \hline
  793. 0.01 & 0.1 & 1.408035e+000 & 3.244944e+001 & 1.408035e+000 & 4.171762e-001\\
  794. \hline
  795. 0.001 & 0.1 & 1.434434e+000 & 3.320048e+002 & 1.425572e+000 & 4.200893e-001\\
  796. \hline
  797. 0.1 & 0.01 & 5.092824e-001 & 9.990538e-001 & 5.092824e-001 & 1.391219e-001\\
  798. \hline
  799. 0.01 & 0.01 & 1.170820e+000 & 4.532761e+000 & 1.170820e+000 & 6.085286e-002\\
  800. \hline
  801. 0.001 & 0.01 & 1.332313e+000 & 6.443278e+001 & 1.332313e+000 & 8.249002e-002\\
  802. \hline
  803. 0.1 & 0.001 & 6.209493e-001 & 8.908931e-001 & 6.209493e-001 & 1.229899e-001\\
  804. \hline
  805. 0.01 & 0.001 & 5.092783e-001 & 2.824843e+000 & 5.092783e-001 & 3.781270e-002\\
  806. \hline
  807. 0.001 & 0.001 & 1.170820e+000 & 1.006486e+001 & 1.170820e+000 & 1.343125e-002\\
  808. \hline
  809. \end{tabular}
  810. {\center Нормы погрешности расчетов: неявная схема, параметр $\omega = 0$}
  811. \end{center}
  812. %\end{table}
  813.  
  814. \end{document}
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement
Public Pastes
Advertisement
We use cookies for various purposes including analytics. By continuing to use Pastebin, you agree to our use of cookies as described in the Cookies Policy.  OK, I Understand
Not a member of Pastebin yet?
Sign Up, it unlocks many cool features!