Модели и методы поддержки сбалансированного развития региональных экономических систем (110
..pdfЛАБОРАТОРНЫЕ РАБОТЫ
Лабораторная работа 1
Рассчитать траекторию развития некоторой РЭС на три года, используя модель поиска оптимальной траектории РЭС и метод Соболя, чтобы суммарный объем валовых выпусков 7 представленных ВЭД с учетом коэффициентов важности был максимальный. Данные для расчетов приведены в таблице ниже.
Величина дополнительного финансового ресурса (в млн руб):
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Период |
|
|
|
|
|
K (t) |
|
|
|
|
|
|
|
|
|
L(t) |
|
|
||||
1 |
|
|
|
|
|
8708 |
|
|
|
|
|
|
|
|
|
|
a |
|
|
|||
2 |
|
|
|
|
|
|
|
b |
|
|
|
|
|
|
|
|
|
6240 |
|
|
||
3 |
|
|
|
|
|
|
|
c |
|
|
|
|
|
|
|
|
|
d |
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
5 |
|
|
|
|
|
|
|
№ |
|
1 |
|
|
|
2 |
|
3 |
|
4 |
|
|
|
|
|
6 |
|
7 |
||||
λi |
|
e |
|
|
f |
|
j |
|
0,67 |
|
0,1 |
|
|
|
0,48 |
|
0,3 |
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
Доля выбытия |
|
Коэффициент |
|
Минимальная |
||||||||||||
|
|
|
|
|
|
основных |
|
|
доля выпуска, |
|||||||||||||
Элемент РЭС |
|
фондоемкости, |
|
|||||||||||||||||||
|
|
|
|
|
|
производственных |
|
|
|
|
|
ϕ j |
|
|
|
идущая на |
||||||
|
|
|
|
|
|
фондов, d j |
|
|
|
|
|
|
|
|
потребление, gi |
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
Промышленность |
0,011 |
|
|
|
|
1,121 |
|
|
|
|
0,2 |
|||||||||||
С/Х |
|
|
|
|
0,007 |
|
|
|
|
0,476 |
|
|
|
|
0,49 |
|||||||
Строительство |
0,025 |
|
|
|
|
0,41 |
|
|
|
|
0,4 |
|||||||||||
Транспорт |
|
|
|
|
0,009 |
|
|
|
|
0,754 |
|
|
|
|
0,45 |
|||||||
Торговая и |
|
|
|
|
0,004 |
|
|
|
|
0,576 |
|
|
|
|
0,62 |
|||||||
комерческая |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
деятельность |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
Отрасли, |
|
|
|
|
0,002 |
|
|
|
|
0,345 |
|
|
|
|
0,4 |
|||||||
оказывающие |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
нерыночные услуги |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
Другие отрасли |
0,003 |
|
|
|
|
0,457 |
|
|
|
|
0,45 |
|||||||||||
Матрица коэффициентов прямых затрат: |
|
|
|
|
|
|
||||||||||||||||
0,1 |
|
0,03 |
|
|
|
0,003 |
|
|
0,03 |
|
|
|
|
0,02 |
0,1 |
|
|
0,03 |
||||
0,01 |
|
0,02 |
|
|
|
0,001 |
|
|
0,01 |
|
|
|
|
0,001 |
0,01 |
|
0 |
|||||
0,0001 |
|
0,0025 |
|
|
0 |
|
|
0,004 |
|
|
|
0,0001 |
0 |
|
|
|
0 |
|||||
0,002 |
|
0,001 |
|
|
|
0,002 |
|
|
0,004 |
|
|
|
0,006 |
0,006 |
0,005 |
|||||||
0,003 |
|
0,0005 |
|
|
0,000025 |
|
0,0008 |
|
|
|
0,0006 |
0,0004 |
0,0005 |
|||||||||
0 |
|
0,001 |
|
|
|
0,0001 |
|
0,0002 |
|
|
|
0,002 |
0,002 |
0,006 |
||||||||
0 |
|
0 |
|
|
|
0 |
|
|
0 |
|
|
|
|
0 |
|
|
0,1 |
|
|
0,01 |
31
Для расчетов для всех элементов РЭС в качестве производственной функции была взята функция Кобба – Дугласа со следующими параметрами:
ВЭД |
1 |
2 |
3 |
4 |
5 |
6 |
7 |
Ai |
0,15 |
0,15 |
0,19 |
0,14 |
106,5 |
5,4 |
0,16 |
αi |
0,9 |
0,8 |
0,9 |
0,98 |
0,38 |
0,23 |
0,78 |
Величины основного капитала, фондов оплаты труда и валового выпуска в начальный период времени:
ВЭД |
1 |
2 |
3 |
4 |
5 |
6 |
7 |
Xi0 |
47320,5 |
4786,1 |
767,3 |
7980,32 |
589321,15 |
15386,12 |
5900,1 |
Ki0 |
774033 |
91336 |
5851 |
65176 |
97640 |
91897 |
126694 |
L0i |
150 |
707 |
369 |
553 |
996 |
s |
775 |
Расчет необходимо сделать при следующих значениях параметров:
Вариант |
a |
b |
c |
d |
e |
f |
j |
1 |
3732 |
14569 |
18080 |
7750 |
0,54 |
0,8 |
0,45 |
2 |
4521 |
10657 |
10233 |
10256 |
0,65 |
0,5 |
0,23 |
3 |
6123 |
9800 |
12005 |
11500 |
0,7 |
0,4 |
0,5 |
4 |
2630 |
16254 |
17385 |
5700 |
0,1 |
0,5 |
0,32 |
5 |
6230 |
8000 |
14678 |
8690 |
0,4 |
0,9 |
0,5 |
6 |
1754 |
12234 |
20127 |
5456 |
0,6 |
0,85 |
0,6 |
Лабораторная работа № 2
На основе данных из первой лабораторной работы и данных для 6 вариантов сделайте новый расчет, используя две функции цели:
F1 |
(t) = λit Pi |
(t) → max ; F2 (t ) = min Li (t ) → max. |
||
|
n |
|
|
|
|
i=1 |
i |
Xi (t ) |
|
|
|
|
|
Лабораторная работа № 3
Используя данные из задачи 4 (глава 1) и задачи 6 (глава 3), сделайте расчет для Воронежской области, применяя модель поиска оптимальной траектории РЭС и метод Соболя, чтобы суммарный объем валовых выпусков представленных ВЭД был максимальный.
32
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
Литература
Основная
1.Баева Н.Б. Математические методы оценки и наращивания экономического потенциала региона : монография / Н.Б. Баева, Д.В. Ворогушина // Воронежскийгосударственныйуниверситет. – Воронеж: ИПЦВГУ, 2012. – 192 с.
2.Белобродский А.В. Методические указания по применению производственных функций в экономико-математическом моделировании // А.В. Белобродский, В.В. Давнис, И.Н. Щепина. – Воронеж, 1996. – 19 с.
3.Гранберг А.Г. Основы региональной экономики / А.Г. Гранберг. –
М. : ГУ ВШЭ, 2000. – 495 с.
4.Жак С.В. Математическая модель менеджмента и маркетинга / С.В. Жак. – Ростов н/Д. : ЛаПо, 1997. – 307 с.
5.Коссов В.В. Межотраслевой баланс / В.В. Коссов – М. : Экономика, 1979. – 270 с.
6.Леонтьев В. Межотраслевая экономика / В. Леонтьев. – М. : Эконо-
мика, 1997. – 250 с.
7.Острейковский В.А. Теория систем : учеб. для студ. вузов / В.А. Острейковский. – М. : Высш.шк., 1997. – 239 с.
8.Экономико-математические методы и модели / под ред. В.В. Федо-
сеева. – М. : Юнити, 2002. – 391 с.
Дополнительная
1.Андронникова Н.Г. Модели и методы оптимизации региональных программ развития / Н.Г. Андронникова [и др.]. – М. : ИПУ РАН, 2001. – 60 с.
2.Баева Н. Б. Математические методы поддержки процесса перехода региональных экономических систем в режим устойчивого развития : монография / Н. Б. Баева, Е. В. Куркин ; Воронежский государственный университет. – Воронеж : Издательский дом ВГУ, 2015. – 180 с.
3.Баева Н.Б. Модели и методы опережающего развития региональной экономики / Н.Б. Баева, Д.В. Ворогушина, Е.А.Сергеева (Е.А. Пронина) // Современная экономика : проблемы и решения. – 2013. – № 1. – С. 169–179.
4.Бурков В.Н. Модели и методы оптимизации региональных программ развития / В.Н. Бурков [и др.]. – М. : ИПУ РАН, 2001. – 60 с.
5.Дрогобыцкий И.Н. Системный анализ в экономике / И.Н. Дрогобыцкий. – М. : Финансы и статистика, 2007. – 512 с.
6.Лексин В.Н. Государство и регион. Теория и практика государственного регулирования территориального развития / В.Н. Лексин, А.Н. Шве-
цов. – М. : ЛКИ, 2007. – 368 с.
7.Пчелинцев О.C. Региональная экономика в системе устойчивого развития / О.С. Пчелинцев. – М. : Наука, 2004. – 260 с.
8.Шелобаев С.И. Математические методы и модели в экономике, финансах, бизнесе : учеб. пособие / С.И. Шелобаев. – М. : ЮНИТИ, 2000. – 249 с.
33
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
П р и л о ж е н и е 1 |
||
|
|
|
|
|
|
|
|
|
|
Сводная таблица производственных функций |
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
Степень |
Постоянство |
|
|
|
|
||
Название, |
|
|
одно- |
эластич- |
|
|
|
|
||||||||||
|
|
ности заме- |
Особенности |
|
||||||||||||||
общий вид |
|
|
родно- |
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
сти* |
ны факто- |
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
ров** |
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Функция |
|
|
|
|
|
|
|
1 |
+ |
|
|
Предназначена для |
моделирования |
|||||
Леонтьева, |
|
|
|
|
|
|
(σ = 0) |
|
строго детерминированных техно- |
|||||||||
|
|
x |
x |
2 |
|
|
|
|
|
|
логий. Обычно используется для |
|||||||
y = min |
|
1 |
, |
|
|
|
|
|
|
|
|
описания мелкомасштабных |
или |
|||||
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
a1 |
a2 |
|
|
|
|
|
полностью |
автоматизированных |
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
производственных объектов |
|
||
Функция |
|
|
|
Коб- |
α1 + α2 |
+ |
|
|
Используется для описания средне- |
|||||||||
ба – Дугласа, |
|
|
|
(σ = 1) |
|
масштабных |
хозяйственных объек- |
|||||||||||
y = a xa1 xa2 |
|
|
|
|
|
|
|
|
|
тов. Вовлечение новой единицы ре- |
||||||||
0 |
|
1 |
2 |
|
|
|
|
|
|
|
|
|
|
|
сурса приносит эффект, пропорцио- |
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
нальный средней производительно- |
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
сти имеющегося ресурса |
|
||
Линейная |
|
|
|
|
|
|
|
1 |
+ |
|
|
Применяется |
для |
моделирования |
||||
функция, |
|
|
|
|
|
|
|
|
(σ = ∞ ) |
|
крупномасштабных систем. Особую |
|||||||
y = a1 x1 + a2 x2 |
|
|
|
|
|
|
роль играет предпосылка о постоян- |
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
стве предельных производительно- |
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
стей факторов или об их неограни- |
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ченной замещаемости |
|
||
Функция |
|
|
|
|
|
|
|
2 |
- |
|
|
Обычно используется для описания |
||||||
Аллена, |
|
|
|
|
|
|
|
|
|
|
|
|
мелкомасштабных |
производствен- |
||||
y = a |
x x |
2 |
− a x2 |
− a |
2 |
|
|
|
ных систем |
с ограниченными |
воз- |
|||||||
0 |
|
1 |
|
|
|
1 |
1 |
|
|
|
|
можностями переработки ресурсов, |
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
в которых чрезмерный рост любого |
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
из факторов оказывает отрицатель- |
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ное воздействие на объем выпуска |
|||
Функция CES, |
|
α3 α4 |
+ |
|
|
Применяется в случаях, когда отсут- |
||||||||||||
y = (a xa3 |
+ a |
2 |
xa3 |
)a4 |
|
1 |
|
ствует точная информация об уров- |
||||||||||
1 |
1 |
|
|
|
|
|
2 |
|
|
(σ = |
|
) |
не взаимозаменяемости производст- |
|||||
|
|
|
|
|
|
|
|
|
|
|
|
1 − a3 |
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
венных факторов и вместе с тем есть |
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
основания предполагать, что этот |
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
уровень существенно не изменится |
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
при изменении объемов вовлекае- |
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
мых ресурсов. Может быть исполь- |
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
зована для |
моделирования систем |
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
любого уровня |
|
|
34
Функция LES, |
|
|
|
|
α3 α0 |
– |
Рекомендуется для описания произ- |
|||||||||||||||||||||
y = xa0 (a x |
|
|
|
+ a |
2 |
x |
2 |
)a3 |
|
|
водственных процессов, |
у которых |
||||||||||||||||
|
1 |
|
|
1 |
1 |
|
|
|
|
|
|
|
|
|
|
|
возможность |
замещения |
вовлекае- |
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
мых факторов существенно зависит |
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
от их пропорций |
|
|
||
Функция Солоу, – |
– |
Может использоваться примерно в |
||||||||||||||||||||||||||
y = (a xa3 |
|
+ a |
2 |
xa4 )a5 |
|
|
тех же случаях, что и функция CES, |
|||||||||||||||||||||
|
1 |
|
1 |
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
но не требуется предположения об |
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
однородности. Может моделировать |
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
системы любого масштаба |
||||
Ограниченная |
|
|
|
|
– |
– |
Предназначена для описания двух- |
|||||||||||||||||||||
функция CES, |
|
|
|
|
|
|
режимного производственного про- |
|||||||||||||||||||||
|
x1 |
|
|
|
|
x2 |
|
|
|
|
|
|
|
a5 |
|
|
|
|
|
цесса, в котором один из режимов |
||||||||
y = min |
|
|
|
, |
|
|
|
,(a3 x1 |
|
+ a4 |
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
a2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
характеризуется отсутствием заме- |
||||||
|
a1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
няемости факторов, другой – нену- |
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
левой постоянной (но не известной |
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
заранее) |
величиной |
эластичности |
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
замены |
|
|
|
|
Многорежимная |
|
k |
+ |
Используется при описании процес- |
||||||||||||||||||||||||
функция, |
|
|
|
|
|
|
|
|
|
|
|
|
α0 αi |
|
сов, в которых уровень отдачи каж- |
|||||||||||||
y = (a xa0 |
+ a |
|
|
xa0 )a |
1 |
|
дой новой единицы ресурса скачко- |
|||||||||||||||||||||
21 |
|
|
||||||||||||||||||||||||||
|
11 |
|
1 |
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
образно меняется в зависимости от |
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
...(a |
xa0 |
|
+ a |
2k |
xa0 )ak |
|
|
|
соотношения |
факторов. |
Функцию |
|||||||||||||||||
1k |
1 |
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
целесообразно применять при нали- |
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
чии априорной информации о числе |
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
режимов |
|
|
|
|
Функция |
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
– |
Имеет смысл использовать в тех |
||||||||||||
линейного |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
случаях, когда выпуск продукции |
|||||||||||||
программирова- |
|
|
|
является |
результатом |
одновремен- |
||||||||||||||||||||||
ния, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ного функционирования k фиксиро- |
||||
y = min |
x1 |
|
, |
|
|
x2 |
|
+ |
|
|
ванных технологий, использующих |
|||||||||||||||||
|
|
|
|
|
|
одни и те же ресурсы |
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
a11 |
|
|
a12 |
|
|
|
|
|
|
|
|
||||||||||||||
... + min |
|
|
x |
|
|
|
, |
|
|
x |
2 |
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
ak1 |
|
|
ak 2 |
|
|
|
|
|
|
|
|
*указывается, если функция обладает свойством однородности; ** в скобках приведены значения эластичности замены по любому фактору произ-
водства.
35
П р и л о ж е н и е 2
Метод Соболя
Предполагаем, что задана математическая модель исследуемой или проектированной системы, и модель эта зависит от параметров a1,…,an .
Слова «задана математическая модель» означают, что имеются формулы (или готовые программы), позволяющие по заданному набору a1,…,an вы-
числить любые интересующие нас характеристики системы. Пространство параметров. Пространством параметров называется
n-мерное пространство, состоящее из точек A с декартовыми координатами (α1, . . . , αn). Таким oбразом, каждой точке A пространства параметров соответствует конкретный набор параметров (α1, . . . , αn) и наоборот.
Как правило, проектировщики могут указать разумные пределы изменения каждого из параметров, которые называются параметрическими ограничениями (2.1)
aminj ≤ aj ≤ amaxj , j = |
|
. |
(2.1) |
1,n |
Ограничения (2.1) выделяют в пространстве параметров параллелепипед П = {A | (2.1)}, объем которого (n-мерный объем) равен произведению
VП = (a1max − a1min )…(anmax − anmin ). |
(2.2) |
В дальнейшем нас будут интересовать только точки A , принадлежащие П: им и только им соответствуют системы, параметры которых удовлетворяют ограничениям (2.1).
Так как метод Соболя основан на зондировании параллелепипеда П
конечным числом пробных точек, то без необходимости |
расширять гра- |
ницы не рекомендуется: при этом объем возрастает, и |
для просмотра мо- |
жет потребоваться больше точек. |
|
Функциональные ограничения. Кроме параметрических ограниче-
ний обычно в условия задачи включаются функциональные |
ограничения |
cminj ≤ f j ( A) ≤ cmaxj , j = 1..t, |
(2.3) |
здесь f j ( A) – некоторые функции от параметров А =(α1, . . . , αn). Они могут
быть заданы явно. Но если, например, функционирование системы описывается дифференциальными уравнениями, то f j ( A) часто представляют со-
бой функционалы, зависящие от интегральных кривых этих уравнений. Предполагается, что все функции f j ( A) непрерывны в П.
Обозначим через G подмножество параллелепипеда П, состоящее из точек А, удовлетворяющих ограничениям ( 2.3 ):
G = {A |(2.2,2.3)}.
36
Множество G, вообще говоря, может быть любым замкнутым множеством. Единственное ограничение: объем G должен быть положительным (VG
> 0).
Можно сказать, что требование VG > 0 исключает из рассмотрения задачи с функциональными ограничениями в форме равенств, например, f(А) = с.
Идея метода Соболя состоит в представлении полученных согласно (2.1) ограничений параллелепипеда П конечным множеством дискретных точек A с набором параметров (а1, а2 … аn). Затем полученные точки A проверяются на принадлежность допустимому множеству с помощью функциональных ограничений, в случае если А удовлетворяет им, считается значение целевой функции, иначе отбрасываем ее. Среди полученных допустимых точек ищется максимальное значение целевой функции.
Для нахождения точек Аi используют следующую формулу:
|
a |
|
= amin + (amax − amin ) * q , q ≤ 1, j = |
|
, |
(2.4) |
|||
j |
1,n |
||||||||
|
|
j |
j |
j |
ij ij |
|
|||
где Qi = (qi1,…, qin ), i = |
|
– точки ЛПτ − последовательностей, которые являются |
|||||||
1, N |
|||||||||
наиболее равномерно распределенными в |
K n (единичный многогранник) |
||||||||
среди всех известных последовательностей. |
|
|
|
|
|||||
Согласно Лемме 1 из [10] точки |
Аi (полученные по 2.4) |
также будут |
образовывать равномерно распределенную последовательность в параллелепипеде П.
Рассмотрим способы получения таких последовательностей точек. В
таблице Д1[10] приведены числители |
rj(l ) направляющих чисел: |
|||||||||
|
|
|
V (l) = r(l) * 2−l , j = |
|
,i = |
|
|
(2.5) |
||
|
|
|
1,51 |
1,20 . |
||||||
|
|
|
j |
j |
|
|
|
|
|
|
Таким образом, по таблице числителей можно генерировать точки Qi с |
||||||||||
номерами 0 ≤ i ≤ 220 |
в кубе K n |
с максимальной размерностью n=51. |
||||||||
1. Номер i |
записываем в |
|
двоичной форме |
i = emem−1 …e1 , |
||||||
где m = 1 + [ln i / ln 2], [z] − целая частьz |
по таблице числителей; |
|
||||||||
2. |
Рассчитываем Vj(l). |
|
|
|
|
|
|
|
||
3. Qi = (qi1 ,…, qin ) |
считаем |
|
|
следующим |
образом: |
|||||
|
1 |
2 |
m |
|
|
|
|
|
|
|
q = eV( ) *eV( ) *…*e V( ). |
|
|
|
|
|
|
|
|||
ij 1 j |
2 |
j |
m j |
|
|
|
|
|
|
|
* – «исключающее ИЛИ».
Однако при больших n удобнее и эффективнее с точки зрения экономии машинного времени использовать способ сверхбыстрого расчета ЛПτ − последовательностей. Для этого порядок следования точекQi меня-
ется так, чтобы каждая следующая точка Qi легко вычислялась по предыдущей точке Qi'−1 и чтобы при этом двоичные участки новой последователь-
37
ности Q0' ,Q1' ,…,Qi' ,… совпадали с двоичными участками последовательностей Q1,Q2 ,…,Qi ,… . Последнее условие гарантирует сохранение всех основных и
дополнительных свойств равномерности.
Пусть Г(i) – так называемый код Грея, соответствующий номеру. По определению:
Г (i) = i *[i / 2],[ z] − целаячастьz.
Два соседних кода Грея Г(i) и Г(i-1) отличаются всегда в одном и том же разряде l = l (i), номер которого можно вычислить по формуле:
|
|
l = 1 + log |
2 |
Г(i)* Г(i − 1) . |
(2.7) |
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Поэтому, полагая Q' = Q |
Г(i) |
, получаем алгоритм для расчета Qi' : |
|
||||||||
|
|
i |
|
|
|
|
|
|
|
|
||
q' |
= … = q' |
= 0. |
|
|
|
|
|
|
|
|
|
|
01 |
0n |
|
|
|
*V (l) , |
|
|
|
(2.8) |
|||
|
|
q' |
= q' |
|
j = |
|
. |
|||||
|
|
|
1,n |
|||||||||
|
|
ij |
i−1, j |
|
|
|
j |
|
|
|
|
Для вычисления последней формулы и формулы (2.6) требуются только логические операции.
38
Учебное издание
МОДЕЛИ И МЕТОДЫ ПОДДЕРЖКИ СБАЛАНСИРОВАННОГО РАЗВИТИЯ РЕГИОНАЛЬНЫХ ЭКОНОМИЧЕСКИХ СИСТЕМ
Учебно-методическое пособие
Составители:
Баева Нина Борисовна, Сергеева Екатерина Александровна
Редактор М.С. Исаева Компьютерная верстка О.В. Шкуратько
Подписано в печать 07.09.2016. Формат 60×84/16. Уч.-изд. л. 2,5. Усл. п. л. 2,3. Тираж 30 экз. Заказ 483.
Издательский дом ВГУ. 394018 Воронеж, пл. Ленина, 10
Отпечатано в типографии Издательского дома ВГУ. 394018 Воронеж, ул. Пушкинская, 3
39