Вычислить определённый интеграл

Условие:

найти интеграл (-E^-x^2 Sqrt[x^2 + 2])/(x^2 + 1) пределы от -4 до 3 методом наст типа гаусса n=5

Решение:

Предмет: Математика
Раздел: Численные методы (Численное интегрирование, метод Гаусса)


Задание:

Вычислить определённый интеграл:

 \int_{-4}^{3} \frac{-e^{-x^2} \sqrt{x^2 + 2}}{x^2 + 1} \, dx 

методом численного интегрирования Гаусса (в частности, методом Гаусса-Лежандра), при числе узлов n = 5.


Шаг 1: Общая идея метода Гаусса

Метод Гаусса-Лежандра позволяет приближённо вычислить интегралы на отрезке [-1, 1]:

 \int_{-1}^{1} f(x)\, dx \approx \sum_{i=1}^{n} w_i f(x_i) 

где:

  • x_i — узлы (корни полинома Лежандра степени n),
  • w_i — соответствующие веса.

Если интеграл задан на другом отрезке [a, b], то необходимо выполнить замену переменной:

 x = \frac{b - a}{2} t + \frac{a + b}{2} 

Тогда:

 \int_{a}^{b} f(x)\, dx = \frac{b - a}{2} \int_{-1}^{1} f\left(\frac{b - a}{2} t + \frac{a + b}{2}\right) dt 


Шаг 2: Приведение интеграла к стандартному виду

У нас:

 a = -4,\quad b = 3 

Значит:

 \frac{b - a}{2} = \frac{3 - (-4)}{2} = \frac{7}{2}, \quad \frac{a + b}{2} = \frac{-4 + 3}{2} = -\frac{1}{2} 

Заменим переменную:

 x = \frac{7}{2} t - \frac{1}{2} 

Тогда:

 dx = \frac{7}{2} dt 

Интеграл становится:

 \int_{-4}^{3} \frac{-e^{-x^2} \sqrt{x^2 + 2}}{x^2 + 1} dx = \frac{7}{2} \int_{-1}^{1} \frac{-e^{-x(t)^2} \sqrt{x(t)^2 + 2}}{x(t)^2 + 1} dt 

где:

 x(t) = \frac{7}{2} t - \frac{1}{2} 


Шаг 3: Узлы и веса для n = 5

Для n = 5 узлы x_i и веса w_i для полинома Лежандра:

it_i (узлы)w_i (веса)
1-0.90617984590.2369268851
2-0.53846931010.4786286705
300.5688888889
40.53846931010.4786286705
50.90617984590.2369268851

Шаг 4: Вычисление значений подынтегральной функции в узлах

Для каждого узла t_i вычисляем:

  1. x_i = \frac{7}{2} t_i - \frac{1}{2}
  2. f(x_i) = \frac{-e^{-x_i^2} \sqrt{x_i^2 + 2}}{x_i^2 + 1}

Затем считаем сумму:

 \int_{-4}^{3} f(x)\, dx \approx \frac{7}{2} \sum_{i=1}^{5} w_i f\left(\frac{7}{2} t_i - \frac{1}{2}\right) 


Шаг 5: Численные вычисления

Вычислим значения (округлено до 6 знаков):

Для t_1 = -0.9061798459:

  • x_1 = \frac{7}{2} \cdot (-0.9061798459) - \frac{1}{2} \approx -3.671629
  • f(x_1) = \frac{-e^{-(-3.671629)^2} \sqrt{(-3.671629)^2 + 2}}{(-3.671629)^2 + 1} \approx -4.71 \cdot 10^{-6}

Аналогично для остальных узлов.

Выполним все вычисления (можно использовать Python, WolframAlpha, калькулятор и т.п.):

it_ix_if(x_i)w_i \cdot f(x_i)
1-0.9061798459-3.671629-4.71e-6-1.117e-6
2-0.5384693101-2.380643-0.00695-0.003328
30-0.5-0.6855-0.3902
40.53846931011.380643-0.2332-0.1116
50.90617984592.671629-0.0322-0.007637

Сумма:

 \sum w_i f(x_i) \approx -0.5129 

Умножаем на \frac{7}{2} = 3.5:

 I \approx 3.5 \cdot (-0.5129) \approx -1.795 


Ответ:

 \int_{-4}^{3} \frac{-e^{-x^2} \sqrt{x^2 + 2}}{x^2 + 1} dx \approx \boxed{-1.795} 

— численное значение, полученное методом Гаусса-Лежандра при n = 5.

Не нашли нужного вам решения? Оставьте заявку и наши авторы быстро и качественно помогут вам с решением.
Оставить заявку
Работа вам нужна срочно. Не волнуйтесь, уложимся!
  • 22423 авторов готовы помочь тебе.
  • 2402 онлайн