85459 (612492), страница 2
Текст из файла (страница 2)
Нетрудно убедится, что эти формулы в точности соответствуют формулам (11).
3.1.1.2. Примеры.
1) Дана некоторая диаграмма касательных усилий (на пальце кривошипа) для некоторой паровой машины. В связи с вопросом о крутильных колебаниях вала представляет интерес выделить гармонические составляющие касательного усилия Т как функции от угла
поворота кривошипа. Сняв с графика двенадцать равноотстоящих ординат, произведем гармонический анализ по указанной схеме:
| T | -7200 | -300 | 7000 | 4300 | 0 | -5200 | -7400 |
| 250 | 4500 | 7600 | 3850 | -2250 | |||
| U | -7200 | -50 | 11500 | 11900 | 3850 | -7450 | -7400 |
| V | -550 | 2500 | -3300 | -3850 | -2950 |
| u | -7200 | -50 | 11500 | 11900 | |||
| -7400 | -7450 | 3850 | |||||
| s | -14600 | -7500 | 15350 | 11900 | |||
| d | 200 | 7400 | 7650 | ||||
| V | -550 | 2500 | -3300 | ||||
| -2950 | -3850 | ||||||
| | -3500 | -1350 | -3300 | ||||
| | 2400 | 6350 | |||||
Теперь по формулам (12):
Таким образом,
Объединим члены, содержащие косинус и синус одного и того же угла:
Мы видим, что наиболее сильное влияние здесь оказывает вторая гармоника.
2) Для того чтобы дать себе отчет в том, с какой примерно точностью получаются коэффициенты Фурье функции по двенадцати ординатам ее графика, мы приложим изложенный метод к некоторым аналитически заданным функциям и сравним приближенные результаты с точными.
Сначала рассмотрим функцию
, которую в промежутке
задается формулой
,
А для остальных значений x определяется по закону периодичности
.
Вычислим табличку:
| x | 0 | | | | | | | | | | | | 2 |
| y | 0 | 0.400 | 0.582 | 0.589 | 0.465 | 0.255 | 0 | -0.255 | -0.465 | -0.589 | -0.582 | -0.400 | 0 |
При этом можно использовать легко проверяемое тождество:
По схеме Рунге по этим значениям yнайдем:
b1=0.608; b2=0.076; b3=0.022;
все числа
, а с ними и все коэффициенты
оказываются нулями.
В то же время формулы (10) непосредственно дают (с помощью трехкратного интегрирования по частям):
,
Так что
;
;
.
Совпадение превосходное!
3) Однако далеко не всегда получается столь точный результат. В виде второго примера мы возьмем функцию с периодом
, которая в промежутке
определяется так:
.
Пользуясь тождеством:
,
составим таблицу:
| x | 0 | | | | | | | | | | | | 2 |
| y | 1 | 0,694 | 0,444 | 0,250 | 0,111 | 0,028 | 0 | 0,028 | 0,111 | 0,250 | 0,444 | 0,694 | 1 |
Тогда по схеме Рунге
числа же
и коэффициенты
- на этот раз нули. Точные значения коэффициентов будут:
в частности,
;
;
.
Таким образом, если для первых двух коэффициентов относительная погрешность не превосходит 1,5-2 %, то для последующих она достигнет10% и даже 20%! Ясно, что для повышения этой точности нужно брать больше ординат.
3.1.1.3. Схема для двадцати четырех ординат.
Положим теперь, что даны или сняты с графика двадцать четыре ординаты:
,
отвечающие значениям аргумента:
,
или
.
На этот раз все множители, на которые при приближенном вычислении коэффициентов Фурье приходится умножать ординаты, сведутся к таким:
Не вдаваясь в подробности (ввиду полной аналогии с предыдущим), приведем сразу схему вычислений, также предложенную Рунге. Ввиду вышеизложенного опыта следующая схема идет без пояснений. Вот она:
| ординаты | | | | | | | | | | | | | |
| | | | | | | | | | | | |||
| Суммы | | | | | | | | | | | | | |
| разности | | | | | | | | | | | |
| Суммы | | | | | | | |
| | | | | | | ||
| Суммы | | | | | | | |
| разности | | | | | | |
| Суммы | | | | | | |
| | | | | | ||
| Суммы | | | | | | |
| разности | | | | | |
| Суммы | | | | | |
| | | | |||
| Суммы | | | | | |
| разности | | | | ||
| разности | | | | ||
| | | ||||
| суммы | | | | ||
| разности | | | |||
Отметим, что с величинами q и r нет надобности проделывать сложения и вычитания.
Теперь через полученные указанным путем величины k, l, m, n, q и r коэффициентов Фурье выразятся следующим образом:














