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

Рассмотрим пример на нахождение координат вектора.

Задача. Даны векторы $ {\overrightarrow {OA}={\bf a}}$ , $ {\overrightarrow {OB}={\bf b}}$ . Вектор $ {\overrightarrow {OC}={\bf c}}$ -- медиана треугольника $ OAB$ . Найдите координаты вектора a в базисе b, c.

Решение. Сначала рассмотрим геометрическое решение (рис. 10.13).




Рис.10.13.Геометрическое разложение вектора


Проведем через конец вектора a прямую параллельно вектору b до пересечения с продолжением вектора c. Получим точку пересечения $ D$ . Легко видеть, что $ {\overrightarrow {OD}=2{\bf c}}$ , $ {\overrightarrow {AD}={\bf b}}$ . Проведем через точку $ A$ прямую параллельно вектору c до пересечения с продолжением вектора b. Получим точку $ F$ . Очевидно, что $ {\vert OF\vert=\vert AD\vert}$ , то есть $ {\overrightarrow {OF}=-{\bf b}}$ . Таким образом, $ {\bf a}=\overrightarrow {OD}+\overrightarrow {OF}=2{\bf c}+(-{\bf b})=(-1){\bf b}+2{\bf c}$ . Получим $ {{\bf a}=(-1;2)}$ .

Аналитическое решение. Получим какое-нибудь уравнение, связывающее векторы a, b, c. Для этого достроим треугольник $ OAB$ до параллелограмма (рис. 10.14).




Рис.10.14.


Тогда $ \overrightarrow {OD}=2{\bf c}$ , $ \overrightarrow {OD}={\bf a}+{\bf b}$ . Получим равенство $ {2{\bf c}={\bf a}+{\bf b}}$ . Откуда $ {{\bf a}=-{\bf b}+2{\bf c}}$ , то есть $ {{\bf a}=(-1;2)}$ .
Ответ:$ {\bf a}=(-1;2)$

Дальнейшее повышение порядка точности формул, подобных квадратурным формулам центральных прямоугольников, трапеций, парабол, можно получить, применяя на отрезках разбиения интерполяцию функции $ f(x)$ многочленами более высокой степени (ниже мы использовали линейные функции, то есть многочлены степени 1, и квадратные трёхчлены -- многочлены степени 2).

Например, если использовать кубическое интерполирование, то есть приближать функцию $ f(x)$ многочленами степени 3, то получится формула, называемая кквадратурной формулой $ \frac{3}{8}$ (три восьмых)":

 

$\displaystyle I\approx I_{\frac{3}{8}}=\frac{1}{8}\sum_{i=1}^n \bigl(f(x_i) +3f(x_{i-\frac{1}{3}}) +3f(x_{i-\frac{2}{3}}) +f(x_{i-1})\bigr)\cdot h_i,$

где $ h_i=x_i-x_{i-1}$ , $ x_{i-\frac{1}{3}}=x_i-\frac{1}{3}h_i$ и $ x_{i-\frac{2}{3}}=x_i-\frac{2}{3}h_i$ , то есть точки $ x_{i-\frac{2}{3}}$ и $ x_{i-\frac{1}{3}}$ делят отрезок разбиения $ [x_{i-1};x_i]$ на три равных части. Число $ \frac{3}{8}$ в названии формулы связано с тем, что если положить $ h=h_i=\frac{b-a}{n}$ постоянным и ввести обозначение $ {\delta}=\frac{h}{3}$ , то формула получает вид

$\displaystyle I\approx I_{\frac{3}{8}}=\frac{3}{8}{\delta}\sum_{i=1}^n \bigl(f(x_i) +3f(x_{i-\frac{1}{3}}) +3f(x_{i-\frac{2}{3}}) +f(x_{i-1})\bigr).$

Если же использовать для интерполяции многочлены шестой степени $ P_i(x)$ и заменять интеграл ог $ f(x)$ на каждом из отрезков $ [x_{i-1};x_i]$ на интеграл от $ P_i(x)$ , то получится квадратурная формула, называемая формулой Уэддля:

$\displaystyle I\approx I_W=\frac{1}{20}\sum_{i=1}^n \bigl(f_{i0}+5f_{i1}+f_{i2}+6f_{i3}+f_{i4}+5f_{i5}+f_{i6}\bigr)\cdot h_i,$

где $ h_i=x_i-x_{i-1}$ и $ f_{ij}=f(x_{i-1}+j\cdot\frac{h_i}{6}.$ Таким образом, при применении формулы Уэддля на каждом очередном отрезке нужно вычислить 6 новых значений функции $ f_{i1},\dots,f_{i6}$ , а значение $ f_{i0}=f_{(i-1)6}$ вычислять заново не нужно, оно было уже вычислено на предыдущем шаге.

Для ошибки формулы Уэддля $ {\varepsilon}_W=I-I_W$ с постоянным шагом $ h_i=h=\frac{b-a}{n}$ имеется такая оценка:

$\displaystyle \vert{\varepsilon}_W\vert\leqslant \frac{b-a}{39\,191\,040}h^6\Bigl(M_6+\frac{h^2}{40}M_8\Bigr),$

где $ \vert f^{(6)}(x)\vert\leqslant M_6$ и $ \vert f^{(8)}(x)\vert\leqslant M_8$ при всех $ x\in[a;b]$ . При этом предполагается, что восьмая производная $ f^{(8)}(x)$ непрерывна на отрезке $ [a;b]$ .

Таким образом, формула Уэддля является квадраткрной формулой шестого порядка точности. На практике формулы более высокого порядка точности, чем формула Уэддля, не используются. Формула "$ \frac{3}{8}$ " используется редко. Если не устраивает формула Симпсона, то сразу переходят к формуле Уэддля.

        Замечание 5.2   Мы не затрагиваем здесь вопросов, связанных с квадратурными формулами Гаусса (в которых, для повышения точности, точки $ x_i\in[a;b]$ , в которых вычисляются значения функции $ f(x)$ , выбираются специальным образом), а также другие вопросы приближённого нахождения интегралов, например, нахождение несобственных интегралов, осциллирующих интегралов и других интегралов специальных типов. Интересующихся этими вопросами мы отсылаем к книгам
Бахвалов Н. С., Жидков Н. П., Кобельков Г. М. Численные методы. -- М.: Наука, 1987,
Амосов А. А., Дубинский Ю. А., Копченова Н. В. Вычислительные методы для инженеров. -- М.: Высш. шк., 1994,
Рябенький В. С. Введение в вычислительную математику. -- М.: Наука, Физматлит, 1994.     

Простейшие интегралы, содержащие квадратный трехчлен

 Рассмотрим интегралы вида

I. .

II. .

III. .

IV. .

 Выделяя из квадратного трехчлена ах2 +bx + c полный квадрат запишем его в виде ах2 +bx + c = а(х +b)2 + q. Если в интегралах I, II, III сделать замену х + b = z, то получим интегралы

I¢. .

II¢. .

III¢. .

 Вычисление этих интегралов в зависимости от знака числа а сводится к вычислению интегралов вида

, , ,

,  ,

каждый из которых представляет собой комбинацию двух интегралов, один из которых табличный, а другой сводится к табличному, применяя равенство d(z2 ± a2) = 2z dz. Интегралы  и  не входят в таблицу (см. таблицу простейших интегралов), но они уже были вычислены ранее. Так как интегралы такого вида часто встречаются в приложениях, а вычисление их технически сложно, то предлагается соответствующие первообразные просто запомнить. Поэтому эти интегралы также называют табличными:

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