中国科学技术大学中国科学院无线光电通信重点实验室, 合肥 230027
2017年6月23日 收稿; 2017年10月23日 收修改稿
基金项目: 国家自然科学基金(61672484)资助
通信作者: 邱玲,E-mail:lqiu@ustc.edu.cn
摘要: 为减少射频链路开销的同时满足系统高容量需求,提出使用有限射频链路的毫米波混合数模预编码,以实现硬件成本和系统性能的折中。在多用户毫米波MIMO系统下,首次提出基于射频链路选择的高能效混合预编码算法。首先利用预设的模拟预编码码本将原三元耦合变量联合优化问题转变为稀疏数字预编码优化问题。由于问题非凸非线性,接着利用顺序凸近似与分数规划理论,设计一种迭代求解算法。仿真结果表明所提算法在能效上十分逼近于穷举法并显著优于其他算法。
关键词: 毫米波MIMO通信混合预编码射频链路选择能量效率
Energy-efficient hybrid precoding design based on RF chains selection for millimeter wave MIMO systems
SUN Jihan, QIU Ling
Key Laboratory of Wireless-Optical Communications of Chinese Academy of Sciences, University of Science and Technology of China, Hefei 230027, China
Abstract: In order to reduce the radio frequency (RF) chains cost and satisfy the high capacity demand for communication systems, millimeter wave hybrid precoding with limited RF chains is proposed to make a tradeoff between hardware cost and system performance. In this work, we study the energy-efficient hybrid precoding design problem based on RF chains selection for a millimeter wave multi-user multiple-input multiple-output (MIMO) system. First, we use the predesigned analog precoding codebook to reformulate the original problem as a sparse digital precoding problem. Then we solve it iteratively using successive convex approximation (SCA) method combined with fractional programming method. The simulation results show that our proposed algorithm gives the performance very close to that of the exhaustive search and better than those of the other existed algorithms.
Keywords: millimeter wave MIMO communicationshybrid precodingradio frequency chains selectionenergy efficiency
随着移动通信技术的迅猛发展,智能终端的快速普及,移动数据容量需求正以惊人的速度增加。毫米波由于具备大量未经使用的频段,被视为第五代移动通信系统中解决容量需求的关键技术之一[1]。毫米波由于波长较短,基站端可以以较小的物理阵列尺寸配置较多的天线。在传统的纯数字基带预编码方案中,每根天线都有对应的基带和射频链路结构,这些射频链路不仅造价昂贵而且功耗较大,在配置较多天线的毫米波通信系统中采用这种预编码方案是不切实际的。为解决上述问题,学术界提出在毫米波MIMO系统中采用混合数/模预编码结构[2]。混合数/模预编码在发射端将数据流经基带数字预编码处理后映射到各个射频链路上,然后通过恒模相移器调整各个射频链路上信号的相位完成模拟预编码。在此结构上,射频链路的数目远远小于天线的数目,从而降低通信系统对硬件的需求,同时也没有对系统的性能造成明显的损失[3-4]。
近年来,由于能源短缺以及温室效应的影响,通信系统的能耗问题也受到广泛关注。能量效率作为权衡系统容量和系统能耗的性能指标,成为未来无线通信研究的热点之一[5]。目前,有大量文献对微波MIMO系统下能效优化问题进行了广泛研究,如Pa Tervo等[6]在多用户MISO场景下提出一种能效最优的波束成形方案;Shi等[7]在干扰广播信道下设计一种优化能效的迭代算法。
然而毫米波通信系统下新型混合预编码结构的提出为能效优化问题带来了更多新的难点:1)模拟预编码器的恒模限制为原目标问题带来了非凸的限制;2)射频链路数对系统能效有很大影响[8],但是由于其数值与模拟预编码矩阵和数字预编码矩阵的维度直接相关,难以通过数值分析得到其最优解.虽然目前已有少部分文献关注毫米波混合预编码系统下的能效优化问题,如文献[9]在给定射频链路数目的条件下通过将毫米波混合预编码的能效优化问题转化为求解欧式距离最小的问题,利用正交匹配追踪算法得出原问题的近似最优值;文献[10]同样利用正交匹配追踪法,在遍历每个可能的射频链路数目后,获得系统能效的最优值。然而上述文献都忽略了难点2),而使用预先设定好的射频链路数目,从而降低了求解难度。这样做一方面忽略了射频链路数对系统能效的影响;另一方面当天线数目较多的时候,穷举搜索每个可能的射频链路数目将会十分耗时。
基于上述研究现状,本文在多用户毫米波MIMO系统下,提出一种基于射频链路选择的能效优化算法。由于原问题难以直接求解,首先引入一个预设的模拟预编码码本将问题等价转换为求解稀疏数字预编码[11-12],而模拟预编码则为从码本中选择出来的NRF个码字,其中NRF为最优的射频链路数目;接着由于转化后的问题仍然为一个非凸非线性问题,我们利用顺序凸近似理论结合Dinkelbach’s理论将问题变为凸问题并进行迭代求解。仿真结果表明,所提算法性能十分接近于穷举法性能,并且远远高于等增益传输(EGT)算法[13]性能。
1 系统模型1.1 信道模型考虑毫米波单小区下行链路场景,如图 1所示。该系统由K个单天线用户和一个配置Nt根天线的基站组成。基站端射频链路的数目为NRF,其取值范围为[K, Nt]。基站采用全连接的混合数\模预编码结构,包含一个NRF×K的基带数字预编码器WBB和一个Nt×NRF的由恒模移相器构成的模拟预编码器WRF。第k个用户接收到的信号可以表示为
${y_k} = \mathit{\boldsymbol{h}}_k^{\rm{H}}{\mathit{\boldsymbol{W}}_{{\rm{RF}}}}{\mathit{\boldsymbol{W}}_{{\rm{BB}}}}\mathit{\boldsymbol{s}} + {n_k}.$ | (1) |
Download: JPG larger image | |
图 1 系统模型 Fig. 1 System model 图 1 系统模型 Fig. 1 System model --> |
式中:s=[s1, s2, …, sK]T;sk~CN(0, 1)代表传输给第k个用户的信号;n~CN(0, σ2IK)是独立同分布的加性高斯白噪声,其均值为0,方差为σ2;基站到K个用户的信道为H=[h1, h2, …, hK]H,其中hk表示从基站到第k个用户的下行信道.本文采用基于角度扩展的Saleh-Valenzuela模型刻画毫米波信道[3],其目前被广泛应用于毫米波混合预编码的研究。hk可以表示为
${\mathit{\boldsymbol{h}}_k} = \sqrt {\frac{{{N_{\rm{t}}}{\rho _k}}}{{{N_{{\rm{ray}}}}}}} \sum\limits_{i = 1}^{{N_{{\rm{ray}}}}} {{\alpha _{ki}}\mathit{\boldsymbol{u}}\left( {{\varphi _i},{\theta _i}} \right)} .$ | (2) |
$\begin{array}{*{20}{c}}{\mathit{\boldsymbol{u}}\left( {{\varphi _i},{\theta _i}} \right) = \frac{1}{{\sqrt {{N_t}} }}\left[ {1, \cdots ,{{\rm{e}}^{{\rm{j}}\frac{{2{\rm{ \mathsf{ π} }}}}{\lambda }d\left( {p\sin {\varphi _i}\sin {\theta _i} + q\cos {\theta _i}} \right)}}} \right.}\\{\left. {, \cdots ,{{\rm{e}}^{{\rm{j}}\frac{{2{\rm{ \mathsf{ π} }}}}{\lambda }d\left( {\left( {\sqrt {{N_t}} - 1} \right)\sin {\varphi _i}\sin {\theta _i} + \left( {\sqrt {{N_t}} - 1} \right)\cos {\theta _i}} \right)}}} \right].}\end{array}$ | (3) |
1.2 能耗模型由于在移动通信系统当中,基站端占主要的功率消耗,因此,本文不考虑用户的功率消耗。基站的总功率消耗通常包括信号传输功耗以及电路功耗,所以毫米波通信系统的一般功耗模型[2]为
${P_{{\rm{total}}}} = \frac{1}{\varepsilon }{P_{\rm{t}}} + {N_{{\rm{RF}}}}{P_{{\rm{RF}}}} + {P_{\rm{c}}}.$ | (4) |
2 基于射频链路选择的高能效混合预编码设计2.1 问题形成上述毫米波系统模型下的能效优化问题可以建模成如下形式:
$\begin{array}{l}\mathop {\max }\limits_{{\mathit{\boldsymbol{W}}_{{\rm{RF}}}},{\mathit{\boldsymbol{W}}_{{\rm{BB}}}},{N_{{\rm{RF}}}}} \frac{{{R_{{\rm{sum}}}}}}{{{P_{{\rm{total}}}}}}\\{\rm{s}}.\;{\rm{t}}.\;{\left| {{\mathit{\boldsymbol{W}}_{{\rm{RF}}}}\left( {i,j} \right)} \right|^2} = \frac{1}{{{N_{\rm{t}}}}}\\\;\;\;\;\;{R_k} \ge {\gamma _k},\forall k = 1, \cdots ,K\\\;\;\;\;\;\sum\limits_{k = 1}^K {{{\left\| {{\mathit{\boldsymbol{W}}_{{\rm{RF}}}}{\mathit{\boldsymbol{w}}_{{\rm{BB}},k}}} \right\|}^2}} \le {P_{\max }}\\\;\;\;\;\;{N_{{\rm{RF}}}} \ge K.\end{array}$ | (5) |
${R_k} = {\log _2}\left( {1 + \frac{{{\mathit{\boldsymbol{h}}_k}{\mathit{\boldsymbol{W}}_{{\rm{RF}}}}{\mathit{\boldsymbol{w}}_{{\rm{BB}},k}}\mathit{\boldsymbol{w}}_{{\rm{BB}},k}^{\rm{H}}\mathit{\boldsymbol{W}}_{{\rm{RF}}}^{\rm{H}}\mathit{\boldsymbol{h}}_k^{\rm{H}}}}{{\sum\limits_{i = 1,i \ne k}^K {{\mathit{\boldsymbol{h}}_k}{\mathit{\boldsymbol{W}}_{{\rm{RF}}}}{\mathit{\boldsymbol{w}}_{{\rm{BB}},i}}\mathit{\boldsymbol{w}}_{{\rm{BB}},i}^{\rm{H}}\mathit{\boldsymbol{W}}_{{\rm{RF}}}^{\rm{H}}\mathit{\boldsymbol{h}}_k^{\rm{H}}} + {\sigma ^2}}}} \right).$ | (6) |
2.2 问题模型转换为了最大化系统能效,需要同时优化式(5)中的3个变量:WRF,WBB以及NRF。由于WRF和WBB的大小与NRF直接相关,并且目标问题非凸非线性,式(5)变得十分复杂而难以直接求解。尽管文献[10]通过穷举法搜索每个可能的NRF下的系统能效从而获得了最优值,但当天线数目较多时,这种方法耗时太多,复杂度太高。为避免穷举搜索并让问题变得可解,下面将对原问题作进一步转换,使原三元耦合变量优化问题变成一个只包含一元变量的稀疏数字预编码优化问题。
考虑模拟预编码矩阵WRF是从一个预设的码本中选出的NRF个码字组成的,这里用符号
$\frac{1}{{\sqrt {{N_{\rm{t}}}} }}{\left[ {1,{{\rm{e}}^{{\rm{j}}\frac{{2{\rm{ \mathsf{ π} }}}}{{{N_{\rm{t}}}}}\left( {k - 1} \right)}}, \cdots ,{{\rm{e}}^{{\rm{j}}\frac{{2{\rm{ \mathsf{ π} }}}}{{{N_{\rm{t}}}}}\left( {k - 1} \right)\left( {{N_t} - 1} \right)}}} \right]^{\rm{T}}},k = 1,2, \cdots ,{N_{\rm{t}}}.$ | (7) |
因此,模拟预编码的设计可以看成是从码本
${\mathit{\boldsymbol{W}}_{{\rm{RF}}}}{\mathit{\boldsymbol{W}}_{{\rm{BB}}}} = {{\mathit{\boldsymbol{\tilde W}}}_{{\rm{RF}}}}{{\mathit{\boldsymbol{\tilde W}}}_{{\rm{BB}}}} = {\mathscr{W}_{{\rm{RF}}}}\mathit{\boldsymbol{Q}}{{\mathit{\boldsymbol{\tilde W}}}_{{\rm{BB}}}} = {\mathscr{W}_{{\rm{RF}}}}{{\mathit{\boldsymbol{\tilde W}}}_{{\rm{BB}}}}.$ | (8) |
利用等式(8),式(5)可以等价转化为
(9) |
通过上述转换, 式(9)只包含一个未知变量,即稀疏的数字预编码矩阵
2.3 算法设计式(9)中的问题是一个经典的分式规划问题,利用Dinkelbach’s理论[16-17],通过引入参数η,将分式规划问题转化为等效的线性规划问题,从而最优单一预编码矩阵可以通过求解J(η)=0而得到。其中J(η)为
$J\left( \eta \right) = \mathop {\max }\limits_{{{\tilde W}_{{\rm{BB}}}}} \sum\limits_{k = 1}^K {{R_k}} - \eta {P_{{\rm{total}}}}.$ | (10) |
首先,由于DFT矩阵是一个酉矩阵,因此有
接下来,引入几个辅助变量,结合式(10)与式(9)中的约束条件, 子问题可以重写成如下形式:
(11) |
对于非凸约束
其次,考虑‖diag(
综合上述所有结果,给定η下的子问题(11)的求解被转化为一个凸问题,问题的数学描述如式(12)所示
$\begin{array}{l}\max \tau \\{\rm{s}}.\;{\rm{t}}.\;\sum\limits_{k = 1}^K {{{\log }_2}\left( {{\beta _k} + 1} \right)} \ge \tau + \eta \left( {\frac{1}{\varepsilon }\sum\limits_{i = 1}^{{N_t}} {{f_i}} + \sum\limits_{i = 1}^{{N_t}} {{x_i}} \times {P_{{\rm{RF}}}} + {P_{\rm{c}}}} \right)\\\;\;\;\;{\left\| {\mathit{\boldsymbol{\bar w}}_i^{\rm{T}},\frac{1}{2}\left( {{f_i} - {x_i}} \right)} \right\|^2} < \frac{1}{2}\left( {{f_i} + {x_i}} \right)\\\;\;\;\;{{\mathit{\boldsymbol{\tilde h}}}_k}{{\mathit{\boldsymbol{\tilde W}}}_{{\rm{BB}},k}} \ge \frac{{{\phi _k}}}{2}z_k^2 + \frac{1}{{2{\phi _k}}}{\beta _k},\forall k\\\;\;\;\;\sum\limits_{i = 1}^{{N_t}} {{f_i}} \le {P_{\max }};0 \le {x_i} \le 1,\forall i;\sum\limits_{i = 1}^{{N_t}} {{x_i}} \ge K\\\;\;\;\;{z_k} \ge {\left( {{\sigma ^2} + \sum\limits_{i \ne k}^K {{{\left\| {{{\mathit{\boldsymbol{\tilde h}}}_k}{{\mathit{\boldsymbol{\tilde w}}}_{{\rm{BB}},i}}} \right\|}^2}} } \right)^{1/2}}\\\;\;\;\;\frac{1}{{\sqrt {{{\bar \gamma }_k}} }}{{\mathit{\boldsymbol{\tilde h}}}_k}{{\mathit{\boldsymbol{\tilde w}}}_{{\rm{BB}},k}} \ge {\left( {{\sigma ^2} + \sum\limits_{i \ne k}^K {{{\left\| {{{\mathit{\boldsymbol{\tilde h}}}_k}{{\mathit{\boldsymbol{\tilde w}}}_{{\rm{BB}},i}}} \right\|}^2}} } \right)^{1/2}}\\\;\;\;\;{\mathop{\rm Im}\nolimits} \left( {{{\mathit{\boldsymbol{\tilde h}}}_k}{{\mathit{\boldsymbol{\tilde w}}}_{{\rm{BB}},k}}} \right) = 0.\end{array}$ | (12) |
Table 1
表 1 基于射频链路选择的毫米波高能效混合预编码算法(RFEE)流程Table 1 Proposed mmWave energy-efficient hybrid precoding algorithm based on RF chains selection
| 表 1 基于射频链路选择的毫米波高能效混合预编码算法(RFEE)流程Table 1 Proposed mmWave energy-efficient hybrid precoding algorithm based on RF chains selection |
在式(12)中,由于xi被松弛为一个[0, 1]上的连续变量,采用一个简单的匹配原则:对于xi>1-ζ,令xi=1;否则,令xi=0。这里ζ是一个十分小的数。下节的仿真结果中表明,这种匹配算法对性能的影响几乎可以忽略不计,这是因为求解得到的xi大部分都十分接近于1或者0。通过匹配后的xi,可以找到码本中被选择的码字,从而组成模拟预编码矩阵WRF。用
3 仿真结果与分析本节将对上述算法进行仿真性能验证,仿真中所用部分参数[2-10]如表 2所示。在不作特殊说明的情况下,用户数K=4,发射天线数Nt=64,最大发射功率Pmax=30 dBm。
Table 2
表 2 仿真参数Table 2 Simulation parameters
| 表 2 仿真参数Table 2 Simulation parameters |
图 2验证了算法1内层循环的收敛性,容错误差设置为gap=10-6。经观察发现,内层循环在几个迭代后就可以收敛至局部最优点。这是因为,在非凸约束条件
Fig. 2
Download: JPG larger image | |
图 2 算法1内层收敛性 Fig. 2 Convergence of inner loop of algorithm 1 图 2 算法1内层收敛性 Fig. 2 Convergence of inner loop of algorithm 1 --> |
为验证算法的有效性,从多个角度对所提算法进行比较。为评估算法1中xi映射为0或1对系统造成的性能损失,将xi映射前采用算法1(记为RFEE-Relax)得到的结果与xi映射后再次采用算法1得到的结果(记为RFEE-Mapping)进行比较。为评估算法1射频链路选择对系统能效的优化,将所提RFEE-Mapping算法与穷举法以及EGT[13, 15]算法进行对比。在EGT算法中,射频链路的数目必须与用户数目保持一致,模拟预编码矩阵WRF中元素值为
图 3比较几种算法的能效随最大传输功率Pmax的变化情况.从仿真结果可以看出,对于每种算法,系统的能效都是随着Pmax的增加,先增加后持平。比较RFEE-Relax算法以及RFEE-Mapping算法,可以看出,RFEE-Relax算法只略微优于RFEE-Mapping算法,这是因为RFEE-Relax算法得到的xi值大部分都十分接近于0、1,所以xi采用简单匹配规则对性能造成的损失几乎可以忽略不计。图 3(a)中,在较小发射天线数的情况下,对比采用穷举法得到最优射频链路数下的系统能效与采用RFEE-Mapping算法得到系统能效的性能差别。可以很直观地看出,所提RFEE-Mapping算法只稍逊于穷举法。而在图 3(b)中,由于天线数目增多,如果采用穷举法,需要重复进行
Fig. 3
Download: JPG larger image | |
图 3 能效与最大传输功率的关系 Fig. 3 Energy efficiency vs. the total transmit power 图 3 能效与最大传输功率的关系 Fig. 3 Energy efficiency vs. the total transmit power --> |
图 4展示不同发射天线数目对系统能效造成的影响。由于本文采用的天线阵列是方阵,因此图 4的横坐标取整数的平方。由仿真结果可以看出,天线数目越多,系统能效越高,但是随着天线数目的增加,系统能效增加速度趋于缓慢,并且RFEE-Mapping算法与EGT算法的性能差异也逐渐变小.这是因为虽然增加天线数目可以增加系统的频效但是也会额外增加系统的功耗,当天线数目大到一定程度的时候,天线数对频效的影响趋小而对能耗的影响趋大,从而使得能效增加趋于缓慢。
Fig. 4
Download: JPG larger image | |
图 4 系统能效与发射天线数目的关系 Fig. 4 Energy efficiency vs. the number of transmit antennas 图 4 系统能效与发射天线数目的关系 Fig. 4 Energy efficiency vs. the number of transmit antennas --> |
4 结束语在毫米波混合预编码结构下,系统能效与射频链路数目的优化问题十分具有挑战性。本文提出一种基于射频链路选择的高能效混合预编码算法.首先利用预设的模拟预编码码本,将原问题等价转换为稀疏数字预编码优化问题,使得原问题的3个耦合变量转化为1个未知变量。随后利用Dinkelbach’s理论结合顺序凸近似设计一种迭代求解算法。结果显示,本文提出的算法可以在避免穷举搜索的情况下优化射频链路数目、有效提升系统能效,其结果与穷举法所获性能十分逼近,并显著高于其他常用算法。
参考文献
[1] | Rappaport T S, Sun S, Mayzus R, et al. Millimeter wave mobile communications for 5G cellular:It will work![J]. IEEE Access, 2013, 1: 335-349. DOI:10.1109/ACCESS.2013.2260813 |
[2] | Méndez-Rial R, Rusu C, González-Prelcic N, et al. Hybrid MIMO architectures for millimeter wave communications:phase shifters or switches?[J]. IEEE Access, 2016, 4: 247-267. DOI:10.1109/ACCESS.2015.2514261 |
[3] | Ayach O E, Rajagopal S, Abu-Surra S, et al. Spatially sparse precoding in millimeter wave MIMO systems[J]. IEEE Transactions on Wireless Communications, 2013, 13(3): 1499-1513. |
[4] | Liang L, Xu W, Dong X. Low-complexity Hybrid precoding in massive multiuser MIMO systems[J]. IEEE Wireless Communications Letters, 2014, 3(6): 653-656. DOI:10.1109/LWC.2014.2363831 |
[5] | Li G Y, Xu Z, Xiong C, et al. Energy-efficient wireless communications:tutorial, survey, and open issues[J]. IEEE Wireless Communications, 2011, 18(6): 28-35. DOI:10.1109/MWC.2011.6108331 |
[6] | Pa Tervo O, Tran L N, Juntti M. Optimal energy-efficient transmit beamforming for multi-user MISO downlink[J]. IEEE Transactions on Signal Processing, 2015, 63(20): 5574-5588. DOI:10.1109/TSP.2015.2453134 |
[7] | Shi Q, Razaviyayn M, Luo Z Q, et al. An iteratively weighted MMSE approach to distributed sum-utility maximization for a MIMO interfering broadcast channel[J]. IEEE Transactions on Signal Processing, 2011, 59(9): 4331-4340. DOI:10.1109/TSP.2011.2147784 |
[8] | Garcia-Rodriguez A, Venkateswaran V, Rulikowski P, et al. Hybrid analog-digital precoding revisited under realistic RF modeling[J]. IEEE Wireless Communications Letters, 2016, 5(5): 528-531. DOI:10.1109/LWC.2016.2598777 |
[9] | Ma C, Shi J, Huang N, et al. Energy-efficient hybrid precoding for millimeter wave systems in MIMO interference channels[C]//2016 IEEE 83rd Vehicular Technology Conference (VTC Spring), Nanjing, 2016: 1-5. |
[10] | Zi R, Ge X, Thompson J, et al. Energy efficiency optimization of 5G radio frequency chain systems[J]. IEEE Journal on Selected Areas in Communications, 2016, 34(4): 758-771. DOI:10.1109/JSAC.2016.2544579 |
[11] | Shi Y, Zhang J, Letaief K B. Group sparse beamforming for green cloud-RAN[J]. IEEE Transactions on Wireless Communications, 2014, 13(5): 2809-2823. DOI:10.1109/TWC.2014.040214.131770 |
[12] | Dai B, Yu W. Sparse beamforming and user-centric clustering for downlink cloud radio access network[J]. IEEE Access, 2014, 2: 1326-1339. DOI:10.1109/ACCESS.2014.2362860 |
[13] | Liang L, Xu W, Dong X. Low-complexity hybrid precoding in massive multiuser MIMO systems[J]. IEEE Wireless Communications Letters, 2014, 3(6): 653-656. DOI:10.1109/LWC.2014.2363831 |
[14] | MacCartney G R, Samimi M K, Rappaport T S. Omnidirectional path loss models in New York city at 28 GHz and 73 GHz[C]//2014 IEEE 25th Annual International Symposium on Personal, Indoor, and Mobile Radio Communication (PIMRC), Washington DC, 2014: 227-231. http://www.researchgate.net/publication/281784747_Omnidirectional_Path_Loss_Models_in_New_York_City_at_28_GHz_and_73_GHz |
[15] | Love D J, Heath R W. Equal gain transmission in multiple-input multiple-output wireless systems[J]. IEEE Transactions on Communications, 2003, 51(7): 1102-1110. DOI:10.1109/TCOMM.2003.814195 |
[16] | Dong Y, Huang Y, Qiu L. Energy-efficient sparse beamforming for multiuser MIMO systems with nonideal power amplifiers[J]. IEEE Transactions on Vehicular Technology, 2017, 66(1): 134-145. |
[17] | Schaible S, Ibaraki T. Fractional programming[J]. European Journal of Operational Research, 1983, 12(4): 325-338. DOI:10.1016/0377-2217(83)90153-4 |
[18] | Beck A, Ben-Tal A, Tetruashvili L. A sequential parametric convex approximation method with applications to nonconvex truss topology design problems[J]. Journal of Global Optimization, 2010, 47(1): 29-51. DOI:10.1007/s10898-009-9456-5 |