Метод Симпсона и численные методы расчета объемов и сигналов

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

Численное интегрирование по формуле Симпсона. Метод Симпсона основывается на аппроксимации подинтегральной функции квадратичными параболами на элементарных отрезках разбиения. Пусть геометрический объект нерегулярной формы ограничен двумя параллельными плоскостями и разделен на равные интервалы шириной d. Если известен ряд площадей поперечных сечений A₁, A₂, A₃, ..., A_n, находящихся на одинаковом расстоянии d друг от друга, то суммарный объем геометрического тела определяется выражением:

V = (d / 3) × [(A₁ + A₇) + 4 × (A₂ + A₄ + A₆) + 2 × (A₃ + A₅)]

Рис. 2.32. Геометрическая схема тела с переменным поперечным сечением, разделенного на равные интервалы d.

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

Практический пример: Оценка объема ствола дерева. Рассмотрим прикладную задачу из лесной таксации и биометрии. Ствол дерева длиной 12 метров имеет переменное поперечное сечение. Измерения площадей сечений с шагом d = 2 м показали следующие значения: 0.52 м², 0.55 м², 0.59 м², 0.63 м², 0.72 м², 0.84 м² и 0.97 м².

Подставив имеющиеся эмпирические данные в формулу Симпсона, получаем:

V = (2 / 3) × [(0.52 + 0.97) + 4 × (0.55 + 0.63 + 0.84) + 2 × (0.59 + 0.72)]

V = (2 / 3) × [1.49 + 4 × (2.02) + 2 × (1.31)] = (2 / 3) × [1.49 + 8.08 + 2.62] = (2 / 3) × 12.19 = 8.13 м³.

Таким образом, объем ствола дерева составляет приблизительно 8.13 кубических метров.

Правило призм для определения объемов правильных и неправильных тел. Правило призм представляет собой частный случай метода Симпсона, адаптированный для тел длиной x, разделенных на два равных интервала тремя равноудаленными поверхностями или поперечными сечениями (A₁, A₂, A₃).

Рис. 2.33. Схема применения правила призм к телу длиной x со срединным сечением A₂.

Согласно данной математической модели, объем определяется формулой:

V = (x / 6) × [A₁ + 4 × A₂ + A₃]

Правило призм позволяет определять абсолютно точные значения объемов для правильных стереометрических фигур, таких как пирамида, конус, сфера, усеченный конус и призмоид.

Практический пример: Расчет емкости конического сосуда. Емкость имеет форму усеченного конуса высотой x = 24 см, диаметр дна составляет 18 см (радиус r₃ = 9 см), а диаметр верхней части — 30 см (радиус r₁ = 15 см). В силу линейного изменения радиуса, на расстоянии 12 см от дна средний радиус равен r₂ = (9 + 15) / 2 = 12 см.

Рис. 2.34. Геометрические параметры ведра в форме усеченного конуса.

Площади поперечных сечений рассчитываются по формуле площади круга A = π × r²:

· A₁ = π × (15)² см²

· A₂ = π × (12)² см²

· A₃ = π × (9)² см²

Подставляя значения в формулу правила призм:

V = (24 / 6) × [π × (15)² + 4 × π × (12)² + π × (9)²]

V = 4 × [706.86 + 1809.56 + 254.47] = 4 × 2770.89 = 11080 см³.

Переводя полученное значение в литры (1 л = 1000 см³), получаем V = 11080 / 1000 = 11 л (с точностью до целого литра).

Расчет средней величины сигнала и прямоугольное интегрирование. В электротехнике, радиотехнике и теории сигналов ключевую роль играет определение среднего значения функциональной зависимости y(t) на интервале длины b. Аналитически среднее значение определяется как отношение площади под кривой к длине основания:

y = (площадь под кривой) / (длина основания b)

Рис. 2.35. Дискретизация функции методом средних ординат.

При использовании правила средних ординат для определения площади под произвольной кривой формула принимает следующий вид:

y = (сумма средних ординат) / (количество средних ординат) = (y₁ + y₂ + y₃ + y₄ + y₅ + y₆ + y₇) / 7

Анализ средних значений синусоидальных волн. Далее подробно будет об этом разъяснение:

Рис. 2.36. Формы синусоидальных сигналов: а) полный период, б) двухполупериодное выпрямление, в) однополупериодное выпрямление.

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

- Полупериод синусоиды: среднее значение составляет 0.637 × Vm (или (2 / π) × Vm, где V — амплитудное значение).
- Двухполупериодный выпрямленный сигнал: равен 0.637 × Vm (Рис. 2.36б).
- Однополупериодный выпрямленный сигнал: равен 0.318 × mV (или (1 / π) × Vm) (Рис. 2.36в).

Вычисление параметров сложных периодических сигналов. Рассмотрим практическое вычисление средних значений для напряжений и токов различной формы волны, представленных на Рис. 2.37.

Рис. 2.37. Осциллограммы сигналов различной формы: а) треугольный сигнал, б) ступенчатый импульсный ток, в) узкий треугольный импульс.

Треугольный сигнал (Рис. 2.37а). Площадь треугольной волны за половину периода (2 мс) рассчитывается как:

S = (1 / 2) × (основание) × (высота) = (1 / 2) × (2 × 10⁻³ с) × (20 В) = 20 × 10⁻³ В·с.

Средняя величина сигнала равна:

U_{ср} = (20 × 10⁻³ В·с) / (2 × 10⁻³ с) = 10 В.

Ступенчатый ток (Рис. 2.37б). Площадь под импульсным сигналом за полупериод длительностью 3 секунды состоит из суммы прямоугольных областей:

S = (1 с × 1 А) + (3 А × 2 с) = 1 + 6 = 7 А·с.

Средний ток составляет:

I_{ср} = (7 А·с) / (3 с) = 2.33 А.

Импульсное напряжение (Рис. 2.37в). Половина периода сигнала составляет 4 мс (от 2 до 6 мс на графике). Треугольный импульс высотой 10 В существует на интервале шириной 3 - 1 = 2 мс.

S = (1 / 2) × [(3 - 1) × 10⁻³ с] × (10 В) = 10 × 10⁻³ В·с.

Средняя величина сигнала за полупериод равен:

U_{ср} = (10 × 10⁻³ В·с) / (4 × 10⁻³ с) = 2.5 В.

Расчет индикаторных диаграмм термодинамических систем. Численные методы интегрирования незаменимы при анализе работы поршневых машин, двигателей внутреннего сгорания и паровых машин по их индикаторным диаграммам.

Рис. 2.38. Индикаторная диаграмма парового двигателя с 7 ординатами.

Пусть основание индикаторной диаграммы длиной 12.0 см разделено на 6 равных интервалов (ширина каждого d = 12.0 / 6 = 2.0 см). Ординаты высоты равны соответственно: 3.6, 4.0, 3.5, 2.9, 2.2, 1.7 и 1.6 см.

Применяя формулу Симпсона, найдем площадь индикаторной диаграммы:

S = (1 / 3) × (2.0) × [(3.6 + 1.6) + 4 × (4.0 + 2.9 + 1.7) + 2 × (3.5 + 2.2)]

S = (2 / 3) × [(5.2) + 4 × (8.6) + 2 × (5.7)] = (2 / 3) × [5.2 + 34.4 + 11.4] = (2 / 3) × 51.0 = 34 см².

Для определения среднего давления в цилиндре сначала находится средняя высота ординат:

h_{ср} = (34 см²) / (12 см) = 2.83 см.

Если масштабом диаграммы установлено, что 1 см соответствует 100 кПа, то среднее индикаторное давление в цилиндре составит:

P_{ср} = 2.83 см × 100 кПа/см = 283 кПа.

Применение численных методик геометрического и функционального анализа гарантирует минимальную погрешность и является стандартом инженерных расчетов в гидродинамике, теплотехнике и электронике.

Сведения об авторах и источниках:

Авторы: Джон Берд

Источник: Инженерная математика

Данные публикации будут полезны инженерно-техническим работникам, студентам и абитуриентам технических вузов и колледжей, а также всем, кому требуется практический математический аппарат для решения инженерных задач.