Разложение функций в ряд Фурье – это математический прием, который можно наблюдать и в природе, если использовать прибор, чувствующий синусоидальные функции.
Данный процесс происходит, когда человек слышит какой-либо звук. Ухо человека устроено таким образом, что может чувствовать отдельные синусоидальные колебания давления воздуха разной частоты, что, в свою очередь, позволяет человеку распознавать речь, слушать музыку.
Ухо человека воспринимает звук не целиком, а через составляющие его ряда Фурье. Струны музыкального инструмента производит звуки, представляющие собой синусоидальные колебания различных частот. Действительность разложения света в ряд Фурье представляет радуга. Зрение человека воспринимает свет через некоторые его составляющие разных частот электромагнитных колебаний.
Преобразованием Фурье является функция, которая описывает фазу и амплитуду синусоид, определенной частоты. Это преобразование используют для решения уравнений, описывающих динамические процессы, которые возникают под действием энергии. Ряды Фурье решают задачу выделения постоянных составляющих в сложных колебательных сигналах, что позволило правильно трактовать полученные данные экспериментов, наблюдений в медицине, химии и астрономии .
Открытие данного преобразования принадлежит французскому математику Жан Батисту Жозефу Фурье. В честь, которого впоследствии было и названо рядом Фурье. Первоначально ученый нашел применение своего метода при изучении и объяснении механизмов теплопроводности. Было предположено, что изначальное нерегулярное распределение тепла можно представить в виде простейших синусоид. Для каждой, из которых будет определен температурный минимум, максимум и фаза. Функция, описывающая верхние и нижние пики кривой, фазу каждой гармоники называется преобразованием Фурье от выражения распределения температуры. Автор преобразования предложил способ разложения сложной функции в виде суммы периодических функций косинуса, синуса .
Целью курсовой работы является изучение ряда Фурье и актуальности практического применения данного преобразования.
Для достижения поставленной цели были сформулированы следующие задачи:
1) дать понятие тригонометрического ряда Фурье;
2) определить условия разложимости функции в ряд Фурье;
3) рассмотреть разложение в ряд Фурье четных и нечетных функций;
4) рассмотреть разложение в ряд Фурье непериодической функции;
5) раскрыть практическое применение ряда Фурье.
Объект исследования: разложение функций в ряд Фурье.
Предмет исследования: ряды Фурье.
Методы исследования: анализ, синтез, сравнение, аксиоматический метод.
1.5. Ряд Фурье для четных и нечетных функций
Рассмотрим симметричный интеграл
где непрерывная или кусочно-непрерывная на. Сделаем замену в первом интеграле. Полагаем. Тогда
Следовательно, если четная функция, то (т.е. график четной функции симметричен относительно оси и
Если - нечетная функция, то (т.е. график нечетной функции симметричен относительно начала координат) и
Т.е. симметричный интеграл от четной функции равен удвоенному интегралу по половинному промежутку интегрирования, а симметричный интеграл от нечетной функции равен нулю.
Отметим следующие два свойства четных и нечетных функций:
1) произведение четной функции на нечетную есть функция нечетная;
2) произведение двух четных (нечетных) функций есть функция четная.
Пусть - четная функция, заданная на и разлагающаяся на этом отрезке в тригонометрический ряд Фурье. Используя полученные выше результаты, получим, что коэффициенты этого ряда будут иметь вид:
Если - нечетная функция, заданная на отрезке и разлагающаяся на этом отрезке в тригонометрический ряд Фурье, то коэффициенты этого ряда будут иметь вид:
Следовательно, тригонометрический ряд Фурье на отрезке будет иметь вид
для четной функции:
(16)
для нечетной функции:
Ряд (16) не содержит синусов кратных углов, то есть в ряд Фурье четной функции входят только четные функции и свободный член. Ряд (17) не содержит косинусов кратных углов, то есть в ряд Фурье нечетной функции входят только нечетные функции .
Определение.
Ряды
являются частями полного ряда Фурье и называются неполными
тригонометрическими рядами Фурье.
Если функция разлагается в неполный тригонометрический ряд (16) (или (17)), то говорят, что она разлагается в тригонометрический ряд Фурье по косинусам (или по синусам).
1.6. Разложение в ряд Фурье непериодической функции
1.6.1. Разложение в ряд Фурье функций на
Пусть функция задана на отрезке и удовлетворяет на этом отрезке условиям теоремы Дирихле. Выполним замену переменной. Пусть, где подберем так, чтобы получившаяся функция аргумента была определена на. Следовательно, считаем, что
Получившуюся в результате замены функцию можно разложить на в ряд Фурье:
где
Сделаем обратную замену ⇒ Получим
где
(19)
Ряд (18) – ряд Фурье по основной тригонометрической системе функций
Таким образом, получили, что если функция задана на отрезке и удовлетворяет на этом отрезке условиям теоремы Дирихле, то она может быть разложена в тригонометрический ряд Фурье (18) по тригонометрической системе функций (20) .
Тригонометрический ряд Фурье для четной функции, заданной на, будет иметь вид
где
для нечетной функции
где
Замечание! В некоторых задачах требуется разложить функцию в тригонометрический ряд Фурье по системе функций (20) не на отрезке, а на отрезке. В этом случае необходимо просто изменить пределы интегрирования в формулах (19) ((15), если, то есть в этом случае
(23)
или, если
(24)
Сумма тригонометрического ряда Фурье периодическая функция с периодом, являющаяся периодическим продолжением заданной функции. А для периодической функции справедливо равенство (4).
1.6.2. Разложение в ряд Фурье функций на
Пусть функция задана на и удовлетворяет на этом отрезке условиям теоремы Дирихле. Такую функцию также можно разложить в ряд Фурье. Для этого функцию нужно доопределить на промежуток и полученную функцию разложить в ряд Фурье на отрезке. При этом полученный ряд следует рассматривать только на отрезке, на котором функция задана. Для удобства вычислений доопределим функцию четным и нечетным образом.
1) Продолжим функцию на промежуток четным образом, то есть построим новую четную функцию, совпадающую на отрезке с функцией. Следовательно, график этой функции симметричен относительно оси и на отрезке совпадает с графиком. По формулам (21) найдем коэффициенты ряда Фурье для функции и запишем сам ряд Фурье. Сумма ряда Фурье для – периодическая функция, с периодом. Она будет совпадать с функцией на во всех точках непрерывности.
2) Доопределим функцию на промежуток нечетным образом, то есть построим новую нечетную функцию, совпадающую на с функцией. График такой функции симметричен относительно начала координат и на отрезке совпадает с графиком. По формулам (22) найдем коэффициенты ряда Фурье для функции и запишем сам ряд Фурье. Сумма ряда Фурье для – периодическая функция с периодом. Она будет совпадать с функцией на во всех точках непрерывности.
Замечания!
1) Аналогично можно разложить в ряд Фурье функцию, заданную на отрезке
2) Так как разложение функции на отрезке предполагает ее продолжение на отрезок произвольным образом, то и ряд Фурье для функции не будет единственным .
1.6.3. Разложение в ряд Фурье функций на
Пусть функция задана на произвольном отрезке длины и удовлетворяет на нем условиям теоремы Дирихле.
Тогда эта функция может быть разложена в ряд Фурье. Для этого функцию нужно периодически (с периодом) продолжить на всю числовую прямую и полученную функцию разложить в ряд Фурье, который следует рассматривать только на отрезке. В силу свойства (3) периодических функций имеем
Поэтому коэффициенты Фурье для полученного продолжения функции можно найти по формулам
(25)
2. Практическое применение рядов Фурье
2.1. Задачи на разложение функций в ряд Фурье и их решение
В тригонометрический ряд Фурье требуется разложить функцию, являющуюся периодическим продолжением заданной на отрезке функции. Для этого необходимо пользоваться алгоритмом разложения периодической функции в ряд Фурье.
Алгоритм разложения периодической функции в ряд Фурье:
1) Построить график заданной функции и ее периодического продолжения;
2) Установить период заданной функции;
3) Определить функция четная, нечетная или общего вида;
4) Проверить выполнимость условий теоремы Дирихле;
5) Составить формальную запись ряда Фурье, порожденного данной функцией;
6) Вычислить коэффициенты Фурье;
7) Записать ряд Фурье для заданной функции, используя коэффициенты ряда Фурье (п.4).
Пример 1. Функцию разложить в ряд Фурье на промежутке.
Решение:
1) Построим график заданной функции и его периодическое продолжение.
2) Период разложения функции.
3) Функция - нечетная.
4) Функция - непрерывная и монотонна на, т.е. функция удовлетворяет условиям Дирихле.
5) Вычислим коэффициенты ряда Фурье.
6) Запишем ряд Фурье, подставив коэффициенты Фурье в формулу
Ответ:
Пример 2. Разложим функцию с произвольным периодом в ряд Фурье .
Решение: функция определена на полуинтервале (-3;3]. Период разложения функции, полупериод. Разложим функцию в ряд Фурье
В начале координат функция разрывная, поэтому каждый коэффициент Фурье будем представлять в виде суммы двух интегралов.
Запишем ряд Фурье, подставив найденные коэффициенты ряда Фурье в формулу.
Пример 3. Разложить функцию на промежутке в ряд Фурье по косинусам. Построить график суммы ряда.
Решение: продолжим функцию на промежуток четным образом, то есть построим новую четную функцию, совпадающую на отрезке с функцией. Найдем коэффициенты ряда Фурье для функции и запишем ряд Фурье. Сумма ряда Фурье для - периодическая функция, с периодом. Она будет совпадать с функцией на во всех точках непрерывности.
Тригонометрический ряд Фурье для функции будет иметь вид
Найдем коэффициенты ряда Фурье
Таким образом, когда найдены коэффициенты, можно записать ряд Фурье
Построим график суммы ряда
Пример 4. Дана функция, определенная на отрезке . Выяснить, можно ли разложить функцию в ряд Фурье. Записать разложение функции в ряд Фурье .
Решение:
1) построим график функции на .
2) функция непрерывна и монотонна на , то есть по теореме Дирихле разлагается в тригонометрический ряд Фурье.
3) вычислим коэффициенты Фурье по формулам (1.19).
4) запишем ряд Фурье, используя найденные коэффициенты.
2.2. Примеры применения рядов Фурье в различных областях деятельности человека
Математика – одна из наук, которая имеет широкое применение на практике. Любой производственно-технологический процесс основан на математических закономерностях. Применение различных инструментов математического аппарата позволяет конструировать устройства и автоматизированные агрегаты, способные выполнять операции, сложные расчеты и вычисления при проектировании зданий, сооружений.
Ряды Фурье применяются математиками в геометрии при решении задач в сферической геометрии; в м атематической физике при решении задач о малых колебаниях упругих сред. Но кроме математики ряды Фурье нашли свое применение и в других областях наук.
Ежедневно люди пользуются различными устройствами. И зачастую эти устройства работают неисправно. Например, звук плохо различим из-за больших шумов или изображение, полученное по факсу, нечеткое. Причину неисправности человек может определить по звуку. Компьютер также может провести диагностику повреждения устройства. Лишние шумы можно убрать с помощью компьютерной обработки сигналов. Сигнал представляют в виде последовательности цифровых значений, которые затем вводят в компьютер. Выполнив определенные вычисления, получают коэффициенты ряда Фурье.
Изменение спектра сигнала позволяет очищать запись от шумов, компенсировать искажения сигнала различными устройствами звукозаписи, менять тембры инструментов, акцентировать внимание слушателей на отдельных партиях.
В цифровой обработке изображений применение рядов Фурье позволяет проводить следующие эффекты: размытие, подчеркивание границ, восстановление изображений, художественные эффекты (тиснение)
Разложение в ряд Фурье применяется в архитектуре при исследовании колебательных процессов. Например, при создании проекта различного вида конструкций рассчитывают прочность, жесткость и устойчивость элементов конструкций.
В медицине для проведения медицинского обследования с помощью кардиограмм, аппарата УЗИ пользуются математическим аппаратом, в основе которого лежит теория рядов Фурье.
Объемные вычислительные задачи оценки статистических характеристик сигналов, фильтрации шумов возникают при регистрации и обработке данных морского непрерывного дна. При постановке измерений, их регистрации перспективны голографические методы, использующие ряды Фурье. То есть ряды Фурье применяются и в такой науке как океанология.
Элементы математики встречаются на производстве практически на каждом шагу, поэтому специалистам важно знать и блестяще ориентироваться в области применения тех или иных инструментов анализа и расчета .
Заключение
Тема курсовой работы посвящена изучению ряда Фурье. Произвольную функцию можно разложить на более простые, то есть можно разложить в ряд Фурье. Объем курсовой работы не позволяет подробно раскрыть все аспекты разложения функции в ряд. Однако, из поставленных задач, представилось возможным раскрыть основную теорию о рядах Фурье.
В курсовой работе раскрыто понятие тригонометрического ряда Фурье. Определены условия разложимости функции в ряд Фурье. Рассмотрены разложения в ряд Фурье четных и нечетных функций; непериодических функций.
Во второй главе приведены лишь некоторые примеры разложения функций, заданных на различных промежутках, в ряд Фурье. Описаны те области наук, где используется данное преобразование.
Существует также комплексная форма представления ряда Фурье, которую не удалось рассмотреть, так как не позволяет объем курсовой работы. Комплексная форма ряда алгебраически проста. Поэтому часто используется в физике и прикладных расчетах.
Важность темы курсовой работы обусловлена тем, что находит широкое применение не только в математике, но в других науках: физике, механике, медицине, химии и многих других.
Список литературы
1. Бари, Н.К. Тригонометрические ряды. [текст]/ Н.К. Бари. - Москва, 1961 . - 936 с .
2. Бермант, А.Ф. Краткий курс математического анализа: учебник для вузов [текст] / А.Ф. Бермант, И.Г. Араманович. – 11-е изд., стер. – СПб.: Издательство «Лань», 2005. – 736 с.
3. Бугров, Я. С. Высшая математика: Учебник для вузов: В 3 т. [текст] / Я. С. Бугров, С. М. Никольский; Под ред. В. А. Садовничего. - 6-е изд., стереотип. - М.: Дрофа, 2004. -512 с.
4. Виноградова, И. А. Задачи и упражнения по математическому анализу: пособие для университетов, пед. вузов: В 2 ч. [текст] / И. А. Виноградова, С. Н. Олехник, В.А. Садовничий; под ред. В.А. Садовничего. – 3-е изд., испр. – М.: Дрофа, 2001. – 712 с.
5. Гусак, А.А. Высшая математика. В 2-х т. Т. 2. Учебник для студентов вузов. [текст] / А. А. Гусак. – 5-е изд. – Минск: ТетраСистемс, 2004.
6. Данко, П.Е. Высшая математика в упражнениях и задачах: учебное пособие для вузов: 2 ч. [текст] / П.Е. Данко, А.Г. Попов, Т.Я. Кожевникова. Москва: ОНИКС: Мир и образование, 2003. – 306 с.
7. Лукин, А. Введение в цифровую обработку сигналов(математические основы) [текст]/ А. Лукин. - М., 2007. - 54 с.
8. Пискунов, Н. С. Дифференциальное и интегральное исчисления для втузов, т.2: Учебное пособие для втузов. [текст] / Н. С. Пискунов. - 13-е изд.- М.: Наука, 1985. - 432 с.
9. Рудин, У. Основы математического анализа. [текст] / У. Рудин. - 2-е изд., Пер. с англ. .- М.: Мир, 1976 .- 206 с.
10. Фихтенгольц, Г. М. Основы математического анализа. Часть 2. [текст] / Г. М. Фихтенгольц. - 6-е изд., стер. - СПб.: Издательство «Лань», 2005. – 464 с.
Оренбург, 2015 г.
В главе 10 описывалось применение рядов Фурье к исследованию упругих колебаний струны. В данной главе мы рассмотрим некоторые вопросы упругого изгиба балок.
Использование рядов Фурье для решения задач статики упругих тел производится по следующей схеме.
Прежде всего из физических соображений выводится соотношение, которое связывает функцию, описывающую геометрическое состояние деформированного тела, с приложенными к телу нагрузками. Это соотношение, вообще говоря, содержит, помимо самой функции состояния, еще и ее производные, а также некоторые интегральные характеристики.
Затем, исходя из геометрических очертаний тела и кинематических условий, ограничивающих его перемещения, выбирается ортогональная система функций, по которой указанная функция состояния разлагается в ряд Фурье.
Подстановка этого ряда Фурье в выведенное соотношение приводит к тождественному равенству двух рядов Фурье, от которого, пользуясь теоремой 2 § 14 главы 9, можно перейти к равенству коэффициентов при одинаковых функциях. Из этих последних равенств можно вычислить значения коэффициентов Фурье и тем самым описать состояние деформированного тела.
Этот процесс подстановки ряда Фурье в характеризующее изгиб соотношение следует осуществлять достаточно осмотрительно, ибо в ходе его приходится несколько раз почленно дифференцировать ряды Фурье, коэффициенты которых вычисляются лишь впоследствии. Убедиться в правомерности этого дифференцирования, т. е. (см. § 10 главы 5) в равномерной сходимости ряда, составленного
из производных членов дифференцируемого ряда, априори довольно затруднительно. Поэтому при решении каждой конкретной задачи мы будем рассуждать примерно следующим образом.
Сначала мы будем предполагать, что написанный с неизвестными пока коэффициентами ряд Фурье можно (в смысле теоремы § 10 главы 5) почленно дифференцировать нужное число раз. Выписывая производные и решая получающиеся уравнения, мы будем находить интересующие нас коэффициенты Фурье. Это будет означать, что если ряд Фурье поддается почленному дифференцированию (и притом столько раз, сколько это требуется), то он является вполне определенным, найденным нами рядом. Если теперь из рассмотрения полученных коэффициентов будет видно, что этот построенный, вполне определенный ряд действительно почленно дифференцируем, то все операции, проделанные фактически именно над этим рядом, были законными, и найденные коэффициенты Фурье - искомые. Если же окажется, что получился недифференцируемый ряд, то это значит, что проделанные с ним ранее действия были математически некорректными, а полученный на их основе результат - необоснованным, хотя, возможно, и верным. Далее мы познакомимся с примерами исходов обоих типов.
Ряды Фурье - это представление произвольно взятой функции с конкретным периодом в виде ряда. В общем виде данное решение называют разложением элемента по ортогональному базису. Разложение функций в ряд Фурье является довольно мощным инструментарием при решении разнообразных задач благодаря свойствам данного преобразования при интегрировании, дифференцировании, а также сдвиге выражения по аргументу и свертке.
Человек, не знакомый с высшей математикой, а также с трудами французского ученого Фурье, скорее всего, не поймет, что это за «ряды» и для чего они нужны. А между тем данное преобразование довольно плотно вошло в нашу жизнь. Им пользуются не только математики, но и физики, химики, медики, астрономы, сейсмологи, океанографы и многие другие. Давайте и мы поближе познакомимся с трудами великого французского ученого, сделавшего открытие, опередившее свое время.
Человек и преобразование Фурье
Ряды Фурье являются одним из методов (наряду с анализом и другими) Данный процесс происходит каждый раз, когда человек слышит какой-либо звук. Наше ухо в автоматическом режиме производит преобразование элементарных частиц в упругой среде раскладываются в ряды (по спектру) последовательных значений уровня громкости для тонов разной высоты. Далее мозг превращает эти данные в привычные для нас звуки. Все это происходит помимо нашего желания или сознания, само по себе, а вот для того чтобы понять эти процессы, понадобится несколько лет изучать высшую математику.
Подробнее о преобразовании Фурье
Преобразование Фурье можно проводить аналитическими, числительными и другими методами. Ряды Фурье относятся к числительному способу разложения любых колебательных процессов - от океанских приливов и световых волн до циклов солнечной (и других астрономических объектов) активности. Используя эти математические приемы, можно разбирать функции, представляя любые колебательные процессы в качестве ряда синусоидальных составляющих, которые переходят от минимума к максимуму и обратно. Преобразование Фурье является функцией, описывающей фазу и амплитуду синусоид, соответствующих определенной частоте. Данный процесс можно использовать для решения весьма сложных уравнений, которые описывают динамические процессы, возникающие под действием тепловой, световой или электрической энергии. Также ряды Фурье позволяют выделять постоянные составляющие в сложных колебательных сигналах, благодаря чему стало возможным правильно интерпретировать полученные экспериментальные наблюдения в медицине, химии и астрономии.
Историческая справка
Отцом-основателем этой теории является французский математик Жан Батист Жозеф Фурье. Его именем впоследствии и было названо данное преобразование. Изначально ученый применил свой метод для изучения и объяснения механизмов теплопроводности - распространения тепла в твердых телах. Фурье предположил, что изначальное нерегулярное распределение можно разложить на простейшие синусоиды, каждая из которых будет иметь свой температурный минимум и максимум, а также свою фазу. При этом каждая такая компонента будет измеряться от минимума к максимуму и обратно. Математическая функция, которая описывает верхние и нижние пики кривой, а также фазу каждой из гармоник, назвали преобразованием Фурье от выражения распределения температуры. Автор теории свел общую функцию распределения, которая трудно поддается математическому описанию, к весьма удобному в обращении ряду косинуса и синуса, в сумме дающих исходное распределение.
Принцип преобразования и взгляды современников
Современники ученого - ведущие математики начала девятнадцатого века - не приняли данную теорию. Основным возражением послужило утверждение Фурье о том, что разрывную функцию, описывающую прямую линию или разрывающуюся кривую, можно представить в виде суммы синусоидальных выражений, которые являются непрерывными. В качестве примера можно рассмотреть «ступеньку» Хевисайда: ее значение равно нулю слева от разрыва и единице справа. Данная функция описывает зависимость электрического тока от временной переменной при замыкании цепи. Современники теории на тот момент никогда не сталкивались с подобной ситуацией, когда разрывное выражение описывалось бы комбинацией непрерывных, обычных функций, таких как экспонента, синусоида, линейная или квадратичная.
Что смущало французских математиков в теории Фурье?
Ведь если математик был в прав в своих утверждениях, то, суммируя бесконечный тригонометрический ряд Фурье, можно получить точное представление ступенчатого выражения даже в том случае, если оно имеет множество подобных ступеней. В начале девятнадцатого века подобное утверждение казалось абсурдным. Но несмотря на все сомнения, многие математики расширили сферу изучения данного феномена, выведя его за пределы исследований теплопроводности. Однако большинство ученых продолжали мучиться вопросом: "Может ли сумма синусоидального ряда сходиться к точному значению разрывной функции?"
Сходимость рядов Фурье: пример
Вопрос о сходимости поднимается всякий раз при необходимости суммирования бесконечных рядов чисел. Для понимания данного феномена рассмотрим классический пример. Сможете ли вы когда-либо достигнуть стены, если каждый последующий шаг будет вдвое меньше предыдущего? Предположим, что вы находитесь в двух метрах от цели, первый же шаг приближает к отметке на половине пути, следующий - к отметке в три четверти, а после пятого вы преодолеете почти 97 процентов пути. Однако сколько бы вы шагов ни сделали, намеченной цели вы не достигните в строгом математическом смысле. Используя числовые расчеты, можно доказать, что в конце концов можно приблизиться на сколь угодно малое заданное расстояние. Данное доказательство является эквивалентным демонстрации того, что суммарное значение одной второй, одной четвертой и т. д. будет стремиться к единице.
Вопрос сходимости: второе пришествие, или Прибор лорда Кельвина
Повторно данный вопрос поднялся в конце девятнадцатого века, когда ряды Фурье попробовали применить для предсказания интенсивности отливов и приливов. В это время лордом Кельвином был изобретен прибор, представляющий собой аналоговое вычислительное устройство, которое позволяло морякам военного и торгового флота отслеживать это природное явление. Данный механизм определял наборы фаз и амплитуд по таблице высоты приливов и соответствующих им временных моментов, тщательно замеренных в данной гавани в течение года. Каждый параметр представлял собой синусоидальную компоненту выражения высоты прилива и являлся одной из регулярных составляющих. Результаты измерений вводились в вычислительный прибор лорда Кельвина, синтезирующий кривую, которая предсказывала высоту воды как временную функцию на следующий год. Очень скоро подобные кривые были составлены для всех гаваней мира.
А если процесс будет нарушен разрывной функцией?
В то время представлялось очевидным, что прибор, предсказывающий приливную волну, с большим количеством элементов счета может вычислить большое количество фаз и амплитуд и так обеспечить более точные предсказания. Тем не менее оказалось, что данная закономерность не соблюдается в тех случаях, когда приливное выражение, которое следует синтезировать, содержало резкий скачок, то есть являлось разрывным. В том случае, если в устройство вводятся данные из таблицы временных моментов, то оно производит вычисления нескольких коэффициентов Фурье. Исходная функция восстанавливается благодаря синусоидальным компонентам (в соответствии с найденными коэффициентами). Расхождение между исходным и восстановленным выражением можно измерять в любой точке. При проведении повторных вычислений и сравнений видно, что значение наибольшей ошибки не уменьшается. Однако они локализируются в области, соответствующей точке разрыва, а в любой иной точке стремятся к нулю. В 1899 году этот результат был теоретически подтвержден Джошуа Уиллардом Гиббсом из Йельского университета.
Сходимость рядов Фурье и развитие математики в целом
Анализ Фурье неприменим к выражениям, содержащим бесконечное количество всплесков на определенном интервале. В общем и целом ряды Фурье, если изначальная функция представлена результатом реального физического измерения, всегда сходятся. Вопросы сходимости данного процесса для конкретных классов функций привели к появлению новых разделов в математике, например теории обобщенных функций. Она связана с такими именами, как Л. Шварц, Дж. Микусинский и Дж. Темпл. В рамках данной теории была создана четкая и точная теоретическая основа под такие выражения, как дельта-функция Дирака (она описывает область единой площади, сконцентрированной в бесконечно малой окрестности точки) и «ступень» Хевисайда. Благодаря этой работе ряды Фурье стали применимы для решения уравнений и задач, в которых фигурируют интуитивные понятия: точечный заряд, точечная масса, магнитные диполи, а также сосредоточенная нагрузка на балке.
Метод Фурье
Ряды Фурье, в соответствии с принципами интерференции, начинаются с разложения сложных форм на более простые. Например, изменение теплового потока объясняется его прохождением сквозь различные препятствия из теплоизолирующего материала неправильной формы или изменением поверхности земли - землетрясением, изменением орбиты небесного тела - влиянием планет. Как правило, подобные уравнения, описывающие простые классические системы, элементарно решаются для каждой отдельной волны. Фурье показал, что простые решения также можно суммировать для получения решения более сложных задач. Выражаясь языком математики, ряды Фурье - это методика представления выражения суммой гармоник - косинусоид и синусоид. Поэтому данный анализ известен также под именем «гармонический анализ».
Ряд Фурье - идеальная методика до «компьютерной эпохи»
До создания компьютерной техники методика Фурье являлась лучшим оружием в арсенале ученых при работе с волновой природой нашего мира. Ряд Фурье в комплексной форме позволяет решать не только простые задачи, которые поддаются прямому применению законов механики Ньютона, но и фундаментальные уравнения. Большинство открытий ньютоновской науки девятнадцатого века стали возможны только благодаря методике Фурье.
Ряды Фурье сегодня
С развитием компьютеров преобразования Фурье поднялись на качественно новый уровень. Данная методика прочно закрепилась практически во всех сферах науки и техники. В качестве примера можно привести цифровой аудио- и видеосигнал. Его реализация стала возможной только благодаря теории, разработанной французским математиком в начале девятнадцатого века. Так, ряд Фурье в комплексной форме позволил совершить прорыв в изучении космического пространства. Кроме того, это повлияло на изучение физики полупроводниковых материалов и плазмы, микроволновой акустики, океанографии, радиолокации, сейсмологии.
Тригонометрический ряд Фурье
В математике ряд Фурье является способом представления произвольных сложных функций суммой более простых. В общих случаях количество таких выражений может быть бесконечным. При этом чем больше их число учитывается при расчете, тем точнее получается конечный результат. Чаще всего в качестве простейших используют тригонометрические функции косинуса или синуса. В таком случае ряды Фурье называют тригонометрическими, а решение таких выражений - разложением гармоники. Этот метод играет важную роль в математике. Прежде всего, тригонометрический ряд дает средства для изображения, а также изучения функций, он является основным аппаратом теории. Кроме того, он позволяет решать ряд задач математической физики. Наконец, данная теория способствовала развитию вызвала к жизни целый ряд весьма важных разделов математической науки (теорию интегралов, теорию периодических функций). Кроме того, послужила отправным пунктом для развития следующих функций действительного переменного, а также положила начало гармоническому анализу.
Этот ряд может быть также записан в виде:
(2),
где , k-я комплексная амплитуда.
Связь между коэффициентами (1) и (3) выражается следующими формулами:
Отметим, что все эти три представления ряда Фурье совершенно равнозначны. Иногда при работе с рядами Фурье бывает удобнее использовать вместо синусов и косинусов экспоненты мнимого аргумента, то есть использовать преобразование Фурье в комплексной форме. Но нам удобно использовать формулу (1), где ряд Фурье представлен в виде суммы косинусоид с соответствующими амплитудами и фазами. В любом случае неправильно говорить, что результатом преобразования Фурье действительного сигнала будут комплексные амплитуды гармоник. Как правильно говорится в Вики «Преобразование Фурье (?) - операция, сопоставляющая одной функции вещественной переменной другую функцию, также вещественной переменной.»
Итого:
Математической основой спектрального анализа сигналов является преобразование Фурье.
Преобразование Фурье позволяет представить непрерывную функцию f(x) (сигнал), определенную на отрезке {0, T} в виде суммы бесконечного числа (бесконечного ряда) тригонометрических функций (синусоид и\или косинусоид) с определёнными амплитудами и фазами, также рассматриваемых на отрезке {0, T}. Такой ряд называется рядом Фурье.
Отметим еще некоторые моменты, понимание которых требуется для правильного применения преобразования Фурье к анализу сигналов. Если рассмотреть ряд Фурье (сумму синусоид) на всей оси Х, то можно увидеть, что вне отрезка {0, T} функция представленная рядом Фурье будет будет периодически повторять нашу функцию.
Например, на графике рис.7 исходная функция определена на отрезке {-T\2, +T\2}, а ряд Фурье представляет периодическую функцию, определенную на всей оси х.
Это происходит потому, что синусоиды сами являются периодическими функциями, соответственно и их сумма будет периодической функцией.
рис.7 Представление непериодической исходной функции рядом Фурье
Таким образом:
Наша исходная функция - непрерывная, непериодическая, определена на некотором отрезке длиной T.
Спектр этой функции - дискретный, то есть представлен в виде бесконечного ряда гармонических составляющих - ряда Фурье.
По факту, рядом Фурье определяется некоторая периодическая функция, совпадающая с нашей на отрезке {0, T}, но для нас эта периодичность не существенна.
Периоды гармонических составляющих кратны величине отрезка {0, T}, на котором определена исходная функция f(x). Другими словами, периоды гармоник кратны длительности измерения сигнала. Например, период первой гармоники ряда Фурье равен интервалу Т, на котором определена функция f(x). Период второй гармоники ряда Фурье равен интервалу Т/2. И так далее (см. рис. 8).
рис.8 Периоды (частоты) гармонических составляющих ряда Фурье (здесь Т=2?)
Соответственно, частоты гармонических составляющих кратны величине 1/Т. То есть частоты гармонических составляющих Fk равны Fk= к\Т, где к пробегает значения от 0 до?, например к=0 F0=0; к=1 F1=1\T; к=2 F2=2\T; к=3 F3=3\T;… Fk= к\Т (при нулевой частоте - постоянная составляющая).
Пусть наша исходная функция, представляет собой сигнал, записанный в течение Т=1 сек. Тогда период первой гармоники будет равен длительности нашего сигнала Т1=Т=1 сек и частота гармоники равна 1 Гц. Период второй гармоники будет равен длительности сигнала, деленной на 2 (Т2=Т/2=0,5 сек) и частота равна 2 Гц. Для третьей гармоники Т3=Т/3 сек и частота равна 3 Гц. И так далее.
Шаг между гармониками в этом случае равен 1 Гц.
Таким образом сигнал длительностью 1 сек можно разложить на гармонические составляющие (получить спектр) с разрешением по частоте 1 Гц.
Чтобы увеличить разрешение в 2 раза до 0,5 Гц - надо увеличить длительность измерения в 2 раза - до 2 сек. Сигнал длительностью 10 сек можно разложить на гармонические составляющие (получить спектр) с разрешением по частоте 0,1 Гц. Других способов увеличить разрешение по частоте нет.
Существует способ искусственного увеличения длительности сигнала путем добавления нулей к массиву отсчетов. Но реальную разрешающую способность по частоте он не увеличивает.
3. Дискретные сигналы и дискретное преобразование Фурье
С развитием цифровой техники изменились и способы хранения данных измерений (сигналов). Если раньше сигнал мог записываться на магнитофон и храниться на ленте в аналоговом виде, то сейчас сигналы оцифровываются и хранятся в файлах в памяти компьютера в виде набора чисел (отсчетов).Обычная схема измерения и оцифровки сигнала выглядит следующим образом.
рис.9 Схема измерительного канала
Сигнал с измерительного преобразователя поступает на АЦП в течение периода времени Т. Полученные за время Т отсчеты сигнала (выборка) передаются в компьютер и сохраняются в памяти.
рис.10 Оцифрованный сигнал - N отсчетов полученных за время Т
Какие требования выдвигаются к параметрам оцифровки сигнала? Устройство, преобразующее входной аналоговый сигнал в дискретный код (цифровой сигнал) называется аналого-цифровой преобразователь (АЦП, англ. Analog-to-digital converter, ADC) (Wiki).
Одним из основных параметров АЦП является максимальная частота дискретизации (или частота семплирования, англ. sample rate) - частота взятия отсчетов непрерывного во времени сигнала при его дискретизации. Измеряется в герцах. ((Wiki))
Согласно теореме Котельникова, если непрерывный сигнал имеет спектр, ограниченный частотой Fмакс, то он может быть полностью и однозначно восстановлен по его дискретным отсчетам, взятым через интервалы времени 
, т.е. с частотой Fd ? 2*Fмакс, где Fd - частота дискретизации; Fмакс - максимальная частота спектра сигнала. Другими слова частота оцифровки сигнала (частота дискретизации АЦП) должна как минимум в 2 раза превышать максимальную частоту сигнала, который мы хотим измерить.
А что будет, если мы будем брать отсчеты с меньшей частотой, чем требуется по теореме Котельникова?
В этом случае возникает эффект «алиасинга» (он же стробоскопический эффект, муаровый эффект), при котором сигнал высокой частоты после оцифровки превращается в сигнал низкой частоты, которого на самом деле не существует. На рис. 5 красная синусоида высокой частоты - это реальный сигнал. Синяя синусоида более низкой частоты - фиктивный сигнал, возникающий вследствие того, за время взятия отсчета успевает пройти больше, чем пол-периода высокочастотного сигнала.
Рис. 11. Появление ложного сигнала низкой частоты при недостаточно высокой частоте дискретизации
Чтобы избежать эффекта алиасинга перед АЦП ставят специальный антиалиасинговый фильтр - ФНЧ (фильтр нижних частот), который пропускает частоты ниже половины частоты дискретизации АЦП, а более высокие частоты зарезает.
Для того, чтобы вычислить спектр сигнала по его дискретным отсчетам используется дискретное преобразование Фурье (ДПФ). Отметим еще раз, что спектр дискретного сигнала «по определению» ограничен частотой Fмакс, меньшей половине частоты дискретизации Fd. Поэтому спектр дискретного сигнала может быть представлен суммой конечного числа гармоник, в отличие от бесконечной суммы для ряда Фурье непрерывного сигнала, спектр которого может быть неограничен. Согласно теореме Котельникова максимальная частота гармоники должна быть такой, чтобы на нее приходилось как минимум два отсчета, поэтому число гармоник равно половине числа отсчетов дискретного сигнала. То есть если в выборке имется N отсчетов, то число гармоник в спектре будет равно N/2.
Рассмотрим теперь дискретное преобразование Фурье (ДПФ).
Сравнивая с рядом Фурье
Видим, что они совпадают, за исключением того, что время в ДПФ имеет дискретный характер и число гармоник ограничено величиной N/2 - половиной числа отсчетов.
Формулы ДПФ записываются в безразмерных целых переменных k, s, где k – номера отсчетов сигнала, s – номера спектральных составляющих.
Величина s показывает количество полных колебаний гармоники на периоде Т (длительности измерения сигнала). Дискретное преобразование Фурье используется для нахождения амплитуд и фаз гармоник численным методом, т.е. «на компьютере»
Возвращаясь к результатам, полученным в начале. Как уже было сказано выше, при разложении в ряд Фурье непериодической функции (нашего сигнала), полученный ряд Фурье фактически соответствует периодической функции с периодом Т. (рис.12).
рис.12 Периодическая функция f(x) с периодом Т0, с периодом измерения Т>T0
Как видно на рис.12 функция f(x) периодическая с периодом Т0. Однако из-за того, что длительность измерительной выборки Т не совпадает с периодом функции Т0, функция, получаемая как ряд Фурье, имеет разрыв в точке Т. В результате спектр данной функции будет содержать большое количество высокочастотных гармоник. Если бы длительность измерительной выборки Т совпадала с периодом функции Т0, то в полученном после преобразования Фурье спектре присутствовала бы только первая гармоника (синусоида с периодом равным длительности выборки), поскольку функция f(x) представляет собой синусоиду.
Другими словами, программа ДПФ «не знает», что наш сигнал представляет собой «кусок синусоиды», а пытается представить в виде ряда периодическую функцию, которая имеет разрыв из-за нестыковки отдельных кусков синусоиды.
В результате в спектре появляются гармоники, которые должны в сумме изобразить форму функции, включая этот разрыв.
Таким образом, чтобы получить «правильный» спектр сигнала, являющегося суммой нескольких синусоид с разными периодами, необходимо чтобы на периоде измерения сигнала укладывалось целое число периодов каждой синусоиды. На практике это условие можно выполнить при достаточно большой длительности измерения сигнала.
Рис.13 Пример функции и спектра сигнала кинематической погрешности редуктора
При меньшей длительности картина будет выглядеть «хуже»:
Рис.14 Пример функции и спектра сигнала вибрации ротора
На практике бывает сложно понять, где «реальные составляющие», а где «артефакты», вызванные некратностью периодов составляющих и длительности выборки сигнала или «скачками и разрывами» формы сигнала. Конечно слова «реальные составляющие» и «артефакты» не зря взяты в кавычки. Наличие на графике спектра множества гармоник не означает, что наш сигнал в реальности из них «состоит». Это все равно что считать, будто число 7 «состоит» из чисел 3 и 4. Число 7 можно представить в виде суммы чисел 3 и 4 - это правильно.
Так и наш сигнал… а вернее даже не «наш сигнал», а периодическую функцию, составленную путем повторения нашего сигнала (выборки) можно представить в виде суммы гармоник (синусоид) с определенными амплитудами и фазами. Но во многих важных для практики случаях (см. рисунки выше) действительно можно связать полученные в спектре гармоники и с реальными процессами, имеющими циклический характер и вносящими значительный вклад в форму сигнала.
Некоторые итоги
1. Реальный измеренный сигнал, длительностью T сек, оцифрованный АЦП, то есть представленный набором дискретных отсчетов (N штук), имеет дискретный непериодический спектр, представленный набором гармоник (N/2 штук).2. Сигнал представлен набором действительных значений и его спектр представлен набором действительных значений. Частоты гармоник положительны. То, что математикам бывает удобнее представить спектр в комплексной форме с использованием отрицательных частот не значит, что «так правильно» и «так всегда надо делать».
3. Сигнал, измеренный на отрезке времени Т определен только на отрезке времени Т. Что было до того, как мы начали измерять сигнал, и что будет после того - науке это неизвестно. И в нашем случае - неинтересно. ДПФ ограниченного во времени сигнала дает его «настоящий» спектр, в том смысле, что при определенных условиях позволяет вычислить амплитуду и частоту его составляющих.
Использованные материалы и другие полезные материалы.
ПРИМЕНЕНИЕ РЯДОВ ФУРЬЕ ДЛЯ ПРОГНОЗИРОВАНИЯ И ОПТИМИЗАЦИИ ПОСТАВОК ПРЕДПРИЯТИЯ ОПТОВОЙ ТОРГОВЛИ В АСПЕКТЕ УПРАВЛЕНИЯ СОБСТВЕННЫМ И АРЕНДУЕМЫМ ТРАНСПОРТОМ
Горлач Борис Алексеевич 1 , Шигаева Наталья Валерьевна 2
1 Самарский государственный аэрокосмический университет имени академика С.П. Королева (НИУ), д.т.н, профессор
2 Самарский государственный аэрокосмический университет имени академика С.П. Королева (НИУ)
Аннотация
В работе рассмотрен механизм моделирования случайного процесса (для статистических данных о предприятии) с использованием аппарата гармонического анализа. Решена задача рационального распределения объемов поставок сырья между собственным и арендуемым транспортом с целью сокращения затрат на хранение продукции.
THE FOURIER SERIES APPLICATION FOR PREDICTION AND OPTIMIZATION OF DELIVERY COSTS
Gorlach Boris Alekseevich 1 , Shigaeva Nathalie Valerievna 2
1 Samara State Aerospace University, doctor of technical Sciences, Professor
2 Samara State Aerospace University
Abstract
Тhe mechanism of simulation of a random process is considered (for the enterprise data). Harmonic analysis is widely adopted in modeling of enterprise costs. The problem of rational distribution of the raw materials deliveries between own transport and rented transport is solved.
Библиографическая ссылка на статью:
Горлач Б.А., Шигаева Н.В. Применение рядов фурье для прогнозирования и оптимизации поставок предприятия оптовой торговли в аспекте управления собственным и арендуемым транспортом // Экономика и менеджмент инновационных технологий. 2014. № 7 [Электронный ресурс]..02.2019).
Введение. Издержки предприятия на создание системы хранения товаров создают необходимость рационального распределения поставок. Решение задачи управления поставками связано с изменением потребностей предприятия в сырье. Для разработки модели рационального распределения проделана обработка статистических данных предприятия о спросе на сырье.
Статья состоит из следующих частей: построение модели случайного процесса, оптимизация поставок на примере упрощенной модели и на примере реальных данных.
Часть первая. Построение математической модели случайного процесса.
В ретроспективном периоде статистические данные о хранении ресурса на складе выглядят следующим образом (Таблица 1). Предполагается, что задана совокупность статистических данных Y i =Y(t i) в виде временного ряда.
Таблица 1 – Статистические данные спроса на ресурс
Как правило, математические модели временных рядов экономических процессов представляются в виде совокупности 4 компонент: сезонной S, циклической C, случайной ξ и тренда U. Данные компоненты образуют аддитивную модель статистических данных .
Составляющая U – тренд – подбирается таким образом, чтобы она не противоречила основной тенденции изменения исследуемой функции и не затрудняла ее анализ . В данной работе подбор тренда осуществляется с помощью функций Excel, а также вручную методом «нормальных уравнений».
После выполнения процедуры подбора наиболее адекватного тренда, выполняется нормализация функции, позволяющая обеспечить моделирование колебательной составляющей. В данном исследовании колебательная составляющая подбирается с использованием модели, представляющей собой тригонометрический ряд Фурье:
.
Коэффициенты ряда Фурье определяются следующим образом:
После проведения поиска в 6 итераций с помощью средств Excel была выявлена следующая функция колебательной составляющей:
S(t) = -0.215sinπt/6 – 0.077cos πt/6 -0.085sin πt/3-0.013cos πt/3+0.001 sin πt/2+0.023cosπt/2-0.035 sin2πt/3+0.055cos 2πt/3+0.003 sin 5πt/6+0.054cos 5πt/6+0.056cos πt
Динамика поставок и хранения ресурса на складе, а также функциональная зависимость объема ресурса после осуществления нормализации представлены на рисунке 1.
Рисунок 1 – Колебательная составляющая для реальных данных
Вычислим коэффициент детерминации для полученной функции.
Коэффициент детерминации для полученной функции равен 0,75. Следовательно, тренд описывает статистические данные на 75 процентов, а вероятность несоответствия полученной функции реальным статистическим значениям равна 0,25.
Часть вторая. Оптимизация процесса поставок
При составлении пропорции в поставках сырья следует учитывать несколько факторов, влияющих на показатель экономической эффективности поставок:
Своевременность и частота поставок
Стоимость поставок
Допустимые сроки хранения сырья
Обеспеченность предприятия складскими помещениями
Другие факторы .
Рассмотрим процесс оптимизации поставок на упрощенном графике. Выделим в нормализованном тренде одну гармонику (одно слагаемое гармонического ряда) и ограничимся рассмотрением одного периода. Получится следующая упрощенная функция поставок:
Рассмотрим в данной работе три варианта поставок.
1. Поставки обеспечиваются только собственным транспортом на уровне y=1, которому соответствует значение s(t)=0.
Процесс накопления ресурсов в первом полугодии и расхода во втором полугодии определяется формулой интеграла функции на рассматриваемом участке.
Накопленные ресурсы полностью расходуются в следующем полугодии. Проблема состоит в том, что объем хранения на складе имеет слишком большой разброс во времени и нуждается в оптимизации.
2. Собственный транспорт обеспечивает поставки, соответствующие минимальной интенсивности расходования ресурсов. Этот вариант подходит фирме, если предприятие имеет меньше капитала и в силу других причин не может позволить себе транспорта больше, чем минимальный уровень потребности в ресурсах, выглядит это следующим образом. Предприятие недополучает ресурсы в размере, равном площади интеграла между s(t) и прямой, характеризующей минимальный уровень поставок.
Предположим, что предприятие решило арендовать транспорт на уровне максимальной потребности в ресурсах в первом полугодии, тогда накопления полностью расходуются во втором полугодии.
3. Собственный транспорт обеспечивает поставки на уровне -h. Недостаток ресурсов компенсируется арендой транспорта.
Вычисляем уровень поставок h из условия равенства площадей накопления и расхода :
При полученном значении h недостаток ресурсов без аренды выглядит следующим образом:
Обобщая полученные результаты, составлен общий график накопления/расхода, который показывает, насколько оптимальный план отличается минимальным количеством складских ресурсов (Рисунок 2).
Рисунок 2 – Минимизация складских ресурсов
Исходя из графика, привлечение арендованного транспорта при осуществлении оптимизации хранения на складе позволяет сократить удельный объем хранения на складе до 10 раз, так как амплитуда значений функции накопления уменьшилась с 10 единиц до 1.
Часть 3. Оптимизация поставок на примере реальных данных
Выполнение оптимизации поставок начинается с выделения периода колебательной составляющей (в нашем примере t i ϵ 11..23) и поиска точек пересечения функции s(t) с осью Ox.
Иллюстрация варианта динамики поступления и расхода ресурса на предприятии, в котором транспортная аренда не предусмотрена, представлена на рисунке 3.
Рисунок 3 – Накопление/расход для реальных данных без аренды
Функция колебательной составляющей выглядит следующим образом:
S(t) = -0.215 sin πt/6-0.077cos πt/6 -0.085 sin πt/3-0.013cos πt/3+0.001 sin πt/2+0.023cos πt/2-0.035 sin 2πt/3+0.055cos 2πt/3+0.003 sin 5πt/6+0.054cos 5πt/6+0.056cos πt
Функция накопления:
Q = ∫S = (1/π)(0.215 *6* cos (πt/6)-0.077*6*sin (πt/6) +0.085*3*cos πt/3 – 0.013*3*sin πt/3 – 0.0013*2*cos πt/2+(0.023*2*sin πt/2+0.0349*6/4 cos 2πt/3+(0.0552*6/4)sin 2πt/3 – (0.0032*6/5) cos 5πt/6 + (0.0538*6/5)sin 5πt/6 + (0.0559*sin π t)
Определим максимальные площади запаса и расхода для функции поставок, при условии равенства нулю интенсивности поставок s(t).
Таблица 2 – Определение площадей запаса и расхода ресурса
Таким образом, Q max =0,9078 – это максимально возможное количество хранимых на складе ресурсов. Ресурсы, накопленные в первом полугодии, полностью расходуются во втором, т.к. тригонометрические функции имеют свойство симметричности.
Оптимизация с привлечением арендованного транспорта – эффективный способ снижения издержек на хранение ресурса на складе. Уровень поставок предприятия собственным транспортом задается величиной Y(t)=1-h , либо S(t)=-h из условия равенства площадей накопления и расхода по полугодиям (рисунок 4).
Рисунок 4 – Определение уровня поставок арендованным транспортом
В этом случае останется потребность в ресурсе в объеме, определяемом площадью прямоугольника с высотой h и основанием, составляющим весь интервал рассмотрения, равного (из свойств симметричности) площади интеграла циклической составляющей над прямой уровня поставок собственным транспортом. Предприятие арендует транспорт на части рассмотренного интервала . Уровень поставок арендуемым транспортом определится из равенства площадей недостатка ресурса (2) и объема аренды (1), изображенных на рисунке 4.
Поиск уровней h осуществляется итерационно. В варианте привлечения арендованных транспортных средств максимальный уровень хранения запасов на складе равен:
Верхний уровень h* находим из условия равенства площадей невосполненного спроса (1) на ресурсы и объема поставок (2), обозначенных на рисунке 4. Уровень аренды определен значением h*=0,144.
После проведения оптимизации была найдена площадь расхода и запаса:
Итоговая площадь запасов сократилась с 0,9 до 0,5:
Q max2 =0,2016+ 0,3137=0,515
Таким образом, оптимизация процесса поставок с помощью арендуемого транспорта дала сокращение складских расходов на 44%, что указывает на успешное выполнение задачи оптимизации.
Результаты и выводы . Предложенный алгоритм рационального распределения поставок между собственным транспортом предприятия и арендуемым в ходе моделирования функции издержек рядом Фурье опирается на характерные особенности нормализованного графика тренда, учитывает ограничения складских площадей, сроки хранения сырья, обеспечивает сокращение складских расходов (уровня хранения ресурсов на складе) до 50% раз для рассмотренных данных функции поставок. Таким образом, привлечение арендованного транспорта является эффективным способом сокращения складских расходов и издержек на хранение при высокой стоимости аренды и содержания складских помещений.
Библиографический список
- Савельев Г.Л. Задача оптимизации ресурсов предприятия в условиях циклического изменения потребности. – Самара: СГАУ, 2010. – 30 стр.
- Чуйкова Ю.С. Оптимизация материального потока в задаче управления запасами предприятия /Сборник научных статей «Управление организационно-экономическими системами». – Самара: СГАУ, 2009. – с. 25-30.
- Rardin R.L. Optimization in Operations Research. Prentice Hall, 1998.
