
Ин мақола алгоритми қаблан аз ҷиҳати назариявӣ пешниҳодшудаи 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 objective | arcLP objective | Mehrotra бақия | arcLP бақия |
|---|---|---|---|---|---|---|
| Adlittle | 15 | 16 | 2.2549e+05 | 2.2549e+05 | 3.4e–08 | 3.0e–11 |
| Afiro | 9 | 9 | –464.7531 | –464.7531 | 8.0e–12 | 6.2e–13 |
| Agg | 22 | 20 | –3.5992e+07 | –3.5992e+07 | 5.2e–05 | 3.7e–06 |
| Agg2 | 20 | 21 | –2.0239e+07 | –2.0239e+07 | 5.2e–07 | 3.1e–08 |
| Agg3 | 18 | 20 | 1.0312e+07 | 1.0312e+07 | 8.8e–09 | 1.5e–08 |
| Bandm | 22 | 20 | –158.6280 | –158.6280 | 8.3e–10 | 3.6e–11 |
| Beaconfd | 11 | 11 | 3.3592e+04 | 3.3592e+04 | 1.4e–10 | 1.8e–12 |
| Blend | 14 | 14 | –30.8122 | –30.8122 | 4.9e–11 | 1.6e–12 |
| Bnl1 | 35 | 34 | 1.9776e+03 | 1.9776e+03 | 3.4e–09 | 2.9e–09 |
| Bnl2+ | 38 | 35 | 1.8112e+03 | 1.8112e+03 | 9.3e–07 | 3.5e–06 |
| Brandy | 19 | 24 | 1.5185e+03 | 1.5185e+03 | 6.2e–08 | 2.4e–06 |
| Degen2+ | 17 | 19 | –1.4352e+03 | –1.4352e+03 | 2.0e–10 | 5.9e–10 |
| Degen3* | 22 | 35 | –9.8729e+02 | –9.8729e+02 | 1.2e–09 | 8.6e–08 |
| fffff800 | 31 | 28 | 5.5568e+05 | 5.5568e+05 | 7.7e–04 | 3.7e–09 |
| Israel | 29 | 27 | –8.9665e+05 | –8.9664e+05 | 1.8e–08 | 3.4e–08 |
| Lotfi | 18 | 16 | –25.2647 | –25.2646 | 2.7e–07 | 7.8e–09 |
| Maros_r7 | 21 | 20 | 1.4972e+06 | 1.4972e+06 | 6.4e–09 | 1.7e–09 |
| Osa_07+ | 35 | 32 | 5.3578e+05 | 5.3578e+05 | 1.5e–07 | 8.4e–10 |
| Osa_14 | 37 | 42 | 1.1065e+06 | 1.1065e+06 | 3.0e–08 | 5.2e–09 |
| Osa_30 | 36 | 42 | 2.1421e+06 | 2.1421e+06 | 1.3e–08 | 1.3e–08 |
| Qap12 | 24 | 23 | 5.2289e+02 | 5.2289e+02 | 6.2e–09 | 2.9e–10 |
| Qap15+ | 44 | 28 | 1.0410e+03 | 1.0410e+03 | 1.5e–05 | 8.4e–08 |
| Qap8+ | 13 | 12 | 2.0350e+02 | 2.0350e+02 | 7.1e–09 | 6.2e–11 |
| Sc105 | 11 | 11 | –52.2021 | –52.2021 | 9.8e–11 | 2.2e–12 |
| Sc205 | 12 | 12 | –52.2021 | –52.2021 | 8.8e–11 | 4.4e–11 |
| Sc50a | 9 | 10 | –64.5751 | –64.5751 | 8.3e–08 | 8.5e–13 |
| Sc50b | 8 | 10 | –70.0000 | –70.0000 | 9.1e–07 | 3.6e–12 |
| Scagr25 | 18 | 19 | –1.4753e+07 | –1.4753e+07 | 4.6e–09 | 1.7e–08 |
| Scagr7 | 17 | 17 | –2.3314e+06 | –2.3314e+06 | 1.1e–07 | 7.0e–10 |
| Scfxm1+ | 22 | 21 | 1.8417e+04 | 1.8417e+04 | 1.6e–08 | 3.3e–05 |
| Scfxm2 | 26 | 24 | 3.6660e+04 | 3.6660e+04 | 2.6e–08 | 4.8e–05 |
| Scfxm3+ | 23 | 23 | 5.4901e+04 | 5.4901e+04 | 9.8e–08 | 1.2e–04 |
| Scrs8 | 30 | 28 | 9.0430e+02 | 9.0430e+02 | 1.8e–10 | 1.0e–10 |
| Scsd1 | 13 | 11 | 8.6666 | 8.6666 | 8.7e–14 | 3.3e–15 |
| Scsd6 | 16 | 16 | 50.5000 | 50.5000 | 8.6e–15 | 2.6e–13 |
| Scsd8 | 14 | 15 | 9.0500e+02 | 9.0500e+02 | 1.3e–10 | 2.6e–13 |
| Sctap1 | 27 | 20 | 1.4123e+03 | 1.4123e+03 | 0.0031 | 1.4e–11 |
| Sctap2 | 21 | 22 | 1.7248e+03 | 1.7248e+03 | 4.4e–07 | 1.4e–12 |
| Sctap3 | 22 | 21 | 1.4240e+03 | 1.4240e+03 | 5.9e–07 | 1.9e–12 |
| Share1b | 25 | 26 | –7.6589e+04 | –7.6589e+04 | 1.5e–06 | 1.9e–07 |
| Share2b | 15 | 15 | –4.1573e+02 | –4.1573e+02 | 7.9e–10 | 1.4e–10 |
| Ship04l | 18 | 19 | 1.7933e+06 | 1.7933e+06 | 2.9e–11 | 1.3e–10 |
| Ship04s | 20 | 19 | 1.7987e+06 | 1.7987e+06 | 4.5e–09 | 3.1e–10 |
| Ship08l | 22 | 20 | 1.9091e+06 | 1.9090e+06 | 1.0e–10 | 1.8e–11 |
| Ship08s | 20 | 19 | 1.9201e+06 | 1.9201e+06 | 4.5e–12 | 1.7e–09 |
| Ship12l | 21 | 21 | 1.4702e+06 | 1.4702e+06 | 1.0e–08 | 3.0e–10 |
| Ship12s | 19 | 21 | 1.4892e+06 | 1.4892e+06 | 2.1e–13 | 5.0e–11 |
| Stocfor1+ | 14 | 13 | –4.1132e+04 | –4.1132e+04 | 1.1e–10 | 8.6890e–11 |
| Stocfor2 | 22 | 22 | –3.9024e+04 | –3.9024e+04 | 1.6e–09 | 4.3e–09 |
| Stocfor3 | 38 | 37 | –3.9976e+04 | –3.9977e+04 | 6.4e–08 | 7.7e–08 |
| Truss | 26 | 24 | 4.5882e+05 | 4.5882e+05 | 9.5e–06 | 5.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 дар ин мақола дубора оварда нашудаанд.

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