Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
1_teoriya-i-metodyi-stat-prog-ya.pdf
Скачиваний:
215
Добавлен:
16.02.2016
Размер:
2.34 Mб
Скачать

98

10

11 7. Математические модели процессов эволюционного развития техники

7.1. Математическое моделирование процессов развития техники

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

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

составление по известным законам природы теоретических причинноследственных математических моделей, определяющих прогнозируемые процессы развития и оценки параметров этих процессов по предыстории и состоянию в настоящий момент;

изготовления и испытания макетов и экспериментальных образцов, а также физическое моделирование;

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

99

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

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

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

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

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

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

Прогнозирование по математической модели выполняется по этапам:

построение математической модели;

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

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

100

– разработка прогноза и анализ его точности и достоверности. При математическом моделировании процессов как развития, так и

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

 

 

Математическая

 

 

 

 

 

 

модель части

 

 

 

 

 

 

 

 

 

 

 

Прогнозируемый

 

 

ЭВМ

Решение

процесс

 

 

 

 

(прогноз)

 

 

 

 

 

 

 

 

 

 

 

 

 

Физическая модель

 

 

 

 

 

 

части процесса

 

 

 

Рис. 7.1. Сочетание математического и физического моделирования

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

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

Система уравнений, определяющая точность полета, может быть разделена на следующие две системы:

уравнения динамики самолета, описывающие движение при заданных командах управления;

уравнения, определяющие величину команд.

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

t < t;

101

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

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

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

7.2.Прогнозная математическая модель динамики замещения

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

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

Пусть n(t) – некоторая мера распространения новой техники в момент времени t (объем распространения);

N (t) – верхняя граница объема распространения новой техники на момент времени t , выраженная в тех же единицах, что и n(t).

При этом процесс замещения нового оборудования в некоторый момент времени t = tк значения N (tк )

0 n(t)< N(t)< ∞; n(t)n(t) при

N (t)

102

n(tк )= N (tк ), n(t)N(tк ) при t tк tк (0,).

Переменную n(t) можно назвать кумулятивной, или интегральной,

функцией внедрения новой техники, значение которой в начальный момент времени t = 0 зависит от того, что считать за начало внедрения: начало выпуска первых образцов новой техники или начальное поступление этих систем. В первом случае n(0)= 0 , в другом n(0)> 0 . Наиболее естественно первое предположение, хотя в ряде моделей может оказаться необходимым второе допущение.

В общем случае границы области распространения новой техники

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

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

m(t)= Nn((tt)).

При N (t )= const динамические характеристики m(t) совпадают с аналогичными характеристиками n(t). Можно считать, что относительная интегральная функция внедрения m(t) – монотонно возрастающая функция времени, изменяющаяся в интервале [0, 1]:

m(0)= 0 ; m(tк )=1.

H (t)=

103

Без существенных искажений реального процесса внедрения можно считать функцию m(t) непрерывной и дифференцируемой. Тогда дифференциальная функция внедрения

 

dn

 

 

 

 

 

 

 

n (t)= dt

 

и дифференциальная относительная функция внедрения

 

dm

 

 

 

dt

 

 

m (t)=

 

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

к концу периода внедрения их

значение монотонно убывает:

 

 

 

 

 

при

t →∞.

n (t),

m (t)0

Поскольку процесс внедрения имеет верхний предел и есть процесс насыщения потребностей, то целесообразно ввести еще одну переменную, характеризующую скорость приближения процесса к концу – интенсивность внедрения. Эта переменная определяет темп приближения процесса распространения нового к уровню полного насыщения в нем. Она определяется следующим образом:

n(t)

N (t)n(t)

или

h(t)=

 

 

m(t)

 

=

 

1

 

 

dm

.

1

m(t)

1m(t)

 

 

 

 

dt

Функция насыщения H (t) или h(t) играет в моделировании динамики замещения центральную роль, так как она определяет конкретный вид функции

n(t) или m(t).

Несколько изменив последнее выражение, получим m(t)= dmdt = h(t)[1m(t)].

104

Проинтегрировав это уравнение при ранее введенных допущениях, получим

t

h(t )dt

m(t)=1e 0

при m(0)= 0

и

 

t

 

h(t)dt → ∞ при

t →∞, так как m()=1.

0

 

Интенсивность замещения (функция насыщения) в общем случае может

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

Поэтому функцию h можно записать в следующем виде: h = h{t , m (t ),[xi ]},

где [xi ] – множество экзогенных факторов, определяющих конкретный процесс замещения.

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

Идентификацию этой функции можно провести, учитывая ее линейное приближение относительно m(t):

h = h0 + h1m(t),

где h0 – уровень насыщения при t = 0 , h1 – скорость роста функции.

Здесь предполагается, что h0 и h1 не зависят от t, но могут зависеть от экзогенных переменных [xi ].

Интегрирование дифференциального уравнения dmdt = [h0 + h1m(t )][1 m(t )]

105

дает следующую интегральную функцию замещения при m(0)= 0 :

 

h0

[1e

(h0 +h1 )t

 

m(t)=

 

(h

 

+h )t].

(7.1)

h

+ h e

0

 

 

1

 

0

1

 

 

 

 

 

Приняв h0 = 0 или h1 = 0 , можно получить два еще более простых варианта модели, которые чаще всего используются для описания динамики замещения.

При h0 = 0 имеем

dmdt =[h1m(t)][1m(t)],

m(t)=

 

1

,

(7.2)

 

+e

k h t

1

1

 

 

 

 

 

 

 

где k = ln1m( ()0), то есть получаем для интегральной функции m 0

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

При h1 =1 выражение (7.1) будет иметь вид

 

m(t)=1eh0t ,

(7.3)

то есть процесс замещения следует функции экспоненциального распределения вероятностей.

Из соотношений (7.1) – (7.3) можно получить также аналитические выражения m(t) и h(t) как функции времени.

Например, для h(t) имеем

при h = h0 + h1m(t)

при h = h0

при h = h1m(t)

h(t)=

 

h0 (h0 + h1 )

 

;

h

+ h e

(h

0

+h

)t

 

 

1

 

0

1

 

 

 

 

 

 

 

h(t)= h0 = const ;

 

 

h(t)=

 

h1

 

 

.

 

 

1+e

k h t

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

106

Графически функции (7.1) – (7.3) имеют вид s-образных кривых, что вполне согласуется с подавляющим большинством наблюдаемых данных. При этом интегральные функции замещения имеют точки перегиба tп, в которых скорость внедрения новой техники имеет максимальное значение.

Для функции (7.1) имеем

t

п

=

k +1

ln

h1

;

m(t

п

)=

1

h0

,

h

+ h

h

2

2h

 

 

 

 

 

 

 

 

 

 

 

0

1

0

 

 

 

 

 

 

1

 

для логистической кривой (7.2)

 

k

 

 

1

 

h1

 

 

h1

 

tп = h

;

m(tп )= 2 ;

 

;

h(tп )= 2 .

m (t)= 4

1

 

 

 

 

 

 

 

 

 

s-образный вид кривых позволяет аппроксимировать массив данных, характеризующий реальный процесс распространения новой техники во времени с помощью различных известных функций распределения вероятностей.

Так, например, обобщив функцию насыщения в выражении (7.3) в виде

h0 = h(t)= bct c1 ,

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

t c

m(t)=1e b ,

которое имеет также точку перегиба

1

tп = bc 11c c .

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

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

107

военной техники (обозначенное на рис. 7.2 ступенчатой линией) на период с 1980 по 2000 г. (всего рассматривалось 156 новых образцов). В данном случае процесс замещения с достаточной точностью аппроксимируется s-образной функцией вида

m (t )= L / {a + exp [b exp (βt )]}.

После определения параметров модели получаем зависимость

m(t )= 109 ,3 / {0,093 + exp [6,997 exp (0,376 (t 1980 ))]},

позволяющую (см. рис. 7.2) прогнозировать динамику замещения, высокая сходимость s-образной кривой с фактическими результатами свидетельствует о том, что предлагаемая модель пригодна для практического использования.

m(t),%

80

60

40

20

1990

1995

2000

2005

2010 t, год

108

12

13 8. Экспертные методы прогнозирования. Морфологический анализ. Прогнозирование технического облика образца изделия

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

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

Наиболее общая классификация методов генерирования идей и поиска новых технических решений выделяет два основных класса таких методов:

трансформационные методы;

морфологические методы.

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