Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

A_S_Safronov_Modelirovanie_innovatsionnogo_ros

.pdf
Скачиваний:
21
Добавлен:
16.03.2015
Размер:
566.36 Кб
Скачать

Таким образом, инновационная модель Вольтерра – Лотка имеет

вид:

x/

= (a bx )x ,

(8.1.4)

 

1

2

1

 

 

= (c + dx )x

,

x/

 

2

1

2

 

где a,b,c,d >0 .

Графически процесс слияний и поглощений выглядит следующим образом (рис.8, 9):

4 4

X1<2 >

 

 

2

X1<3 >

 

 

 

 

0 0

0

5

10

0

X1<1>

10

Рис. 8. Стационарное состояние модели роста фирм Вольтерра – Лотка при слияниях и поглощениях

6

 

 

 

 

5

 

 

X4<2>

 

 

 

X4<3>

 

 

 

0

0

 

 

 

0

5

10

 

0

X4<1>

10

Рис. 9. Динамическое состояние модели роста фирм Вольтерра – Лотка при слияниях и поглощениях

51

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

x = c

 

1

d

(8.1.5)

 

вызывает снижение количества фирм «хищников».

Сокращение количества фирм «хищников» и фирм «жертв» происходит до тех пор, пока число фирм «хищников» не достигнет величины, равной

x2

= a .

(8.1.6)

 

b

 

С этого момента начинается рост фирм «жертв», а через некоторое время начинается рост фирм «хищников». На графике четко видна цикличность процесса слияния и поглощения, причем количество

фирм «хищников» и фирм «жертв» колеблется возле величин x1 = dc

и x2 = ba . Периодичность процесса видна на фазовой плоскости: фа-

зовая кривая (x1 (t), x2 (t)) – замкнутая линия.

Самая левая точка этой кривой x2 = ba – это точка, в которой

число фирм «жертв» достигает наименьшего значения, самая правая точка – точка максимального количества фирм «жертв». Между этими точками количество фирм «хищников» сначала убывает до нижней

фазовой точки кривой x1 = dc , где достигает своего наименьшего зна-

чения, а затем растет до верхней точки фазовой кривой. Фазовая плоскость представлена на графике (рис.10).

52

3.15719

4

 

 

 

 

 

 

 

 

 

 

 

 

 

X1<3>

3

 

 

 

 

 

 

X2<3

>

 

 

 

 

 

 

 

<3

>

2

 

 

 

 

 

 

X3

 

 

 

 

 

 

 

 

X4<3

>

1

 

 

 

 

 

 

0.664972

0 0

1

2

3

4

5

6

 

 

 

 

0.534142

 

X1<2>, X2<2> ,X3<2>,X4<2>

5.011845

Рис. 10. Фазовая плоскость модели роста Вольтерра – Лотка при слияниях и поглощениях

На языке дифференциальных уравнений это означает, что система инновационного роста имеет стационарное состояние

x/

= 0, x/

= 0 , которое достигается в точке

c

, a

. Если система

 

 

1

2

 

 

 

 

 

 

d

b

 

инновационного роста в начальный момент времени находится в стационарном состоянии, то решения x1 (t) и x2 (t) не будут изме-

няться во времени, останутся постоянными.

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

§ 8.2. Модель Вольтерра – Лотка с логистической поправкой

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

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

53

ведение потребителей, сезонность спроса и предложения на инновационную продукцию и др.

Модель Вольтерра – Лотка с логистической поправкой имеет

вид:

x1/ = (a bx2 )x1 −αx12 ,

x2/ = (c + dx1)x2 −αx22.

(8.2.1)

В этом случае поведение решений в окрестности стационарной точки меняется в зависимости от величины и знака логистической поправки α .

Ниже на графике (рис.11,12) представлены динамическое состояние инновационной системы и ее фазовый портрет на плоскости.

6

6

 

 

 

 

 

X4<2 >

4

 

 

 

 

 

X4<3 >

2

 

 

 

 

 

0

0

 

 

 

0

5

10

 

0

X4<1 >

10

Рис. 11. Динамическое состояние модели инновационного роста фирм модели Вольтерра – Лотка с учетом логистической поправки

3.068 144

4

 

 

 

 

 

 

 

 

 

 

 

X1<3 >

3

 

 

 

 

 

X2<3 >

 

 

 

 

 

 

X3<3 >

2

 

 

 

 

 

 

 

 

 

 

 

X4<3 >

1

 

 

 

 

 

0.718 673

0

 

 

 

 

 

 

 

 

 

 

 

 

0

1

2

3

4

5

 

0.6 8568

X1<2> , X2<2 > , X3<2 > ,X 4<2 >

4.316 99

Рис. 12. Фазовый портрет модели инновационного роста фирм модели Вольтерра – Лотка с учетом логистической поправки

54

Из приведенных ниже графиков следует, что стационарная точка превращается в устойчивый фокус, а решение – в затухающие колебания. При любом начальном состоянии через некоторое время состояние системы инновационного роста становится близким к стационарному и стремится к нему приt →∞ .

§ 8.3. Модель Холлинга – Тэннера

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

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

Модель Вольтерра – Лотка неустойчива относительно возмущений, поскольку ее стационарное состояние – центр.

Существует другой вид моделей инновационного роста, в которых возникают незатухающие колебания, – это модели, имеющие на фазовых портретах предельные циклы. Такая модель применяется в системах конкурентной среды – это модель Холлинга – Тэннера.

Скорость роста количества фирм «жертв» x1/ в этой модели равна сумме трех величин:

-скорости роста в отсутствие фирм «хищников» rx1 ;

-влияния конкуренции между фирмами «хищниками» за поглощения и слияния с фирмами «жертвами» при ограниченном коли-

честве последних rx

x1

;

 

 

 

 

 

 

 

1 K

 

 

x1

 

- влияния фирм «хищников» wx

2

в предположении,

D + x

 

 

 

 

 

 

 

 

1

 

что у фирм «хищников» пропадает интерес к поглощениям и слияниям с конкурентами в случае их насыщения.

Скорость роста количества фирм «хищников» x2/ строится также, как и в модели Вольтерра – Лотка, в предположении, что фирмы

55

«жертвы» встречаются редко. Если для поддержания активного состояния одной фирмы «хищника» требуется J число фирм «жертв»,

то количество фирм «жертв», равных x1 , сможет обеспечить актив-

ность xJ1 фирм «хищников». Модель роста численности фирм «хищ-

ников», в которой их число не может превысить эту критическую величину, имеет вид:

 

 

x2/ = x2 (s

sJ

x2 ) .

 

(8.3.1)

 

 

 

 

 

 

 

 

 

 

 

 

x1

 

 

 

Таким образом, модель Холлинга – Тэннера имеет следующий

вид:

 

 

 

x1

 

 

 

 

 

x1

 

 

x/

= r(1

 

 

)x

wx

,

 

 

 

 

 

 

1

 

 

K

1

2

D + x1

 

(8.3.2)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

J

 

 

 

 

 

 

 

x/

= s(1

)x

,

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

x1

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где r,s, K, D, J >0 .

Приs <

r

 

K D 2

на фазовом портрете системы будет ус-

K

1+ D

 

 

 

тойчивый предельный цикл.

56

ТЕМА 9. Многокритериальная модель выбора оптимальной модели инновационного развития

§ 9.1. Постановка задачи многокритериального выбора

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

ление u0 из множества наихудшему критерию:

~

 

 

U , которое доставляет наилучшее значение

u0

= arg max min Rk [u].

(9.1.1)

 

u U k K

 

Задачу многокритериального выбора представим в форме минимакса; при этом нормализованные критерии Rk [u]= 1 Rk [u] ми-

нимизируются (нормализация

осуществляется по формуле

 

 

k [u]=

Rk [u]

 

,k K , где

Rk*

максимальное значение k-го крите-

R

 

Rk*

 

 

 

 

 

 

 

рия), что соответствует максимизации исходных критериев, а принцип

минимакса записывается в форме

 

k [

 

]

.

(9.1.2)

 

 

[

u

0 ]

 

 

u

 

 

 

 

= min max

 

 

 

R

 

 

R

 

 

 

 

 

 

u U

k K

 

 

 

 

 

 

Рис. 13. Формирование гипербол, аппроксимирующих множество Парето

57

u*n .

Для выбора управления u0 необходимо определить К векторов управления, обеспечивающих такие сочетания критериев, при которых значения (К-1) критериев фиксированы, а один критерий достигает минимума. Далее определяются коэффициенты аппроксимирующей поверхности Г (рис. 13), после чего вычисляются координаты центра С аппроксимирующей поверхности. Сочетание критериев в центре аппроксимирующей поверхности и соответствующий вектор управления представляют собой приближенное решение многокритериальной задачи. Уточнение приближенного решения выполняется с помощью итерационной процедуры.

§ 9.2. Метод выбора управления на графе Парето-оптимальных управлений

При этом многокритериальный выбор предусматривает сопоставление вершин графа по критерию Ωm , являющемуся количественной характеристикой относительной предпочтительности управления u*m по сравнению с другими Парето-оптимальными управлениями:

K

K

 

K

 

K K

 

 

Ωm =Sjm =(K1)

 

j um*

 

 

m u*j

 

∑∑

 

 

u*j ,

(9.2.1)

R

R

Ri

j=1

j=1

 

j=1

 

j=1 i=1

 

 

jm

 

 

 

 

jm

 

jmim

 

 

где Snm (рис. 14) представляет собой векторную характеристику перехода от управления u*n к управлению u*m – при Snm>0, управление u*m является более предпочтительным по векторному критерию, чем управление

Параметр Ωm может использоваться в качестве интегрального критерия выбора компромиссного управления из множества Парето:

u0

= arg max Ωm [u].

(9.2.2)

 

m K

 

58

Рис. 14. Граф управлений (а) и циклические подграфы (б)

Управление u0 является компромиссным в том смысле, что при переходе к нему от других управлений относительные приросты критериев максимально превышают относительные потери критериев.

§ 9.3. Алгоритм выбора пути инновационного развития

Шаг 1. Определение моделей инновационного развития с оптимальными критериями их эффективности.

Формирование набора X k* , k = 1,...,K моделей инновационного роста, оптимизирующих каждый из критериев эффективности. Расчет значений критериев оптимизации f1 (X ), f2 (X ),..., fn (X ) для каждой модели X k , k K , полученной в результате расчета оптимальных критериев без учета остальных критериев оптимизации.

Выделение моделей инновационного роста{X1, X2 ,..., Xn Xk },

которым соответствует максимальное либо минимальное значение соответствующего критерия оптимизации:

X1{f1max (X ), f21 (X ),..., fn1 (X )},

 

X 2 {f12 (X ), f2min (X ),..., fn2 (X )},

(9.3.1)

................................................,

 

X n {f1n (X ), f2n (X ),..., fnmax (X )},

 

59

где

f1max (X ) = max{f1k (X )}, k =1..K ,

f2min (X ) = min{f2k (X )}, k =1..K ,

......................................................., fnmax (X ) = max{fnk (X )}, k =1..K .

Шаг 2. Производится нормализация критериев. Исходя из того, что критерии оптимизации {f1k (X ), f2k (X ), ..., fnk (X )} имеют раз-

личную размерность и экономический смысл, необходимо произвести их нормализацию:

 

 

f

 

(X ) f min

,k 1..n ,

 

fk (X ) =

k

 

 

k

(9.3.2)

f (X )kmax fkmin

 

 

 

 

где fk (X ) – текущее значение критерия для k - той структуры;

fk (X ) – нормализованное значение k-го критерия;

f (X )mink – минимальное значение k-го критерия, полученное в

результате решения однокритериальной задачи оптимизации без учета остальных критериев;

f (X )maxk – максимальное значение k-го критерия, полученное в

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

 

 

Шаг 3. Определим параметр hi, j

по формуле:

 

 

 

 

 

 

 

k

 

 

 

 

i, j

 

 

f j (X ) fi (X )

 

f j (X ) fi (X )

 

f j (X ) fi(X )

,

(9.3.3)

hk

=

 

 

+

 

 

+...+

 

 

 

f1max (X ) f1min (X )

f2 max (X ) f2 min (X )

fn max (X ) fn min (X )

 

 

 

 

 

 

 

 

где hi, j

– параметр, отражающий долю прироста (потерь) k-го крите-

 

k

 

 

 

 

 

 

 

 

рия относительно отклонения его максимального значения от минимального при переходе от модели инновационного роста i к модели j .

Шаг 4. Построение графа, где вершинам соответствуют модели, оптимальные по каждому критерию оптимизации. Ребра графа представляют собой управление, по которому осуществляется переход от модели j к модели i .

60

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]