Лабораторная работа №1 3 Общие сведения 3 Цель работы 3 План выполнения 3



страница6/7
Дата22.02.2016
Размер1.08 Mb.
ТипЛабораторная работа
1   2   3   4   5   6   7

Литература


  1. Алексинская Т. В. Учебное пособие по решению задач по курсу экономико-математические методы и модели. — Таганрог: Изд-во ТРТУ, 2002. — 153 с.

  2. Таха Х. Введение в исследование операций. — М.: Вильямс, 2005. — 912 с.

  3. Эддоус М., Стэнсфилд Р. Методы принятия решений. — М.: Юнити, 1997. — 587 с.


Лабораторная работа № 6


Методы анализа временных рядов

Общие сведения




Цель работы


Приобрести навыки прогнозирования значений временного ряда, в частности, выделения устойчивых тенденции и учета сезонной составляющей, а также навыки использования средств Пакет Анализа и Поиск решения, входящих в MS Excel.

План выполнения


  1. Изучить теоретическую часть;

  2. Получить задание преподавателя;

  3. Выполнить задания 1 и 2:

    1. Построить график значений временного ряда;

    2. Построить графики прогнозируемых значений (для каждого из полученных прогнозов);

    3. Вычислить среднее абсолютных отклонений или среднее относительных ошибок;

    4. Ответить на вопросы задачи;

  4. Составить отчёт по лабораторной работе. Отчёт должен иметь следующую структуру:

    1. Титульный лист, который должен содержать следующую информацию:

      1. Название университета и кафедры, ответственной за дисциплину;

      2. Заголовок — номер и название лабораторной работы;

      3. Подзаголовок — номер варианта и номера задач;

      4. ФИО и должности студента и преподавателя;

      5. «г. Санкт-Петербург, 2012 год»;

    2. Отчёт о решении заданиях 1 и 2, содержащий следующее информационное наполнение:

      1. Формулировка индивидуального задания;

      2. Графики значений временного ряда и прогноза;

      3. При необходимости, снимки экрана монитора, содержащие основные моменты решения задачи;

      4. Ответы на вопросы задания;

      5. Результаты решения и выводы.



Теоретическая часть

Общие сведения



Временным рядом называется последовательность значений некоторого показателя во времени (например, объемов продаж, см. Рисунок ).


Рисунок . Временный ряд объёма продаж
Анализ временных рядов является способом выявления тенденций прошлого и продления их в будущее. Методы анализа временных рядов осуществляют прогноз путем экстраполяции значений отдельной переменной на основе статистических данных за прошлый временной период. Основное допущение, которое при этом делается, заключается в том, что происшедшее в прошлом дает хорошее приближение в оценке будущего.

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





Рисунок . Зависимость объёма продаж от стадии жизненного цикла продукта

Кривые тренда сглаживают временной ряд значений показателя, выделяя общую тенденцию. Именно выбор кривой тенденции во многом определяет результаты прогнозирования.

В большинстве случаев временной ряд, кроме тенденции и случайных отклонений от него, характеризуется еще сезонной составляющей. Сезонная составляющая — это периодические изменения показателя. Обычная продолжительность сезонной составляющей измеряется днями, неделями или месяцами.

Методы без сезонной составляющей

Вначале рассмотрим несколько простейших методов прогнозирования, не учитывающих наличия сезонности во временном ряде. Предположим, что в журнале РБК приведена сводка за последние 12 дней (включая сегодняшний) цен на апельсины, сложившихся на момент закрытия биржи. Используя эти данные, нужно предсказать завтрашнюю цену на какао (также на момент закрытия биржи). Рассмотрим несколько способов сделать это.



  1. Если последнее (сегодняшнее) значение наиболее значимо по сравнению с остальными, то оно является наилучшим прогнозом на завтра.

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

  3. Если все значения существенны, но сегодняшнее 12-е значение наиболее значимо, а предыдущие 11-е, 10-е, 9-е и т.д. имеют все меньшую и меньшую значимость, следует найти взвешенное среднее всех 12 значений. Причем весовые коэффициенты для последних значений должны быть больше, чем для предыдущих, и сумма всех весовых коэффициентов должна равняться 1.

Первый способ называется «наивным» прогнозом и достаточно очевиден. Рассмотрим подробнее остальные способы.

Метод скользящего среднего


Одним из предположений, лежащих в основе данного метода, является то, что более точный прогноз на будущее можно получить, если использовались недавние наблюдения, причем, чем «новее» данные, тем их вес для прогноза должен быть больше. Удивительно, но такой «наивный» подход оказывается чрезвычайно полезным для практики. Например, многие авиакомпании используют частный тип скользящего среднего для создания прогнозов спроса на авиаперелеты, которые, в свою очередь, используются в сложных механизмах управления и оптимизации доходов. Более того, практически все программные пакеты управления запасами содержат модули, выполняющие прогнозы на основе того или иного типа скользящего среднего.

Рассмотрим следующий пример. Маркетологу нужно спрогнозировать спрос на производимые его компанией станки. Данные по объемам продаж за последний год работы компании находятся в файле «ЛР6.Пример 1.Станки.xls».



Простое скользящее среднее. В этом методе среднее фиксированного числа N последних наблюдений используется для оценки следующего значения временно ряда. Например, используя данные о продажах станков за первые три месяца года, менеджер получает для апреля значение, используя формулу, приведённую ниже:



(1).

В случае произвольного числа N узлов расчетная формула обобщается следующим образом:





(2).

Менеджер вычислил объем продаж на основе простого скользящего среднего за 3 и 4 месяца. Однако требуется определить, какое количество узлов даёт более точный прогноз. Для оценки точности прогнозов используются среднее абсолютных отклонений (САО) и среднее относительных ошибок, в процентах (СООП), вычисляемые по формулам (3) и (4).





(3),



(4).

где xi i-ое реальное значение переменной в i-й момент времени, а xi i-ое спрогнозированное значение переменной в i-й момент времени, N — количество прогнозов.

Согласно результатам, полученным на листе «Простое ск. среднее» рабочей книги «ЛР6.Пример 1.Станки.xls» (см. Рисунок ), скользящее среднее за три месяца имеет значение САО равное 12,67 (ячейка D16), тогда как для скользящего среднего за 4 месяца значение САО равно 15,59 (ячейка F16). Тогда можно выдвинуть гипотезу, что использование большего количества статистических данных скорее ухудшает, чем улучшает точность прогноза методом скользящего среднего.





Рисунок . Пример 1 – результаты прогнозирования методом простого скользящего среднего
На графике (см. Рисунок ), построенном по результатам наблюдений и прогнозов с интервалом 3 месяца, можно заметить ряд особенностей, общих для всех применений метода скользящего среднего.
Рисунок . Пример 1 – график кривой прогноза методом простого скользящего среднего и график реального объёма продаж

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

Основной недостаток метода простого скользящего среднего возникает в результате того, что при вычислении прогнозируемого значения самое последнее наблюдение имеет такой же вес (т. е. значимость), как и предыдущие. Это происходит потому, что вес всех N последних наблюдений, участвующих в вычислении скользящего среднего, равен 1/N. Присвоение равного веса противоречит интуитивному представлению о том, что во многих случаях последние данные могут больше сказать о том, что произойдет в ближайшем будущем, чем предыдущие.

Взвешенное скользящее среднее. Вклад различных моментов времени можно учесть, вводя вес для каждого значения показателя в скользящем интервале. В результате получается метод взвешенного скользящего среднего, который математически можно записать так:



(4).

где — вес, с которым используется показатель при расчете.

Вес — это всегда положительное число. В случае, когда все веса одинаковы, вырождается метод простого скользящего среднего.

Теперь маркетолог может использовать метод взвешенного скользящего среднего за 3 месяца. Но прежде требуется понять, как выбрать веса. Используя средство Поиск решения, можно определить оптимальный вес узлов. Чтобы определить вес узлов с помощью средства Поиск решения, при котором значение среднего абсолютных отклонений было бы минимально, выполните следующие действия:


  1. Выберите команду Сервис -> Поиск решения.

  2. В диалоговом окне Поиска решения установите ячейку G16 целевой (см. лист «Веса»), минимизируя её.

  3. Изменяемыми ячейками укажите диапазон В1:В3.

  4. Установите ограничения В4 = 1,0; В1:ВЗ ≥ 0; В1:В3 ≤ 1; B1 ≤ В2 и В2 ≤ В3.

  5. Запустите поиск решения (результат отображает).




Рисунок . Пример 1 – результат поиска весов значений показателей при использовании метода взвешенного скользящего среднего
Полученные результаты показывают, что оптимальное распределение весов таково, что весь вес сосредоточен на самом последнем наблюдении, при этом значение среднего абсолютных отклонений равно 7,56 (см. также Рисунок ). Этот результат подтверждает предположение о том, что более поздние наблюдения должны иметь больший вес.
Рисунок . Пример 1 – график кривой прогноза методом взвешенного скользящего среднего и график реального объёма продаж

Метод экспоненциального сглаживания


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

Очевидно, что в методе взвешенного скользящего среднего существует множество способов задавать значения весов так, чтобы их сумма была равной 1,и экспоненциальное сглаживание – один из таких способов. В этой схеме метода взвешенного среднего для любого t > 1 прогнозируемое значение в момент времени (t+1) представляет собой взвешенную сумму фактического объема продаж , за период времени t и прогнозируемого объема продаж , за период времени t (см. формулу (5)).





(5).

Экспоненциальное сглаживание имеет вычислительные преимущества перед скользящим средним. Здесь, чтобы вычислить , необходимо знать только значения и , (а также значение α), вместо значений показателя ряда во всех узлах, по которым происходит сглаживание. Сохраняя значение α и последний прогноз, мы также неявно сохраняем и все предыдущие прогнозы.

Рассмотрим некоторые свойства модели экспоненциального сглаживания. Для начала заметим, что если t > 2, то в формуле (5) t можно заменить на t–1, т.е.:



(6).

Подставив формулу (6) в первоначальную формулу (5), получим:





(7).

Выполняя последовательно аналогичные подстановки, получим следующее выражение для :





(8).

Поскольку из неравенства 0 < α < 1 следует, что 0 < 1 – α < 1, то α>α(1-α)>α(1-α)2… Другими словами, наблюдение , имеет больший вес, чем наблюдение , которое, в свою очередь, имеет больший вес, чем . Это иллюстрирует основное свойство модели экспоненциального сглаживания — коэффициенты при убывают при уменьшении номера k. Также можно показать, что сумма всех коэффициентов (включая коэффициент при ), равна 1.

Из приведенной формулы видно также, что значение является взвешенной суммой всех предыдущих наблюдений (включая последнее наблюдение ). Последнее слагаемое этой суммы является не статистическим наблюдением, а «предположением» (можно предположить, например, что). Очевидно, что с ростом t влияние , на прогноз уменьшается, и в определенный момент им можно будет пренебречь. Даже если значение α достаточно малое (такое, что (1 – α) приблизительно равно 1), значение будет быстро убывать.

Значение параметра α сильно влияет на функционирование модели прогнозирования, поскольку α представляет собой вес самого последнего наблюдения . Это значит, что следует назначать большее значение α в том случае, когда в модели наиболее прогностическим является именно последнее наблюдение. Если же α близко к 0, это означает практически полное доверие к прошлому прогнозу и игнорирование последнего наблюдения.

Перед маректологом возникает проблема: как наилучшим образом подобрать значение α? В этом поможет средство Поиск решения. Чтобы найти оптимальное значение α (т.е. такое, при котором прогнозная кривая будет менее всего отклоняться от кривой значений временного ряда), выполните следующие действия.


  1. Выберите команду Сервис -> Поиск решения.

  2. В диалоговом окне Поиск решения установите целевую ячейку G16 (см. лист «Экспо»), минимизируя значение.

  3. Укажите изменяемой ячейкой ячейку В1.

  4. Установите ограничения В1 > 0 и B1 < 1

  5. Запустите расчёт (см. результат, Рисунок ).




Рисунок . Пример 1 – результат использования метода экспоненциального сглаживания
Опять, как и в методе взвешенного скользящего среднего, наилучший прогноз будет получен, если назначить весь вес последнему наблюдению (см. Рисунок ). Следовательно, оптимальное значение α равно 1, при этом среднее абсолютных отклонений равно 6,82 (ячейка G16).
Рисунок . Пример 1 – график кривой прогноза методом взвешенного скользящего среднего и график реального объёма продаж
Метод экспоненциального сглаживания хорошо работает в ситуациях, когда переменная имеет стационарный характер, а ее отклонения от постоянного значения вызваны случайными факторами и не носят регулярного характера. Но этим методом, как и методами скользящего среднего не удастся спрогнозировать монотонно возрастающие или монотонно убывающие данные. Прогнозируемы значения будут всегда меньше или больше наблюдаемых, соответственно, а точность данных будет сравнима с точностью «наивного прогноза». Эти методы также не учитывают сезонных изменений показателя ряда.

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


Подбор кривой тренда

В качестве примера, воспользуемся данными объемов продаж из файла «ЛР6.Пример 2.Продажи.xls». Вначале построим точечную диаграмму, отображающие реальные объемы продаж. Чтобы построить по этим данным линию тренда, отражающую тенденцию в изменении объемов продаж, надо выполнить такие действия.



  1. Выберите любую точку выбранного ряда данных. В результате будут выделены все точки ряда.

  2. Вызовите через контекстное меню функцию «Добавить линию тренда».

  3. В диалоговом окне Линия тренда по умолчанию будет выбран линейный тип функции.

  4. Нажмите «Ок».

После этого на графике появится прямая линия тренда (см. Рисунок ).

Рисунок . Пример 2 – линия тренда на графике Excel

Для того чтобы осуществить прогноз, нужно в диалоговом окне Линия тренда отметить интересующий интервал времени в пункте «вперед на» (или «назад на»).

В меню Линия тренда можно также задать параметры подбираемой кривой. Например, он может быть экспоненциальной или полиномом заданной степени.

В рассматриваемом случае кривой тренда является прямая линия с уравнением y = ax + b. Коэффициенты a и b для этой кривой можно также найти с помощью надстройки Пакет Анализа, выбрав средство Регрессия.


Метод Хольта

Метод Хольта представляет собой развитие метода экспоненциального сглаживания, с учетом наличия тренда. Формулировка метода имеет вид





(8),

где



(9),



(10).

Метод Хольта позволяет прогнозировать на k периодов времени вперед. Метод, как видно, использует два параметра α и β, значения которых находятся в пределах от 0 до 1. Переменная L, указывает на долгосрочный уровень значений или базовое значение данных временного ряда. Переменная Т указывает на возможное возрастание или убывание значений за один период, т.е. на присутствие тренда.

Рассмотрим работу этого метода на следующем примере. Светлана работает аналитиком в большой брокерской фирме. На основе имеющихся у нее квартальных отчетов компании RusWind Airlines она хочет спрогнозировать доход этой компании в следующем квартале. Имеющиеся данные и диаграмма, построенная на их основе, находятся в файле «ЛР6.Пример 3.RusWindAirlines.xls» (см. Рисунок ).


Рисунок . Пример 3 - данные о доходе авиакомпании RusWind
Видно, что данные имеют явный тренд (почти монотонно возрастают). Светлана хочет применить метод Хольта, чтобы спрогнозировать значение прибыли на одну акцию на тринадцатый квартал. Для этого необходимо задать начальные значения для L и Т. Есть несколько вариантов выбора:


  1. L равно значению прибыли на одну акцию за первый квартал и T = 0;

  2. L равно среднему значению прибыли на одну акцию за 12 кварталов и T равно среднему изменению за все 12 кварталов.

Существуют и другие варианты начальных значений для L и Т, но Светлана выбрала первый вариант.

Она решила воспользоваться средством Поиск решения, чтобы найти оптимальное значение параметров α и β, при которых значение среднего абсолютных ошибок в процентах было бы минимально. Для этого нужно выполнить следующие действия.



  1. Выбрать команду Сервис -> Поиск решения.

  2. В диалоговом окне Поиск решения задать ячейку F18 целевой, минимизируя её.

  3. Выбрать диапазон ячеек В1:В2 изменяемыми.

  4. Добавить ограничения В1:В2 > 0 и В1:В2 < 1.

  5. Выполнить поиск.

Полученный прогноз показан отображают Рисунок , Рисунок , Рисунок . Как видно, оптимальными оказались значения α = 0,59 и β = 0,42, при этом среднее абсолютных ошибок в процентах равно около 38%.


Рисунок . Пример 3 - результат прогноза с использованием метода Хольта


Рисунок . Пример 3 - результат прогноза с использованием метода Хольта. Формулы

Рисунок . Пример 3 - результат прогноза с использованием метода Хольта. График
Однако метод Хольта, как и рассмотренные ранее простейшие методы прогнозирования, не учитывает наличие во временном ряде сезонных изменений.

Учет сезонных изменений

Спрос на значительное число товаров меняется в течение года. Например, если посмотреть на объемы продаж мороженого по месяцам, то можно увидеть в теплые месяцы (с июня по август в северном полушарии) более высокий уровень продаж, чем зимой, и так каждый год. Здесь сезонные колебания имеют период в 12 месяцев. Другой пример: анализируются еженедельные отчеты о количестве постояльцев, которые оставались на ночь в отеле, расположенном в бизнес-центре города. Предположительно можно сказать, что большое число клиентов ожидается в ночи на вторник, среду и четверг, меньше всего клиентов будет в ночи на субботу и воскресенье, и среднее число постояльцев ожидается в ночи на пятницу и понедельник. Такая структура данных, отображающая количество клиентов в разные дни недели, будет повторяться через каждые семь дней.

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


  1. На основе исходных данных определяется структура сезонных колебаний и период этих колебаний.

  2. Используя численный метод, описанный далее, из исходных данных исключают сезонную составляющую.

  3. На основе данных, из которых исключена сезонная составляющая, делается наилучший возможный прогноз.

  4. К полученному прогнозу добавляется сезонная составляющая.

Проиллюстрируем этот подход на данных об объемах сбыта угля (измеряемого в тысячах тонн) в США на протяжении девяти лет. Пусть некто Фрэнк работает менеджером в компании Gillette Coal Mine, и ему необходимо спрогнозировать спрос на уголь на ближайшие два квартала. Он ввел данные по всей угольной отрасли в рабочую книгу «ЛР6.Пример 4.Уголь.xls» и построил по этим данным график (см. Рисунок ).

Рисунок . Пример 4 - данные о спросе на уголь в компании Gillette Coal Mine

Определение структуры и периода сезонных колебаний. Из графика, приведённого выше видно, что объемы продаж выше среднего уровня в первом и четвертом кварталах (зимнее время года) и ниже среднего во втором и третьем кварталах (весенне-летние месяцы).

Исключение сезонной составляющей. Прежде всего необходимо вычислить среднее значение всех отклонений за один период сезонных изменений. Чтобы исключить сезонную составляющую в пределах одного года, используются данные за четыре периода (квартала). А чтобы исключить сезонную составляющую из всего временного ряда, вычисляется последовательность скользящих средних по T узлам, где T — продолжительность сезонных колебаний. Для выполнения необходимых вычислений Фрэнк использовал столбцы С и D, как показано на рисунке ниже. Столбец С содержит значения скользящего среднего по 4 узлам на основе данных, которые находятся в столбце В.

Далее требуется назначить полученные значения скользящего среднего средним точкам последовательности данных, на основе которых эти значения были вычислены. Эта операция называется центрированием значений. Если T нечетное, то первое значение скользящего среднего (среднее значений от первой до T-й точки) надо присвоить (T + 1)/2 точке (например, если T = 7, то первое скользящее среднее будет назначено четвертой точке). Аналогично среднее значений от второй до (T + 1)-й точки центрируется в (T + 3)/2 точке и т. д. Центр n-го интервала находится в точке (T+(2n-1))/2.

Если T четное, как в рассматриваемом случае, то задача несколько усложняется, поскольку центральные (средние) точки расположены между точками, по которым вычислялось значение скользящего среднего. Поэтому центрированное значение для третьей точки вычисляется как среднее первого и второго значений скользящего среднего. Например, первое число в столбце D отцентрированных средних (см. Рисунок ), слева равняется (1613 + 1594) / 2 = 1603.

Рисунок . Пример 4 – нахождение центрированного среднего по исходным данным


Графики исходных данных и отцентрированных средних показаны ниже (см. Рисунок ).
Рисунок . Пример 4 - график центрированного скользящего среднего
Далее найдём отношения значений точек данных к соответствующим значениям отцентрированных средних (см. Рисунок ). Поскольку точкам в начале и конце последовательности данных нет соответствующих отцентрированных средних (см. первые и последние значения в столбце D), такое действие на эти точки не распространяется. Эти отношения показывают степень отклонения значений данных относительно типового уровня, определяемого отцентрированными средними. Заметим, что значения отношений для третьих кварталов меньше 1, а для четвертых — больше 1.


Рисунок . Пример 4 – отношение данных к центрированным средним скользящим
Эти отношения являются основой для создания сезонных индексов. Для их вычисления группируются вычисленные отношения по кварталам, как показано в столбцах C—K (см. Рисунок ).


Рисунок . Пример 4 - отношения исходных данных к центрированным скользящим средним, сгруппированные по сезонам
Затем требуется найти средние значения отношений по каждому кварталу (см. столбец A, Рисунок ). Например, среднее всех отношений для первого квартала равно 1,108. Это значение является сезонным индексом первого квартала, на основе которого можно сделать вывод, что объем сбыта угля за первый квартал составляет в среднем около 110,8% относительного среднего годового объема сбыта.

Сезонный индекс — это среднее отношение данных, относящихся к одному сезону (в данном случае сезоном является квартал), ко всем данным. Если сезонный индекс больше 1, значит, показатели этого сезона выше средних показателей за год, аналогично, если сезонный индекс ниже 1, то показатели сезона ниже средних показателей за год.

Наконец, чтобы исключить из исходных данных сезонную составляющую, следует поделить значения исходных данных на соответствующий сезонный индекс. Результаты этой операции приведены в столбцах F и G (см. Рисунок ). График данных, которые уже не содержат сезонной составляющей, отражает Рисунок .



Прогнозирование без сезонной составляющей. На основе данных, из которых исключена сезонная составляющая, строится прогноз. Для этого используется соответствующий метод, который учитывает характер поведения данных (например, данные имеют тренд или относительно постоянны). В этом примере прогноз строится с помощью простого экспоненциального сглаживания. Оптимальное значение параметра α находится с помощью средства Поиск решения. Графики прогноза и реальных данных с исключенной сезонной составляющей приведены ниже (см. Рисунок ).


Рисунок . Пример 4 – сезонные индексы и данные без сезонной составляющей


Рисунок . Пример 4 - данные без сезонной составляющей


Рисунок . Пример 4 – графики прогноза и реальных данных с исключенной сезонной составляющей
Учет сезонной структуры. После нахождения прогноза без сезонной составляющей требуется учесть в полученном прогнозе (1726,5) сезонную составляющую. Для этого следует умножить 1726 на сезонный индекс первого квартала 1,108, в результате чего получим значение 1912. Аналогичная операция (умножение 1726 на сезонный индекс 0,784) даст прогноз на второй квартал, равный 1353. Результат добавления сезонной структуры к полученному прогнозу показан ниже (см. Рисунок ).

Рисунок . Пример 4 – итоговый результат прогнозирования с сезонной составляющей




Поделитесь с Вашими друзьями:
1   2   3   4   5   6   7


База данных защищена авторским правом ©psihdocs.ru 2017
обратиться к администрации

    Главная страница