Лекции и конспекты по математике Примеры решения задч


Используя оценку остаточного члена в форме Лагранжа, можно провести анализ погрешности в формулах приближённого дифференцирования, предполагая шаг $ h$ малым.

Пусть функция $ f(x)$ разложена по формуле Тейлора, с остаточным членом в форме Лагранжа, в точке $ x_0$. Положим $ x=x_0+h$, тогда

$\displaystyle f(x_0+h)=f(x_0)+f'(x_0)h+\frac{f''(x_{{\theta}})}{2}h^2.$

Отсюда

$\displaystyle f'(x_0)=\dfrac{f(x_0+h)-f(x_0)}{h}+{\varepsilon}(x_0;h),$

где

$\displaystyle {\varepsilon}(x_0;h)=\frac{f''(x_{{\theta}})}{2}h$ --

погрешность формулы приближённого дифференцирования, получающаяся при замене $ f'(x_0)$ на разностную производную $ \dfrac{f(x_0+h)-f(x_0)}{h}$.

Следовательно,

$\displaystyle \vert{\varepsilon}(x_0;h)\vert\leqslant \frac{m_2}{2}h,$

где

$\displaystyle m_2=\max_{x\in[x_0;x_0+h]}\vert f''(x)\vert.$

[an error occurred while processing this directive]

Как правило, заранее известна более грубая оценка для $ f''$ на некотором отрезке $ [a;b]$, включающем в себя $ [x_0;x_0+h]$:

$\displaystyle M_2=\max_{x\in[a;b]}\vert f''(x)\vert\geqslant m_2,$

и $ M_2$ не зависит от $ x_0$ и $ h$. Тогда

$\displaystyle \vert{\varepsilon}(x_0;h)\vert\leqslant \frac{M_2}{2}h;$

из этой оценки и определяют погрешность вычислений при данном шаге $ h$.

Аналогично, можно получить оценку погрешности для разностной производной вида

$\displaystyle \dfrac{f(x_0+h)-f(x_0-h)}{2h}.$

Ошибку $ {\varepsilon}(x_0;h)$ при замене $ f'(x_0)$ на это отношение можно оценить исходя из разложения $ f(x)$ в точке $ x_0$ по формуле Тейлора с остаточным членом в форме Лагранжа порядка 3:

$\displaystyle f(x_0+h)=f(x_0)+f'(x_0)h+\frac{f''(x_0)}{2}h^2+ \frac{f'''(x_{{\theta}})}{6}h^3,$

где $ x_{{\theta}}\in(x_0;x_0+h)$. Подставляя сюда $ -h$ вместо $ h$, получаем:

$\displaystyle f(x_0-h)=f(x_0)-f'(x_0)h+\frac{f''(x_0)}{2}h^2- \frac{f'''(x_{{\theta}_1})}{6}h^3,$

где $ x_{{\theta}_1}\in(x_0-h;x_0)$. Вычтем из первой формулы вторую:

$\displaystyle f(x_0+h)-f(x_0-h)=f'(x_0)\cdot2h+ \frac{f'''(x_{{\theta}})}{6}h^3- \frac{f'''(x_{{\theta}_1})}{6}h^3.$

Отсюда

$\displaystyle f'(x_0)=\dfrac{f(x_0+h)-f(x_0-h)}{2h}-\dfrac{1}{12}\left( f'''(x_{{\theta}})- f'''(x_{{\theta}_1})\right)h^2.$

Если теперь предположить, что

$\displaystyle \max_{x\in[a;b]}\vert f'''(x)\vert=M_3,$

то оценка погрешности получится такая:

        Упражнение 6.4   Исследуйте приближённую формулу $\displaystyle \vert{\varepsilon}(x_0;h)\vert=\vert f'(x_0)-\dfrac{f(x_0+h)-f(x_... ...)\vert+\vert f'''(x_{{\theta}_1})\vert\right)h^2\leqslant \dfrac{1}{6}M_3h^2.$

$\displaystyle f'(x_0)\approx\dfrac{f(x_0-2h)-8f(x_0-h)+8f(x_0+h)-f(x_0+2h)}{12h}.$

Какая степень приращения $ h$ будет множителем в оценке ошибки $ {\varepsilon}(x_0;h)$? Оценки каких производных войдут в формулу для оценки ошибки?     

Как было сказано выше, минором матрицы порядка s называется определитель матрицы, образованной из элементов исходной матрицы, находящихся на пересечении каких - либо выбранных s строк и s столбцов.  

Определение. В матрице порядка m´n минор порядка r называется базисным, если он не равен нулю, а все миноры порядка r+1 и выше равны нулю, или не существуют вовсе, т.е. r совпадает с меньшим из чисел m или n. Столбцы и строки матрицы, на которых стоит базисный минор, также называются базисными.  В матрице может быть несколько различных базисных миноров, имеющих одинаковый порядок.  

Определение. Порядок базисного минора матрицы называется рангом матрицы и обозначается Rg А.  Очень важным свойством элементарных преобразований матриц является то, что они не изменяют ранг матрицы.

Определение. Матрицы, полученные в результате элементарного преобразования, называются эквивалентными.  Надо отметить, что равные матрицы и эвивалентные матрицы - понятия совершенно различные.  

[an error occurred while processing this directive]

Теорема. Наибольшее число линейно независимых столбцов в матрице равно числу линейно независимых строк.  Т.к. элементарные преобразования не изменяют ранг матрицы, то можно существенно упростить процесс нахождения ранга матрицы.     Пример. Определить ранг матрицы. ~ ~,  RgA = 2.  Пример: Определить ранг матрицы. ~ ~ ~,  Rg = 2. Пример. Определить ранг матрицы. ~, Þ Rg = 2.  Если с помощью элементарных преобразований не удается найти матрицу, эквивалентную исходной, но меньшего размера, то нахождение ранга матрицы следует начинать с вычисления миноров наивысшего возможного порядка. В вышеприведенном примере – это миноры порядка 3. Если хотя бы один из них не равен нулю, то ранг матрицы равен порядку этого минора.

Теорема о базисном миноре.  Теорема. В произвольной матрице А каждый столбец (строка) является линейной комбинацией столбцов (строк), в которых расположен базисный минор. Таким образом, ранг произвольной матрицы А равен максимальному числу линейно независимых строк (столбцов) в матрице.    Если А- квадратная матрица и detA = 0, то по крайней мере один из столбцов – линейная комбинация остальных столбцов. То же самое справедливо и для строк. Данное утверждение следует из свойства линейной зависимости при определителе равном нулю.

Абсолютная и условная сходимость несобственных интегралов по бесконечному промежутку

  До сих пор рассматривались интегралы от знакоположительных (знакопостоянных) функций. Теперь пусть подынтегральная функция таких ограничений не имеет, то есть может быть и знакочередующейся функцией.

 Если наряду с собственным интегралом по бесконечному промежутку сходится и интеграл  по этому же промежутку, то первый интеграл называется абсолютно сходящимся.

 Если интеграл сходится, а интеграл  расходится, то первый интеграл называется условно сходящимся.

  Пример 8. Исследовать на абсолютную сходимость интеграл: .

  Решение. В начале исследуется данный интеграл вообще на сходимость, для чего проведем интегрированние по частям: пусть  тогда , далее . Так как последний интеграл сходится, то по признаку сравнения сходится и интеграл , причем абсолютно. Исходный интеграл  при этом является сходящимся (кстати, сходимость этого можно определить быстрее с помощью признака сходимости Дирихле, который будет рассмотрен позже). Чтобы исследовать исходный интеграл на абсолютную сходимость, надо рассмотреть интеграл: . Так как  при , то имеем: . Интеграл  аналогично исходному интегралу  сходится, а интеграл  расходится; стало быть, и интеграл  является расходящимся. При этом исходный интеграл является условно сходящимся.

 Упражнение 8. Установить условную сходимость интеграла: .

  Пример 9. Исследовать на абсолютную сходимость интеграл: .

Математический анализ Типовые расчеты по математике