markdown:
Самарский научный центр Российской академии наук
В.К. СЕМЁНЫЧЕВ, В.Н. КОЖУХОВА
АНАЛИЗ И ПРЕДЛОЖЕНИЯ МОДЕЛЕЙ ЭКОНОМИЧЕСКОЙ
ДИНАМИКИ С КУМУЛЯТИВНЫМ
ЛОГИСТИЧЕСКИМ
ТРЕНДОМ
Монография
Самара
Издательство «СамНЦ РАН»
2013
1
УДК 330
ББК 65В6
С30
Печатается по решению редакционно-издательского совета СамНЦ РАН Рецензенты:
Хмелева Г.А. – д.э.н., профессор Гераськин М.И. – д.э.н., профессор С 30 Семёнычев, В.К.
Анализ и предложения моделей экономической динамики с ку мулятивным логистическим трендом: монография / В.К. Семё нычев, В.Н. Кожухова. – Самара: Изд-во «СамНЦ РАН», 2013. 156 с.
ISBN 978-5-906605-10- В монографии представлены оригинальные результаты иссле дований моделей экономической динамики с кумулятивным (накоп ленным) логистическим трендом, и в том числе с дополнительными трендовыми и колебательными компонентами.
Рассмотрены логистические модели трендов с фиксированной и с произвольной асимметрией: как решения дифференциальных урав нений роста с разного рода ограничениями, так и феноменологиче ские. Впервые в отечественной литературе приведен «атлас моделей»
в функции параметров, исследована точность идентификации моде лей в зависимости от объема выборки и положения точки перегиба как при растущем тренде, так и при падающем.
Наряду с известными логистическими моделями трендов приве дены их расширения и новые модели, даны многочисленные приме ры известных и новых приложений с высокой точностью моделиро вания и прогнозирования на относительно малых выборках. Ориги нальны предложенная методика оценки точности моделирования и прогнозирования, области приложения рассмотренных моделей и методов их идентификации.
Монография ориентирована на студентов экономических специ альностей вузов, аспирантов и практикующих аналитиков.
ISBN 978-5-906605-10- ББК 65в © Семёнычев В.К., Кожухова В.Н., СОДЕРЖАНИЕ ВВЕДЕНИЕ...................................................................................................................................... ГЛАВА 1. КУМУЛЯТИВНЫЕ МОДЕЛИ ЛОГИСТИЧЕСКОЙ ДИНАМИКИ, ИХ ПРИЛОЖЕНИЯ, ИСПОЛЬЗУЕМЫЙ ИНСТРУМЕНТАРИЙ ИДЕНТИФИКАЦИИ И ЗАДАЧИ ЕГО РАЗВИТИЯ........... 1.1. Использование моделей с логистическим трендом для социально-экономических приложений.......................................................... 1.2. Формирование логистических моделей трендов как решений дифференциальных уравнений роста СЭС.................................. 1.2.1. Логистические модели трендов с фиксированной асимметрией............................................................................................................................... 1.2.2. Логистические модели трендов с произвольной асимметрией............................................................................................................................... 1.3. Феноменологические логистические модели трендов...................... 1.4. Актуальные задачи развития инструментария идентификации динамики с логистическим трендом....................................... ВЫВОДЫ ПО ПЕРВОЙ ГЛАВЕ...................................................................................... ГЛАВА 2. ИССЛЕДОВАНИЕ ТОЧНОСТИ МОДЕЛЕЙ РЯДОВ ДИНАМИКИ С ЛОГИСТИЧЕСКИМ ТРЕНДОМ И ОБОСНОВАНИЕ ВЫБОРА МЕТОДОВ ИХ ИДЕНТИФИКАЦИИ..........................................................
Экономическая система определяется обычно как сложная ве роятностная система, охватывающая процессы производства, обмена, распределения и потребления материальных благ [40].
В СЭС, наряду с материально-вещественными потоками, учиты ваются и производственные отношения людей в обществе, воспроиз веденные и развитые экономической системой производственные отношения, социально-экономические, социально-политические, со циокультурные и демографические процессы. В СЭС рассматривают и вопросы инвестиционной, финансовой, социальной и/или экологи ческой политики, ценообразования на рынке как системы отноше ний купли-продажи между продавцами и покупателями, распредели тельных отношений, лагов и т.д.
СЭС в каждый момент времени характеризуется множеством показателей макро - (страна), мезо - (регион страны, отрасль хо зяйства, муниципальное образование) и микроуровней (конкрет ное предприятие, домохозяйство).
Показателями являются, например, доход, потребление, оборот розничной торговли, уровень безработицы, средняя заработная пла та, объём инвестиций и др. Динамика конкретных СЭС описывается единицами, десятками, а в отдельных случаях и сотнями показателей.
Будем рассматривать динамику СЭС, которая полностью вмеща ет в себя проблемы статики: оптимальное распределение и исполь зование имеющихся производственных ресурсов, максимальное удовлетворение сложившихся общественных и индивидуальных по требностей, анализ структуры и взаимосвязей хозяйства страны, ба лансирование производства и потребления, изучение допустимых и рациональных состояний и др.
В рамках анализа динамики показателей СЭС осуществляется более общий анализ указанных проблем, включая и воспроизводст венный подход к возможным траекториям деградации или развития, к последовательностям состояний и переходам от одних состояний к другим.
Показатель динамики СЭС определим как ее обобщающую коли чественную характеристику, в общем случае вектор Y, в конкретных условиях места и времени.
Показатель называют также траекторией, результирующей, оп ределяемой, объясняемой (explained variable), внутренней или эндо генной переменной СЭС. Эндогенная переменная Y является принципиально случайной (стохастической), так как климатические и природные явления, структура материальных и духовных потребностей членов общества могут быть определены только с некоторой вероятностью. Слож ность и динамичность реальных социально-экономических процес сов приводят к тому, что затраты на производство, экономический эффект, производительность труда, результаты научных исследова ний и разработок, эффективность новой техники и т.д. поддаются предварительному расчёту и экономическому анализу только с тем или иным уровнем достоверности. В общем случае определяющий показатель Y формируется в процессе и внутри функционирования СЭС под воздействием некото рого вектора x : числа i других процессов или явлений X i, называе мых факторными, внешними, экзогенными (exogenous variable) или определяющими переменными, которые (или часть которых) подда ются регистрации и планированию. Каждый из факторов при анализе динамики Y может быть так же функцией физического времени t.
Экзогенные переменные описывают условия функционирова ния СЭС, могут задаваться извне анализируемой системы. Примерами экзогенных переменных могут быть физическое время, налоги, госу дарственные закупки товаров, устанавливаемая цена на благо, рента, налоги, банковские ставки, официальный курс доллара, численность экономически активного населения и др.
Современным и плодотворным методом исследования динами ки показателей СЭС является моделирование (model-building) и про гнозирование (forecasting).
Моделью динамики показателей СЭС будем считать образ про цесса, явления в форме математических соотношений, отражающий существенные свойства моделируемой СЭС и замещающий его в ходе исследования и управления.
Модели могут быть различными в зависимости от формулиров ки цели моделирования и прогнозирования, размера выборки, ие рархического уровня СЭС и т.д. В качестве моделей траектории будут выступать континуальные (непрерывные) и дискретные аналитиче ские выражения.
При декомпозиции динамической траектории, а это один из важнейших подходов в моделировании динамики показателей СЭС, на отдельные компоненты к модели будем относить как аналитиче ские выражения и характеристики (амплитудные, временные) ком понент, так и структуру (аддитивную и/или мультипликативную, смешанную) их взаимодействия.
Декомпозиция представляет динамику траектории как супер позицию медленного детерминированного процесса T (t ) (векового уровня, secular trend, главной тенденции, эволюторной компоненты) тренда и более быстрых процессов – детерминированной колеба тельной C (t ) и стохастической компонент (t ).
Последнюю компоненту называют также нерегулярной, шумом, помехой (irregular component). Именно присутствие в структуре ряда стохастической компоненты делает анализируемую траекторию принципиально случайной.
В анализе динамики СЭС в сравнении с анализом статики по вышаются требования к точности моделирования и прогнозирова ния на коротких интервалах наблюдения (на коротких выборках), так как эволюция процессов и явлений ведёт к нестационарности мо делей показателей (аналитических выражений моделей и/или их па раметров, характера их взаимодействия компонент).
Огромное значение при этом приобретает мониторинг эволю ции моделей, использование их для прогнозного моделирования воз можных траекторий развития, для оценки эффективности принятых или возможных управленческих (технологических, маркетинговых) решений.
Эволюция выражается обычно в высоких темпах спада или рос та показателей (неестественных с точки зрения стабильных эконо мик), появлении или исчезновении колебательной компоненты, эво люции ее параметров, изменении характера взаимодействия компо нент ряда динамики показателей, увеличении мощности и/или по явлении гетероскедастичности (нестационарности дисперсии) поме хи.
Модели динамики СЭС принципиально сложнее моделей стати ки не только за счет введения в них дополнительного параметра – времени. Они, как правило, приобретают нелинейный характер и по переменным, и по параметрам.
Под логистическим законом изменения тренда показателя СЭС понимаем случай, когда он проходит в своем развитии несколько стадий. Первая из них является стадией медленного, близкого к экспоненциальному по своему характеру, роста. Вторая стадия име ет характер почти линейного роста. Третья стадия – стадия замед ляющегося, близкого к гиперболическому, роста, стремящегося к постоянному уровню (уровню насыщения).
Чаще в экономической практике рассматривают растущие ло гисты, однако сложный эволюционный характер изменения может иметь и уменьшающаяся (падающая) логистическая динамика: с начального уровня, асимптотически уменьшаясь, она также прохо дит три стадии до уровня спада.
Графики растущей и падающей логистических траекторий тренда приведены на рисунке В.1.
Tt уровень насыщения точка перегиба уровень спада 0 t Tt I Tt II Рис. В.1. Графики логистических траекторий В качестве аргумента логистической траектории чаще рас сматривается время, но также может быть и другой экономический показатель.
Например, тренд T I на рисунке В.1 может отражать характер динамики продаж товара (услуги) в функции рекламного бюджета.
В этом случае можно говорить о логистической пространственной траектории, когда дополнительные затраты на рекламу практиче ски не будут увеличивать продажу.
Казалось бы, эволюционирующую динамику показателя СЭС можно характеризовать набором трех различных моделей (экспо нентой, линейной функцией и гиперболой) с трудноразличимыми точками трансформаций (преобразований модели).
Однако существуют теоретические возможности и практиче ское подтверждение использования на всех трех стадиях единой аналитической модели. Кривые, которыми описывают подобный характер динамики рядов динамики показателей, называют логи стическими, S-образными кривыми, сигмоидами, логистами [10, 12, 16, 37, 42, 50, 51, 54-66, 80-82].
В растущих логистах первая производная неотрицательна:
растет от нуля (как правило, при нулевом значении аргумента), достигая при уровне насыщения вновь нулевого значения.
В падающих логистах первая производная отрицательна, она изменяет свою величину от нулевого значения вновь до нуля при достижении уровня спада.
Особенностью логистических кривых является наличие точки перегиба, которая соответствует максимуму ее первой производной и нулевому значению второй. В маркетинге точка перегиба назы вается обычно точкой нерасширяемого спроса.
Уровень насыщения логистической кривой характеризует ем кость рынка или потенциал отрасли.
В отечественной эконометрической литературе логистиче ские кривые сравнительно мало изучены, рассматриваются редко и лишь в перечислении других моделей нелинейной регрессии [1, 8, 10, 11, 37, 42, 81, 82].
При этом, как правило, рассматривают лишь две наиболее из вестные (еще с 1825-1838 гг.) модели логистической динамики Ферхюльста (в другой транскрипции Верхулста;
ее называют также моделью Перла-Рида) и Гомпертца (Гомперца).
Для параметрической идентификации этих моделей обычно рекомендуют использовать методы, которые предполагают апри орное экспертное назначение некоторых параметров (обычно уровня насыщения или спада).
Указанные методы идентификации, однако, не удовлетворяют современным требованиям по точности моделирования и прогно зирования, а также по возможностям существующего математиче ского аппарата эконометрики и его компьютерной поддержки.
Кроме того, зачастую интерес для исследователя представля ет определение оценки уровня насыщения (спада) логистической кривой, что существенно снижает возможную область применения этих методов.
Крайне мало в известной научной литературе публикаций, в которых сравниваются возможные методы идентификации и обос новывается их выбор по достигаемой ими точности моделирова ния и прогнозирования в зависимости от положения выборки на ло гистической кривой, а также от длины выборки, соотношения дис персий стохастической компоненты и детерминированной компо ненты (выборки, содержащей логистический тренд).
В наибольшей мере это актуально для обеспечения точности прогнозирования, которая как конечная цель почти любого моде лирования представляет основной интерес для многих приложе ний.
В зарубежной литературе вопросам моделирования логисти ческими функциями посвящено существенно большее количество публикаций, особенно в последние годы. Увеличилось и количество используемых для приложений в технике, экономике, биологии, при моделировании социальной динамики логистических моделей тренда. Но и в данных публикациях недостает результатов для от ражения многообразия логистических трендов, обоснования мето дов идентификации и области их применения.
Актуально решение задачи исследований логистических трендов с произвольной относительно точки перегиба асимметри ей, которые моделируют динамику рынков с преобладанием облас ти расширяемого спроса (асимметрия «влево», как, например, в функции Гомпертца) и нерасширяемого спроса (асимметрия «впра во»).
Если ордината точки перегиба равна половине уровня насы щения кривой, то такая функция является симметричной (как, на пример, функция Ферхюльста).
Большинство разработанных зарубежными учеными моделей логист и методов их идентификации [87, 101, 103, 112] применя лось для отдельных частных случаев приложений [97, 106, 132, 133]. Они конструировались под конкретные задачи моделирова ния (редко прогнозирования) на больших выборках (до 50-100 на блюдений), при малой зашумленности исходных данных (в преде лах 5-10% соотношения дисперсий помехи и детерминированной компоненты) [80, 93] и только лишь при конкретном (определен ном данным приложением) наборе параметров моделей.
Возможность применения тех или иных моделей логистиче ских трендов, методов их идентификации при другом наборе пара метров или при других соотношениях дисперсий поме хи/детерминированной компоненты в известной научной литера туре практически не рассматривалась.
Отметим также малое число исследований по точности анали за (под ним будем понимать моделирование и прогнозирование) падающих логист, что можно объяснить малой точностью извест ных методов их идентификации.
Не анализировалась и возможность идентификации рядов с логистическим трендом при наличии других дополнительных трендов (мультитрендовых рядов) и/или колебательных компо нент, на наличие которых указывает визуальный анализ многих рядов динамики в приложениях.
Считаем, что учет многокомпонентности реальных рядов эко номической динамики позволит увеличить точность идентифика ции как самой логистической динамики (тренда), так и других де терминированных компонент ряда, обеспечивая тем самым боль шую точность прогнозирования траекторий.
Еще один положительный результат рассмотрения мульти трендовых рядов видится в возможности конструирования ими бо лее сложной модели логистической динамики (билогистической или мультилогистической) путем использования алгебры (сложе ния и/или умножения) взаимодействия отдельных трендов.
Можно заключить, что до настоящего времени отсутствует инструментарий (комплекс логистических моделей, методов и приемов их идентификации), допускающий применение для широ кого круга приложений, который бы обеспечивал высокую точ ность моделирования и прогнозирования на относительно малых выборках в различных (широких) диапазонах сочетаний парамет ров моделей, соотношениях дисперсий помехи и детерминирован ной компоненты.
Практически неизвестны и методики исследований точности такого инструментария. В отечественной литературе можно упо мянуть лишь одну попытку такого оценивания [80], проведенную на единственной тестовой выборке для модели Ферхюльста, а в за рубежной литературе исследование в широком диапазоне сочета ний параметров проведено также для той же модели [93].
Теоретическую базу эконометрического моделирования и про гнозирования логистической динамики составили труды таких оте чественных ученых, как С.А. Айвазян, В.Н. Афанасьев, В.В. Давнис, Т.А. Дуброва, И.И. Елисеева, Е.В. Зарова, Н.Ш. Кремер, В.С. Мхитарян, Ю.П. Лукашин, В.П. Носко, А.А. Пересецкий, С.Г. Светуньков, В.К. Семёнычев, В.Л. Сергеев, Н.П. Тихомиров, Г.Р. Хасаев, А.А. Цып лаков, Е.М. Четыркин и др.
Использованы научные результаты исследований и зарубеж ных ученых: Р. Айреса, Ф. Басса, Дж. Бокса, А. Грюблера, Г. Дженкинса, К. Доугерти, Р.А. Кашьяпа, В. Мааджана, П. Майера, Т. Модиса, Н. Накиченовича, Дж. Рамсея, Ф. Ричардса, К. Скиадаса, Г. Тейла и др.
Целью монографии является представление результатов ис следований комплекса моделей с логистическим трендом, предло жения по их развитию;
определение области возможного примене ния;
сравнение методов и приемов их идентификации, обеспечи вающих повышение точности моделирования и прогнозирования на относительно малых выборках.
Для достижения указанной цели были поставлены и решены следующие задачи:
– показать область применения логистических кривых при анализе социально-экономической динамики, определить и срав нить их основные характеристики (начало и завершение стадий роста, точек перегиба, уровня насыщения, симметрии и асиммет рии);
– дополнить известные логисты новыми моделями, обеспе чив для всех них возможность идентификации при включении в логистические модели дополнительных трендовых и колебатель ных компонент;
– разработать для моделей логист методы идентификации на относительно коротких выборках (до 30-35 наблюдений), методику исследования их точности в различных диапазонах изменения па раметров моделей и соотношениях дисперсий детерминированной компоненты и помехи;
– реализовать программные средства для идентификации рассматриваемого комплекса моделей, сравнения и реализации ме тодик исследования точности;
– применить модели и методы на реальных, важных для при ложений, выборках социально-экономической динамики показате лей микро-, мезо - и макроуровня агрегирования.
Социально-экономические процессы эволюции показателей анализировались авторами при помощи 24 эконометрических мо делей с трендами на основе логистических функций и дополни тельных трендовых и/или колебательных компонент.
При этом получены следующие результаты:
- впервые в отечественной научной литературе обобщены и сравнены результаты по этим моделям логистической динамики, которые наглядно представлены атласами моделей в функции па раметров;
- предложены обобщения логистической модели Рамсея на случай смещения вдоль оси абсцисс и ординат, логистической мо дели Гомпертца для случая асимметрии вправо, а также предложе ны новые модели билогистической динамики;
- для моделей с растущим и падающим логистическими трен дами, с фиксированной и произвольной асимметрией обоснованы методы идентификации параметров, дающие высокую точность моделирования и прогнозирования в функции соотношения дис персий помехи и детерминированной компоненты. Показаны воз можности точечной оценки точности параметров и моделей для тестовых выборок, а также точечной и интервальной оценки точ ности на реальных выборках;
- предложенные методы идентификации позволяют анализи ровать многокомпонентные ряды динамики с логистическим трен дом в случаях, когда логистическая тенденция дополнена другими трендами и/или колебательными компонентами.
Практическая ценность полученных результатов заключается в следующем:
- приведены примеры достижения высокой точности модели рования и прогнозирования динамики реальных показателей соци ально-экономических систем для разных уровней их агрегирова ния, при произвольной асимметрии трендов, на относительно ко ротких и многокомпонентных выборках;
- даны иллюстрации точности на сотнях тысяч тестовых вы борок, рекомендации по выбору участков выборок логистических трендов.
ГЛАВА 1. КУМУЛЯТИВНЫЕ МОДЕЛИ ЛОГИСТИЧЕСКОЙ ДИНАМИКИ, ИХ ПРИЛОЖЕНИЯ, ИСПОЛЬЗУЕМЫЙ ИНСТРУМЕНТАРИЙ ИДЕНТИФИКАЦИИ И ЗАДАЧИ ЕГО РАЗВИТИЯ 1.1. Использование моделей с логистическим трендом для социально-экономических приложений Логистическая динамика была изначально обнаружена при рассмотрении модели роста численности населения.
Идея прогнозирования численности популяций по логистиче ской кривой принадлежала Б. Гомпертцу (1825) [98], чуть позже – П.Ф. Ферхюльсту и А. Кетле (1838) [137]. Дальнейшее развитие тео рия получила в трудах американских ученых Р. Перла и Л. Рида (1920) [123, 124].
Логистическими кривыми до сих пор описывают численность населения: например, логистический характер имеет численность населения г. Орландо1 (штат Флорида, США) по десятилетиям с 1890 по 2010 гг.
Динамика этого показателя прошла точку перегиба (2010 г. 200 тыс. чел.), но еще не достигла уровня насыщения.
Еще одна область применения логистических кривых – описа ние процессов, происходящих в развитии технологий. В соответст вии с логистической кривой эффективность технологии растет медленно, затем набирает темп и ускоряется [76].
В конечном итоге рост технологического параметра неизбеж но останавливается. Таким образом, чтобы вовремя перейти к но вой технологии, необходимо знать о технологическом пределе ста рой.
Время перехода от одной технологии к другой называется технологическим разрывом (рис. 1.1).
Технологический разрыв характеризует различие в эффек тивности P (параметре) новой и старой технологий P P2 P, а также разницу в объемах средств K K 2 K1, необходимых для инвестирования в новую технологию.
URL: http:// www.census.gov.
Рис. 1.1. Технологический разрыв Технологический разрыв можно рассматривать и во времени t t2 t1, которое должно пройти для смены технологии. Наступает момент, когда инвестировать средства в совершенствование новой технологии выгоднее, чем в совершенствование старой [6, 43, 69, 78].
Часто при построении нескольких логистических кривых сме ны технологий оказывается, что они составляют огибающую кри вую (рис. 1.2), также изменяющуюся по логистическому закону [85, 117].
показатель развития техноло гии время Рис. 1.2. Огибающая логистическая кривая смены технологий [109] Логистический характер носит и динамика развития инфра структуры и промышленности [86, 94, 106, 116]. На рисунке 1. приведена динамика развития сетей транспорта и коммуникаций в США, подчиняющаяся логистическим закономерностям. Видно, что способы сообщения транспорта и коммуникаций сменяют друг друга с периодом примерно в 55 лет. На том же рисунке видно при сутствие колебательной компоненты в траекториях с логистиче ским трендом.
Отмечено, что множество видов инфраструктуры в странах мира (не только в США) также подчиняются логистическим зако номерностям.
56 лет 55 лет Рис. 1.3. Динамика развития инфраструктуры США В экономике и социологии также существует много примеров динамики показателей, которые в своем развитии проходят после довательно несколько стадий эволюции (развития или деграда ции).
Под эволюцией в данном случае понимаем любое изменение характера динамики показателя, например изменение скорости роста населения какого-либо региона с течением времени, измене ние основной тенденции потребительского спроса на товар после экономического кризиса и т.д.
Так, маркетологи часто сталкиваются со сложным характером динамики продаж при выведении нового товара на рынок и приня тии его потребителями.
Обычно на первом этапе товар пользуется спросом лишь у не большой части потребителей, на втором этапе все большее число потребителей осваивает новинку, а на третьем этапе число потре бителей становится относительно постоянным – определяется так называемая емкость рынка для данного товара, т.е. потенциально возможный объем продаж товара на рынке.
Известен случай применения логистических кривых при опи сании спроса на товары Q в зависимости от дохода потребителей M (рис. 1.4) на основе предположения, что цены и другие факторы фиксированы [39].
Данный подход позволяет оценить изменение коэффициента M эластичности, по определению равного E Q '( M ).
Q Q Q Q Q M Рис. 1.4. Логистические функции спроса на товар в зависимости от дохода в менее обеспеченных (1) и более обеспеченных (2) группах населения При увеличении дохода потребителя спрос на многие товары с какого-то момента времени начинает замедляться и в конце кон цов останавливается на некотором уровне.
Поэтому можно утверждать, что в более обеспеченных груп пах населения, где уровень удовлетворенного спроса более высок, склонность населения увеличивать потребление товара (при одном и том же относительном приросте доходов) будет ниже, чем в ме нее обеспеченных группах населения.
Другими словами, при одинаковом доходе потребителей зна чение коэффициента эластичности E в более обеспеченных группах увеличивается меньше, чем в менее обеспеченных группах: в точке Q2 на рисунке 1.4 коэффициент эластичности ниже, чем в точке Q1.
Кроме того, в различных группах при одинаковом уровне по требления реакция на изменение дохода будет также различной (точки Q2 и Q3). В более обеспеченной группе реакция на рост дохо да может оказаться более сильной, поскольку потребление одина ково.
Но так как у более обеспеченной группы существуют большие возможности увеличивать потребление данного товара, то в точке Q3 коэффициент эластичности выше, чем в точке Q2.
Известна статистическая модель изменения цен на жилье на примере г. Москвы, основанная на логистической закономерности роста цен (в долларовом эквиваленте) при переходе от централи зованной к рыночной экономике и корреляционной связи между приростом месячных темпов роста курса доллара и долларовых цен на жилье [73].
В приведенных примерах определяемыми показателями яв ляются численность населения, уровень спроса на товары или объ ем продаж товаров.
Подобные примеры могут иллюстрировать развитие показа телей не только во времени, но и в пространстве (к примеру, как уже указывалось выше, отражать зависимость исследуемого пока зателя от значений других показателей, например уровень сбыта продукции в функции рекламного бюджета фирмы) [17].
В монографии будет рассматриваться динамика показателей во времени t (временная динамика), хотя ряд результатов может быть применен и для пространственной динамики.
Показатель (или иначе траекторию, результирующую, опреде ляемую, эндогенную, внутреннюю, зависимую переменную) в даль нейшем будем обозначать как Y.
Как правило, значения рассматриваемого показателя форми руются под воздействием множества факторов (называемых внеш ними, экзогенными, определяющими переменными).
Эти факторы далеко не всегда удается учесть, поскольку во многих случаях имеет место сравнительно большое число опреде ляющих переменных, а также вследствие нелинейного характера действия факторов на определяемый показатель, малости исход ных выборок и задачи прогнозирования на малый горизонт про гноза (упреждения).
Большое количество факторов и принимаемый обычно ли нейный характер их связи с определяемым показателем создают зачастую лишь «иллюзию точности» моделирования. Редукция размерности и переход к нелинейным связям сложны.
Однако фактор времени как единственная экзогенная пере менная позволяет отразить совокупное (интегрирующее) влияние всего множества факторов. Именно поэтому на практике широко применяются модели временной динамики, рассматриваемые в данной работе. Зафиксированные в разные моменты времени зна чения показателя, собранные в хронологическую последователь ность, образуют временной ряд, или ряд динамики.
Чаще всего используются эквидистантные временные ряды, когда значения показателя фиксируются через одинаковые перио ды опроса (день, неделя, месяц, квартал, год).
Для таких регулярных наблюдений непрерывная переменная времени t может быть заменена дискретной t k, образуя ряд на блюдений Yk, где k 0,1, 2,..., n 1 – номера наблюдений, n – объем выборки.
Здесь и далее для случаев непрерывного времени будет ис пользоваться переменная времени t, для случаев дискретного вре мени – номер наблюдения k.
Если динамика показателей (технических, природных, соци ально-экономических) отражает процесс перехода от одного ста бильного состояния к другому, сопровождающийся изменением скорости роста или спада показателя, а также имеет S-образную форму, то такая динамика носит название логистической, или сиг моидальной траектории.
Положение точки перегиба логистической кривой определяет, равны ли во времени стадии медленного роста и угасания (в случае симметричного расположения точки перегиба на уровне половины насыщения) либо одна из них больше другой.
Как будет показано в п. 1.2, 1.3 монографии, логистические траектории обладают еще несколькими критическими точками, расчет которых всегда представляет интерес для исследователя и может быть осуществлен либо аналитически, либо численно.
Модели логистической динамики относятся к нелинейным моделям, для них характерно ограничение каким-либо пределом, а скорость роста такой динамики изменяется с течением времени (у функции существует точка перегиба).
Применение логистических кривых «правильнее» некоторых кривых роста (например линейной, экспоненциальной функций) по той причине, что в долгосрочном периоде (т.е. при прогнозирова нии на срок свыше 10 лет) рост ни одного из указанных показате лей не может быть бесконечным в силу реальной ограниченности динамических диапазонов (ресурсов) любых показателей объекта анализа.
Это же относится и к убывающей динамике – она также стре мится к некоторому пределу [104].
Упомянутые выше другие кривые роста (компоненты логи сты) могут быть использованы для прогнозирования лишь в крат косрочном периоде (т.е. при прогнозировании на срок до одного го да или с горизонтом до одного периода дискретизации).
Ограничения на рост (убывание) показателей существуют в живой и неживой природе. Например, на динамику добычи полез ных ископаемых, на рост популяций растений и животных, на чис ленность людей и т.д. влияют природные ограничения.
На социальные показатели влияют происходящие в стране со циально-политические процессы, а также отношение общества к потребляемым товарам и услугам, например к услугам Интернета и сотовой связи.
Последние, в свою очередь, определяются и неживой приро дой: развитие технологии ограничено влиянием конкурентных разработок, наличием материальных ресурсов и т.д.
Обратим внимание и на известную гипотезу о том, что сим метричные логистические тренды в большей мере характерны для явлений и процессов неживой природы и техники, а асимметрия преимущественно относится к живой природе, к социальным явле ниям и процессам (понимаемым в широком смысле, предполагаю щим не только сообщества людей, но и растений, животных).
Перечисленные особенности логистических кривых опреде ляют область их применения: как правило, это глобальная кумуля тивная (накопленная, суммарная) динамика, в том числе динамика макроэкономических показателей.
Можно сказать, что логистическая динамика в макро масштабе отражает кумулятивную эволюцию развития экономики при ограничениях внешней среды.
Отметим также еще одно важное приложение логистических кривых и их производных, заключающееся в том, что они модели руют жизненный цикл продукта (товара, услуги, организации, бренда).
Примером может быть динамика продаж товара/услуги после выведения продукта на рынок. При этом логистическая кривая (S-curve) отражает кумулятивную динамику, т.е. описывает накоп ленные, суммарные продажи к текущему моменту, а первая произ водная логистической функции является импульсной кривой и опи сывает текущие продажи товаров.
Импульсная кривая (bell-shaped curve) также известна своим применением, прежде всего, при описании добычи невозобновляе мых природных ресурсов (нефти, газа, угля, золота и т.д.), при мо делировании многих социальных процессов [25, 29, 66, 91].
Возможны случаи, когда в социально-экономической системе параллельно развиваются несколько логистических процессов. На пример, второй логистический процесс чаще начинается в момент достижения первым логистическим процессом 50% и более про центов роста, или когда они начинаются практически одновремен но, но динамика роста различна.
Билогистическая (в более общем случае - мультилогистиче ская) динамика имеет место, когда одна логистическая функция, достигнув уровня насыщения, переходит в другую, то есть верхняя горизонтальная асимптота первой логисты является нижней гори зонтальной асимптотой для другой.
Приведем некоторые известные примеры билогистических и мультилогистических процессов, представленные в западной науч ной литературе [100, 115]: средний рост мальчиков (второй про цесс запускается после 10-летнего возраста [139]);
рост числа аме риканских университетов (первый импульс начался с 1800 г., вто рой - с 1950 года, когда первый достиг 95% уровня насыщения [140]);
генерация электрической мощности в США (первый процесс начался примерно в 1910 г., пришел к насыщению в 1926 г., второй начался в 1940 г., оказался длиннее и существенно выше по уровню насыщения, достиг его к 2000 году, после чего начался третий про цесс [121]).
Представленные примеры ставят перед разрабатываемым ин струментарием задачу разделения отдельных логистических трен дов, взаимодействующих во времени (мультитрендовых логист).
При этом можно предположить и относительно простое алгебраи ческое взаимодействие отдельных логистических трендов – адди тивное или мультипликативное.
Актуальными сферами применения моделей логистической динамики являются описание аудитории СМИ, изменение числа семей, имеющих радио и телевизоры, численность аудитории сети Интернет.
В общем случае логистическая тенденция может дополняться другими составляющими ряда динамики: дополнительными трен дами и/или (для временных рядов) колебательными компонента ми.
Примером может быть аппроксимация динамики Рунета про стой экспонентой с использованием исторических данных о разви тии в прошлые годы, верна только для первых 2-3 месяцев [14]. По следнее время рост числа пользователей Рунета замедлился. В ка честве примера на рисунке 1.5 приведена динамика аудитории Ру нета г. Москвы1.
Данные приведены помесячно в % от численности населения г. Москвы в возрасте от 12 лет и старше.
URL: http:// www.etarget.ru.
Июл ' Июл ' Июл ' Июл ' Июл ' Мар' Мар ' Мар ' Мар ' Мар ' Мар ' Ноя ' Ноя ' Ноя ' Ноя ' Ноя ' Ноя ' Рис. 1.5. Динамика аудитории Рунета г. Москвы Обратим внимание на то, что падающие логистические кри вые не нашли до настоящего времени столь широкого применения, как растущие, тем не менее приведем ряд характерных примеров.
Поскольку такая логистическая кривая берет свое начало из биологии, то отметим, что в природе популяции имеют не только максимальную численность, определяемую величиной экологиче ской ниши (в экономике – емкость рынка), но и минимальную кри тическую численность (в экономике это можно интерпретировать как минимальный уровень продаж).
Интересен пример описания падающими логистическими кривыми стоимости 1МВт энергии ветра в Великобритании (в це нах 2000 года) [121]. Такой же характер носит динамика производ ства первичной стали и добычи угля в ряде развитых стран мира [133, 127].
Декомпозиционный подход, широко используемый при ана лизе рядов динамики, позволяет разделить исходный ряд на более простые составляющие (компоненты).
Обычно в структуре временного ряда выделяют следующие компоненты [8, 42, 82]:
– тренд Tt (главную тенденцию, вековой уровень);
– циклические колебания Ct ;
– сезонные колебания St ;
– стохастическую компоненту t (шум, помеху, нерегулярную компоненту).
Тренд, циклическая и сезонная компоненты образуют детер минированную (неслучайную) компоненту Dt. Уровни наблюдений сезонной и циклической компонент (их порой объединяют терми ном колебательная компонента) обычно меньше уровней тренда, а динамика (скорость изменения уровней) больше, чем у тренда.
Те или иные компоненты (кроме стохастической) ряда дина мики могут отсутствовать как принципиально, в силу причинно следственных связей, так и в зависимости от длины интервала на блюдения (объёма выборки). Считается, что чем короче интервал наблюдения, тем проще может быть модель, адекватная выборке.
Наблюдения показателя Yt являются случайными из-за нали чия наиболее динамичной и существенно меньшей в сравнении с трендом и колебательной компонентой стохастической компонен ты t, которая имманентно присутствует в наблюдениях, в то время как некоторые другие компоненты ряда могут и отсутствовать.
Следует, однако, помнить, что соотношение уровней сезонной и циклической компонент на практике может изменяться в широ ком диапазоне значений, поэтому оно может не являться критери ем декомпозиции. Декомпозицию можно проводить и исходя из свойства периодичности и цикличности выделяемых компонент.
Под цикличностью понимают повторяемость явления в общих чертах, а под периодичностью – частный случай цикличности, ко гда явление повторяется в деталях на каждом следующем цикле.
Трендовая компонента не обладает свойствами цикличности и/или периодичности. В первом приближении обычно считают, что колебательная (во всяком случае сезонная) компонента периодич на.
Ряды динамики, состоящие из тренда, сезонной и стохастиче ской компонент, принято называть тренд-сезонными, а ряды, со стоящие из тренда, циклической и стохастической компонент, тренд-циклическими.
В структуре ряда динамики можно дополнительно учитывать и другие компоненты:
календарную компоненту, происхождение которой связано с разным количеством дней в месяцах года продолжительностью, национальными и религиозными праздниками;
выбросы: аномальные движения временного ряда, связан ные с редко происходящими событиями, которые резко и кратко временно отклоняют ряд от общего закона движения, например авария на атомной станции в Японии, выброс нефти в Мексикан ском заливе и др.;
структурные сдвиги: аномальные движения временного ря да, связанные с редко происходящими событиями (например с из менением ставок налогов, ставок акцизных сборов, существенным изменением климатических условий), имеющие скачкообразный характер и изменяющие тенденцию;
инфляционную компоненту.
Обычно анализируемые выборки перед моделированием и прогнозированием «очищаются» от выбросов, календарной и ин фляционной компонент [44, 81], поэтому здесь они рассматривать ся не будут. Анализ структурных сдвигов возможен [81, 82].
Декомпозицию можно применить как к определяемым пока зателям динамики СЭС при многокомпонентной структуре ряда, так и к моделям отдельных компонент ряда.
Компоненты Tt, Ct, St, t ряда могут взаимодействовать раз личными способами, образуя разные структуры модели этого ряда при декомпозиции.
В силу того что наблюдаемыми являются только уровни Yt оп ределяемой переменной ряда, а Tt, Ct и St ненаблюдаемы (unob served components), их выбор в общем случае неоднозначен как по моделям компонент, так и по структуре их взаимодействия и явля ется в определенной мере искусством.
Наиболее часто предполагают аддитивную структуру ряда Yt Tt Ct St t или реже – мультипликативную Yt Tt Ct St t.
Эти структуры называют также классическими, или канони ческими.
Аддитивная структура взаимодействия компонент позволяет использовать для них параметрическую (аналитическую) запись всех компонент, отражает их независимость между собой, является наиболее распространенной в силу простоты идентификации и ис пользования для выполнения прогнозирования.
Однако практике отвечают и более сложные структуры рядов:
мультипликативный и аддитивно-мультипликативный (смешан ный) характеры взаимодействия компонент, особенно если речь идет о рядах стоимостных рядов показателей (широко распростра ненных в приложениях) [66].
Приведенное выражение для мультипликативной структуры отражает взаимодействие всех компонент ряда динамики и обла дает в силу этого большей общностью, чем предположение о неза висимости в аддитивной структуре.
Однако каноническая мультипликативная структура распро странена в прикладной эконометрике существенно меньше. Глав ная трудность ее идентификации состоит в том, что речь идет о влиянии стохастической на колебательную компоненту, а затем о влиянии колебательной компоненты на тренд, причем о влияниях в пропорциях (в долях).
Моделирование влияния в пропорциях не получило известно го простого аналитического решения. Кроме того, на колебатель ную и стохастическую компоненты нужно наложить условие, чтобы они не меняли знаки уровней тренда. Стохастическая компонента должна быть неотрицательной и иметь в силу этого ненулевое и положительное математическое ожидание [1], что заставляет вве сти в модель дополнительный постоянный параметр. Параметри ческое (аналитическое выражение) в канонической мультиплика тивной структуре может иметь только модель тренда, что затруд няет реализацию прогнозирования уровней ряда, которое обычно и является конечной целью моделирования.
На коротком горизонте прогноза, используемом обычно для мониторинга эволюции, прогноз в большей мере связывают имен но с колебательной компонентой.
Поставим задачу рассмотрения и других структур взаимодей ствия компонент ряда динамики.
Многим приложениям, отраженным, например, в [5], соответ ствует более сложная пропорционально-мультипликативная форма взаимодействия компонент, для которой можно предложить, на пример, следующую структуру [66]:
Yt Tt 1 St Ct t.
В этом случае новая колебательная компонента Tt St оказыва ется пропорциональной тренду.
Аддитивная часть рассматриваемой пропорционально мультипликативной структуры иллюстрирует независимость ком понент, а мультипликативная, наоборот, демонстрирует наличие пропорциональной связи между колебательной компонентой и трендом, которая часто встречается в экономической практике [5].
Достоинством последней структуры ряда является то, что она допускает аналитическое моделирование всех компонент и в силу этого может широко использоваться на практике.
Опираясь на идею пропорционально-мультипликативного взаимодействия компонент, можно предложить еще ряд структур временных рядов [66], причем применительно к анализу кумуля тивных логистических они в известной отечественной и зарубеж ной эконометрической литературе не обнаружены.
Известны два общих подхода к моделированию и прогнозиро ванию рядов динамики социально-экономических показателей – параметрический и непараметрический.
Непараметрический (или алгоритмический) подход не позво ляет получить аналитического выражения, которое описывало бы исходный временной ряд. Достоинством такого подхода является, как правило, его простота, что обеспечивается, однако, требовани ем больших объемов исходных выборок [41]. При непараметриче ском подходе прогнозирование и мониторинг эволюции моделей крайне затруднены.
При параметрическом (или аналитическом) подходе к моде лированию ряда динамики строится аналитическая модель вре менного ряда и идентифицируются ее параметры. В этом случае результатом является построение зависимости вида Yto f, t, где Yto – модельные значения ряда динамики, – вектор па раметров модели.
Реализация параметрического подхода, как правило, сложнее.
Она требует выбора (обоснования) класса адекватных (соответст вующих цели моделирования) моделей, а также выбора (разработ ки) метода идентификации, обеспечивающего высокую точность моделирования и прогнозирования.
В монографии, несмотря на большие трудности реализации, выбран параметрический подход, так как он, в дополнение к ука занным достоинствам, позволяет осуществить и декомпозицию ря да динамики: выделить логистический тренд и, что часто имеет ме сто на практике, дополнительный тренд и колебательные компо ненты.
При использовании параметрических моделей соответствую щим выбором метода идентификации можно на относительно ма лых выборках обеспечить и прогнозирование рядов с логистиче ским кумулятивным трендом.
1.2. Формирование логистических моделей трендов как решений дифференциальных уравнений роста СЭС Логистические модели трендов чаще всего вводятся как ана литические (или численные) решения соответствующих диффе ренциальных уравнений для показателей моделей роста СЭС с тем или иным видом его ограничения. Число таких логистических трендов, описанных в литературе, можно оценить в несколько де сятков [102].
Особенностью некоторых из предложенных зарубежными учеными моделей является наличие произвольной асимметрии по отношению к точке перегиба, в то время как асимметрия двух наи более часто рассматриваемых в отечественной литературе кривых Ферхюльста и Гомпертца фиксирована, определяется лишь через уровень насыщения функции.
В монографии будут рассматриваться три критические точки логистической кривой, две из которых в известной литературе не определяются.
Первая критическая точка – упомянутая точка перегиба t, ко торая разделяет кривую на две части, не обязательно равные по длительности.
На рисунке 1.6 приведены симметричная и асимметричная ло гистические кривые. Для уровня насыщения принято обозначение A0, а цифрами обозначены:
1 – симметричная логиста, в которой точка перегиба делит кривую на две равные части, а ее ордината соответствует половине уровня насыщения;
2 – асимметричная логиста, в которой точка перегиба смещена влево, ее ордината составляет менее половины уровня насыщения, первая часть кривой короче второй;
3 – асимметричная логиста, в которой точка перегиба смещена вправо, ее ордината составляет более половины уровня насыще ния, а первая часть кривой длиннее второй.
Уровень насыще- A T (t ) ния T3 t A0 / T1 t T2 t 0 t t Рис. 1.6. Точка перегиба логистической кривой Отметим, что логистическая динамика является эволюционирую щей, поскольку представляет собой смену трех стадий роста пока зателя во времени (рис. 1.7).
T (t ) I II III 0 t t II t III Рис. 1.7. Три стадии логистического роста Точки t II и t III, в которых третья производная логистической функции равна нулю, обозначают начало второй (линейного роста) и третьей (замедляющейся, гиперболической) стадий логистиче ского роста соответственно.
На рисунке 1.8 схематично изображены симметричная отно сительно точки перегиба логиста и три ее производные, а также отмечены точки t, t II и t III.
T ''(t ) T (t ) t t III t II t t II t t III t T '''(t ) T '(t ) t t II t t III t II t t III t Рис. 1.8. Логистическая функция и три ее производные Для многих приложений целесообразно подробнее рассмат ривать логистические модели трендов с фиксированной и произ вольной асимметрией и приводить, где это возможно, аналитиче ский расчет всех трех критических точек.
1.2.1. Логистические модели трендов с фиксированной асимметрией Симметричная трехпараметрическая логистическая кривая Ферхюльста [123, 137] задается дифференциальным уравнением (1.1):
dT T 1 T / A0, T (0) A0 / (1 A1 ), (1.1) dt решением которого является функция (1.2):
A T (t ). (1.2) 1 A1et Здесь и далее параметр отражает скорость роста функции.
Параметр A0 является верхней горизонтальной асимптотой и отве чает за уровень насыщения логистической кривой.
В частности, для (1.3):
A A0. (1.3) lim t 1 A e t Некоторые модификации модели Ферхюльста, определяю щиеся, однако, тем же дифференциальным уравнением (1.1), были предложены Менсфилдом (1961), Фишером и Праем (1971) [103].
Например, функцию Фишера и Прая обычно записывают сле дующим образом:
T A et C F1 T t C, или T (t ) 0 t C.
t C e, F1, или ln A0 T 1 e 1 F1 A Дифференциальное уравнение при этом имеет вид (1.1), а точка перегиба определяется из условия равенства нулю второй производной: t C /, T t * A0 / 2. Точка перегиба «классической»
логисты (1.2) также соответствует половине уровня насыщения (1.4):
t* ln 1/ A1 /, T t * A0 / 2. (1.4) Точки t II, t III находятся численным решением следующего уравнения: e2 t A12 4 A1e t 1 0.
Далее будем придерживаться классического определения ло гисты Ферхюльста, без переобозначения параметров.
Вид функции Ферхюльста при изменении параметров пред ставлен на рисунке 1.9. Здесь и далее стрелками показано направ ление увеличения параметров.
6 6 T(t) T(t) A A0 4 A0 0, 4 A A0 0, A0 0, 3 A1 0, 3 A1 0, A0 0, A1 0, 2 A1 0, t t 0 2 4 6 8 10 12 14 16 18 8 10 12 14 16 18 02 4 а) б) 6 T(t) A0 0, A1 0, A0 0, A1 0, t 0 2 4 6 8 10 12 14 16 18 в) Рис. 1.9. Вид логисты Ферхюльста при изменении знаков и величин параметров В настоящее время логистической кривой Ферхюльста в ос новном описывают динамику роста или уменьшения численности популяций в биологии, реже – производство электроэнергии или добычу энергоресурсов, смену технологий и др. [101, 103, 104].
Дифференциальное уравнение для функции Гомпертца [98] имеет вид:
dT T ln A0 / T, T (0) A0e e.
t dt Аналитическим решением уравнения является функция (1.5):
t t T (t ) A0e e. (1.5) Запись функции (1.5) включает в качестве параметра абсциссу точки перегиба t0 (1.6):
t* t0, T (t*) A0 / e. (1.6) Точки t II и t III находятся численным решением уравнения:
e 3e2 e3 0.
t t t t t t 0 0 В отечественной эконометрической литературе логисту Гом пертца традиционно записывают в следующем виде (1.7) [1, 83]:
T (t ) KAB.
t (1.7) В настоящее время логистой Гомпертца пользуются при моде лировании динамики роста опухолей, числа абонентов сотовой те лефонии, численности населения, потребительских товаров дли тельного пользования, инноваций в сельском хозяйстве и др. [101].
Следует отметить, что поскольку функция Гомпертца асим метрична по отношению к половине уровня насыщения A0 / 2, то ее точка перегиба находится левее по отношению к точке перегиба симметричной функции Ферхюльста.
На рисунке 1. показаны функции Гомпертца и Ферхюльста при совпадающих значениях абсциссы точки перегиба t, скорости рос та и уровня насыщения A0.
100 T(t) t* 80 Кривая Ферхюльста 60 A A e Кривая Гомпертца t 0 3.8 7.6 11.4 15.2 19.0 22. Рис. 1.10. Сравнение функций Ферхюльста и Гомпертца при A0 100;
0,5;
t0 ln 1 / A1 11,4;
A1 Предложим расширение модели Гомпертца (1.8): преобразуем модель (1.5) так, чтобы получить кривую с правой асимметрией (точка перегиба будет правее точки половины уровня насыщения):
.
t t0 t T (t ) A0 1 e e, T (0) A0 1 e e (1.8) Точка перегиба определится выражениями: T (t*) A0 (1 1/ e), t* t0. На рисунке 1. представлены асимметричные (влево и вправо относительно точки перегиба) кривые Гомпертца.
В известной литературе не обнаружены модели с фиксирован ной асимметрией, смещенной вправо относительно точки перегиба.
20 T(t) t левая асимметрия A0/ правая асимметрия t 0 5 10 15 Рис. 1.11. Сравнение функций Гомпертца с левой и правой асимметрией при A0 20;
0,5;
t0 Решением дифференциального уравнения (1.9) dT 3T 1 T / A0 3, T (0) A0 1 e t (1.9) dt при A1 0 является четырехпараметрическая логистическая функ ция фон Берталанффи с фиксированной асимметрией, смещенной влево относительно точки перегиба [88]:
, T (t*) 27 A.
t t T (t ) A0 1 e Точки t II и t III находятся численным решением уравнения:
e 8e2 9e3 0.
t t t t t t 0 0 Существуют и другие дифференциальные уравнения, опреде ляющие логистические функции, но не всегда имеющие аналитиче ское решение.
Для точки перегиба в этом классе функций аналитически можно определить только ординату. На практике непрерывный случай заменяется дискретным по схеме Эйлера Tk H (Tk 1, ), где H – нелинейная в общем случае функция, – вектор параметров мо дели.
В такой постановке очевидно, что параметром является и на чальное значение T0, на основе которого осуществляется пересчет дискретного ряда.
Наибольший интерес при этом представляют модели, точка перегиба которых определяется несколькими параметрами, а не только одним, отвечающим за уровень насыщения.
Правая асимметрия достигается в этих случаях за счет исполь зования функций с произвольной асимметрией.
1.2.2. Логистические модели трендов с произвольной асимметрией Рассмотрим более сложные четырех - и пятипараметрические логисты, точки перегиба которых могут располагаться достаточно произвольно по отношению к половине уровня насыщения, по скольку определяются двумя или даже тремя параметрами модели.
Примером четырехпараметрической логисты может служить модель Басса (1969) [83]:
A dT T C 1 T / A0, T (0) 0. (1.10) 1 C dt Уравнение (1.10) имеет аналитическое решение, а с переобо значением параметров модель Басса можно записать и в виде (1.11):
A0 e ( pq ) t, p C, q.
T (t ) (1.11) 1 ( q / p )e ( p q ) t Модель Басса используется в маркетинге и отражает процессы того, как новый товар принимается покупателями, как происходит его продвижение на рынке в результате взаимодействия тех по требителей, которые уже пользуются товаром, и его потенциаль ных потребителей.
При этом коэффициент p носит название коэффициента инно вации, коэффициента внешнего влияния или рекламного эффекта.
Коэффициент q – это коэффициент имитации, внутреннего влия ния или так называемого «маркетинга из уст в уста» или «сарафан ного радио» [113].
Точка перегиба рассчитывается следующим образом:
ln p q A p t, T (t ) 0.
pq 2 2q При p C 0 получим как частный случай симметричную ло гистическую кривую Ферхюльста.
Точки t II и t III находятся численным решением уравнения:
q2 e2t 4 pqe t p 2 0.
p q p q Динамику смены технологий Шариф и Кабир в 1976 г. пред ложили описать четырехпараметрической моделью (1.12), добавив в модель дополнительный параметр [132, 131]:
A0 3 1 dT T 1 T / A, T (0) T0, T (t*). (1.12) 4 dt 1 T 1 / A Ординаты точек t II и t III находятся численным решением уравнения:
6( 1)2 T 4 18 A0 ( 1)T 3 A0 (19 4 )T 2 2 A0 (4 )T A0 0.
2 3 Схожий вид имеет четырехпараметрическая логиста Скиадаса, Джиованиса и Димотикалиса (1985) [97]:
A0 dT T 1 T / A, T (0) T0, T (t*). (1.13) dt 1 T 1 / A0 Ординаты точек t II и t III находятся численным решением уравнения:
( 1)2 T 4 4 A0 ( 1)T 3 6 A0 T 2 2 A0 (2 )T A0 0.
2 3 При 1 модель Скиадаса становится моделью Ферхюльста.
На рисунке 1. представлено сравнение кривых Скиадаса и Шарифа при одинаковых значениях параметров.
50 T(t) Кривая Скиадаса Кривая Шарифа t 0 10 20 30 40 Рис. 1.12. Сравнение кривых Скиадаса и Шарифа при A0 50;
0,9;
6;
T0 0. Представляется возможным использование и пятипараметри ческой логисты более общего вида:
dT T 1 T / A, 1, T (0) T0.
dt 1 T 1 / A Другим пятипараметрическим усложнением модели (1.13) яв ляется дифференциальное уравнение, содержащее параметр g.
Точка перегиба в этом случае определяется уже тремя пара метрами модели:
dT 1 T / A0 T / A0 T g, T (0) T, A02 2 A0 g g 2 A0 g.
T (t*) dt 1 T / A0 T g / A Ф. Ричардс в 1959 г. предложил обобщенную логистическую кривую, добавив в дифференциальное уравнение Ферхюльста (1.1) дополнительный параметр и получив уравнение (1.14) [129]:
, dT TM 1 T / A M t T (0) A0 1 e. (1.14) M dt При этом одной из наиболее употребимых записей решения дифференциального уравнения (1.14) является пятипараметриче ская функция (1.15) с добавлением константы С (для того чтобы обеспечить смещение функции вдоль оси ординат):
t t M T (t ) C A0 1 e, (1.15) которая имеет логистический характер при M 0, A1 0.
Функция Ричардса является действительно общим случаем для нескольких известных логист: при M 1 получаем функцию Ферхюльста T (t*) A0 / 2 ;
при M – функцию Гомпертца T (t*) A0 / e.
Также одним из частных случаев функции Ричардса при M 3, A1 0 является функция фон Берталанффи. Для произвольно го M 0 точка перегиба функции Ричардса определится соотноше ниями:
1 e t M, T (t*) A0 1.
t* ln A1M M Точки t II и t III находятся численным решением уравнения:
3M 1 e tt M 2e2 t t 1 0.
0 На рисунке 1. представлена логиста Ричардса при различных значениях параметра M.
100 T(t) 40 M=–0, M=– M=– t 14. 0 5 9.9 19.9 24. Рис. 1.13. Функция Ричардса при различных значениях параметра М, A0 100;
0,5;
C 0;
t0 5,10, Приведем примеры и других пятипараметрических логист, положение точек перегиба которых зависит уже от четырех пара метров.
Например, логиста Джюланда (1981) [103] определяется диф ференциальным уравнением (1.16):
dT T C 1 T / A0, T (0) T0.
(1.16) dt При 1 получаем логисту Басса, запись которой в общем случае имеет вид T (t*) A0 1 C / / (1 ).
Функция Изингвуда (1981) была предложена для описания динамики розничных продаж товаров, а также инноваций в сфере образования [103]. Она также включает пять параметров (1.17):
T C 1 T / A0, T (0) T0, dT (1.17) dt однако аналитически определить ординату точки перегиба невоз можно. Ордината точки перегиба T * T (t*) есть решение следую щего уравнения:
1 T * A0 T * p 0.
Перейдем к рассмотрению феноменологических моделей логи стического тренда, под которыми будем понимать эмпирически ус танавливаемые закономерности между малым числом входных и выходных переменных по выборкам малого объема [72].
1.3. Феноменологические логистические модели трендов Некоторые из рассматриваемых в этом параграфе феномено логических моделей получены в результате действий сложения и умножения над уже рассмотренными аналитическими решениями дифференциальных уравнений, другие же предложены эмпириче ски [66].
Так, новыми свойствами по отношению к рассмотренным вы ше моделям обладает асимметричная двухпараметрическая логи стическая кривая Рамсея (1.18), ранее использовавшаяся для обес печения робастности оценок регрессии [38, 128]:
T (t ) B0 (1 (1 t )et ). (1.18) Данная логиста отличается тем, что у нее практически отсут ствует стадия медленного (экспоненциального) роста, поскольку функция в этом случае выходит из начала координат.
Логиста Рамсея (1.18) имеет уровень насыщения, равный B0, и точку перегиба с координатами:
t* 1/, T t * B0 1 2 / e, т.е. обладает левой фиксированной асимметрией: 0,264B0 B0 / 2.
Принципиально новым положительным свойством кривой Рамсея является отсутствие нелинейных операций над экспонен циальной функцией: деления (как в функции Ферхюльста), возве дения в степень (как в функции Ричардса) или взятия логарифма (как в функции Гомпертца).
Это позволяет упростить ее идентификацию, поскольку воз можно применение аналитического метода конструирования обобщенных параметрических ARMA-моделей при помощи Z преобразования [61], без использования численных методов, а так же допускается возможность введения в модель с логистическим трендом Рамсея других трендовых и колебательных компонент для идентификации.
Аналогичная функция приведена в монографии [33], где она рассматривалась как модификация модели Монти-Кляйна [107, 119] для регулирования депозитной ставки банка в функции затрат на привлечение средств. Таким образом, модель Рамсея можно счи тать и модификацией модели Монти-Кляйна, а предложенные ниже модели – расширениями модели Рамсея (или Монти-Кляйна).
Целесообразно предложить и обобщение модели (1.19): трех параметрическую модель Рамсея в виде суммы функции (1.18) и параметра C, который введен для того, чтобы в начальный момент времени значение логистической функции было ненулевым:
T (t ) C B0 (1 (1 t )et ). (1.19) Такое требование бывает необходимым, например, при моде лировании жизненных циклов продуктов (товаров, услуг), когда в начальный момент запуска товара на рынок уровень продаж отли чен от нуля.
Уровень насыщения логисты при этом изменится и станет равным сумме B0 C, а координаты точки перегиба логистической кривой (1.19) будут следующими:
t 1/, T (t ) B0 1 2 / e C 0,264 B0 C.
Модель (1.19) асимметрична, так как в точке перегиба орди ната в общем случае не равна половине уровня насыщения:
0,264B0 C 0,5 B0 C.
Поскольку в данной модели тренда практически отсутствует стадия экспоненциального роста, то для нее и ее модификаций, предложенных ниже, будет отсутствовать точка t II. Вид первых двух производных для феноменологической модели Рамсея в срав нении с двумя производными модели Ферхюльста представлен на рисунке 1.. Производные модели Рамсея показаны жирными ли ниями.
T '(t ) T ''(t ) t t 1 t t II t III t t 0 5 10 15 20 0 5 10 15 t III - Рис. 1.14. Вид двух производных для модели Рамсея в сравнении с двумя производными модели Ферхюльста Точка t III для кривой Рамсея может быть вычислена аналити чески:
t III 2 /, T (t III ) B0 1 3 / e2 C.
Вторым предложением по расширению функциональных воз можностей классической модели Рамсея является еще одна, уже трехпараметрическая, модифицированная модель (1.1.20) с уров нем насыщения B0 C :
T (t ) C ( B0 t )et. (1.20) На рисунке 1.5 показано изменение модифицированной функ ции Рамсея (1.20) при вариации параметра B0.
B0 = 0 T(t) C 0, 0 B0 = - B0 = - B B0 = -3 t - 10 -5 0 5 B0 = -4 - Рис. 1.5. Графики модифицированной функции Рамсея при изменении величины B Для (1.20) можно найти точку перегиба и точку t III. Логисти ческая кривая при этом обладает произвольной асимметрией:
3 B 2 B, T t C, t III, T (t III ) C t.
exp(3 B0 ) exp(2 B0 ) Заметим, что при условии B0 2 точка перегиба модифици рованной логисты Рамсея будет находиться в отрицательной об ласти, а сама кривая приобретает характер кривой роста, а не куму лятивной логистической модели.
При условии B0 0 функция имеет характер импульсной мо дели. При B0 2 ордината точки перегиба равна нулю.
Существует еще один вариант возможного расширения моде ли (1.18): обобщенная логиста Рамсея (1.21) с уровнем насыщения B0 C, в которой количество параметров равно четырем, а функ циональные возможности существенно шире:
T (t ) C ( B0 B1t )et. (1.21) В этом случае изменяется характер кумулятивной логистиче ской функции (вместо стадии медленного роста на начальном эта пе происходит кратковременное уменьшение динамики), а коорди ната точки перегиба определяется уже тремя параметрами модели:
B0 B 2 2 B1 B0 3B1 B B, T t C B, T t III C 2B e 3B e t, t III.
B1 B1 Отметим, что дальнейшее увеличение порядка полинома ( B B1t ) в (1.21) до второго ( B0 B1t B2t 2 ) и выше может существен но изменить форму функции.
Таким образом, обобщенная модель Рамсея соответствующи ми изменениями параметров C, B0, B1 может перемещаться и по оси абсцисс, и по оси ординат, что демонстрирует ее гибкость и позво ляет рассчитывать на ее широкое применение.
Известно применение и модели Чантера [61] (1.22):
A T (t ) (1.22), A2 t 1 A1 exp 1 e которая является лишь «перепараметризацией» модели Ферхюль ста, поскольку симметрична.
В отдельных приложениях использовалась еще одна симмет ричная модель – кусочная функция Андрюса [61], определяемая со отношением (1.23):
t t C C 1 cos, T (t ) C (1.23) 2C, t C.
Параметры логисты Андрюса таковы: уровень насыщения ра C, T (t ) C.
вен 2C, координаты точки перегиба t Известна также асимметричная функция Моргана, Мерсера и Флодина (1.24) [120], которая имеет S-образный вид в первой ко ординатной четверти при 1:
t T (t ). (1.24) ( t ) Параметр 0 определяет уровень насыщения, а параметр 0 смещает функцию вдоль оси абсцисс.
Точка перегиба при заданных ограничениях определяется со отношением:
1 1 t, T (t ).
1 Точки t II и t III находятся численным решением уравнения:
2 3 2 t 2 4 1 2 t 2 2 3 2 0.
Известная функция распределения Вейбулла (1.25) [138] так же имеет логистическую форму:
T (t ) exp( t ), (1.25) где параметр 0 также определяет уровень насыщения, а 1.
Последнее требование необходимо, чтобы функция имела ло гистический характер, поскольку при других значениях параметров ее вид может быть отличен от S-образного.
Точка перегиба определится следующими соотношениями:
1 /, T (t ) e t.
Точки t II и t III находятся численным решением уравнения:
2 2t 2 3 1 t 2 3 2 0.
Функции Моргана и Вейбулла нелинейны, и, помимо этого, вид данных функций сильно зависит от значений параметров (как уже было отмечено, кривая необязательно будет иметь S-образную форму), что сильно ограничивает возможности их идентификации.
Тем не менее известно их применение при оценке качества жизни населения, в частности при моделировании динамики пока зателей заболеваемости [84, 96, 108, 111, 141].
К феноменологическим моделям относится и полиномиальная модель Гомпертца (рис. 1.16) с числом параметров m + 2, имеющая t m вид T (t ) A0 exp bi u i du.
0 i 0 20 T(t) t 0 4 8 12 16 Рис. 1.6. Вид полиномиальной функции Гомпертца третьего порядка при различных сочетаниях параметров Данной моделью, до 5 порядка (m = 5) включительно, модели руют, например, средний рост и вес людей, в частности уже упоми навшийся в примерах средний рост мальчиков [125].
Перейдем теперь к рассмотрению моделей би - и мультилоги стической динамики.
Рассмотрим, например, модель в виде суммы двух логистиче ских кривых Ферхюльста, что дает известную билогистическую ад дитивную шестипараметрическую модель (1.26):
A01 A T (t ). (1.26) 1 A11e1t 1 A12 e 2t Другую модель можно предложить путем перемножения логи стических кривых Ферхюльста.
Ее можно, по аналогии с предыдущей моделью, назвать било гистической мультипликативной шестипараметрической моделью (1.27):
A01 A T (t ). (1.27) 1 A11e1t 1 A12 e2t Сравним графики исходных логистических кривых Ферхюль ста, билогистических аддитивной и мультипликативной кривых при одинаковых значениях параметров (рис. 1.7).
T(t) T(t) мультипликатив- мультипликативная ная модель модель аддитивная аддитивная модель модель t t а) б) T(t) мультипликативная модель аддитивная модель t в) Рис. 1.7. Случай взаимодействия растущей и падающей (а), двух растущих (б) и двух падающих (в) логистических кривых Ферхюльста Две точки перегиба образуются от взаимодействия двух логи стических кривых, а абсциссы остаются прежними:
1 1 t1, t2 ln.
ln 2 A 1 A11 Ординаты изменяются и не равны половине уровня насыще ния каждой из логистических кривых.
В случае билогистической аддитивной модели будем иметь:
A02 A A A T (t1 ) 2, T (t2 ) 02.
2 1 A12 A11 A11 A Для билогистической мультипликативной модели получим:
A01 A02 A01 A T (t1 ), T (t2 ).
2 1 2 A12 A11 1 2 A11 A12 В качестве суммируемых или перемножаемых логистических функций могут выступать и другие рассмотренные выше логисти ческие кривые.
Например, суммированием функции Ферхюльста и функции Гомпертца можно добиться того, что в билогистической динамике одна часть будет асимметричной, а другая – симметричной.
Импульсная кривая, как показано на рисунке 1.7 а, формиру ется двумя различными логистическими кривыми.
Импульсную динамику можно описывать и одной кривой (это известные функции Хабберта, Хаммонда-Маккея, Капицы (Коши), логнормального распределения, некоторые дробно-рациональные функции [31, 34, 57, 91, 90, 105]).
Например, одной из наиболее известных и давно используе мых импульсных моделей для моделирования жизненного цикла добычи нефти является симметричная колоколообразная функция Хабберта (1.28) [105]:
2Ymax T (t ), (1.28) 1 ch( (t t1 )) где t1 и Ymax – абсцисса и ордината вершины импульсной кри вой (пика), а параметр регулирует наклон кривой этапов роста и падения относительно ее пика.
Асимметрия кривой жизненного цикла характерна для боль шинства выборок в экономической практике. Задание этой асим метрии также является одной из областей применения логистиче ских кривых: введением в них для управления наклонов этапов роста и спада кривой жизненного цикла логистической функции [90].
Например, можно использовать модели Гомпертца для реали зации такой асимметрии путем изменения параметра по логи стическому закону: от уровня С на этапе роста до уровня C A0 на этапе спада. Скорость этого перехода будет задаваться параметром 1/.
На рисунке 1.8 представлены модели Хабберта с различными значениями параметра :
1. (t ) 1.
t t 2. (t ) C A0e e – функция Гомпертца с левой асимметрией.
t t 3. (t ) C A0 1 e e – функция Гомпертца с правой асим метрией.
4.
T(t) t Рис. 1.8. Описание импульсной динамики при помощи логистических кривых Если цель моделирования иная или если в выборке анализи руемого показателя присутствует только часть данных (только этап роста до пика либо только этап спада), то имеет смысл рас сматривать отдельно две независимых логистических кривых.
Более сложная динамика может формироваться и более чем из двух логистических кривых.
1.4. Актуальные задачи развития инструментария идентификации динамики с логистическим трендом Уникальным свойством логистической кривой является ее способность описывать качественные (структурные) изменения в развитии динамики, т.е. моделировать эволюцию социально экономических показателей с помощью одной модели тренда. Та кое свойство логисты определяется наличием точки перегиба функции, которая означает смену знака ее второй производной.
Предположение о едином тренде для эволюционирующей ди намики означает возможность прогнозировать дальнейшее разви тие, имея лишь некоторый начальный этап траектории, не доходя до уровня насыщения.
При планировании будущих продаж товара, оценки эффек тивности вложенных инвестиций, необходимости дальнейшего сервисного обслуживания проданных товаров и т.п. возникает не обходимость получить оценку емкости рынка (потенциала отрас ли), сделав прогноз уровня насыщения по логистической кривой.
Его оценка важна для производителей товаров (услуг, техно логий), поскольку завышенный прогноз приведет к увеличению за трат на исследования и разработку новых товаров, а заниженный – препятствует выведению новых товаров на рынок.
В большинстве исследований, однако, при рассмотрении логи стических кривых ограничиваются экспертно задаваемыми оцен ками уровней насыщения, что связано в первую очередь с трудно стью идентификации моделей, содержащих логистический тренд.
Возникает задача рассмотрения методов идентификации мо делей логистических трендов, позволяющих рассчитать оценку уровня насыщения при моделировании, определить координаты точек перегиба.
В известной литературе не приводится также расчет двух критических точек логистической кривой, определяющих начало стадий линейного роста и угасания. Тем не менее подобный анализ может быть интересен исследователю, поскольку смена скорости роста в динамике любого социально-экономического показателя всегда связана с какими-либо событиями внутренней и/или внеш ней среды.
Особое внимание следует уделить логистическим функциям с произвольным положением точки перегиба, которая не соответст вует середине интервала по оси ординат между началом логисти ческого процесса и уровнем насыщения.
Сложность идентификации таких моделей состоит в том, что они, как правило, содержат не менее четырех-пяти параметров, входящих в модель нелинейно. Существующие методы линеариза ции, которые будут подробно описаны во второй главе, касаются только простейших моделей Ферхюльста и Гомпертца с фиксиро ванной асимметрией и не могут быть применимы к остальным мо делям.
Помимо этого, подобные методы игнорируют присутствие в эконометрической модели стохастической компоненты.
Для получения достоверных прогнозов уровней насыщения с помощью логистических кривых до сих пор не решена задача опре деления требуемых объемов выборки: достаточно ли для прогно зирования лишь начального этапа динамики, либо показатель в своем развитии должен «пройти» точку перегиба и т.д.
Не изучены и сферы применения симметричных и асиммет ричных логистических кривых. В известной литературе представ лены модели, которые разрабатывались авторами под каждый конкретный случай, без рекомендаций по дальнейшему примене нию. Особенно это характерно для более новых дифференциальных и феноменологических моделей. При этом модели, разработанные еще в XIX веке, применяются повсеместно, поскольку содержат меньшее число параметров и могут считаться более «простыми».
Однако разработанные методы идентификации таких моде лей не были протестированы в различном диапазоне изменения сочетаний параметров и дисперсии помехи. Во всех известных при мерах дисперсия стохастической компоненты (шума) составляет, как правило, не более 5% от дисперсии детерминированной ком поненты [84, 97, 108, 104, 125, 92, 121].
В то же время известной рекомендацией в теории измерений является требование обеспечивать удовлетворительную точность применяемого метода (в нашем случае идентификации модели) в диапазоне соотношений дисперсий помехи/детерминированной компоненты до 20-30%.
В примерах, приведенных выше, показано, что ряды динамики с логистическим трендом могут содержать и колебательную ком поненту, и дополнительный тренд.

Content-Type: text/html; charset=WINDOWS-1251


markdown:
Trang chủ
Nhà đất
Huyện Ân Thi
2
1
Về MuaBanNhaDat | Nhận bản tin | FAQ | Điều khoản thỏa thuận | Bản quyền | Quy chế Sàn GDTMĐT | Tuyên bố từ chối bảo đảm | Liên hệ | Tuyển dụng | Bảng giá dịch vụ | Sitemap
MuaBanNhaDat © 2007 - 2016. Ghi rõ nguồn "MuaBanNhaDat.vn" khi phát hành lại thông tin từ website này.
Giấy phép số: 40/GP-TTDT, do Sở Thông tin và Truyền Thông TP. Hồ Chí Minh cấp ngày 22/04/2015
Giấy Chứng Nhận Đầu Tư Số: 411022000534; Ngày cấp: 20/01/2014 do Sở Kế Hoạch và Đầu Tư Tp.HCM cấp
Địa chỉ: Lầu 3 Tòa nhà Copac, Số 12, Tôn Đản, Phường 13, Quận 4, TP.HCM
Follow us on Facebook Linkedin Google+
Technology powered by - Copyright ©2016 by Ringier Vietnam

WARC-Type: response

Content-Type: application/http; msgtype=response

Cache-Control: private

Content-Type: text/html; charset=utf-8

Server: MBND's Server

X-Frame-Options: SAMEORIGIN

Connection: close


markdown:
Tác giả
KOMEYAMA Shiwo
Tình trạng
Đang tiến hành
Thể loại
Adult - Romance - Manga
Bán thiên thần là loài vật mang đầy tội lỗi từ khi mới sinh ra, họ phải chịu lời nguyền chết chóc vào lúc 18 tuổi. Muốn giải lời người, những bán thiên thần phải mạo hiểm tìm quỷ máu thuần và cạn sạch máu của nó. Tsukimiya là thiên thần lai ma cà rồng, trong lúc chiến đấu tìm quỷ thuần, cô gặp Hinata và quyết định hợp tác với anh để dễ dành xóa bỏ ấn nguyền trên ngực mình. Nhưng Hinata là con người không đáng tin cậy, vậy thực sự anh xem Tsukimiya là gì?


markdown:
Tác giả
đang cập nhập
Tình trạng
Đang tiến hành
Thể loại
Slice of Life - Horror
Truyện ngắn kinh dị, bựa, hài...


markdown:
Thêm bài hát này vào danh sách Playlist
Lời đăng bởi: meegovn
Uhhh ah...
Kỉ niệm còn đâu đây nhớ nơi sân trường ấy
Uhhh ah...
Kỉ niệm học sinh khi vui đùa trong gió chiều
Nhìn theo tà áo dài tan học về trên phố tung tăng
Còn đâu, còn đâu, oh oh...
Thời gian dần trôi cùng với tuổi thơ
Ngày xưa hồn nhiên cùng bao ước mơ
Và khi tôi cách xa tháng năm học trò
Còn nhớ ngày chia tay ấy ... Buồn!
Ngoài sân phượng rơi hồng trang kỉ niệm
Chợt nhớ về những ngày xưa đã xa
Và hôm nay chúng ta đã quay trở về
Bạn ơi, giờ đây tôi muốn nói
Rằng lòng tôi vẫn nhớ
Người thầy ngày xưa dạy chúng tôi không còn
Trở lại nơi chốn ấy
Từng năm tháng chợt hiện trong ký ức
Ngày xưa đó mãi mãi bay xa, xa rời ta về nơi phương trời
Cùng bao cành hoa phượng rơi ngập đầy nắng vàng
Và những chiều mưa, ta cùng bước trên đường xưa
Ngày xưa đó mãi mãi bay xa, xa rời ta về nơi chốn ấy
Bao bạn thân khi xưa bên nhau
Giờ đây còn lại mình ta
Trong nỗi nhớ thiết tha
Uhhh ah...
Kỉ niệm còn đâu đây nhớ nơi sân trường ấy
Uhhh ah...
Kỉ niệm học sinh khi vui đùa trong gió chiều
Nhìn theo tà áo dài tan học về trên phố tung tăng
Còn đâu, còn đâu, oh oh...
Bao nhiêu yêu thương giấc mơ
Đã đi xa về nơi trái tim
Và còn lại trong ta những tháng năm dài
Nhớ về những thơ ấu kỉ niệm đã xa từ lâu
Giờ còn nụ hồng khi xưa ta tặng bạn
Ép trong trang vở lưu bút ấy
Thế gian có đổi thay, tháng năm có phai
Ta sẽ mãi nhớ đến, nhớ đến nhau
Rằng lòng tôi vẫn nhớ
Người thầy ngày xưa dạy chúng tôi không còn
Trở lại nơi chốn ấy
Từng năm tháng chợt hiện trong ký ức
Ngày xưa đó mãi mãi bay xa, xa rời ta về nơi phương trời
Cùng bao cành hoa phượng rơi ngập đầy nắng vàng
Và những chiều mưa, ta cùng bước trên đường xưa
Ngày xưa đó mãi mãi bay xa, xa rời ta về nơi chốn ấy
Bao bạn thân khi xưa bên nhau
Giờ đây còn lại mình ta
Trong nỗi nhớ ấy
Đâu? Đâu? Biết đâu rồi?
Ngày bạn và tôi vẫn thường đến lớp
Giờ tìm nơi đâu?
Còn lại đây chút nắng hồng ngây ngô bằng lăng tím
Vẫn nhớ hoài những ngày thơ ấu
Ta vui chơi dưới sân trường
Nhưng giờ đây xa nhau mất rồi
Bạn thân ơi...Chớ quên
Ngày xưa đó mãi mãi bay xa, xa rời ta về nơi phương trời
Cùng bao cành hoa phượng rơi ngập đầy nắng vàng
Và những chiều mưa, ta cùng bước trên đường xưa
Ngày xưa đó mãi mãi bay xa, xa rời ta về nơi chốn ấy
Bao bạn thân khi xưa bên nhau
Giờ đây còn lại mình ta
Trong nỗi nhớ thiết thaaaaaaa
(Ha hà ha há ha ....)
Uhhh ah...
Kỉ niệm còn đâu đây nhớ nơi sân trường ấy
Uhhh ah...
Kỉ niệm học sinh khi vui đùa trong gió chiều
Nhìn theo tà áo dài tan học về trên phố tung tăng
Còn đâu, còn đâu, oh oh...
Thu gọn
Lynk Lee
Lâm Chấn Huy
V.A
Lynk Lee, Binz
xem thêm
rút gọn
Đơn vị chủ quản: Công ty Cổ phần NCT
Địa chỉ: Tòa nhà HAGL Safomec, 7/1 Thành Thái, P14, Q10, TP.HCM
Người chịu trách nhiệm nội dung: Ông Nhan Thế Luân - Email: support@nct.vn - Tel: (08) 3868 7979
Giấy phép MXH số 499/GP-BTTTT do Bộ Thông Tin và Truyền thông cấp ngày 28/09/2015
## Bình luận

Server: nginx+

Vary: Accept-Encoding, User-Agent

P3P: CP='NOI ADM DEV PSAi COM NAV OUR OTRo STP IND DEM'

Content-Type: text/html; charset=UTF-8


markdown:
![ lvhgt](http://www.oto-hui.com/styles/xenith/xenforo/avatars/avatar_l.png)
- Hoạt động cuối:
- 24/10/13
- Tham gia ngày:
- 28/11/11
- Số km:
- 0
- Được đổ xăng:
- 0
- Mã lực:
[0](members/lvhgt.107384/trophies)
- Xăng dự trữ:
-
[307 lít xăng](bank/transfer?to=lvhgt&rtn=ref) [$307.00](credits/transfer?receiver=lvhgt)
### Chia sẻ trang này
# lvhgt
Tài xế O-H
- lvhgt được nhìn thấy lần cuối:
- 24/10/13
-
- Đang tải...
- Đang tải...
-
### Tương tác

Cache-control: private, max-age=0

Vary: Accept-Encoding,User-Agent


markdown:
Cụ thể, theo một bản báo cáo do hãng tin BBC công bố hôm 20/12/2007, ngành công nghiệp game tại xứ sở sương mù đạt tổng doanh thu 1,52 tỷ bảng (tương đương 3 tỷ USD) trong năm 2007 - số liệu này chưa kể đến doanh thu của 2 tuần cuối cùng tháng 12. Mức tăng trưởng hàng năm đạt 25%, trong đó có trên 78 triệu bản đĩa trò chơi đã được tiêu thụ.
BBC cũng cho biết đã có hơn 11 triệu chiếc console được bán ra trong suốt quý 3/2007, tăng 33% so với cùng kỳ năm ngoái. Hơn thế nữa, số lượng đĩa game chiếm hơn 79% số lượng phần mềm được bán ra tại quốc gia này. Báo cáo của BBC cũng nêu rõ, chỉ tính riêng đĩa trò chơi, doanh thu của tuần cuối vừa qua đạt đến 87,9 triệu bảng (khoảng 174 triệu USD). Theo giới phân tích, nếu tính luôn doanh thu của 2 tuần cuối cùng tháng 12 – mùa mua sắm nhân dịp lễ Noel và cuối năm – thì doanh thu của thị trường game tại Anh sẽ còn cao hơn.
Một tuần trước đó, công ty nghiên cứu thị trường NPD cũng đã công bố báo cáo về doanh thu của thị trường game tại Mỹ trong tháng 11/2007. Theo đó, tổng doanh thu (bao gồm máy chơi game, đĩa game và phụ kiện - không xét cho hệ PC) đạt 2,63 tỷ USD.
Trong lĩnh vực máy chơi game, Wii qua mặt Xbox 360 với tổng số lượng máy bán ra lần lượt đạt 981.000 chiếc, Xbox 360 là 770.000 chiếc. Xem ra thì PlayStation 3 - 466.000 chiếc vẫn chưa đủ sức để thoát khỏi cái bóng của người anh PlayStation 2 - 496.000 chiếc. Trong khi đó, Nintendo DS vẫn tiếp tục là máy chơi game bán chạy nhất với 1,53 triệu chiếc, đối thủ Sony PlayStation Portable chỉ đạt 567.000 chiếc. Theo Nintendo, tính cuối đến cuối tháng 11 vừa qua, hãng đã bán được tổng cộng trên 6 triệu chiếc DS tại Mỹ.
Danh sách 10 game bán chạy nhất trong tháng 11/2007:
1. Call Of Duty 4: Modern Warfare (Xbox 360) - phát hành Activision - 1,57 triệu bản2. Super Mario Galaxy (Wii) - phát hành: Nintendo - 1,12 triệu bản3. Assassin's Creed (360) - phát hành: Ubisoft - 980.000 bản4. Guitar Hero III: Legends of Rock (PS2) - phát hành: Activision - 967.000 bản5. Wii Play (with remote) (Wii) - phát hành: Nintendo - 564.000 bản6. Mass Effect (360) - phát hành: Microsoft - 473.000 bản7. Call Of Duty 4: Modern Warfare (PS3) - phát hành: Activision - 444.000 bản8. Guitar Hero III: Legends of Rock (Wii) - phát hành: Activision - 426.000 bản9. Halo 3 (360) - phát hành: Microsoft - 387.000 bản10. Assassin's Creed (PS3) - phát hành: Ubisoft - 377.000 bản
Văn Tám


markdown:
Thi Viện
×
Tên tác giả
Tên bài thơ
Nội dung bài thơ
Tên nhóm bài thơ
Tên chủ đề diễn đàn
Tìm với
G
o
o
g
l
e
Chú ý:
gõ đúng chữ hoa/thường với các ký tự có dấu trong tiếng Việt!
Toggle navigation
Tác giả
Danh sách tác giả
Tác giả Việt Nam
Tác giả Trung Quốc
Tác giả Nga
Danh sách nước
Danh sách nhóm bài thơ
Thêm tác giả...
Thơ
Các chuyên mục
Tìm thơ...
Thơ Việt Nam
Cổ thi Việt Nam
Thơ Việt Nam hiện đại
Thơ Trung Quốc
Đường thi
Thơ Đường luật
Thêm bài thơ...
Tham gia
Diễn đàn
Các chủ đề mới
Các chủ đề có bài mới
Tìm bài viết...
Thơ thành viên
Danh sách nhóm
Danh sách thơ
Khác
Giới thiệu
Thống kê
Giao diện rút gọn (cho di động)
Đăng ký...
Danh sách thành viên
Từ điển Hán Việt trực tuyến
Đổi mã font tiếng Việt
Đăng nhập
×
Đăng nhập
Tên đăng nhập:
Mật khẩu:
Nhớ đăng nhập
Đăng nhập
Quên mật khẩu?
Đăng nhập bằng Facebook
Đăng ký
Lỗi
Chưa xác định bài thơ
Nội dung liên quan
© 2004-2016 Vanachi
RSS

WARC-Truncated: length

Set-Cookie: PHPSESSID=tud0v9ucf5pv6j7knbpns8s1f1; expires=Sat, 22-Oct-2016 00:35:54 GMT; Max-Age=3600; path=/

Cache-Control: private, must-revalidate

Pragma: no-cache

Vary: Accept-Encoding

Content-Type: text/html


markdown:
Tại đây bạn có thể chơi Zwingo.
Zwingo là một trong những Trò Chơi Thể Thao đã được chúng tôi tuyển chọn đặc biệt.
[Sports Heads: Basketball Championship](http://www.trochoi.net/tr%C3%B2+ch%C6%A1i/sports-heads-basketball-championship.html), [Uphill Rush 2](http://www.trochoi.net/tr%C3%B2+ch%C6%A1i/xe+m%C3%A1y+nh%E1%BA%A3y+d%E1%BB%91c/uphill+rush+2.html), [Bicycle Kick Champ 2014](http://www.trochoi.net/tr%C3%B2+ch%C6%A1i/bicycle+kick+champ+2014.html) và [Puppet Soccer Champions](http://www.trochoi.net/tr%C3%B2+ch%C6%A1i/puppet-soccer-champions.html) hiện cũng được yêu thích trên trang web này.
Báo cáo trò chơi bị lỗi
Hủy báo cáo trò chơi bị lỗi

Connection: Close

Server: nginx

Cache-Control: max-age=900, public

Vary: User-Agent, Accept-Encoding


markdown:
Công ty Cổ phần Vật Giá Việt Nam. Số GCNDT: 011032001615, cấp ngày 21/06/2012, nơi cấp: UBND thành phố Hà Nội
Trụ sở chính: phố Thái Thịnh, phường Trung Liệt, quận Đống Đa, Hà Nội - Chi nhánh Hồ Chí Minh: Đường Lữ Gia, Phường 15, Quận 11, Hồ Chí Minh
Trụ sở chính: phố Thái Thịnh, phường Trung Liệt, quận Đống Đa, Hà Nội - Chi nhánh Hồ Chí Minh: Đường Lữ Gia, Phường 15, Quận 11, Hồ Chí Minh

Set-Cookie: u_nologin_id=589389362; expires=Mon, 30-Jan-2017 00:17:06 GMT; path=/; httponly

Server: Apache


markdown:
Trụ sở chính: phố Thái Thịnh, phường Trung Liệt, quận Đống Đa, Hà Nội - Chi nhánh Hồ Chí Minh: Đường Lữ Gia, Phường 15, Quận 11, Hồ Chí Minh
|
|Giá sản phẩm : |
|
Thỏa thuận
Liên hệ người đăng:
**Nguyen Thanh Bao**
(Thái độ phục vụ khách hàng)
(Tin chính xác như mô tả)
Chưa có người dùng nào đánh giá tin rao vặt này.
Hãy là người đầu tiên đánh giá tin rao vặt này.
[Click để đánh giá](#box_quick_review)
**0919 218 468 – 0908 725 159** **Công ty TNHH ** **1374A QL1A, P. Thới An, Q.12, HCM** **(gần cầu vượt Tân Thới Hệp Q12)** **Email: ** **baont2010@gmail.com** **Yahoo chat: baont2010** **http://dongfeng.mysite.vn/**
[ dongfeng trường giang Xe tải](/3182/xe-tai.html) *Nhà phân phối các dòng xe chuyên dùng: xe tec, xe trộn, xe tải* **Một số sản phẩm do công ty phân phối:** **DONGFENG 3.1 TẤN**
Tổng trọng :6975 kg
-YCD4A22-84
- Dung tích: 2672 (cm3)
- Công suất: 62/3200 (Kw/rpm)
- Kích thước thùng: 5150 x 1990 x 450 mm
- Lốp xe: Trước8.25-16 / Sau (Đôi) 8.25-16
- không có thùng lửng
**DONGFENG 4.9 TẤN**
Tổng trọng : 11000 kg
-YC4D130-20
- Dung tích: 4214 (cm3)
- Công suất: 96/2800 (Kw/rpm)
- Kích thước thùng: 6360x2350x2260
- Lốp xe: Trước 9.00-20 / Sau 9.00-20 (Đôi)
- không có thùng lửng
**DONGFENG 6.8 TẤN**
- Tổng trọng : 13600 kg
- Dung tích: 4257 (cm3)
- Công suất: (100/2800 Kw/rpm)
- Kích thước thùng: 7400x2360x2370
- Lốp xe: Trước. 10.00-20/ Sau 10.00-20 (Đôi)
- không có thùng lửng
**DONGFENG 6.9 TẤN**
- Tổng trọng : 13600 kg
- YC4E135-21
- Dung tích: 4257 (cm3)
- Công suất: (100/2800 Kw/rpm)
- Kích thước thùng: 7400x2360x2370
- Lốp xe: Trước. 11.00-20/ Sau 11.00-20 (Đôi)
- không có thùng lửng
**XE BEN DONGFENG** **BEN DONGFENG 3.45 T**
- Tổng trọng : kg
- YC4E135-21
- Dung tích: (cm3)
- Công suất: Kw/rpm)
- Kích thước thùng:
- Lốp xe: Trước. 10.00-20/ Sau (Đôi)
**BEN DONGFENG 4.95 T**
-Tổng trọng : 10605 kg
- YC4D130-20
- Dung tích: 4214 (cm3)
- Công suất: 96/2800 Kw/rpm)
- Kích thước thùng: 4030 x 2180 x 670 mm
- Lốp xe: Trước. 9.00-20 / Sau 9.00-20 (Đôi)
**BEN DONGFENG 6.98T**
-Tổng trọng : 13275 kg
- YC4E135-21
- Dung tích: 4257 (cm3)
- Công suất: 100/2800 Kw/rpm)
- Kích thước thùng: 4050 x 2170 x 770 mm
- Lốp xe: Trước. 10.00-20/ Sau 10.00-20 (Đôi)
**BEN DONGFENG 6.95T**
-Tổng trọng : 13595 kg
- YC4E135-21
- Dung tích: 4257 (cm3)
- Công suất: 100/2800 Kw/rpm)
- Kích thước thùng: 4180 x 2190 x 810 mm
- Lốp xe: Trước. 11.00-20 / Sau 11.00-20 (Đôi
**XE BEN 2 CẦU** **BEN DONGFENG 4.98T**
-Tổng trọng : 11605 kg
- YC4D130-20
- Dung tích: 4214 (cm3)
- Công suất: 96/2800 Kw/rpm)
- Kích thước thùng: 4030 x 2180 x 670 mm
- Lốp xe: Trước. 90.00-20/ Sau 90.00-20 (Đôi)
**BEN DONGFENG 6.5T**
-Tổng trọng : 13605 kg
- YC4E135-21
- Dung tích: 4257 (cm3)
- Công suất:100/2800 Kw/rpm)
- Kích thước thùng: : 4180 x 2190 x 810 mm
- Lốp xe: Trước. 11.00-20/ Sau 11.00-20 (Đôi)
|Thông tin chung|
|Hãng sản xuất|
Thỏa thuận...
Thỏa thuận...
Thỏa thuận...
Thỏa thuận...
602.000.000 VNĐ
Thỏa thuận...
Thỏa thuận...
Thỏa thuận...
Thỏa thuận...
200.000.000 đ
1.230.000.000 đ
329.500.000 đ
Thỏa thuận...
287.000.000 đ
287.000.000 đ
245.000.000 đ
287.000.000 đ
Thỏa thuận...
287.000.000 đ
287.000.000 đ
329.000.005 đ
216.000.000 đ
245.000.000 đ
245.000.000 đ
155.000.000 đ
Thỏa thuận...
553.000.000 đ
152.500.000 đ
553.000.000 đ
553.000.000 đ
469.000.000 đ

Set-Cookie: u_nologin_id=589387673; expires=Mon, 30-Jan-2017 00:16:12 GMT; path=/; httponly

Cache-Control: no-store, no-cache, must-revalidate, post-check=0, pre-check=0


markdown:
Bóng đá
Châu Âu
Giao thông
Trung Đông
Đời sống
ASEAN
(Vietnam+)
23/06/2015 06:31 GMT+7
Bản in
Châu Mỹ
Châu Âu
Kinh doanh
Đời sống
ASEAN
Châu Phi
Điểm đến

Cache-Control: no-cache


markdown:
Để
**tải phim Hộp Đêm Đồng Tính**về máy, bạn hãy nhập mã xác nhận rồi bấm nút "Lấy link download". Link tải phim sẽ hiện ra cho bạn **download phim Hộp Đêm Đồng Tính (Shortbus)**miễn phí. Hãy chọn server phù hợp mà bạn muốn download để tải phim với tốc độ tốt nhất. [Hoàn tất VietSub HD](http://www.vietsubhd.com/phim/nhuc-bo-doan-3d-1583/)
![ Nhục Bồ Đoàn 3d](https://images2-focus-opensocial.googleusercontent.com/gadgets/proxy?container=focus&gadget=a&no_expand=1&refresh=604800&url=http://www.vietsubhd.com/poster/lh3/-FXVrGETgLxY/VxePlJE926I/AAAAAAAAa-A/OY1xsupefTY/w300-h400/Film.Poster.1583.jpg%3F_iv%3D1475870925)
[Nhục Bồ Đoàn 3d](phim-le/phi-vu-rio_24/)Sex And Zen: Extreme Ecstasy 2011 [Hoàn tất VietSub HD](http://www.vietsubhd.com/phim/mot-thoi-hoang-kim-2410/)
![ Một Thời Hoàng Kim](https://images2-focus-opensocial.googleusercontent.com/gadgets/proxy?container=focus&gadget=a&no_expand=1&refresh=604800&url=http://www.vietsubhd.com/poster/lh3/-wrm5kCxHLlg/VqeR_8t7VmI/AAAAAAAAX5A/8S8zy2NZlqE/w300-h400/Film.Poster.2410.jpg%3F_iv%3D1475870925)
[Một Thời Hoàng Kim](phim-le/phi-vu-rio_24/)Youth Paradise Lost 2011 [Hoàn tất VietSub HD](http://www.vietsubhd.com/phim/tinh-duc-la-chuyen-nho-2-857/)
![ Tình Dục Là Chuyện Nhỏ 2](https://images2-focus-opensocial.googleusercontent.com/gadgets/proxy?container=focus&gadget=a&no_expand=1&refresh=604800&url=http://www.vietsubhd.com/poster/lh3/-MWyCq5C4Et0/Vp-UGbULKaI/AAAAAAAAV7g/GfLajilUQE8/w300-h400/Film.Poster.857.jpg%3F_iv%3D1475870925)
[Tình Dục Là Chuyện Nhỏ 2](phim-le/phi-vu-rio_24/)Sex Is Zero 2 2007 [Hoàn Tất NoSub HD](http://www.vietsubhd.com/phim/cam-giac-toi-loi-5637/)
![ Cảm Giác Tội Lỗi](https://images2-focus-opensocial.googleusercontent.com/gadgets/proxy?container=focus&gadget=a&no_expand=1&refresh=604800&url=http://www.vietsubhd.com/poster/lh3/--Y-AxuHl9_s/Vi-YnD3-8DI/AAAAAAAAIZU/QbxSTUe2UGs/w300-h400/Film.Poster.5637.jpg%3F_iv%3D1475870925)
[Cảm Giác Tội Lỗi](phim-le/phi-vu-rio_24/)Marriage With A Liar 2010 [Hoàn tất VietSub HD](http://www.vietsubhd.com/phim/hop-dem-dong-tinh-1136/)
![ Hộp Đêm Đồng Tính](https://images2-focus-opensocial.googleusercontent.com/gadgets/proxy?container=focus&gadget=a&no_expand=1&refresh=604800&url=http://www.vietsubhd.com/poster/lh3/-MVzNhJZVS24/V0Sfjlc2aRI/AAAAAAAAblw/qxpvtRNZW5I/w300-h400/Film.Poster.1136.jpg%3F_iv%3D1475870925)
[Hộp Đêm Đồng Tính](phim-le/phi-vu-rio_24/)Shortbus 2006

X-Mod-Pagespeed: Powered By VIETSUBHD.COM

Server: Apache/2.2.15 (CentOS)

Cache-Control: max-age=0, no-cache, no-store, must-revalidate, post-check=0, pre-check=0


markdown:
3
![ Điệp Vụ Tam Giác Vàng](https://images2-focus-opensocial.googleusercontent.com/gadgets/proxy?container=focus&gadget=a&no_expand=1&refresh=604800&url=http://www.vietsubhd.com/poster/lh3/-chUxobO5rBc/V_8zlGQQ10I/AAAAAAAAgeU/orQO2B5xyyU/w40-h40/Film.Poster.11429.jpg%3F_iv%3D1475870925)
[Điệp Vụ Tam Giác Vàng](http://www.vietsubhd.com/phim/diep-vu-tam-giac-vang-11429/)Hoàn Tất VietSub HD
4
![ Anh Cần Em](https://images2-focus-opensocial.googleusercontent.com/gadgets/proxy?container=focus&gadget=a&no_expand=1&refresh=604800&url=http://www.vietsubhd.com/poster/lh3/-kKt5X38fTWk/Vhu2449cZYI/AAAAAAAAHeo/abp4_Cep_D0/w40-h40/Film.Poster.5426.jpg%3F_iv%3D1475870925)
[Anh Cần Em](http://www.vietsubhd.com/phim/anh-can-em-5426/)Hoàn Tất VietSub HD
5
![ Vương Triều Nhục Dục](https://images2-focus-opensocial.googleusercontent.com/gadgets/proxy?container=focus&gadget=a&no_expand=1&refresh=604800&url=http://www.vietsubhd.com/poster/lh3/-s91KjtGGGYM/VZ8fYT07fII/AAAAAAAAAFA/zhMb15iksZ0/w40-h40/Film.Poster.3671.jpg%3F_iv%3D1475870925)
[Vương Triều Nhục Dục](http://www.vietsubhd.com/phim/vuong-trieu-nhuc-duc-3671/)Hoàn Tất VietSub HD
6
![ Biệt Đội Cảm Tử](https://images2-focus-opensocial.googleusercontent.com/gadgets/proxy?container=focus&gadget=a&no_expand=1&refresh=604800&url=http://www.vietsubhd.com/poster/lh3/-NEaqWjzzGrA/V74Tz01H25I/AAAAAAAAfY8/Wr6VPPpbJG8/w40-h40/Film.Poster.6504.jpg%3F_iv%3D1475870925)
[Biệt Đội Cảm Tử](http://www.vietsubhd.com/phim/biet-doi-cam-tu-6504/)Hoàn Tất VietSub HD
7
![ Cheeky Gợi Cảm](https://images2-focus-opensocial.googleusercontent.com/gadgets/proxy?container=focus&gadget=a&no_expand=1&refresh=604800&url=http://www.vietsubhd.com/poster/lh3/-eO8mDqkAED8/VcyfOLvZteI/AAAAAAAABWc/WLOygUScyVg/w40-h40/Film.Poster.3997.jpg%3F_iv%3D1475870925)
[Cheeky Gợi Cảm](http://www.vietsubhd.com/phim/cheeky-goi-cam-3997/)Hoàn Tất VietSub HD
8
![ Người Vận Chuyển 4](https://images2-focus-opensocial.googleusercontent.com/gadgets/proxy?container=focus&gadget=a&no_expand=1&refresh=604800&url=http://www.vietsubhd.com/poster/lh3/-ezrWK8W39ug/VZF95iq660I/AAAAAAAADH8/yKpebrLsVXk/w40-h40/Film.Poster.3543.jpg%3F_iv%3D1475870925)
[Người Vận Chuyển 4](http://www.vietsubhd.com/phim/nguoi-van-chuyen-4-3543/)Hoàn Tất VietSub Bản Đẹp
9
![ Vợ Chồng Dân Chơi](https://images2-focus-opensocial.googleusercontent.com/gadgets/proxy?container=focus&gadget=a&no_expand=1&refresh=604800&url=http://www.vietsubhd.com/poster/lh3/-tnHeG4UC8CY/VfP8-4ujwQI/AAAAAAAADhA/lIG8NEaTqEw/w40-h40/Film.Poster.4489.jpg%3F_iv%3D1475870925)
[Vợ Chồng Dân Chơi](http://www.vietsubhd.com/phim/vo-chong-dan-choi-4489/)Hoàn Tất VietSub HD
10
![ Bảy Tay Súng Huyền Thoại](https://images2-focus-opensocial.googleusercontent.com/gadgets/proxy?container=focus&gadget=a&no_expand=1&refresh=604800&url=http://www.vietsubhd.com/poster/lh3/-dnQ0Jga9xQA/WAUbKAI-PfI/AAAAAAAAglM/FTr4WqHPkIA/w40-h40/Film.Poster.11450.jpg%3F_iv%3D1475870925)
[Bảy Tay Súng Huyền Thoại](http://www.vietsubhd.com/phim/bay-tay-sung-huyen-thoai-11450/)Hoàn Tất VietSub CAM