Тадқиқоти академӣ, забони фаҳмо

Verianla | Тадқиқоти академӣ ва илм ба забони тоҷикӣ

27 сентябр 2026, якшанбе
VERİANLAНашри мустақили илмӣ
Кушодан ё бастани меню
...
Саҳифаи асосӣ / Илмҳои амалӣ / MATLAB / ArcLP: татбиқи Matlab-и алгоритми нуқтаи дохилии infeasible бо ҷустуҷӯ рӯи камон ва ҳади мураккабии O(√nL) барои барномасозии хаттӣ
MATLAB

ArcLP: татбиқи Matlab-и алгоритми нуқтаи дохилии infeasible бо ҷустуҷӯ рӯи камон ва ҳади мураккабии O(√nL) барои барномасозии хаттӣ

Yaguang Yang'ın çalışması, doğrusal programlama için daha önce teorik olarak geliştirilen ve polinom iterasyon sınırı O(√nL) olarak verilen yay-arama tabanlı uygunsuz iç-nokta algoritmasını kullanılabilir bir Matlab çözücüsüne dönüştürmektedir.

26/08/2026  Veri Anla 26 боздид
ArcLP: татбиқи Matlab-и алгоритми нуқтаи дохилии infeasible бо ҷустуҷӯ рӯи камон ва ҳади мураккабии O(√nL) барои барномасозии хаттӣ

Ин мақола алгоритми қаблан аз ҷиҳати назариявӣ пешниҳодшудаи interior-point бо ҷустуҷӯи камонӣ барои барномасозии хаттиро ба татбиқи амалии Matlab табдил медиҳад. Барои алгоритм ҳади полиномии O(√nL) зикр шудааст. arcLP масъалаи стандартии \(\min c^Tx\) бо шартҳои \(Ax=b,\;x\ge0\)-ро ҳал карда, пешкоркард, интихоби нуқтаи ибтидоӣ, системаҳои парокандаи Cholesky, кунҷи қадами камонӣ, параметри марказонӣ, коркарди degeneracy, паскоркард ва санҷиши оптималиро якҷо мекунад.

Дар санҷишҳои Netlib барнома бо татбиқи Matlab-и усули Mehrotra predictor-corrector дар шароити якхела муқоиса шудааст. Аз рӯи ҳисоби мустақими 51 масъалаи Ҷадвали 1, arcLP дар 24 ҳолат итератсияи камтар, дар 11 ҳолат баробар ва дар 16 ҳолат итератсияи бештар истифода мекунад. Маҷмӯи итератсияҳо 1084 барои arcLP ва 1099 барои Mehrotra аст. Қимати бақияи constraint дар 34 масъала барои arcLP камтар, дар 1 масъала баробар ва дар 16 масъала бештар мебошад.

Маҳдудияти асосӣ ин аст, ки мақолаи соли 2026 исботи пурраи ҳади O(√nL)-ро дубора пешниҳод намекунад; он ба кори пешина такя мекунад. Саҳми мустақими ин матн бештар муҳандисии нармафзор ва санҷиши ҳисоббарорӣ мебошад.

Масъалаи primal ва dual

\[ \min c^Tx,\qquad Ax=b,\quad x\ge0. \tag{1} \]

\[ \max b^T\lambda,\qquad A^T\lambda+s=c,\quad s\ge0. \tag{2} \]

Алгоритм аз нуқтаи дохилии \((x^0,s^0)>0\) оғоз карда, мусбатии \(x^k\) ва \(s^k\)-ро нигоҳ медорад.

Интихоби нуқтаи ибтидоӣ

\[ \max\left\{ \|Ax^0-b\|, \|A^T\lambda^0+s^0-c\|, \frac{(x^0)^Ts^0}{n} \right\}. \tag{3} \]

Аз ду номзад нуқтае интихоб мешавад, ки ин меъёр барои он хурдтар аст.

Диапазони элементҳои матритса

\[ \frac{\max |A_{i,j}|}{\min\{|A_{k,l}|:A_{k,l}\neq0\}}. \tag{4} \]

Муаллиф мегӯяд, ки scaling-и умумӣ дар санҷишҳои пешина самаранокии ҳисоббарориро пайваста беҳтар накардааст ва аз ин рӯ ба татбиқи асосӣ дохил нашудааст. Аммо ин нисбат барои интихоби яке аз қоидаҳои pre-processing истифода мешавад.

Системаи парокандаи Cholesky

\[ AD^2A^Tu=L\Lambda L^Tu=v. \tag{5} \]

Ҳалли чунин системаҳо қисми асосии арзиши ҳисоббарории алгоритмро ташкил медиҳад.

Кунҷи камонӣ ва параметри марказонӣ

arcLP ба ҷойи ҳаракати хаттӣ ба самти оптимум аз траекторияи камонӣ истифода мекунад. Кунҷи қадам аналитикӣ ҳисоб шуда метавонад ва \(\sigma_k\) бо варианти махсуси Golden Section интихоб мешавад.

\[ \alpha_k=\min\{0.9999\alpha_k,\;0.99\pi/2\}<0.99\pi/2. \]

Муайянкунии infeasible ва unbounded

\[ \max_i|x_i|>10^{10} \]

ҳамчун аломати unbounded ва

\[ \max_i|\lambda_i|>10^{10} \]

ҳамчун яке аз аломатҳои infeasible истифода мешавад. Инҳо дар худи мақола ҳамчун эвристикаҳое тавсиф мешаванд, ки дар санҷишҳои ададӣ муайянкунии қобили қабул нишон додаанд.

Меъёри қатъ

\[ \frac{\|r_b^k\|}{\max\{1,\|b\|\}}+ \frac{\|r_c^k\|}{\max\{1,\|c\|\}}+ \frac{\mu_k}{\max\{1,|c^Tx^k|,|b^T\lambda^k|\}} <10^{-8}, \]

\[ r_b^k=Ax^k-b,\qquad r_c^k=A^T\lambda^k+s^k-c. \]

Натиҷаҳо чӣ мегӯянд?

Натиҷаҳо нишон медиҳанд, ки arcLP solver-и коркунандаи Matlab мебошад ва дар маҷмӯи санҷидашудаи Netlib ҳалли оптималӣ ба даст овардааст. Он аз рӯи итератсия ва бақия бо Mehrotra рақобатпазир аст, аммо дар ҳамаи масъалаҳо бартарӣ надорад.

Аҳамияти эҳтимолӣ барои Тоҷикистон

Тадқиқот бо маълумоти Тоҷикистон анҷом нашудааст. Бо вуҷуди ин, барномасозии хаттӣ дар банақшагирии энергетика, нақлиёт, логистика, истеҳсол, тақсимоти захираҳо ва таҳқиқоти амалиётӣ истифода мешавад. arcLP метавонад барои омӯзиш ва муқоисаи алгоритмҳо дар муассисаҳои илмӣ ва донишгоҳии Тоҷикистон муфид бошад, вале татбиқи амалӣ бояд дар шароити маҳаллии ҳисоббарорӣ ва системаҳои воқеӣ санҷида шавад.

Усул ва Натиҷаҳои Таҳқиқот

[x,obj,kk,infe,lambda,s,exflag]=arcLP(A,b,c,d,tol,iter)

exflagМаъно
0Анҷоми муваффақ
1Тибқи API-и манбаъ instance-и infeasible
2Тибқи API-и манбаъ instance-и unbounded
3Ҳам primal ва ҳам dual infeasible
4Маълумоти A, b ё c нопурра аст

Мисоли оддӣ

\[ \min x_1,\qquad x_1+x_2=5,\quad x_1,x_2\ge0. \]

Барои \(A=[1\;1]\), \(b=5\), \(c=[1\;0]^T\), arcLP пас аз панҷ итератсия \(x^*=[0\;5]^T\)-ро медиҳад.

Ҷадвали пурраи benchmark

МасъалаMehrotra ит.arcLP ит.Mehrotra objectivearcLP objectiveMehrotra бақияarcLP бақия
Adlittle15162.2549e+052.2549e+053.4e–083.0e–11
Afiro99–464.7531–464.75318.0e–126.2e–13
Agg2220–3.5992e+07–3.5992e+075.2e–053.7e–06
Agg22021–2.0239e+07–2.0239e+075.2e–073.1e–08
Agg318201.0312e+071.0312e+078.8e–091.5e–08
Bandm2220–158.6280–158.62808.3e–103.6e–11
Beaconfd11113.3592e+043.3592e+041.4e–101.8e–12
Blend1414–30.8122–30.81224.9e–111.6e–12
Bnl135341.9776e+031.9776e+033.4e–092.9e–09
Bnl2+38351.8112e+031.8112e+039.3e–073.5e–06
Brandy19241.5185e+031.5185e+036.2e–082.4e–06
Degen2+1719–1.4352e+03–1.4352e+032.0e–105.9e–10
Degen3*2235–9.8729e+02–9.8729e+021.2e–098.6e–08
fffff80031285.5568e+055.5568e+057.7e–043.7e–09
Israel2927–8.9665e+05–8.9664e+051.8e–083.4e–08
Lotfi1816–25.2647–25.26462.7e–077.8e–09
Maros_r721201.4972e+061.4972e+066.4e–091.7e–09
Osa_07+35325.3578e+055.3578e+051.5e–078.4e–10
Osa_1437421.1065e+061.1065e+063.0e–085.2e–09
Osa_3036422.1421e+062.1421e+061.3e–081.3e–08
Qap1224235.2289e+025.2289e+026.2e–092.9e–10
Qap15+44281.0410e+031.0410e+031.5e–058.4e–08
Qap8+13122.0350e+022.0350e+027.1e–096.2e–11
Sc1051111–52.2021–52.20219.8e–112.2e–12
Sc2051212–52.2021–52.20218.8e–114.4e–11
Sc50a910–64.5751–64.57518.3e–088.5e–13
Sc50b810–70.0000–70.00009.1e–073.6e–12
Scagr251819–1.4753e+07–1.4753e+074.6e–091.7e–08
Scagr71717–2.3314e+06–2.3314e+061.1e–077.0e–10
Scfxm1+22211.8417e+041.8417e+041.6e–083.3e–05
Scfxm226243.6660e+043.6660e+042.6e–084.8e–05
Scfxm3+23235.4901e+045.4901e+049.8e–081.2e–04
Scrs830289.0430e+029.0430e+021.8e–101.0e–10
Scsd113118.66668.66668.7e–143.3e–15
Scsd6161650.500050.50008.6e–152.6e–13
Scsd814159.0500e+029.0500e+021.3e–102.6e–13
Sctap127201.4123e+031.4123e+030.00311.4e–11
Sctap221221.7248e+031.7248e+034.4e–071.4e–12
Sctap322211.4240e+031.4240e+035.9e–071.9e–12
Share1b2526–7.6589e+04–7.6589e+041.5e–061.9e–07
Share2b1515–4.1573e+02–4.1573e+027.9e–101.4e–10
Ship04l18191.7933e+061.7933e+062.9e–111.3e–10
Ship04s20191.7987e+061.7987e+064.5e–093.1e–10
Ship08l22201.9091e+061.9090e+061.0e–101.8e–11
Ship08s20191.9201e+061.9201e+064.5e–121.7e–09
Ship12l21211.4702e+061.4702e+061.0e–083.0e–10
Ship12s19211.4892e+061.4892e+062.1e–135.0e–11
Stocfor1+1413–4.1132e+04–4.1132e+041.1e–108.6890e–11
Stocfor22222–3.9024e+04–3.9024e+041.6e–094.3e–09
Stocfor33837–3.9976e+04–3.9977e+046.4e–087.7e–08
Truss26244.5882e+054.5882e+059.5e–065.2e–07

Verianla Live: ҷараёни arcLP

Ҷараён танҳо марҳилаҳои дар манбаъ тавсифшударо нишон медиҳад.

МарҳилаШарҳSource
1. ВурудA, b, c ва параметрҳо қабул мешаванд.Main function
2. Pre-processПешкоркард ва санҷиши ибтидоии статус.Pre-process
3. Нуқтаи оғозНомзади беҳтар бо формулаи (3) интихоб мешавад.Equation 3
4. Arc-searchИтератсия рӯи камон бо x ва s мусбат.Arc-search
5. Системаи парокандаСистемаи Cholesky ҳал мешавад.Equation 5
6. ҚатъБақия, μ, қадам ва шартҳои статус санҷида мешаванд.Termination criteria
7. Post-processСанҷиши ниҳоӣ ва баргардондани натиҷаҳо.Main function
 

Ёддошт оид ба Манбаъ ва Усул

Муаллиф: Yaguang Yang, Independent Researcher, US.

DOI: 10.5334/jors.674.

Маҷалла: Journal of Open Research Software 14(1):57, Ubiquity Press.

Нашр: 12 августи 2026.

Статус: Software Metapaper-и аз peer review гузашта.

Литсензияи мақола: CC BY 4.0. Литсензияи барнома: BSD 3-Clause.

Версия: 1.0; санаи барнома дар манбаъ “08/02/2026” аст ва бинобар норавшании формат тафсир карда намешавад.

Маблағгузорӣ: маҷалла аз ҳаққи нашр даст кашидааст; гранти алоҳидаи тадқиқотӣ нишон дода нашудааст.

Бархӯрди манфиатҳо: муаллиф набудани онро эълом кардааст.

Маҳдудият: исботи муфассали O(√nL) ва performance profile дар ин мақола дубора оварда нашудаанд.


Мубодила:

Шарҳҳо пас аз баррасӣ нашр мешаванд.Шарҳи шумо ба раванди тасдиқ фиристода шуда, пас аз пазируфта шудан намоён мегардад.

Шарҳ гузоред

Нишонии почтаи электронии шумо нашр намешавад. Майдонҳои ҳатмӣ бо * нишон дода шудаанд

Иҷозат додан ба кукиҳо таҷрибаи шуморо дар ин сомона беҳтар мекунад. Сиёсати кукиҳо