Бул изилдөөдө сунушталган Эбола вирус оорусунун модели эпидемиянын динамикасын Капутонун бөлчөк тартиптеги туундулары менен көрсөтүлгөн эс-тутум таасирлери, инкубациялык кечигүү, адамдын жүрүм-турумуна байланышкан каныккан жугуу жана социалдык-маданий сөөк коюу практикалары менен бирге эске алган кечигүүсү бар компартменттик модель болуп саналат. Модель сезгич адамдарды, клиникалык жактан жуккан адамдарды, карантинге алынган адамдарды, айыккандарды жана ар түрдүү сөөк коюу шарттарындагы өлүмдүн эки классын бири-биринен бөлөт. Изилдөөнүн негизги максаты чыныгы Эбола эпидемиясын так алдын ала айтуу эмес; бул механизмдердин математикалык туруктуулукка, эпидемиялык босогого жана термелме динамикага тийгизген таасирин теориялык жана сандык жактан изилдөө болуп саналат.
Моделде негизги репродукция саны \(R_0\) негизги босоголук чоңдук болуп саналат. Булакта \(R_0<1\) болгондо оорусуз тең салмактын локалдык асимптотикалык туруктуу экени; \(R_0>1\) болгондо болсо эндемикалык тең салмак пайда болору көрсөтүлөт. Эндемикалык тең салмак үчүн инкубациялык кечигүү \(\tau\) белгилүү бир критикалык мааниге жеткенде Хопф бифуркациясы пайда болушу мүмкүн жана инфекциянын деңгээли туруктуу тең салмакка жакындагандын ордуна үзгүлтүксүз термелүүлөрдү көрсөтүшү мүмкүн.
Сандык симуляциялар ошондой эле кечигүү убактысынын, бөлчөк тартиби \(q\) нун жана көзөмөлгө алынганга чейинки өлүмдүн үлүшүн көрсөткөн \(p_1\) дин эпидемиялык ийри сызыктарды олуттуу өзгөртө аларын көрсөтөт. Бирок бул натыйжалар моделдин чыгыш натыйжалары; алар клиникалык же эпидемиологиялык талаа байкоолору катары чечмеленбеши керек. Изилдөөдө чыныгы эпидемиялык маалыматтар менен параметрлерди баалоо же көз карандысыз божомолду текшерүү жүргүзүлгөн эмес.
Булактагы сандык натыйжалардын айрым бөлүктөрүнүн ортосунда параметрлер жана чоңдуктар боюнча дал келбестиктер бар. Айрыкча, кыскача мазмун менен 3–6-сүрөттөрдүн тексттери жана аталыштары айрым \(\tau\), \(p_1\) жана инфекциянын чокулук маанилери боюнча бирдей маалымат бербейт. Ошондуктан Verianla баяны бул айырмачылыктарды бириктирип, жаңы “оңдолгон” натыйжа чыгарбайт.
Изилдөө караган негизги маселе
Эбола вирус оорусунун эпидемияларында жугуу активдүү жуккан адамдар менен болгон байланыш менен гана чектелбеши мүмкүн. Жуккан адам ооруну алып жүргөн, бирок али жугуштуу деп эсептелбеген инкубациялык мезгил, ооруга байланышкан коркуу жана маалымдуулуктан улам адамдардын жүрүм-турумун өзгөртүшү жана жуккан адамдардын сөөктөрү менен жакын байланышты камтыган айрым салттуу сөөк коюу практикалары эпидемиянын динамикасын өзгөртө алат.
Изилдөө адабиятта көргөн негизги боштук — бул таасирлердин маанилүү бөлүгү мурда өз-өзүнчө моделдештирилгенине карабастан, эс-тутум таасири + убакыт кечигүүсү + жүрүм-турумдук каныгуу + социалдык-маданий сөөк коюу жолдору бир эле бөлчөк тартиптеги алкакта бирге жетиштүү деңгээлде каралбагандыгы.
Бул жердеги “эс-тутум” адамдын биологиялык эс-тутумун билдирбейт. Математикалык жактан бул системанын учурдагы өзгөрүү ылдамдыгы ошол учурдагы абалдан гана эмес, өткөн абалдардын салмакталган суммасынан да таасир ала турганын билдирет. Бөлчөк туундулардын классикалык биринчи тартиптеги туундулардан негизги айырмасы ушунда.
Капутонун бөлчөк туундусу бул Эбола моделинде эмнени билдирет?
Капутонун бөлчөк туундусу моделдин бүгүнкү эпидемиялык динамикасынын өткөн абалдарга көз карандылыгын математикалык түрдө көрсөткөн эс-тутумдуу туунду оператору болуп саналат; \(q=1\) классикалык биринчи тартиптеги моделге жакындаса, \(0<q<1\) өткөн абалдар моделге ташылып келинген бөлчөк динамиканы билдирет.
Капутонун бөлчөк туундусунун аныктамасы
Формула эмнени эсептейт? \(f(t)\) функциясынын \(q\) тартибиндеги Капутонун бөлчөк туундусун эсептейт.
- \(q\): бөлчөк туундунун тартиби. Изилдөөнүн эпидемиялык моделинде \(0<q\le1\) аралыгы колдонулат.
- \(\Gamma(\cdot)\): гамма-функция; факториал түшүнүгүнүн чыныгы жана комплекстүү сандарга кеңейтилген түрү.
- \(\xi\): өткөн убакыт боюнча интеграл алууда колдонулуучу жардамчы убакыт өзгөрмөсү.
- \((t-\xi)^{q+1-n}\): өткөн абалдардын бүгүнкү туундуга кандай салмак менен салым кошорун аныктоочу ядронун бөлүгү.
- \(f^{(n)}(\xi)\): функциянын \(n\)-классикалык туундусу.
Бөлчөк туундунун интегралды камтышы өзгөчө маанилүү. Классикалык дифференциалдык теңдемеде өзгөрүү көбүнчө системанын ошол учурдагы абалына көз каранды. Бул жерде болсо интеграл \(a\) дан \(t\) ге чейинки өткөн убакытты камтыгандыктан, система өзүнүн тарыхын математикалык түрдө алып жүрөт.
Изилдөөнүн түшүндүрмөсүнө ылайык, \(q\) маанисинин 1 ге жакындашы эс-тутум таасиринин азайышына туура келет. Сандык симуляцияларда \(q\) нун 1 ге жакындашы инфекциянын чокулук маанилерин жогорулатып, чоку убакыттарын өзгөртө ала турганы байкалган. Бул натыйжа колдонулган модель жана параметрлер үчүн гана жарактуу.
Моделдеги алты эпидемиологиялык класс
| Белги | Биологиялык мааниси |
|---|---|
| \(S_H(t)\) | \(t\) убакыттагы сезгич адамдар |
| \(I_H(t)\) | Клиникалык жактан жуккан адамдар |
| \(Q_H(t)\) | Аныкталган жана карантинге алынган жуккан адамдар |
| \(R_H(t)\) | Айыккан адамдар |
| \(D_1(t)\) | Салттуу же диний практикалардын алкагында дене менен байланыш болушу мүмкүн болгон шарттарда коюлган каза болгон адамдар |
| \(D_2(t)\) | Саламаттык сактоо органдары көзөмөлдөгөн жана моделде жаңы жугуу жаратпайт деп кабыл алынган коопсуз сөөк коюу классы |
Негизги параметрлер
| Параметр | Мааниси |
|---|---|
| \(A\) | Популяцияга кошулуу / төрөлүү ылдамдыгы |
| \(\beta_1\) | Жуккан адамдардан сезгич адамдарга жугуу коэффициенти |
| \(\beta_2\) | \(D_1\) классындагы каза болгон адамдардан сезгичтерге жугуу коэффициенти |
| \(\mu\) | Табигый өлүм ылдамдыгы |
| \(\sigma\) | Жүрүм-турумдук/каныгуу таасирин көрсөткөн параметр |
| \(\tau\) | Таасирге кабылуу менен жугуштуулуктун башталышынын ортосундагы инкубациялык кечигүү |
| \(\alpha\) | Жуккан класстан чыгуу ылдамдыгы; \(1/\alpha\) орточо жугуштуу мезгилди билдирет |
| \(p_1\) | Аныкталганга же табигый айыкканга чейин оорудан каза болгон жуккан адамдардын үлүшү |
| \(p_2\) | Аныкталып, карантинге алынган жуккан адамдардын үлүшү |
| \(\epsilon_1\) | \(D_1\) классында өлүмдөн сөөк коюу процессине өтүү ылдамдыгы |
| \(\epsilon_2\) | \(D_2\) классында өлүмдөн сөөк коюу процессине өтүү ылдамдыгы |
| \(\delta\) | Карантинден чыгуу ылдамдыгы; \(1/\delta\) орточо карантин мезгилин билдирет |
| \(f\) | Карантинге алынган адамдар арасындагы ооруга байланыштуу өлүмдүн үлүшү |
| \(q\) | Бөлчөк туундунун тартиби / эс-тутум таасиринин параметри |
Алты теңдемеден турган негизги Эбола модели
Булакта модель (3.1) төмөнкүдөй берилет:
Биринчи теңдеме сезгич адамдардын өзгөрүшүн сүрөттөйт. \(A\) жаңы адамдарды кошот; жуккан адамдар жана \(D_1\) классындагы жугуштуу сөөктөр менен байланыш сезгич класстан адамдарды чыгарат; табигый өлүм да \(\mu S_H\) өлчөмүндө азайууну жаратат.
Экинчи теңдеме жаңы клиникалык инфекциялардын пайда болушун сүрөттөйт. Ошол эле эки жугуу каналы бул жолу оң мүчө болуп саналат; анткени сезгич класстан чыккан жаңы инфекциялар \(I_H\) классына кирет. \((\alpha+\mu)I_H\) болсо оору процесси же табигый өлүмдүн айынан жуккан класстан чыккан адамдарды билдирет.
\(D_1\), \(Q_H\), \(D_2\) жана \(R_H\) теңдемелери жуккан адамдардын өлүм, аныктоо-карантин жана айыгуу жолдору боюнча бөлүнүшүн көзөмөлдөйт. Булактын маанилүү божомолдорунун бири — \(D_2\) классындагы көзөмөлдөнгөн сөөк коюулар жаңы жугууга алып келбейт.
Каныккан жугуу мүчөсү эмне үчүн колдонулат?
Бөлчөктүн алымындагы \(\beta S\times I\) же \(\beta S\times D_1\) түзүлүшү мүмкүн болгон байланыштын көлөмүн көрсөтөт. Бөлүмүндөгү \(1+\sigma S_H(t)\) болсо популяция көбөйгөн сайын байланыштын чексиз сызыктуу өсүшүнө жол бербейт. Булак муну маалымдуулук, коркуу, коргоочу жүрүм-турум, чектелген байланыш мүмкүнчүлүктөрү жана коомдук саламаттык сактоо чаралары жаратышы мүмкүн болгон каныгуу таасири менен байланыштырат.
\(\sigma\to0\) болгондо бөлүм болжол менен 1 ге барабар болуп, жүрүм-турумдук каныгуу алсырайт. Булак бул абалды коргоочу жүрүм-турум азыраак натыйжалуу болгон моделдөө чеги катары түшүндүрөт.
Эмне үчүн моделди үч теңдемеге чейин кыскартууга болот?
Авторлор \(S_H\), \(I_H\) жана \(D_1\) теңдемелери \(Q_H\), \(D_2\) жана \(R_H\) теңдемелерине көз каранды эмес экенин белгилешет. Ошондуктан негизги динамика жана туруктуулук анализи үчүн алгачкы үч теңдемеден турган подсистема колдонулат:
Бул кыскартуу \(Q_H\), \(D_2\) жана \(R_H\) класстары биологиялык жактан маанилүү эмес дегенди билдирбейт. Бул үч класс кайтарым байланыш аркылуу алгачкы үч дифференциалдык теңдеменин оң жагына кирбейт дегенди гана билдирет.
Оң маанилүүлүк жана чектелгендик эмне үчүн текшерилет?
Эпидемиологиялык модель адамдардын терс санын чыгарбоого тийиш. Ошондуктан авторлор модель оң баштапкы шарттардан башталганда анын чечимдери терс аймакка өтпөй турганын жана чектелген бойдон каларын көрсөтүшөт.
Бул туюнтма кыскартылган системанын биологиялык жактан уруксат берилген абалдар мейкиндигин аныктайт. \(A/\mu\) жогорку чеги популяцияга кирүү менен табигый өлүмдүн ортосундагы тең салмак масштабынан келип чыгат.
Бул жердеги \(E_q\) — Миттаг–Леффлер функциясы. Ал бөлчөк дифференциалдык теңдемелерде классикалык экспоненциалдык функция аткарган ролго окшош ролду аткарат. Бул теңсиздик баштапкы маани ылайыктуу аймакта болсо, жалпы өзгөрмө \(A/\mu\) чегинен сыртка көзөмөлсүз өсүп кетпей турганын көрсөтүү үчүн колдонулат.
Негизги репродукция саны R0 бул моделде кандай чечмеленет?
\(R_0\) — толугу менен сезгич популяцияда бир жуккан адам өзүнүн жугуштуу мезгили ичинде түз жана салттуу сөөк коюуга байланышкан кыйыр жолдор аркылуу орточо канча жаңы инфекция жарата аларын көрсөткөн босоголук параметр; \(R_0<1\) жугуунун басаңдашы, ал эми \(R_0>1\) моделде эндемикалык тең салмактын болушу жагынан критикалык чек болуп саналат.
Бул формула эки негизги салымдан турат.
- Түз жуккан адамдын салымы: \(\beta_1\) камтылган биринчи мүчө.
- Салттуу сөөк коюу классынан келип чыккан салым: \(\beta_2\), \(p_1\) жана \(\epsilon_1\) камтылган экинчи мүчө.
\(\beta_1\) же \(\beta_2\) көбөйгөндө башка параметрлер туруктуу болсо, \(R_0\) жогорулайт. \(\sigma\) бөлүмдө болушу күчтүүрөөк каныгуу/жүрүм-турум таасири жугуу босогосун азайтуу багытында математикалык салым кошорун көрсөтөт. \(\epsilon_1\) экинчи мүчөнүн бөлүмүндө болушу болсо \(D_1\) классында өткөрүлгөн убакыт сөөк аркылуу жугуу жолунда маанилүү экенин көрсөтөт.
Оорусуз жана эндемикалык тең салмак
Оорусуз тең салмак:
Бул чекитте сезгич популяция бар, бирок \(I_H\) жана \(D_1\) нөлгө барабар.
Булакта \(R_0>1\) болгондо оң эндемикалык тең салмак \(\mathcal{E}_1=(S_H^\ast,I_H^\ast,D_1^\ast)\) аныкталат. Булактагы (3.3) теңдеме төмөнкү түзүлүштө:
\((R_0-1)\) көбөйткүчүнүн \(I_H^\ast\) жана \(D_1^\ast\) туюнтмаларында көрүнүшү босоголук логиканы ачык көрсөтөт: эндемикалык тең салмак үчүн модель \(R_0>1\) шартын талап кылат.
Локалдык туруктуулук жана мүнөздөмөлүк теңдеме
Тең салмак чекиттеринин туруктуулугун изилдөө үчүн сызыктуу түргө келтирилген системанын мүнөздөмөлүк тамырлары каралат. Оорусуз тең салмак үчүн булак төмөнкү мүнөздөмөлүк түзүлүштү берет:
Кечигүү болгондо коэффициенттердин ичинде \(e^{-s\tau}\) мүчөлөрү орун алат. Бул экспоненциалдык кечигүү мүчөсү мүнөздөмөлүк теңдемени классикалык полином болуудан чыгарып, трансценденттик түзүлүшкө айландырат. Ошондуктан кечигүү тамырлардын комплекстүү тегиздиктеги жайгашуусун жана демек туруктуулукту өзгөртө алат.
Булакта көңүл бурулушу керек болгон математикалык эскертүү: (3.9) теңдеменин астында берилген \(a_0\) туюнтмасынын удаалаш жазылган эки теңдигинин ортосунда масштаб/өлчөм жагынан ачык кайра уюштуруу кадамы көрсөтүлгөн эмес. Verianla бул туюнтманы булактын ордуна кайра чыгарып “оңдобойт”; туруктуулукту түшүндүрүүдө авторлор берген \(R_0\) босогосу жана теоремалар негиз катары алынат.
Кечигүү эмне үчүн Хопф бифуркациясына алып келиши мүмкүн?
Кечигүү \(\tau\) мүнөздөмөлүк теңдемеде \(e^{-s\tau}\) көбөйткүчтөрүн пайда кылып, айрым өздүк маанилердин комплекстүү тегиздиктеги жоромол окту кесип өтүшүнө мүмкүндүк бере алат; критикалык \(\tau^\ast\) маанисинде бул өтүү болгондо туруктуу тең салмак туруктуулугун жоготуп, үзгүлтүксүз термелүүгө өтүшү мүмкүн жана булак муну Хопф бифуркациясы катары талдайт.
Анализде \(s=i\omega\) деп алынып, мүнөздөмөлүк теңдеме чыныгы жана жоромол бөлүктөргө ажыратылат. Андан кийин кечигүүнүн талапкер маанилери:
түрүндө берилет жана критикалык кечигүү:
деп аныкталат.
Бул жердеги \(\omega\) термелүүнүн мүнөздүү бурчтук жыштыгына байланышкан оң маани. \(G_1(\omega)\) жана \(G_2(\omega)\) — мүнөздөмөлүк теңдеменин чыныгы жана жоромол компоненттеринен алынган функциялар. Максат — комплекстүү тамырлардын жуптары жоромол огуна жеткен эң кичине кечигүүнү табуу.
Булактагы Теорема 3.3, \(q\in(0,1]\), \(R_0>1\), \(\tau>0\) жана берилген трансверсалдуулук шарты аткарылганда эндемикалык тең салмак \(0\le\tau<\tau^\ast\) аралыгында асимптотикалык туруктуу экенин жана \(\tau=\tau^\ast\) чекитинде Хопф бифуркациясы пайда болорун билдирет.
Routh–Hurwitz критерийинин ролу
Кечигүүсүз учурда \(\tau=0\) деп алынып, \(s^q=\lambda\) алмаштыруусу жүргүзүлөт. Ошентип, мүнөздөмөлүк теңдемени классикалык полином формасына келтирүүгө болот. Эндемикалык тең салмак үчүн булак:
полиномун алат жана
Routh–Hurwitz шарттарын колдонот. Бул шарттар мүнөздөмөлүк тамырлардын чыныгы бөлүктөрүн терс бойдон сактоого жана ошону менен локалдык асимптотикалык туруктуулукту камсыз кылууга багытталган.
Lyapunov ыкмасы жана глобалдык эндемикалык туруктуулук
Булак \(R_0>1\) үчүн эндемикалык тең салмактын белгилүү кошумча шарттардагы глобалдык асимптотикалык туруктуулугун изилдөө максатында Lyapunov функционалын түзөт:
Биринчи компонент:
Ушуга окшош логарифмдик түзүлүштөр \(I_H\) жана \(D_1\) үчүн да түзүлөт; кечигип келген абалдарды эске алуу үчүн Риман–Лиувилл бөлчөк интегралдары кошулат.
Далилдөөнүн негизги математикалык куралы:
теңсиздиги. Теңдик \(x=1\) болгондо гана аткарылат. Авторлор ылайыктуу \(H_1-H_4\) монотондуулук шарттарында \(D^qV\le0\) натыйжасына жетип, эндемикалык тең салмак үчүн глобалдык асимптотикалык туруктуулук жыйынтыгын алышат.
Изилдөөнүн ыкмасы жана жыйынтыктары
Сандык чечим үчүн баштапкы шарттар
Сандык симуляцияларда кичинекей жамаат болжолдонуп, баштапкы маанилер:
деп алынган.
Adams–Bashforth–Moulton predictor–corrector ыкмасы
Бөлчөк тартиптеги кечигүүсү бар дифференциалдык система түздөн-түз жабык түрдө чечилбегендиктен, авторлор Adams–Bashforth–Moulton (ABM) predictor–corrector ыкмасын колдонушат.
Жалпы бөлчөк баштапкы маани маселеси:
Бул жерде \(y(t)\) — абал вектору, \(f\) — системанын оң жагы, ал эми \(\psi(t)\) — кечигүү аралыгындагы өткөн баштапкы функция.
Булак жалгыз чечим үчүн Lipschitz шартын колдонот:
Андан кийин бөлчөк дифференциалдык теңдеме ага эквиваленттүү Volterra интегралдык теңдемесине өзгөртүлөт:
Бул өзгөртүү маанилүү; анткени ABM алгоритми бөлчөк эс-тутум таасирин интегралдык салмактар аркылуу эсептейт.
Убакыт тору үчүн булак \(h=\tau/m\) кадам өлчөмүн колдонот жана чечимди \(y_h(t_j)\approx y(t_j)\) түрүндө дискреттештирет.
Predictor кадамы
Predictor кийинки убакыт чекитиндеги чечим үчүн алгачкы болжолду чыгарат.
Corrector кадамы
Андан соң трапециялык продукт квадратурасы колдонулуп, оңдолгон маани эсептелет. Булактагы негизги corrector түзүлүшү:
түрүндө. Predictor аркылуу алынган болжол оң жактагы белгисиз \(y_h(t_{n+1})\) ордуна колдонулуп, corrector кадамы аяктайт.
Симуляцияларда колдонулган негизги параметрлер топтому
| Параметр | Маани |
|---|---|
| \(A\) | 189 |
| \(\mu\) | 0.000044 |
| \(\sigma\) | 0.09 |
| \(f\) | 0.1 |
| \(p_2\) | 0.3 |
| \(\epsilon_1\) | 0.5 |
| \(\epsilon_2\) | 0.5 |
| \(\alpha\) | 0.21 |
| \(\delta\) | 0.135 |
| \(q\) | \((0,1]\) аралыгында өзгөртүлгөн |
Булак \(\beta_1\), \(\beta_2\) жана \(p_1\) параметрлерин тиешелүү симуляцияга жараша \((0,1)\) аралыгынан тандайт. Параметрлердин маанилүү бөлүгү мурдагы изилдөөдөн алынган; бул изилдөө аларды жаңы талаа маалыматтары менен кайра баалабайт.
R0 < 1: оорусуз тең салмакка жакындашуу
2-сүрөттөгү сценарийде \(\tau=0\), \(\beta_1=0.00985\), \(\beta_2=0.0075\) тандалып, \(R_0=0.5717\) алынган. \(q\) 0.6 дан 1 ге чейин 0.1 кадам менен өзгөртүлгөн. Симуляциялар чечим траекториялары оорусуз тең салмакка жакындай турганын көрсөтөт.
Бул табылганын мааниси моделдин ичинде: \(R_0<1\) тандалган параметрлер топтомунда инфекция туруктуу уланбайт. Бул жыйынтык чыныгы коом үчүн автоматтык түрдө “Эбола жок болот” деген божомол эмес; чыныгы \(R_0\) жана башка параметрлер өзүнчө өлчөнүшү керек.
R0 > 1: эндемикалык тең салмак
Булак тексти бул тажрыйба үчүн \(q=0.8\), \(\tau=0\) деп айтса да, 3-сүрөттүн аталышында \(\tau=1\) деп жазылган. Сүрөттүн аталышы ошондой эле \(\beta_1=0.00985\), \(\beta_2=0.00755\), \(\sigma=0.09\) жана \(R_0=1.33\) маанилерин берет. Булак ичиндеги бул дал келбестиктен улам \(\tau\) мааниси Verianla тарабынан тандалып, оңдолгон эмес.
Жалпы илимий билдирүү: \(R_0>1\) шартында жана тиешелүү аз кечигүү сценарийинде модель траекториялары эндемикалык тең салмакка жакындай алат жана узак мөөнөттүү, үзгүлтүксүз чоң термелүүлөрдү көрсөтпөшү мүмкүн.
Критикалык кечигүү жана уланган термелүүлөр
Сандык натыйжалардын тексти \(q=0.8\) жана \(R_0>1\) шартында \(\tau=2.5\) үчүн болжол менен 1000 күндүк симуляция жүргүзүлгөнүн билдирет. Ошол эле беттеги 4-сүрөттүн аталышында болсо \(\tau=2\) деп жазылган. Текстке ылайык, жуккандардын саны баштапкы 300 учурдан болжол менен 47 күндүн ичинде болжол менен 45.873 ке чейин көтөрүлүп, андан кийин болжол менен 650 деңгээлине түшүп, болжол менен 600–700 аралыгында уланып турган термелүүлөр пайда болгон.
Авторлор бул жүрүм-турумду критикалык кечигүү босогосунун ашылышы жана Хопф бифуркациясы менен байланыштырышат. Бирок \(\tau=2\) менен \(\tau=2.5\) ортосундагы булак ичиндеги аталыш/текст дал келбестигинен улам бул сандык натыйжа “так гана \(\tau=2.5\) тин натыйжасы” деп кайра белгиленбеши керек.
τ = 2.5 жана τ = 3 салыштыруусу
5-сүрөт \(q=0.8\), \(\beta_1=0.00985\), \(\beta_2=0.00755\), \(\sigma=0.09\) жана \(R_0=1.33\) шартында \(\tau=2.5\) менен \(\tau=3.0\) сценарийлерин салыштырат.
Булак текстине ылайык, \(\tau=3.0\) болгондо 300 жуккан адам менен башталган сценарий болжол менен 61-күнү болжол менен 101.173 жуккан адамга жетет. Андан кийин болжол менен 960 деңгээлине төмөндөп, уланып жаткан термелүүлөрдө болжол менен 65.845 чокулук маанилери пайда болот.
Авторлор бул салыштырууда кечигүүнүн өсүшүн чокунун чоңдугун көбөйтүү менен гана түшүндүрүүгө болбой турганын; айрыкча эпидемия чокусунун убактысын кечиктирерин белгилешет.
Социалдык-маданий сөөк коюу практикаларын көрсөткөн p1 таасири
\(p_1\) жуккан адамдардын аныкталганга же айыкканга чейин оорудан каза болуп, \(D_1\) жолуна өтүү үлүшүн билдирет. Модель жагынан \(p_1\) дин өсүшү салттуу сөөк коюу классына көбүрөөк адам өтүшүнө алып келиши мүмкүн жана бул класс \(\beta_2\) аркылуу жаңы инфекция жаратуучу механизмге байланышкандыктан, жугуу динамикасын күчөтүшү мүмкүн.
6-сүрөттүн графиктик легендасында \(p_1=0.1\) жана \(p_1=0.25\) көрүнсө, сүрөттүн аталышында кошумча түрдө \(p_1=0.35\) үчүн \(R_0=1.3614\) туюнтмасы берилген. Кыскача мазмун болсо \(p_1>0.4\) үчүн эпидемиянын көлөмү %35 көбөйгөнүн билдирет. Булактагы бул үч башка берүү бир мааниге түшүрүлгөн эмес.
Ишенимдүү жалпы жыйынтык мындай: моделде көзөмөлсүз/эрте өлүмдүн жогорку үлүшү жана ага байланышкан байланыштуу сөөк коюу жолу чоңураак же көбүрөөк кечиккен эпидемиялык толкундар менен байланыштуу. Бул белгилүү бир чыныгы коомдо салттуу сөөк коюу практикалары так ошол эле сандык өсүштү жаратат дегенди далилдебейт.
Бөлчөк тартиби q нун таасири
7-сүрөттө \(\tau=3\) туруктуу кармалып, \(q\) өзгөртүлгөн. Графиктерде \(q=0.8\), \(q=0.9\) жана \(q=1.0\) сценарийлери салыштырылат. \(q\) нун 1 ге жакындашы менен инфекция чокуларынын чоңойгону жана алардын убактысы өзгөргөнү көрүнөт.
Булак муну эс-тутум таасиринин алсыроосу менен байланыштырат. Кыскача мазмун \(q\) нун 1 ге жакындашы инфекциянын чокулук деңгээлдеринде болжол менен %20 өсүш жаратарын билдирет; графиктердин чийки сандык чекиттери берилбегендиктен, Verianla бул пайызды кайра эсептебейт.
Булакта берилген сандык сценарийлерди бирге окуу
| Сценарий | Булакта берилген шарт | Негизги билдирүү | Түшүндүрмөнүн чеги |
|---|---|---|---|
| 2-сүрөт | \(\tau=0\), \(R_0=0.5717\), \(q=0.6\)–1 | Оорусуз тең салмакка жакындашуу | Симуляциянын натыйжасы |
| 3-сүрөт | Текст: \(\tau=0\); аталыш: \(\tau=1\), \(R_0=1.33\) | Эндемикалык тең салмакка жакындашуу | Булакта \(\tau\) дал келбестиги бар |
| 4-сүрөт | Текст: \(\tau=2.5\); аталыш: \(\tau=2\) | Уланган термелүүлөр / Хопф түшүндүрмөсү | Булакта \(\tau\) дал келбестиги бар |
| 5-сүрөт | \(\tau=2.5\) жана 3, \(q=0.8\), \(R_0=1.33\) | Узагыраак кечигүүдө чоку убактысы жана термелүүлөр өзгөрөт | Чийки убакыт катары берилген эмес |
| 6-сүрөт | \(p_1\) өзгөртүлгөн | Көзөмөлсүз өлүм/байланыштуу сөөк коюу жолунун таасири өсөт | Легенда, сүрөт аталышы жана кыскача мазмун ортосунда \(p_1\) мааниси дал келбейт |
| 7-сүрөт | \(\tau=3\), \(q=0.8,0.9,1.0\) | \(q\to1\) болгондо жогорураак жана убактысы өзгөргөн толкундар | Модель натыйжасы; талаа өлчөмү эмес |
Изилдөө колдогон жыйынтыктар
- Бөлчөк тартиптеги жана кечигүүсү бар Эбола моделин математикалык түрдө түзүүгө болорун.
- Моделдин чечимдери ылайыктуу баштапкы шарттарда оң жана чектелген бойдон кала аларын.
- \(R_0\) негизги босоголук параметр экенин.
- Булактагы Теорема 3.2 чегинде \(R_0<1\) үчүн оорусуз тең салмактын локалдык асимптотикалык туруктуулугун.
- \(R_0>1\) үчүн оң эндемикалык тең салмактын бар экенин.
- Ылайыктуу шарттарда кечигүүнүн критикалык мааниден ашышы Хопф бифуркациясын жана термелме динамиканы жарата аларын.
- Моделдин ичинде \(p_1\), \(q\) жана \(\tau\) өзгөрүүлөрү эпидемиялык траекторияларды өзгөртөрүн.
- Белгилүү \(H_1-H_4\) шарттарында эндемикалык тең салмактын Lyapunov ыкмасы менен глобалдык асимптотикалык туруктуулугун.
Изилдөө колдобогон жыйынтыктар
- Бул модель чыныгы өлкөдөгү келечектеги Эбола учурларынын санын ырасталган түрдө божомолдобойт.
- 101.173, 45.873 же башка симуляциялык маанилер чыныгы эпидемияда байкалган учурлардын саны эмес.
- Моделдеги \(\tau^\ast\) босогосу ар бир Эбола эпидемиясы үчүн универсалдуу биологиялык туруктуу эмес.
- Модель маданий сөөк коюу практикалары бардык коомдордо бирдей сандык таасир жаратат деп көрсөтпөйт.
- Изилдөө эч кандай дарынын, вакцинанын же клиникалык дарылоонун натыйжалуулугун текшербейт.
- Симуляция натыйжалары чыныгы эпидемиялык маалыматтар менен калибрленген клиникалык же коомдук саламаттык сактоо божомолу эмес.
Моделдин чектөөлөрү жана келечектеги изилдөөлөр
Авторлор биринчи маанилүү чектөө катары чыныгы эпидемиялык маалыматтардын чектелүү болушунан улам параметрлер колдонуудагы адабияттан алынганын белгилешет. Чыныгы эпидемиялык маалыматтар менен эмпирикалык калибрлөө моделдин белгилүү эпидемияларга колдонулушун жогорулата алат.
Экинчиден, модель мейкиндиктик гетерогендүүлүктү жана адамдардын мобилдүүлүгүн ачык түрдө камтыбайт. Келечекте географиялык айырмачылыктарды, саякат кыймылдарын жана аймактар аралык жугууну кошуу сунушталат.
Үчүнчүдөн, модель детерминисттик. Авторлор жүрүм-турумда, айлана-чөйрөдө жана чакан популяцияларда пайда болгон кокустукту көрсөтүү үчүн келечекте стохастикалык компоненттерди кошууну сунушташат.
Бул чектөөлөр өзгөчө маанилүү: математикалык жактан татаалыраак модель автоматтык түрдө чыныгы эпидемияларды такыраак божомолдогон модель дегенди билдирбейт. Чыныгы дүйнөдөгү тактык үчүн параметрлерди калибрлөө, көз карандысыз маалымат менен текшерүү жана контекстке жараша баалоо талап кылынат.
Булак жана ыкма жөнүндө эскертүү
Түпнуска изилдөө: A delayed fractional-order model of Ebola virus disease with human behavioral dynamics
Авторлор: Eva Lusekelo; Mlyashimbi Helikumi; Salamida Daudi; Steady Mushayabasa.
Жооптуу автор: Steady Mushayabasa.
Тең салым: Булакта тең биринчи авторлук же тең салым жөнүндө билдирүү аныкталган эмес.
Мекемелер: Department of Mathematics and Statistics, College of Natural and Mathematical Sciences, University of Dodoma, Dodoma, Tanzania; Mbeya University of Science and Technology, Department of Mathematics and Statistics, College of Science and Technical Education, Mbeya, Tanzania; Department of Education and Mathematics, National Institute of Transport, Dar-es-Salaam, Tanzania; Department of Mathematics and Computational Sciences, University of Zimbabwe, Mount Pleasant, Harare, Zimbabwe.
Журнал: Results in Control and Optimization.
Басмакана: Elsevier B.V.
Платформа: ScienceDirect.
Том / макаланын номери: Volume 23 (2026), Article 100738.
DOI: 10.1016/j.rico.2026.100738.
Алынган датасы: 14 март 2026.
Кайра кароо датасы: 7 май 2026.
Кабыл алынган датасы: 20 май 2026.
Онлайн жарыяланган датасы: 28 май 2026.
Булактын түрү: Рецензияланган илимий журнал макаласы; теориялык математикалык эпидемиология жана сандык симуляция изилдөөсү.
Лицензия: Creative Commons Attribution-NonCommercial 4.0 International (CC BY-NC 4.0).
Илимий ыкма: Капутонун бөлчөк тартиптеги кечигүүсү бар дифференциалдык теңдемелери; каныккан инциденттүүлүк функциясы; оң маанилүүлүк жана чектелгендик анализи; оорусуз жана эндемикалык тең салмак; негизги репродукция саны; мүнөздөмөлүк теңдеме жана өздүк маанилер анализи; Routh–Hurwitz критерийи; Хопф бифуркациясы; Lyapunov функционалы; Adams–Bashforth–Moulton predictor–corrector сандык чечими.
Маалыматтардын жеткиликтүүлүгү: Булак изилдөөдө маалымат колдонулбаганын ачык билдирет. Сандык натыйжалар математикалык модель жана адабияттан алынган параметрлер менен түзүлгөн.
Каржылоо: Авторлор изилдөө мамлекеттик, коммерциялык же коммерциялык эмес каржылоо уюмунан атайын грант албаганын билдиришет.
Кызыкчылыктардын кагылышы: Авторлор изилдөөгө таасир этиши мүмкүн болгон белгилүү каржылык кызыкчылык же жеке мамиле жок экенин билдиришет.
Жасалма интеллект жөнүндө билдирүү: Авторлор эмгекти даярдоодо тилди жана окулушун жакшыртуу максатында жасалма интеллект колдонушканын, андан кийин мазмунду карап чыгып түзөтүшкөнүн жана жарыянын мазмуну үчүн толук жоопкерчиликти өздөрүнө алышканын билдиришет.
CRediT / автордук салымдар: Булакта өзүнчө CRediT authorship contribution statement аныкталган эмес; ошондуктан авторлордун милдеттери Verianla тарабынан божомолдонгон эмес.
Булак ичиндеги маанилүү сандык дал келбестик: Кыскача мазмун, 3–6-сүрөттөрдүн аталыштары жана тиешелүү жыйынтык абзацтарында айрым \(\tau\), \(p_1\) жана инфекция чокуларынын маанилери бирдей эмес. Айрыкча 3-сүрөт үчүн текст \(\tau=0\), сүрөт аталышы \(\tau=1\); 4-сүрөт үчүн текст \(\tau=2.5\), сүрөт аталышы \(\tau=2\) деп берет. Кыскача мазмундагы “1200 дөн 10.000 ден жогору” деген билдирүү да негизги текст/сүрөттөрдөгү айрым чоңураак чокулук маанилер менен түздөн-түз дал келбейт. Бул айырмачылыктар унчукпай оңдолгон эмес.
Туруктуулук билдирүүсүнө байланыштуу булак ичиндеги чек: Теорема 3.2, \(R_0<1\) үчүн оорусуз тең салмакты “локалдык асимптотикалык туруктуу” деп билдирет. Жыйынтык бөлүмү болсо ошол эле абалды “глобалдык асимптотикалык туруктуу” деп мүнөздөйт.
Пикир калтырыңыз
E-mail дарегиңиз жарыяланбайт. Милдеттүү талаалар * менен белгиленген