CN114337751B - 一种时间反转ofdm多用户通信系统的功率分配方法 - Google Patents

Abstract

本发明公开了一种时间反转OFDM多用户通信系统的功率分配方法。针对时间反转OFDM多用户通信系统计算系统和速率，以系统和速率最大化为目标进行用户信号发送功率的优化分配方法。首先给出在发送总功率约束条件下多个用户功率联合分配的数学模型，该模型中的多目标联合优化是非凸优化，然后给出一种改进的迭代注水功率分配算法来联合优化所有用户的发送功率。相比较平均功率分配方法，本发明给出的方法能获得更高的系统和速率。

Description

一种时间反转OFDM多用户通信系统的功率分配方法

技术领域

本发明涉及信息通信领域，具体是一种时间反转OFDM多用户通信系统中用户信号发送功率的分配方法。

背景技术

时间反转(Time Reversal,TR)可以利用丰富的多径无线电传播环境来产生空间-时间共振效应，即所谓的聚焦效应，从而在减少干扰的同时提高接收信号强度。在TR通信过程中，接收端先向发送端发送导频脉冲，接收端根据收到的信号估计信道的脉冲响应(Channel Impulse Response,CIR)；然后发送端对估计的信道脉冲响应做时间反转和共轭处理得到时间反转波形，最后将其与发送信号进行卷积之后经过信道发送给接收端。根据信道互异性，TR本质上是把多径信道作为匹配滤波器，文献[Wang Beibei,Wu Yongle,HanFeng,et al.Green wireless communications:a time-reversal paradigm[J].IEEEJournal on Selected Areas in Communications,2011,29(8):1698-1710.doi:10.1109/JSAC.2011.110918.]证明了由于CIR的固有特性，经过多径传输后信号的能量聚焦于特定的时域和空域，即TR可以利用多径传播从周围环境中以极低的复杂度收集信号的能量。

在无线通信系统中，符号间干扰(Inter Symbol Interference,ISI)和用户间干扰(Inter User Interference,IUI)恶化了每个用户的服务质量(Quality of Service,QoS)。由于TR具有时间聚焦特性，减少了符号间的能量泄露，因此可以显著降低ISI的影响，这样可以使得接收机结构简单且不需要高复杂度的均衡器。除此之外，TR的空间聚焦效应抑制将大部分信号能量聚焦在预期用户处，这意味着可以提高接收机的信干噪比(Signalto Interference plus Noise Ratio,SINR)。文献[Emami M,Vu M,Hansen J,etal.Matched filtering with rate back-off for low complexity communications invery large delay spread channels[C]//Conference Record of the Thirty-EighthAsilomar Conference on Signals,Systems and Computers,2004.Pacific Grove,CA,USA:IEEE Press 2004:218-222.doi:10.1109/ACSSC.2004.1399123.]证明TR系统中，可以通过提高速率回退因子(Rate Back-off Factor,BOF)来改善接收信干噪比(Signal toInterference plus Noise Ratio,SINR)，BOF定义为采样速率与符号速率的比值。这里采样速率即离散信号的样值速率，符号速率为在信道上承载信息的符号的传输速率。最大样值速率由奈奎斯特准则决定，通常一个样值对应一个符号，此时符号速率与样值速率相同。若不是每个样值都对应一个符号，则符号速率就低于样值速率，即BOF大于1的情况。BOF大于1时，符号速率低于带宽限制下的最高速率，频带利用率降低，但可减轻符号间干扰。

文献[Lei Weijia,Yao Li.Performance analysis of time reversalcommunication systems[J].IEEE Communications Letters,2019,23(4):680-683.]推导了TR系统中接收端信噪比(Signal-to-noise ratio，SNR)的概率密度函数和累积分布函数，分析了遍历容量、中断概率和BPSK调制下误码率，并说明了ISI随着上采样因子的增加而减少。文献[Viteri-Mera C A,Teixeira F L.Equalized time reversal beamformingfor frequency-selective indoor MISO channels[J].IEEE Access,2017,5:3944-3957.doi:10.1109/ACCESS.2017.2682160.]推导了时间反转系统中的ISI功率的近似闭式表达式，分析了时间反转波束成形技术的空间聚焦特性和时间压缩性能，仿真结果显示用均衡时间反转技术消除ISI，系统可以获得更好的误比特率性能。在多用户系统中，由于不同用户的路径信道特性不同，时间反转传输可以利用多径作为区分不同用户的方式，进而可以将时间反转引入多址接入系统。文献[Han F,Yang Y,Wang B,et al.Time-reversaldivision multiple access over multipath channels[J].IEEE Transactions onCommunications,2012,60(7):1953-1965.doi:10.1109/TCOMM.2012.051012.110531.]研究了基于时间反转的多用户通信系统，提出了一种新的无线信道接入方式即时间反转时分多址接入(Time Reversal Division Multiple Access,TRDMA)，在大时延扩展信道下的多用户下行系统中使用TR结构，不同用户的信号通过多径区分，通过时间反转技术，利用多径信道实现下行多用户传输，可以减少ISI和IUI，提高SINR，并研究了TRDMA系统的SINR、可实现系统速率以及可实现的中断速率等系统性能指标。文献[Nguyen H T,Kovcs I Z,EggersP C F.A time reversal transmission approach for multiuser UWB communications[J].IEEE Transactions on Antennas and Propagation,2006,54:3216-3224.doi:10.1109/TAP.2006.883959.]通过实测的方式评估了多用户时间反转通信系统的性能，实验结果表明，时间反转传输的空间聚焦效应可以减小系统的BER，并降低对信号发射功率的要求。文献[Chen Y,Yang Y,Han F,et al.Time-reversal wideband communications[J].IEEE Signal Processing Letters,2013,20(12):1219-1222.doi:10.1109/LSP.2013.2285467.]研究了TRDMA系统的可达速率和计算复杂度，结果表明，当带宽足够大时，TRDMA系统可以获得比OFDM系统更高的可达速率，并且计算复杂度更低。文献[Han Y,Chen Y,Wang B,et,al.Time-reversal massive multipath effect:a single-antenna“massive MIMO”solution[J].IEEE Transactions on Communications,2016,64(8):3382-3394.doi:10.1109/TCOMM.2016.2584051.]进一步对TRDMA系统进行研究，结果表明采用时间反转技术可以将无线信道中的多径分量转换为虚拟的发送天线，从而在发送端仅配备单根天线的情况下，获得类似于MIMO传输的特性。文献[Yang Y H,Wang B,Lin,etal.Near-optimal waveform design for sum rate optimization in time-reversalmultiuser downlink systems[J].IEEE Transactions on Communications,2013,12(1):346-357.doi:10.1109/TWC.2012.120312.120572.]研究了时域TR多用户下行系统中使和速率最大的功率分配和波形设计优化问题，联合优化传输波形和功率分配，可以获得更好的系统性能。文献[Yang Y H,Liu K J R.Waveform design with interference pre-cancellation beyond time-reversal systems[J].IEEE Transactions on WirelessCommunications,2016,15(5):3643-3654.doi:10.1109/TWC.2016.2524526.]研究了时域TR多用户下行系统中使接收信号的总均方根误差最小的功波形设计优化问题，提出了一种带干扰预消的波形设计，仿真结果表明该方法比传统的波形设计方法有显著的性能改善。

TR传输可以在时域实现，也可以在频域实现，目前大多数关于TR的研究都是针对在时域实现的TR系统。文献[Dubois T,Helard M,Crussiere M,et al.Performance oftime reversal precoding technique for MISO-OFDM systems[J].Eurasip Journal onWireless Communications&Networking,2013,2013(1):1-16.doi:10.1109/VTCSpring.2015.7146002.]证明了通过在频域对信号进行预处理，也可以实现与时域TR预处理相同的效果。正交频分复用(Orthogonal Frequency Division Multiplexing,OFDM)是第4代、第5代移动通信系统中的关键技术，也将是新一代移动通信的关键技术，将TR技术应用到OFDM系统中具有明确的实际意义，而在频域中进行TR预处理更适合于OFDM系统。频域TR预处理过程是各子信道传输信号与子信道频率响应的共轭的乘积，运算复杂度低，根据傅里叶变换性质，这个预处理与在时域上将发送信号与TR预处理滤波器的脉冲响应进行卷积是等效的。文献[Nguyen T,Monfared S,Determe J,et al.Performance analysis offrequency domain precoding time-reversal MISO OFDM systems[J].IEEECommunications Letters,2019,24(1):48-51.doi:10.1109/LCOMM.2019.2949556.]推导了频域TR预处理系统接收信号的均方根误差的闭式表达式，并分析了BOF变化时均方根误差的变化情况，理论分析和仿真结果表明，与时域TR系统类似，增大BOF也能降低信号均方根误差，提高系统的信噪比(Signal to Noise Ratio,SNR)。文献[Golstein S,Nguyen T,Horlin F,et al.Physical layer security in frequency-domain time-reversal SISOOFDM communication[C]//2020International Conference on Computing,Networkingand Communications(ICNC),Big Island,HI,USA:IEEE Press 2020:222-227.doi:10.1109/ICNC47757.2020.9049811.]研究了频域TR预编码系统中的物理层安全，为了使系统的保密速率更大，设计了一种适用于此系统的人工噪声，结果表明所提出的方案可以提高物理层安全。对于TR-OFDM系统，当BOF大于1时，TR的聚焦效应能提高期望接收机的SINR，降低在非期望接收机处的信号泄露，因此在多用户系统中才能更好地发挥TR-OFDM的优势。在频域TR多用户系统中，可采用平均分配的方法为用户各信号分配功率，但平均功率分配往往不能获得最好的系统性能，因此可以对用户信号的功率分配进行优化，从而获得更好的系统性能。

发明内容

本发明的目的在于给出一种时间反转OFDM多用户通信系统中，在满足各用户发送总功率约束下，增强系统和速率性能的用户信号发送功率优化分配方法，该方法根据信道状态信息分配功率，实现系统和速率的最大化。

为了实现上述目的本发明采用如下技术方案：在常规时间反转(Time Reversal，TR)OFDM多用户通信系统的基础上，推导系统和速率公式，并构造以最大化系统和速率为目标的多个用户功率联合分配的优化问题，然后利用KKT条件得到用户的功率分配的解，最后利用迭代算法来联合优化所有用户的发送功率。

(1)对信号的传输过程进行分析，构建通信系统模型，根据构建的系统模型得到用户接收信干噪比，进而得到系统和速率；所述通信系统模型在常规时间反转OFDM多用户通信系统的基础上，发送端和接收端均配备单根天线，发送端发送给多个用户的信号先经过TR预编码矩阵，然后被分配到OFDM子信道上传输。

(2)以和速率最大化为优化目标，给出在发送总功率约束条件下多个用户功率联合分配的数学模型。

接收信干噪比为：

其中，K为用户数量，p_k,m为发送端分配给用户k的符号m的功率，D为速率回退因子，为信道的噪声方差，H _k,m为对角矩阵H _k的第m个对角元素，H _k为频域信道矩阵，表示为

其中，Q为OFDM系统子载波数，对角元素H _k,q,q＝1,2,...,Q为频域信道系数，表示为h_k,l为用户k的多径信道的第l径的系数，l＝1,2,...,L，这里L为信道脉冲响应的长度；

矩阵为H _k的共轭形式，即

系统和速率为：

式中，M表示向每个用户发送的符号数。

(3)利用Karush-Kuhn-Tucker(KKT)条件得到用户的功率分配的解。

(4)利用迭代算法来联合优化所有用户的发送功率。

本发明在对用户信号发送功率进行优化分配的过程中，根据信道的信道状态信息分配功率。由于只优化用户信号发送功率分配，本发明的优化复杂度较低。并且仿真实验表明，本发明给出的方法明显能获得比平均功率分配方法更高的系统和速率。本发明在较低优化复杂度下，提升了系统的和速率性能。

附图说明

图1为本发明的通信系统模型；

图2为用户k的子载波利用示意图；

图3为本发明的等效通信系统模型；

图4为本发明与平均功率分配及传统迭代注水方法的系统和速率的对比；

图5为本发明及传统迭代注水方法的系统和速率与迭代次数的关系；

图6为速率回退因子对系统和速率的影响；

图7为用户数对系统和速率的影响。

具体实施方式

对本发明中涉及到的符号有以下注释：带下划线的字母表示频域符号，没有下划线的表示时域符号；小写粗体字母表示矢量，大写粗体字母表示矩阵；白体字母表示标量；|·|、||·||、(·)^*和(·)^H分别表示求绝对值或模值、2范数、共轭及共轭转置。

TR-OFDM多用户通信系统的模型如图1所示，系统有K个用户，发送端和接收端均配备单根天线。发送端同时向K个用户发送相互独立的数据矢量，其中向用户k发送的数据矢量记为x _k，k＝1,2,...,K。x _k由M个数据符号x _k,m组成，m＝1,2,...,M，即x _k＝[x _k,1,x _k,2,...,x _k,M]。首先数据符号被大小为Q×M的扩展矩阵S加扩，加扩的目的是将数据符号分配到OFDM子载波上。Q为OFDM系统子载波数，Q＝DM。加扩后，每个数据符号以间隔M分配到D个不同的子载波上，这里的D即速率回退因子，相当于时域TR系统中的上采样因子。S由D个M×M的对角矩阵串联而成，其对角元素等概取值这是为了避免将相同的数据符号分配到不同OFDM子载波上，从而造成高峰均功率比，并保证加扩过程不改变信号功率，即S ^H S＝I _M，I _M表示M×M的单位阵。S表示为

TR-OFDM系统中通过扩展矩阵S加扩，相当于时域TR系统中的上采样。加扩后，用户k的子载波利用示意图如图2所示。加扩后的符号在发送前需要与预编码矩阵相乘进行预编码，相当于时域TR系统中的预滤波。预编码后的符号通过快速傅里叶逆变换(Inverse FastFourier Transform,IFFT)得到时域信号，加入循环前缀后发送出去。在接收端进行快速傅里叶变换(Fast Fourier Transform,FFT)再转化为频域信号。发送端与用户k间的多径信道是一个线性系统，在离散域中期脉冲响应可以表示为其中h_k,l为第l径的系数，也就是脉冲响应中的l个脉冲的强度，l＝0,1,...,L-1，这里L为信道脉冲响应的长度，也就是多径信道的路径数。接收机对接收信号进行FFT变换，得到频域符号。可以将IFFT、多径信道、FFT合在一起，等效为Q个频域并行信道，信道系数为信道输出的频域符号为发送的频域符号与对角矩阵H _k相乘。H _k的第q个对角元素为H _k,q，即

相应，系统模型可以等效为图3。其中v _k＝[v _k(1),v _k(2),...,v _k(Q)]为频域加性高斯白噪声(Additive White Gaussian Noise,AWGN)序列，为时域AWGN噪声序列v_k的傅里叶变换，v _k和v_k的元素具有相同的方差时域TR系统中的预滤波器为信道脉冲响应的共轭反转，发送信号经过该预滤波器对应频域中为频域符号与预编码矩阵/>相乘。/>为

用户接收到的序列需要与S ^H相乘进行解扩，用户k解扩后的信号为

用户k的接收SINR为用户接收的有用信号功率与用户间干扰功率和信道噪声功率之和的比值，即

其中，p_k,m为发送端分配给用户k的符号m的功率，D为速率回退因子大小，为信道的噪声方差，H _k,m为对角矩阵H _k的第m个对角元素。和速率是评价系统传输性能常用的指标，本发明以最大化系统和速率为目标对信号功率分配进行优化。系统和速率为

功率分配的优化过程中，应保证分配给用户各信号的功率之和小于等于用户的发送总功率，即

式中P_k为用户k的总发送功率。优化问题构造为

上述多目标联合优化是非凸优化，但对于用户k，若其余用户此时的功率分配是固定的，则用户k的功率分配最优解满足KKT条件。构造拉格朗日函数为

其中，λ＝[λ₁,λ₂,...,λ_K]^T是非负的拉格朗日乘子矢量，令Λ(p,λ)对p_k,m求导并等于0，则最优p和λ满足以下等式

其中，t_k,m包含了用户k对其他用户造成的干扰，定义为

解得用户k的功率分配为

式中，[x]⁺表示取0和x两者中的最大值，x为上式中括号内的内容。且因满足功率约束条件，有

令

发送端为信道条件更好的子信道上的信号分配更高的功率，当对其他用户造成过度的干扰时，t_k,m会降低分配给该信号的功率。λ_k可用二分法找到，当满足功率约束等式的λ_k非正值时，则λ_k为0。采用迭代算法求解所有用户的功率分配。为各用户各信号分配一个初始功率值，基于已赋的初值，可得到t，t为t_k,m的集合，即t＝{t_k,m}，k＝1,2,…,K，m＝1,2,…,M。首先从用户1开始，由于t和其他用户的功率分配都是固定的，因此可以得到用户1的功率分配；然后优化用户2的功率分配，求解时使用已经更新过的p_1,1,p_1,2,...,p_1,M值，而t和从用户3开始的功率分配仍为之前的值；余下其他用户的功率分配也进行类似的更新，得到更新后的p，p为p_k,m的集合，即p＝{p_k,m}，k＝1,2,…,K，m＝1,2,…,M。p更新完后，更新t，重复上述过程，迭代更新各用户的功率分配，直到连续两次迭代得到的功率分配的相对差值小于预先设定的值或者达到预定的迭代次数为止，最后得到各用户各信号的功率分配最优解p^*。具体的迭代算法如表1所示，其中n表示迭代次数，τ₁为控制迭代结束的相对差值，即收敛因子。

表1迭代算法

下面将结合附图，对本发明做进一步的详细描述。每个图中给出的系统和速率的仿真结果是1万组信道实现下的平均值。除非特别指明，仿真中的参数设置如下：采用一个SISO OFDM系统，其子载波数Q＝1024。信道为Rayleigh衰落信道，信道总带宽B＝15.36MHz，信道可分辨路径数L＝10，信道脉冲响应系数服从均值为零的复高斯随机分布。信道的方差为其中σ_T＝10/B为路径的均方根延迟，T_s＝1/B为采样周期。η_k＝η₀(d_k/d₀)^-c为信道的大尺度衰落系数。其中，c＝4为路径损耗指数，η₀为参考距离处的传输损耗，d₀＝10m为参考距离，d_k为发送端与用户k的距离，假设d_k是随机的，在100m到2000m之间。设定η₀＝10^-5，信道的噪声功率设为1×10^-11W。求解功率分配的迭代算法中最大迭代次数设为100，收敛因子τ₁＝1×10^-5。系统总发送功率为P，每个用户的发送总功率P_k＝P/K。

图4是是用户数K＝4，BOF D＝1，路径数L＝10时，系统和速率随发送功率增加的变化的仿真结果。图中改进迭代注水算法的结果即为本发明给出的用户信号发送功率的优化分配方法的结果，同时给出传统迭代注水功率分配方法及平均功率分配的系统和速率的仿真结果。传统迭代注水算法的区别在于t_k,m为0，即各信号分配功率时不管对其他用户的干扰，只需使自己的速率最大化，而平均功率分配则是不管信道条件，为所有用户所有信号分配相同的功率，即p_k,m＝P_k/M,m＝1,2,...,M。图中和速率是所有信道状态下的瞬时和速率的平均值。从图中可以看出，随着发送功率增加，系统和速率增大，然后逐渐趋于一个上限。这是因为当发送功率较小时，干扰功率相对于信道噪声功率较小，发送总功率增大，和速率相应增大。当发送功率增大到一定程度后，干扰功率将明显高于信道噪声功率。由于信号功率与干扰功率随发送功率同步增长，此时，发送功率增加时和速率不再有明显的增加。从图中可以看出，平均功率分配的系统和速率是最低的，这是因为采用平均功率分配时，所有信号分配到的功率是一样的，而不同的信号被分配到不同的子载波上，历经不同的子信道增益，若为信道增益低的和信道增益高的子信道上的信号分配相同的功率，则系统性能肯定不能达到最好。而采用传统的迭代注水算法时，由于各用户优化自己的信号功率分配时，只考虑使自己的速率最高，而忽略了给其他用户带来的干扰，因此往往给其他用户造成了较大的干扰，因此也不能实现最好的和速率性能。本发明给出的改进迭代注水功率分配算法既使自己的速率最大，同时t_k,m的存在使带给其他用户的干扰也减小，因此能最有效地提高系统和速率。如仿真结果所示，传统迭代注水算法下的和速率只能提高了8％，而改进迭代注水算法的和速率提高了13.8％，证明了本发明给出的改进迭代注水算法的性能更好。

图5是用户数K＝4，BOF D＝1，路径数L＝10，发送功率为30dBm时，本发明给出的功率分配算法及传统迭代注水算法的系统和速率与迭代次数的关系。可以从图中看出两种算法的收敛性，本发明给出的改进迭代注水算法收敛速度更快并且收敛时有更好的解。传统迭代注水算法在迭代次数达到6次才能收敛，而迭代注水算法在迭代次数达到4、5次时收敛，证明了改进迭代注水算法的收敛性更好，因此在实际应用中不需要等到完全收敛，只需迭代4-6次即可。

图6是用户数K＝4，路径数L＝10时，随发送功率增大，不同BOF下系统和速率的仿真结果。仿真结果显示，随着BOF增大，系统的和速率在减小，这是因为增大BOF虽然可以减小IUI，但同时BOF的增大使得能够传输的数据符号数减少了，也就导致符号速率的降低，且IUI的减小不足以弥补符号速率的大幅度降低，因此系统和速率减小。图7(a)和(b)分别是是BOF D＝1，路径数L分别为10和20时，随发送功率增大，不同用户数下的系统和速率。可以看到当发送功率较小时，随着用户数增多，系统和速率增大，这说明了频域TR预编码技术能够明显降低多用户系统中的IUI，但当发送功率较大时，噪声功率相对很小，可以忽略不计，干扰功率与发送功率同步增长，较多用户的干扰功率会更大，因此多用户的SINR会更小，系统和速率也就更低。即如图中所示，随发送功率增大，K＝2和K＝4的曲线会相交，并且用户数少的和速率会超过用户数多的。并且，从图7(a)和(b)可以看出，路径数增大，系统和速率会增大，这说明了频域TR适用于丰富的多径环境，能有效提高多径环境中的通信系统性能，且路径数越多，获得的聚焦增益越大。

Claims (6)

1.一种时间反转OFDM多用户通信系统的功率分配方法，其特征在于：包括以下步骤：

(1)构建通信系统模型，根据构建的系统模型得到用户接收信干噪比，进而得到系统和速率；

接收信干噪比为：

其中，K为用户数量，p_k,m为发送端分配给用户k的符号m的功率，D为速率回退因子，为信道的噪声方差，H _k,m为对角矩阵H _k的第m个对角元素，H _k为频域信道矩阵，表示为

其中，Q为OFDM系统子载波数，对角元素H _k,q,q＝1,2,...,Q为频域信道系数，表示为h_k,l为用户k的多径信道的第l径的系数，l＝1,2,...,L，这里L为信道脉冲响应的长度；

矩阵为H_k的共轭形式，即

系统和速率为：

式中，M表示向每个用户发送的符号数；

(2)以和速率最大化为优化目标，给出在发送总功率约束条件下多个用户功率联合分配的数学模型；

(3)利用KKT条件得到用户的功率分配的解；

(4)利用迭代算法来联合优化所有用户的发送功率。

2.根据权利要求1所述一种时间反转OFDM多用户通信系统的功率分配方法，其特征在于：所述通信系统模型在常规时间反转OFDM多用户通信系统的基础上，发送端和接收端均配备单根天线，发送端发送给多个用户的信号先经过TR预编码矩阵，然后被分配到OFDM子信道上传输。

3.根据权利要求1所述一种时间反转OFDM多用户通信系统的功率分配方法，其特征在于：步骤(2)所述数学模型以分配给用户各信号的功率之和小于等于用户的发送总功率为约束条件，以系统和速率为优化目标，具体为：

发送功率约束条件为：

式中P_k为用户k的总发送功率；

优化问题构造为：

4.根据权利要求3所述一种时间反转OFDM多用户通信系统的功率分配方法，其特征在于：步骤(3)所述利用KKT条件得到给定用户的功率分配的解，具体包括：优化用户k的功率分配时，由于其他用户的功率分配是固定的，因此用户k的功率分配满足KKT条件，构造满足功率约束的优化问题的拉格朗日函数为

其中，λ＝[λ₁,λ₂,...,λ_K]^T是非负的拉格朗日乘子矢量，令Λ(p,λ)对p_k,m求导并等于0，则最优p和λ满足以下等式

其中，t_k,m包含了用户k对其他用户造成的干扰，定义为

功率分配解为

式中，[x]⁺表示取0和x两者中的最大值。

5.根据权利要求1所述一种时间反转OFDM多用户通信系统的功率分配方法，其特征在于：步骤(4)所述迭代算法具体为：为各用户各信号分配一个初始功率值，基于已赋的初值，可得到t，首先从用户1开始，由于t和其他用户的功率分配都是固定的，因此可以得到用户1的功率分配；然后优化用户2的功率分配，求解时使用已经更新过的p_1,1,p_1,2,...,p_1,M值，而t和从用户3开始的功率分配仍为之前的值；余下其他用户的功率分配也进行类似的更新，p更新完后，更新t，重复上述过程，迭代更新各用户的功率分配，直到达到预定的迭代次数为止，最后得到各用户各信号的功率分配最优解p^*。

6.一种计算机可读存储介质，其存储有计算机程序，其特征在于：所述计算机程序被执行时，可实现权利要求1-5任一项所述的功率分配方法。