Формирование фрактальной модели магнитоплазменного электродинамического ускорителя Текст научной статьи по специальности «Математика»

и

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

А. А. Михайлов, С. А. Базуева Южно-Российский государственный политехнический университет, Новочеркасск

Аннотация: В статье рассмотрены проблемы построения фрактальной модели магнитоплазменного электродинамического ускорителя. Моделирование такой сложной системы определяется как процесс вычисления обобщенного решения уравнения Фоккера-Планка-Колмогорова, описывающих реальные физические процессы в эредитарных (hereditarity) системах. Исследованы решения данного уравнения, получаемые после обобщения на дробные производные Римана-Лиувилля.

Ключевые слова: магнитоплазменный электродинамический ускоритель (ЭДУ), эредитарная ^^^г^) модель системы, уравнение Фоккера-Планка-Колмогорова, дробные производные Римана-Лиувилля.

Общий анализ объекта исследования

Современные системы высокоскоростного разгона на базе магнитоплазменных электродинамических ускорителей (ЭДУ) [1] представляют собой сложную динамическую многоуровневую систему [2] функционально связанных подсистем, состояния которых характеризуются большим количеством технологических параметров. Примером данных систем является магнитоплазменный электродинамический ускоритель (ЭДУ) с управляемым разгоном объекта управления (рис. 1), который состоит из импульсного источника энергии (ИИЭ), устройство предварительного ускорения (УПУ), рельсотрона (РК) с движущимся плазменным поршнем (ПП) и толкаемым им телом [3].

Рис. 1. Структура ИИУС магнитоплазменного ЭДУ с управляемым разгоном

системы (ИИУС), система автоматического управления содержит шунтирующий ключ (ШК), измерительный преобразователь скорости (ИПС), вычислительное устройство (ВУ), устройство управления (УУ) с каналом управления РК и каналом управления коммутаторами ИИЭ. Управление разгоном ПП в ЭДУ заключается в том, что по направляющим РК (рельсам) протекает электрический ток от ИИЭ через подвижную проводящую перемычку в виде ПП, который начально является частью цепи. Созданное током I, идущим через рельсы, магнитное поле В между ними, перпендикулярно току, проходящему через ПП и смежный рельс взаимодействует с током в ПП и порождает электромагнитную силу Лоренца ¥„ разгоняющую электропроводный ПП вдоль рельсов. Ограничение скорости ПП на заданном уровне осуществляется за счет прерывания в определенный момент процесса передачи энергии от ИИЭ к РК путем шунтирования входа рельсотрона с помощью ШК, срабатывающего по сигналу УУ.

Современный инженерный подход для описания состояния динамических систем обычно использует оценки для гауссовского распределения [4], что определяет возможность применения гарантированного (робастного) подхода при формировании модели движения плазмы в ЭДУ в виде системы N обыкновенных дифференциальных уравнений [1]:

^ = / (г, х (г), и, х)

ш

с начальными условиями:

х(г0) = х0 + Ах0,

V 0 / 0,ном 0 ?

где х (г)-вектор состояния, и-управление, х = Хном +Ах -вектор возмущаемых параметров, хном и Ах -векторы номинальных значений и возмущений параметров процесса ускорения, Х0ном -вектор номинальных начальных условий, Ахп -вектор их возмущений.

Однако обычно сумма случайных величин сходится не к гауссовским, а к классу устойчивых распределений «Леви-Парето» с тяжелыми хвостами (фрактальными распределениями, со степенными распределениями или паретианами), выборочные средние которых неустойчивы и малоинформативны [5]. Поэтому подмена распределений «Леви-Парето» при формировании модели динамических систем нормальным законом распределения приводит к тому, что модель в виде винеровского процесса, соответствующее нормальному закону распределения, остается непрерывной, в то время модель, например, в виде процесса Коши, соответствующая распределению Коши из класса распределений «Леви-Парето», время от времени терпит разрывы. В связи с этим приходится использовать математическую модель [5], включающую параметры (координаты) системы, определяющие ее состояние, и закон пространственно-временной эволюции состояний в общем виде, который на основе знаний о состоянии в начальный момент времени г0 в точке пространства х0 определяет состояние системы в любой момент времени г>г0. Причем 1111 в модели системы задается фрактальной структурой (в том числе и динамической) [6], что определяет переход к эредитарной (НегеШИагИу) модели системы [7], для которой при заданной топологии множества состояний часть состояний от их общего числа необратимо теряется в процессе эволюции и становятся недоступными для системы.

Целью статьи является формирование модели ЭДУ в виде эредитарной модели с «остаточной» памятью, описывающей скрытый порядок в системах данного класса.

Выбор концептуальной модели системы

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

з (0 = $ к(I - т)/(х)ах,

(1)

который преобразует входной сигнал АА) в выходной сигнал 3(1) с помощью импульсной функции (функции памяти) к(А). Если

ступенчатая функция, то выражение (1) будет интегралом первого порядка.

В настоящее время для описания функционирования динамических систем различных классов [8-10] широкое распространение получили марковские модели процессов. Система с марковским процессом (с полным

состояния хг-1 к состоянию х, которая определяется уравнением Чепмена-Колмогорова-Смолуховского (ЧКС) [11], в процессе эволюции теряет все свои состояния кроме одного с бесконечно большой плотностью. При этом к(?)=8(?)-функция Дирака, а преобразование (1) в виде

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

Для эредитарных систем [7], занимающих промежуточное место между марковскими системами с полным отсутствием памяти и простыми системами с полной памятью, функция памяти для интеграла свертки (1) интерполируется дробной функцией к(?)=ки(?) между 8-функцией (интегрирование нулевого порядка) и ступенчатой функцией (интегрирование первого порядка). Причем модель эредитарной системы

отсутствием памяти), заданная вероятностью перехода Р(х{, от

А (0 = $ щ - т) А (т^т

(2)

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

Оценка размерности множества состояний системы

Для исследования структуры множества £ состояний системы {50, £1,..., £п} (0=х0<х1<.. ,<хп) мощности /=|х| воспользуемся ступенчатой функцией

^ если х е[х, х21 (3)

П(х1 <х < Х2) ЧЛ г , (3)

[0, если х вне [ х1, х2],

лапласовский образ которой с параметром преобразования Лапласа р равен

П(х1 < х < х2) ^ — рх1]{1 - exp[-р(х2 - х1)]}, Р

а для ступенчатой функции (3) на к-й стадии имеет вид

П(хкп) < х < хк+1) ^ - exp[-рхк") ]{1 - exp[-рАи ]}, Р

„„„ Д _ у(п) _ у(п) ^ДС л к+1 Лк .

Если повторить этот процесс для (п+1)-го этапа, то получим

П(хкп) < х < хкп) + Ап+1) + п(хкп) + Ап - Ап-1 < х < хкп) + Ап) ^

^ — exp(- рх(п))[1 - exp(- рАп)] +-exp(- рхкп)) exp[-р(Ап - А„+1] х Р р

х [1 - exp(- рА п+1)] = 1 - ^рА п+1 [1 + е -р (А»exp(- рх(п))] .

р( п+1)

Используя связь хк(п)=хк+1(п+1) и вновь повторяя эту процедуру, имеем

2«+1 1 п

IП (( < х < хк++1)) - [1 - exp(- рА п+1 )]П{1 + exp[- р(Ак-1 - Ак)]}. (4)

к=1 р к=1

Воспользовавшись при определении Ап для функции Дх) (с образом Лапласа ^(р)) последовательным делением с коэффициентом деления мощности множества состояний /, имеем Ап=£,п/, а из (2) для совокупности ступенчатых функций

(29

для которого, используя (4), имеем

1 t 2 N

J(x)=Twi I < x < XN )f (x):

J(x)^(p)=l ^ p/£N ) Qn (ps(1 - £))) (p), (5)

где

N-1

- N

Qn (z) = 2-N П[1 + exp(-z£n)] с z=p/(l-£). (6)

Для относительно больших N (N>>1), p/£N|<<1 из (5) следует

0(p)=qn(z)f(p). (7)

Из выражения (6) следует

Q (z) [1 + exp(-z£N-1)] Q Qn (z) =-2-Qn-l( z),

причем QNz), определенное выражением (6), удовлетворяет также уравнению [12]:

Qn (z / £ ) = 1/2qn-!(z). (8)

При N—ю существует предел QN(z). Используя неравенство 0<|exp(-

p/£n(1-£))|<1, можно принять, что для любогоp/ (0<p/|<ro)

0<limqn(z) = Q(z) < 1.

N —

Поэтому в пределе уравнение (8) преобразуется к виду

Q(z / £ )= 1/2Q( z),

решение которого [12] имеет вид Q(z)= Avz-, откуда Av (z/£)-v = 2 Avzи после

логарифмирования имеем v=ln2/ln(l/£) (9)-размерность множества состояний [12]. Причем если параметры случайные и имеют вид £.=£+8. со

n . Ill

случайными отклонениями Ъ. от среднего £, то £ = n- , £ - £. = |5г | << £ < 1,

i

то [12] в формулах необходимо сделать замену £ — £, где

n=0

5, <5">=п 115" (¿=1, 2,..., т)-средние значения

г

множества {8г}г.

При этом минимальное значение предельного интервала состояния Аг=Ах^г при оценивании размерности А множества £ состояний {£0, 51,..., £п} не превышает заданное значение погрешности 8г- оценивания параметра

Определение условий для функции плотности вероятности перехода

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

3 (г) = |кя(г-г) / (т)Шт.

0

Данное выражение можно аналогично (7) преобразовать к виду

Ф(р)=Р^х)Е(р\ (10)

где Рд{х)- функция переходов по предшествующим состояниям.

Определим структуру переходной функции Р(г, х, 5) (г>0, 5еВ) на произвольном фазовом пространстве (Е, В), т.е. вероятность того, что объект с состоянием х, через время г попадает в множество 5, и для нее выполнены условия [13]:

а. При фиксированном г и х функция Р(г, х, 5) является мерой на а-алгебре В.

б. При фиксированном г и 5 Р(г, х, 5) есть В-измеримая функция точки х.

е. Р(г, х, Е)<1.

г. Р(0, х, Ех)=0.

д. Р("+г, х, 5)= |Р(", х, Шу)Р(г, у, 5) (¿, г>0)

Е

Переходная функция определяет стандартный процесс, т.е. стационарный процесс X с мерой на а-алгебре В (не обязательно конечная), такая что: а. если ц(5)=0, то Е\5 всюду плотно в Е;

5 =

exp

М

5

м

252

+.

Ь. для всех г>0, 5еВ Р(г, х, 5)= | Ж (г, х, у)\1(<$у), причем при г>0 функция

Г

Ж(г, х, у) непрерывна по х и ограничена и для любых х, у еЕ Ж(г, х, у)=Ж(г, у, х). Если Ж(г, х; у) - переходная плотность, то формула

¡Ж(г, х, у)йу при г > 0

Р (г, х, 5) =

(хеЕ, 5еВ)

X 5(х) при г = 0

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

Для марковского процесса £,(г, ю), геТ условная вероятность

Р(х(^)1 x(w-l), •^^(хл^хлг-^Х т.е. при фиксированном состоянии процесса в настоящий момент времени гп-1 будущее для момента времени гп не зависит от прошлого в моменты гп-2,..., г1.

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

дРШг) < х|£(г,) = х ,} , , ч

—--п=± = Ж (х\хп 1)

дх п

которая удовлетворяет условиям неотрицательности Ж(хп|хп-1)>0 и

ад

нормировки |Ж(xn\xn-1)dxn = 1. При этом цепочка состояний, описывающих

-ад

эволюцию непрерывного марковского процесса, на множестве состояний в дискретные моменты времени

^={хюе5 0,., х^еБп} (0=г0< г1<...<гп=с, 51,., ^еВ) задается через его двумерный закон распределения

P(x(1),., х(^))=Р(х(^)1 x(W-1),., x(1))P(x(1),., х(^-1))=Р(х(^)1 x(N-1))P(x(1),., х(^-1))=. =Р(х(^)1 х(^-1))Р(х(^-1)|х(^-2)). Р(х(2)1 x(1))P(x(1)),

т.е. при начальной плотности р(х^ и переходной вероятности P(xi\xi-1) можно определить многомерную плотность вероятности вектора {хь..хп} [13, 14]:

П

p(Xl,..., Xn)=p(Xl) П P(Х^ X_!) .

1=2

Исследуем свойства переходной вероятности Р(х1\х1-1) при фиксированном ? и х функции Р(1, х, 5)

Пр( х,к_1),

1=2

где 5к={хкп) < х < х(п) + Ап+1}, А1=Ах£1 - состояния, задающие эволюцию

непрерывного марковского процесса, на множестве его состояний в дискретные моменты времени. Рассмотрим один из сомножителей данного произведения и учтем, что функции Р(х{\х{-1) подобны друг другу, а поскольку р(х1)=Р(х1 |х0)р(х0), то

Р(х1)

р(х1|х0)-

Р( х0)

(11)

Функция Р(х1\х;_1) при различных значениях 1 соответствует различным уровням рассмотрения исследуемой системы и является масштабно-инвариантной, поскольку функции Р(х1\х;_1) подобна! друг другу, т. е. если Р(х1|х0)

зависит не от х1 и х0 в отдельности, а только от их безразмерной комбинации

г \ X

V х0 У

Р(Х1|Х0) =Р Используя 11 и 12, получаем

с \

х

V х0 У

(12)

Р(х1) _ Р( Х0)

Р

Г \

х

V х0 У

(13)

р(х) р(1)

Из (13) при х0=1 имеем =Р (х1), а при х1=1 получаем =Р

Р(1) Р( х0)

V Х0 У

и в результате

(х Л

Р(х,)Р(1/х0)=Р — . (13')

V х0 У

Прологарифмировав (12') и осуществив преобразование 1пР(х1)=^(г1), 21=1п х1, z0=1п х0 получим

и +и (-zo) =и (^). (14)

Продифференцировав (14) по z1 и переходя к пределу z0^z1, получаем

и (z) =и (0) =сопб11, что возможно, если только функция U(z) линейна:

и =a+bz или 1пP(x)=a+b1пx, откуда для постоянных положительных А=еа= и а=_Ь выражение для переходной вероятности

Р(—11—0) =р(—1)/р(—0) = . (15)

—1

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

Из выражения (15) следует, что

а =--V ^ ,-и =1ПР(П/1П(1/£).

1жг-1 - 1п х. 1п^-1) - 1п(^) 7

Усреднив параметр а по совместной плотности распределения двух соседних состояний системы, и учтя, что Р(хг\хг_1)=р(х.)/р(хг_1) получаем

а = V { 1п1(, хг-1)_Н(х.)}/1п2, (16)

где Н(х.)-дифференциальная энтропия .-го состояния системы, т^х., хм)-интеграция между х. и х._1 состояниями системы [15].

Применив преобразование Лапласа к (15) и учтя из (6), что z=p/(1_£), из (10) имеем

Ф (р)-Аа(1_£Г(р/)_аР(р), что соответствует представлению ) в форме дробного интеграла [12, 1618] (Приложение 1)

3 (г) =А а[г (1-^)]"а[Г (а)]-

1 г А г

1|(г - т)а-1 /(т)^т =-а-1(1 - и)а-1 /(иг0<а<1, 0<т<г.

0 (1 - ^)аГ(а) 0

Таким образом, порядок производной а определяется связностью множества состояний, которая задается размерностью топологии множества состояний V, а также энтропией Н(х1) состояния х^ системы и значением

интеграции т1:(х, х{-1) между состояниями системы х^ и хг-1.

Анализ полученных результатов

Марковская модель, определяющая переходы системы на мультифрактальных множествах ее состояний вдоль траектории разгона ПП, согласно проведенным в работе исследованиям, определяется фрактальной размерностью (9), учитывающей ее временную и пространственную память, и математическим аппаратом обобщенных дробных производных Римана-Лиувилля (Приложение 1), путем замены в уравнении Фоккера-Планка-Колмогорова (ФПК) целочисленных производных на обобщенные производные дробного порядка задаваемых параметром а по выражению (16).

Примеры решения задачи формирования уравнения ФПК, описывающих эредитарные системы, приведены в приложении 2. Эти решения сводятся к определению порядка а дифференциального оператора, который характеризует топологию множества состояний системы, обусловленной параметрами среды. При этом размерность множества состояний, определяемая долей сохранившихся состояний множества состояний, сходится к показателю диферинтегралов V (0<у<1) и охватывает случаи полностью замкнутой (у=1) и марковской (у=0) систем, когда все состояния вырождаются в одно (два) с бесконечно большой плотностью.

Переход к производной дробного порядка по времени позволяет учитывать эффекты памяти системы [16], процессы в которой

классифицируются как процессы с «остаточной» памятью [12, 16], часть которых сохраняется, а другая часть соответствует необратимым потерям. Так расщепление дифференциального оператора - а д/д 2 на

произведение диферинтегралов половинного порядка

~4а ддх)(д //д/12 ддх) интерпретируется как два эквивалентных

решения [16]. Одно из них соответствует влиянию в прямом направлении (х>0), а второе решение соответствует обратному влиянию (х<0). Для полубесконечного пространства плотность состояний для прямого процесса становится преобладающей и поэтому половина состояний теряется.

Анализ выражения (9) показывает, параметр определяющий топологию множества состояний, задает интервал регистрации контролируемого параметра через неопределенность (погрешность) его измерения и равен интервалу оценивания вероятности перехода в уравнении ФПК, т. е. интервалу управления, что согласуется с принципами цели и двойственности для системы измерения и управления [19].

Выводы

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

2. Особенности топологии множества состояний системы высокоскоростного разгона на базе магнитоплазменных ЭДУ определяют введение в уравнение ФПК, определяющее функцию плотности

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

3. Уравнение с производной половинного порядка определяется как решение, соответствующее влиянию в системе в прямом направлении (х>0), а система, которая описывается данным уравнением, содержит каналы, входящие в состав ветвящейся структуры, порождаемой внешней средой.

4. Порядок дробной производной а в уравнении ФПК определяется связностью множества состояний, которая задается его размерностью V, характеризующейся топологией множества состояний (долей сохранившихся каналов), а также энтропией Н(х ) состояния х системы и интеграцией т1;(х, хг1) между х^ и хг-1 состояниями системы.

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

Приложение 1

Общий анализ дробных операторов

О. Гёльдер [17] показал, что если функция/х) не дифференцируема, то может существовать предел /а)(х)= Шп Д//Дха, который является производной

с дробной размерностью а (показатель Гёльдера а). При этом функция /х) должна удовлетворять условию Гёльдера с показателем а на отрезке [а, Ь], для которой существуют действительное число А>0 такое, что

| ^(х1) /х2) |<А | х1х2 I а, 0<а<1, х1х2е [а, Ь].

Обобщение формулы Коши на нецелые порядки интегро-дифференциальных операторов определяет диферинтегралы дробного порядка [11, 16-18]:

/) = — и в\/(г) = 01} /(т) ¿т,

а Г(Р( а (г - т)1-р Г(П - в) йгпа (г - т)р-и+1

где в, аеЯ, и-1<Р<и, /арг -интегральный оператор порядка в, Варг -

дифференциальный оператор порядка р.

Дробный дифференциал удовлетворяет правилу композиции [18]

в1в1у(х) = в;;ру(х), (П1.1)

где В0х, В0рх и -операторы дробного дифференцирования Римана-

Лиувилляпорядка а, в и а+в, действующие на интервале от 0 дох.С помощью (П1.1) целочисленная производная представляется в виде двух дробных производных

^ = в-Дг у( х). ах

Приложение 2

Учет дробного порядка в обобщенном уравнении Фокера-Планка

-Колмогорова

Обобщенное уравнение ФПК дробного порядка получается на основе уравнения ЧКС для условной плотности вероятности Ж(у1г1 |у2г2) имеет вид [20]

Ж(У1 |у2; г+Дг)= | йуУ (у^; г)Ж(уъ\у2; Дг). (П2.1)

Для вывода обобщенного уравнения ФПК исходим из выражения

Ьу)®^, (П2.2)

-ш дг

где 0<а<1, производная дробного порядка определена соотношением [18]

дау (у| у.,г а, у(у|у„т)дт, (П2.3)

дга Г(1 - а) зг (г - т)а

где Г(г)-гамма-функция Эйлера.

Учтя соотношения (П2.2), (П2.3) уравнение ЧКС (П2.1) приводится к

виду

7б(У)дЖ^*) dy = !dydy2dt2Q(у)Ж(у2|у„*2)Ж(а)(у2|у;* - *2), (П2.4)

-ш д*

1 д * Ж(у|у • z - * )

где введена функция Ж(а)(у2у, * - *2) = —---!-.

1(1 - а) д* 0 (* - z)

Воспользуемся разложением функции Q(y) в обобщенный ряд Тейлора [21]

ш ( у - у )п+в дР+"

Q(y)= У (у у2)—-д—Q(y)\ ,

^ ¿0 Г(п + в +1) дув+п^Ау=у2'

где 0<Р<1. В результате уравнение (П2.2) можно привести к виду

+ ш даЖ ( у|у • *) ш Q (П+в)( у )

! Q(у) —д^^у = У! dy2dt2 -О—в+^АП+в (у2; * - *2)Ж (у21уп t2),

-ш да* «=0 Г( п + в +1) в

где А;+в(у2; * - *2) = !dy(y - у2)п+вЖ(а)(у2|у; * - *2) (П2.5)_обобщенные моменты.

При а=1 и в=1 (П2.5) переходит в обычное выражение для момента. Проинтегрировав (П2.5) по частям, получим дифференциальное уравнение дробного порядка

. ^ВуМ = у (-1)П Л!![(у;* -,)Ж(у|у,;о]. (П2.6)

дГ п=0 Г(п + в +1) ду Л п+в ^ У }

Уравнение (П2.6) представляет собой обобщенное уравнение ФПК дробного порядка и определяет класс стохастических процессов с обобщенными моментами

АП+в(у2; * - *2) = А;+в§(* -*2).

В результате получим уравнение

дЖ(уу1; *) = -_!_ (ж (у|у.; *))+_+_ Л1 (Ж (у|ур *)), (П2.7)

дГ Г(1 + в) дув Ух ' Г(2 + в) ду

которое при а=1 и в=1 переходит в уравнение ФПК [19].

I. Решение задачи Коши, состоящей из дифференциального уравнения

да/(х, *) - в д2/(х, *) = 0 дГ дх2 '

при 0<а<1, *>0, |х| < ш и с начальным условием к уравнению (П2.1) в виде

д / (х,0) = А( х),

*

а-1

после прямого и обратного преобразования Фурье по переменной х известно [18]

^ а 1 да да

f(x, t)=— Jdk Jdxexp(ik(x -x')A(x')Eaa(-Dk21а), (П2.8)

2 я а,а

да —

где Еа„ (-£а) = X (-1)"-- функция Миттаг-Леффлера [17].

-=0 Г(ап + в)

Для А(х)=5(х) для а=1 и E1,1(-z)=exp(-z) решение (П2.8) имеет известный вид

1

' x2 ^

/х, г)= , exp--.. (П2.9)

л/2лВг Я 4Вг) к '

В случае, например а=1/2, 1/4 функция Еа,а(--а) принимает вид

Е1/2,1/2(-^г )= п-1/2 - V- exp(z)[1-erf( V-)],

1 V— 71/4

Е1/4,1/4(-£1/4)=Г(1/4)1^1(1; 1/4; -)+^1^(1; 3/4; z)-z3/4exp(z)[1+erf(V-)]-•

II. Для решения уравнения (6) при начальном условии /(х, г=0)=/(х, 0) преобразуем его к виду [22]

да/(х,г) /(х,0) . ВдМ=0. (П2.10)

dtа Г(1 - a)tа dx2 Второе слагаемое в (П2.10) результат действия дробной производной на начальное условие. После преобразования Фурье по x и Лапласа по t имеем

да

fx, t)= J dk exp(-ikx) f (k ,0) Ea1(-Dk 2ta)

—да

Для f(x, t=0)=S(x) имеем

да

f(x, t)= J dk cos(kx )Ea1 (-Dk 2ta).. (П2.11)

0

При a=1 решения (П2.7) и (П2.9) совпадают, а для а=1/2, 1/4 функции Миттаг-Леффлера имеют вид:

Ву2 1 (-z1/2)=exp(z)[1-erf(fz)], Вш 1(-z1/4)=exp(z)[1+erf(Vi)-y(1/4, z) - y(3/4, z)^

Г(1/4) Г(3/4)

где y(a; z) - неполная Гамма-функция.

III. Решение задачи

dа f (x, t)/ dtа - Ddа f (x, t )/dxа = 0 для ограниченной области 0<x<l, с краевыми условиями fx, 0)=9(x),

— t

■ д/ (х, г) / дх + к0 / (х, г) = 0 при х = 0 и д/ (х, г) / дх + к0 / (х, г) = 0 при х = I, полученная методом разделения переменных [23] имеет вид

к

( h Л

fx, 0= jC cos(Xnx) + ^srn(Xnx) Вал(- X Dt)

n=0

где Сг

2X

1 2 i 1,2

i(X2 + ho2) + ho + h X + ho 0

j>( x)

h0

cos(X nx) + —Lsin(X nx)

X

(П2.12) dx, а корни

X2 + h2

уравнения tg(Xl)

_X(h0 + hn )

X2 - ho h

Решение (П2.12) при а=1 совпадает с известным решением, а для а=1/2, имеет вид

Ax, t)= j ся

n=0

h0

cos(X nx) + —0Lsin(X nx) nXn

exp(- X A2 [1 - erf(XA )]).

Решение уравнения (П2.11) при t^-ro имеет степенной характер [24]

exp(z)[1-erf (Vz )]=exp(z)erfc (Vz)?

1

1 + j (-1)"

1 • 3 • ...(2" -1)

"=1

(2z)"

л/п • z

и отличается от экспоненциального асимптотического решения (П2.9).

Литература

1. Azanov I. B., Aleksandrov V. A., Obydennikov S. S., Tjutin V. K., Khrustalev M. M., Judas V. I. Macroparticle launch velocity control in rail acceleration // IEEE Transactions on Magnetics. 1997. V. 33, № 1. pp. 213 - 219.

2. Анищенко В. С. Динамические системы // Соросовский образовательный журнал. 1997. №11. С. 77 - 84.

3. Кириевский Е. В. Методы и средства измерительного преобразования скорости движения плазмы для информационно-измерительных и управляющих систем электродинамических ускорителей: Автореф. дис. докт. техн. наук: 05.11.16.-Новочеркасск, 2009.-34 с.

4. Михайлов А. А. Робастные устройства контроля скорости в системах управления специализированными электрофизическими установками: Автореф. дис. канд. техн. наук: 05.13.05.-Новочеркасск, 1994.-18 с.

5. Никитин А.П., Чернавская О.Д., Чернавский Д.С. Распределение Парето в динамических системах, находящихся в шумовом поле//Труды института общей физики РАН им. А.М. Прохорова. 2009, Том 65. C. 107-123.

6. Кобелев Я. Л. Феноменологические модели описания больших систем с фрактальными структурами: Автореф. дис. канд. техн. наук: 01.04.07.-Екатеринбург, 2001.-24 с.

7. Вольтерра В. Теория функционалов и интегральных и интегро-дифференциальных уравнений. М.: Наука. 1982. 304 с.

8. Соседко В. В., Янишевская А. Г. Математическая модель единой системы конструкторско-технологической подготовки и производства на промышленном предприятии//Инженерный вестник Дона, 2012, №4 URL: ivdon .ru/magazine/archive/49N4t2y12/1404.

9. Михайлов А. А., Базуева С. А. Использование скрытой марковской модели при синтезе стохастического алгоритма решения задачи//Инженерный вестник Дона, 2015, №2, ч.2 URL: ivdon.ru/ru/magazine/archive/n2p2y2015/3073

10. Mikhaylov A. A., Bazuyeva S. A. Probabilistic Approach to the Synthesis of Algorithm for Solving Problems//Modern Applied Science/Canadian Center of Science and Education. 2015. V. 9, №5. pp. 125-132

11. Волков И. К., Зуев С. М., Цветкова Г. М. Случайные процессы. М.: Изд-во МГТУ им. Н. Э. Баумана. 1999. 448 с.

12. Нигматулин Р. И. Дробный интеграл и его физическая интерпретация //ТМФ. 1992. Т.90, №3. С. 354-368

13. Дынкин Е. Б. Марковские процессы. М.: Физматгиз. 1963. 860 с

14. Городецкий А. Я. Информационные системы. Вероятностные модели и статистические решения. СПб: Изд-во СПбГПУ. 2003. 326 c.

15. Михайлов А. А., Михайлова С. А. Оценка структурной интеграции информационно-измерительных систем//Известия высших учебных заведений. Северо-Кавказский регион. Технические науки. 2007. №6. С. 22-26.

16. Бабенко Ю. И. Тепломассообмен. Метод расчета тепловых и диффузионных потоков. Л.: Химия. 1986. 144 с.

17. Нахушев А.М. Элементы дробного исчисления и их применение. Нальчик: Изд-во КБНЦ РАН. 2000. 299 с.

18. Самко C. Г., Килбас Ф. Ф., Маричев О. И. Интегралы и производные дробного порядка и некоторые приложения. Минск: Наука и техника. 1987. 688 с.

19. Перегудов Ф. И., Тарасенко Ф.П. Введение в системный анализ. М.: Высшая школа. 1989. 367 с.

20. Гарднер К. В. Стохастические методы в естественных науках. М.: Мир. 1986. 528 с.

21. Градштейн И. С., Рыжик И. М. Таблицы интегралов, сумм, рядов и произведений. М.: Физико-математическая литература. 1963. 1108 с.

22. Кочубей А. Н. Диффузия дробного порядка//Дифференциальные уравнения. 1990. Т.26, № 4. С. 660-670.

23. Геккиева С. Х. Краевая задача для обобщенного уравнения переноса с дробной производной по времени//Доклады Адыгской АН. 1994. Т.1, № 1. С. 17-18.

24. Нахушев А. М. Уравнения математической биологии. М.: Высшая школа. 1995. 301 с.

References

1. Azanov I. B., Aleksandrov V. A., Obydennikov S. S., Tjutin V. K., Khrustalev M. M., Judas V. I. Macroparticle launch velocity control in rail acceleration // IEEE Transactions on Magnetics. 1997. V. 33, № 1 pp. 213 - 219.

2. Anishhenko V. S. Sorosovskij obrazovatel'nyj zhurnal. 1997. №11. pp. 77

- 84.

3. Kirievskij E. V. Metody i sredstva izmeritel'nogo preobrazovanija skorosti dvizhenija plazmy dlja informacionno-izmeritel'nyh i upravljajushhih sistem jelek-trodinamicheskih uskoritelej. [Methods and tools for measurement conversion speed plasma for information-measuring and control systems electrodynamic accelerators]. Novocherkassk, 2009. 34 p.

4. Mihajlov A. A. Robastnye ustrojstva kontrolja skorosti v sistemah upravlenija specializirovannymi jelektrofizicheskimi ustanovkami. [Robust device for speed monitoring in control systems specialised electrophysical installations]. Novocherkassk, 1994. 18 p.

5. Nikitin A.P., Chernavskaja O.D., Chernavskij D.S. Trudy instituta obsh-hej fiziki RAN im. A.M. Prohorova. 2009. Tom 65. pp. 107-123.

6. Kobelev Ja. L. Fenomenologicheskie modeli opisanija bol'shih sistem s fraktal'nymi strukturami. [Phenomenological models describe large systems with fractal structures]. Ekaterinburg, 2001. 24 p.

7. Vol'terra V. Teorija funkcionalov i integral'nyh i integro-differencial'nyh uravnenij. [Theory of functionals and of integral and integro-differential equations]. M.: Nauka. 1982. 304 p.

8. Sosedko V.V., Janishevskaja A.G. Inzenernyj vestnik Dona (Rus), 2012, №4 URL: ivdon.ru/magazine/archive/49N4t2y12/1404.

9. Mihajlov A. A., Bazueva S. A. Inzenernyj vestnik Dona (Rus), 2015, №2, ch.2 URL: ivdon.ru/ru/magazine/archive/n2p2y2015/3073.

10. Mikhaylov A. A., Bazuyeva S. A. Probabilistic Approach to the Synthesis of Algorithm for Solving Problems//Modern Applied Science/Canadian Center of Science and Education. 2015. V. 9, №5. pp. 125-132.

11. Volkov I. K., Zuev S. M., Cvetkova G. M. Sluchajnye processy. [Random processes]. M.: Izd-vo MGTU im. N. Je. Baumana. 1999. 448 p.

12. Nigmatulin R. I. TMF. 1992. T.90, №3. pp. 354-368

13. Dynkin E. B. Markovskie processy. [Markov processes]. M.: Fizmatgiz. 1963. 860 p.

14. Gorodeckij A. Ja. Informacionnye sistemy. Verojatnostnye modeli i sta-tisticheskie reshenija. [Information system. Probabilistic models and statistical decision]. SPb: Izd-vo SPbGPU, 2003. 326 p.

15. Mihajlov A. A., Mihajlova S. A. Izvestija vysshih uchebnyh zavedenij. Severo-Kavkazskij region. Tehnicheskie nauki. 2007. №6. pp. 22-26.

16. Babenko Ju. I. Teplomassoobmen. Metod rascheta teplovyh i dif-fuzionnyh potokov. [Heat and mass transfer. Method of calculation of heat and diffusion fluxes]. L.: Himija. 1986. 144 p.

17. Nahushev A.M. Jelementy drobnogo ischislenija i ih primenenie. [Elements of fractional calculus and their applications]. Nal'chik: Izd-vo KBNC RAN. 2000. 299 p.

18. Samko C. G., Kilbas F. F., Marichev O. I. Integraly i proizvodnye drobnogo porjadka i nekotorye prilozhenija. [Integrals and derivatives of fractional order and some applications]. Minsk: Nauka i tehnika. 1987. 688 p.

19. Peregudov F.I., Tarasenko F.P. Vvedenie v sistemnyj analiz. [Introduction to system analysis]. M.: Vysshaja shkola, 1989. 367 p.

20. Gardner K. V. Stohasticheskie metody v estestvennyh naukah. [Stochastic methods in natural Sciences]. M.: Mir. 1986. 528 p.

21. Gradshtejn I. S., Ryzhik I. M. Tablicy integralov, summ, rjadov i proiz-vedenij. [Tables of integrals, sums, series, and products]. M.: Fiziko-matematicheskoj literatury. 1963. 1108 p.

22. Kochubej A. N. Differencial'nye uravnenija. 1990. T.26, № 4. pp. 660-670.

23. Gekkieva S. H. Doklady Adygskoj AN. 1994. T.1, № 1. pp. 17-18.

24. Nahushev A. M. Uravnenija matematicheskoj biologii. [Equations of mathematical biology]. M.: Vysshaja shkola. 1995. 301 p.