Введение Проблема выбора принципа (назовем его принципом развития), определяющего направление и интенсивность развития конкретной системы («направление притяжения») является одной из наиболее важных для любых систем. Это направление притяжения может быть обозначено в виде определенной цели или - более сложно - совокупности целей. Цели могут быть формализованы до уровня определенных целевых функций. В свою очередь, процесс достижения целевых значений этих функций часто сводят к процессу контроля некоторого набора (целевых) показателей, которые и определяют конкретное содержание действий всех субъектов, связанных с данной системой. Таким образом, содержание описанной пирамиды, всех происходящих в ней процессов, а также куда (в какое состояние) придет система в результате своего развития в существенной степени определяются выбранным принципом развития. Некоторые примеры принципов развития: тяготение к идеальному состоянию, к некоторому (заданному) уровню (или уровням) развития, максимальное удаление от некоторых нежелательных состояний, обеспечение стабильности. Данная работа опирается на принцип тяготения к некоторому набору состояний. Этот принцип частично охватывает также принцип идеального состояния для случая, когда это идеальное состояние может быть описано некоторым набором уровней. Приведем примеры социально-экономических систем, опирающихся на принцип тяготения. А. Профессиональный спорт: основными целями подготовки спортсмена, обусловливающими его спортивные достижения, являются: 1) физическая подготовка; 2) волевая подготовка; 3) быстрота реакции; 4) уровень спортивного мышления; 5) спортивные параметры субъекта; 6) умение контролировать себя в экстремальных спортивных ситуациях; 7) нацеленность на высокие достижения. Каждое из перечисленных качеств имеет свои специфические методы развития и совершенствования. Таким образом, в сфере спортивной подготовки базовым является принцип тяготения к некоторым уровням по перечисленному набору из семи показателей. Б. Рыночная экономика: основной принцип, которым руководствуется хозяйствующий субъект в рыночной экономике, - получение максимальной прибыли. Однако указанная цель в процессе практической реализации в современных условиях в развитых странах дополняется совокупностью ряда требований и стимулов, среди которых выделим: 1) соблюдение правовых и нормативных требований; 2) поддержание должного имиджа; 3) стабильность государственной и социально-политической среды хозяйственной деятельности субъекта; 4) уровень криминальной опасности. Каждое из перечисленных качеств (требований и стимулов) имеет свой уровень достижения. Более того, низкий или недостаточно высокий уровень хотя бы по одному из этих качеств может явиться причиной сворачивания бизнеса. Вследствие этого рыночную экономику можно описать моделью развития, основанной на принципе тяготения, где целевыми функциями являются степень отклонения (точнее, «недостижения») определенного уровня по каждому из перечисленных пяти качеств (включая уровень прибыли). Аналогичным образом можно описать моделями притяжения сферы образовательной деятельности, медицинского обслуживания, туристической деятельности и др. Ранее одним из авторов настоящего исследования были получены дифференциальные уравнения, описывающие поведение во времени ряда характеристик указанных моделей при наличии заданных фиксированных уровней притяжения по каждому качеству (показателю) [1]. Представляет интерес анализ решений полученных уравнений и выявление оптимальных уровней по каждому показателю с учетом возможных затрат, связанных с достижением указанных уровней. Указанная задача и рассматривается в данной работе. Таким образом, целью работы является анализ решений и оптимальный выбор уровней развития на основе исследования построенной ранее математической модели развития системы, функционирующей на основе принципа тяготения к заданному набору состояний. Исследований по построению формализованных моделей функционирования систем на основе определенных принципов или целевых функций крайне мало (см. [2]). Наиболее близкими являются работы [1, 2]. Представляют интерес по теме исследования работы [3-5]. Основные предположения Приведем формализованное описание рассматриваемой в работе модели. Рассматривается система, состоящая из некоторой совокупности субъектов. Каждый субъект характеризуется определенным набором свойств и качеств Q1, Q2, …, QK, в котором количество контролируемых качеств фиксировано. Примеры формирования указанных качеств для различных сфер приведены в [1]. В модели предполагается, что среда деятельности субъекта нацелена на развитие каждого из качеств Qi. В рамках рассматриваемой модели понятие развития сводится к оценке отклонения текущего уровня качества Qi от заданного (желаемого) с помощью некоторой функции, и тогда понятие «развитие качества Qi» означает уменьшение значения целевой функции по данному качеству. Процесс изменения качеств у каждого субъекта вызывается определенным набором воздействий на него. Все множество возможных воздействий разделено на две группы: внешние воздействия (со стороны других субъектов, государственных и иных внешних структур, текущих жизненных обстоятельств) и внутренние воздействия, порождаемые непосредственно внутри субъекта на основе определенных стимулов, воззрений, рассуждений. Возможный потенциал развития качества Qi у каждого субъекта характеризуется некоторым максимальным уровнем развития данного качества. Обозначим как vi(t) уровень развития по качеству Qi у субъекта; очевидно, что для любого t > 0. Тогда предполагается выполнение следующих предположений относительно процесса изменения vi(t) в условиях наличия воздействий на субъект, суммарная величина которых в момент t равна Zi(t) [1]. Для каждого i изменение уровня развития i-го качества Qi за малый промежуток времени прямо пропорционально следующим величинам: 1. Достигнутому уровню vi(t) развития i-го качества: , или, в векторной форме, , (1) где знак указывает на пропорциональную (линейную) зависимость левой части от правой, а - константы, не зависящие от времени t. 2. Оставшимся потенциальным возможностям субъекта по каждому из качеств Qi. 3. Величине суммарного воздействия Zi(t). 4. Длительности воздействия . Анализ и обсуждение данных предположений проведены в [1]. При сделанных предположениях в [1] была получена следующая система дифференциальных уравнений: , (2) Для анализа динамики развития качеств в системе в целом в [1] была введена следующая характеристика: (,) есть число субъектов, достигших в момент t некоторого уровня развития такого, что для всех i и имеющих потенциальные возможности ( для всех i). Здесь индекс i у переменных , xi, yi, zi указывает на то, что эти переменные связаны с качеством Qi. При этом функция удовлетворяет следующему граничному условию: для любого t > 0. Предполагаются выполненными следующие условия: 1) в начальный момент t = 0 все субъекты имеют уровни развития по каждому признаку Qi, пропорциональные их полным возможностям , с некоторым универсальным коэффициентом пропорциональности , что влечет за собой равенство (3) где (;) - константы; - заданная функция, описывающая распределение (количество) «новорожденных» субъектов (хозяйствующие субъекты, спортсмены и др.) по потенциальным возможностям по качествам Qi; 2) доля субъектов, воздействие на которые имеет «эффект воздействия (КПД)», лежащий в промежутке по качеству Qi (;), равна , где и - заданная функция, описывающая распределение всех субъектов по степени воспринимаемого воздействия по качествам Qi (); 3) для каждого качества Qi задан некоторый уровень bi, на который ориентируются все субъекты; это выражается в том, что величина воздействия Zi(t) в момент t на субъект с уровнем развития vi(t) равна: , где - некоторая абсолютная константа, зависящая от суммарной средней величины внешнего воздействия и различающаяся в случаях bi > vi(t) и . Уровень bi может находиться внутри области возможных значений качества Qi, что для субъектов с уровнем качества Qi большим, чем bi, приводит к деградации качества Qi у этих субъектов. Элементы совокупности {bi; } являются управляющими параметрами модели, с помощью которых можно контролировать процесс развития качеств Qi в рассматриваемом коллективе субъектов. В [1] было получено следующее уравнение в частных производных для функции : (4) где - средняя величина воздействия по качеству Qi. Таким образом, решение уравнения в частных производных (4) при начальном условии (3) и граничном условии для любого t > 0 позволяет найти функцию , где функции {} находятся как решения системы (2). Поиск решения уравнения для Для решения системы (4) необходимо прежде всего найти первые интегралы следующей системы дифференциальных уравнений [6]: , (5) Разложим выражение в правой части равенства (5) на простейшие дроби. Для того чтобы не загромождать преобразование несущественными деталями, ниже индекс i опускается и вводится обозначение . Имеем: (6) где Ci (i = 1, 2, 3) - неизвестные константы. Для нахождения C1 умножим обе части (6) на выражение , получим откуда, полагая , выводим: (7) Аналогично, путем последовательного умножения обеих частей (6) на и , получаем: , (8) С учетом (7) и (8) соотношение (5) можно записать в следующем виде (индекс i опущен): (9) Отсюда, после интегрирования, получаем , где G0 - некоторая константа. Умножив обе части последнего соотношения на и пропотенцировав, имеем где - некоторая константа. Наконец, возвратив индекс i, окончательно получаем первые интегралы уравнения (5) (): (10) Тогда общее решение уравнения (5) запишется в виде , где - любая абсолютно непрерывная функция от K переменных. Для нахождения частного решения рассматриваемой задачи воспользуемся начальным условием 1, и, ввиду (10), выводим: (11) Поскольку функция f( ) монотонно возрастает по всем переменным zi, то последнее соотношение можно переписать в виде где в условии при записи функции Ti все остальные переменные опущены. Из последнего соотношения, после замены переменных ui = Ti(z1, …, zK), , получаем: (12) где и {, } есть решение системы уравнений , . (13) Отметим, что значения правой части (13) при zi = bi и находятся по непрерывности. Решение уравнение (13) для K > 1 пока не представляется возможным. Для K = 1 решение приведено в (1). Однако численное решение системы (13) при K< 10 15 может быть получено на основе математических программных систем MathCad, MathLab. Соотношения (11)-(13) позволяют найти искомую функцию . Для этого вначале на основе (12), (13) находится функция f( ), а затем из соотношения (11) - функция f( ). Выбор оптимальных уровней достижения Полученные выше соотношения для функции X( ) позволяют поставить задачу выбора оптимального набора уровней достижения bi, . Предположим, что известна некоторая весовая функция ( ), определяющая вес каждого из состояний системы развития. Состояние в каждый текущий момент времени t потенциально зависит от заданных уровней достижения bi, причем указанные характеристики являются управляемыми параметрами системы, а также зависит от достигнутого уровня развития X, т. е. определяется набором (bi, ; ). Таким образом, . Отметим, что функция ( ) является возрастающей (точнее, неубывающей) по переменным bi. Далее, задан некоторый регламентный промежуток времени [Tнач; Tкон], в конце которого необходимо достичь: I) либо заданного конечного состояния системы; II) либо максимальной суммарной величины уровня развития на всем промежутке - эта суммарная величина определяет потенциал развития на регламентном периоде. При этом для простоты предполагаем, что на промежутке [Tнач; Tкон] величины уровней развития bi и приложенные усилия Zi не изменяются. Тогда можно записать следующие оптимизационные постановки задач развития для случаев I и II: ;(14) (15) где максимум берется по управляемым параметрам {bi}. Постановка I представляет интерес при оперативном управлении, когда управляющие воздействия формируются на основе последовательных воздействий на малых интервал управления. Постановка II соответствует стратегической линии развития требуемых качеств и определяется выбранной политикой в развитии заданных свойств. Например, если рассматривается модель развития спорта, описанная выше в примере А, то предпочтение может отдаваться развитию наиболее одаренных в спортивном плане лиц как основных локомотивов инновационного развития - в этом случае весовая функция ( ) по переменным при больших значениях этих переменных может иметь более высокий порядок роста (вплоть до экспоненциального), чем при малых значениях. Обратное поведение функция будет иметь в случае, когда предпочтение отдается развитию прежде всего массового спорта. Далее, часто функция ( ) считается независимой на регламентном промежутке от достигнутого текущего уровня развития . Задача (15) представляет значительно больший, чем задача (14), поскольку все основные проблемы развития заданных качеств (в частности, подготовки спортсменов, обучения и воспитания) решаются в течение длительных промежутков времени, т. е. решаются на основе стратегического управления. Ввиду соотношений (11)-(13) задачи (14) и (15) относятся к задачам классической оптимизации функций от многих переменным. Точное их решение затруднительно из-за сложного и громоздкого вида функции. Однако полученные соотношения позволяют провести численное моделирование. Пример для случая одного параметра Для случая одного параметра (т. е. K=1) соотношения (11)-(13) могут быть разрешены, что приводит к соотношениям, полученным в [2]: - в случае (16) - (17) где В качестве примера рассмотрим показатель уровня образования, т. е. описывает среднее число субъектов в момент времени t, имеющих потенциал (т. е. объем максимально возможных знаний, которые они могут воспринять) и реальный уровень образования y < z. Примем, что функция , задающая долю субъектов с уровнем развития , является экспоненциальной: , где e - константа, описывающая эластичность уменьшения доли субъектов с заданным потенциалом. Далее примем, что политика образования предполагает равное обучение всех субъектов независимо от их потенциала и уровня образования, т. е. весовая функция не изменяется во времени. Тогда целевая функция (15) перепишется в виде Возьмем значения параметров равными: Tнач = 0, Tкон = 11; e = 0,01; = 0,1. На рисунке приведен график функции. Отметим, что для вычисления функции была написана программа в среде MathCad для решения задачи Б (см. (15)) с использованием (16), (17), когда . График функции при b = 4, t = b/5 Полученное решение b = 4,71. Таким образом, разработанные процедуры нахождения оптимального уровня позволяют численно находить оптимальный уровень притяжения. Заключение В работе приведено решение полученного ранее уравнения в частных производных для распределения субъектов по уровню развития их различных качеств. Поставлены две оптимизационные задачи по выбору оптимального набора уровней развития, в первой из которых достигается максимальный уровень развития в конечный момент заданного регламентного промежутка, а во второй - максимальное значение суммы уровней по всем субъектам на всем регламентном промежутке. Для случая одного показателя качества продемонстрирована возможность решения второй оптимизационной задачи. Результат решения проиллюстрирован графически. Решение поставленной в работе задачи открывает возможности для постановки оптимизационных задач, связанных с выбором оптимальных управлений, обеспечивающих наибольший суммарный объем требуемых качеств в любом из возможных аспектов, в частности на определенный момент времени; по всей или определенной части совокупности субъектов; по суммарным характеристикам.