本申请分别要求于 2006年 11月 27日、 2007年 1月 16日提交中国专利 局、 申请号为 200610149051.6、 200710000242.0、 发明名称分别为"一种动态 功率频谱管理方法及频语优化系统及客户端设备 "的中国专利申请的优先权, 其全部内容通过引用结合在本申请中。
背景技术
数字用户线( DSL, Digital Subscriber Line )技术是一种通过电话双绞线, 即无屏蔽双绞线(UTP , Unshielded Twist Pair )进行数据传输的高速传输技 术, 包括非对称数字用户线( ADSL, Asymmetrical Digital Subscriber Line ), 甚高速数字用户线( VDSL, Very-high-bit-rate Digital Subscriber Line )、基于综 合业务数字网( ISDN, Integrated Services Digital Network )的用户数字线( IDSL, ISDN Digital Subscriber Line )和单线对高速数字用户线( SHDSL, Single-pair High-bit-rate Digital Subscriber Line )等。
随着各种数字用户线技术(xDSL )使用的频带的提高, 串扰, 尤其是高 频段的串扰问题表现得日益突出。 如图 1所示, 包括接入复用器及终端(图中 以两个终端为例), 所述接入复用器包括端口 1和端口 2。 由于 xDSL上下行 信道采用频分复用, 近端串扰(NEXT )对系统的性能不产生太大的危害; 但 远端串扰(FEXT )会严重影响线路的传输性能。 当一捆电缆内有多路用户都 要求开通 xDSL业务时, 会因为远端串扰使一些线路速率低、 性能不稳定、 甚 至不能开通等, 最终导致 DSLAM的出线率比较低。
在一个采用离散多音频调制( DMT, Discrete Multi-Tone Modulation )的 K
个用户, N个 TONE的通信模型中, 各个 TONE上信号传输可独立地表示为: yn = Ηηχη + ση
作为噪声, 则第 k个用户在第 η个 tone上可达到的数据速率可用香农信道容 量公式计算:
由上式可以看出, 串扰严重影响了线路的传输容量,也就是降低了线路速 动态功率频谱管理 ( DSM, Dynamic Spectrum Management ) 可以用于减 少串扰的影响。 具体地, DSM就是自动调整网络中各个 modem上的传输功率 来达到减小串扰的目的。 特别是在 CO/RT混合应用的情况下短线对长线的串 扰影响较大。 DSM的目标就是通过调整发射功率使每个 modem在达到自身速 率最大化和减少对其它 modem的串扰影响之间达到一个平衡。
DSM的目标为在每一个用户发送的总功率不超过门限时, 通过调整所有 用户分别在每一个子频带上的发送功率值, 来最大化所有用户的加权速率之 和。 因此, DSM问题的数学形式描述如下:
Κ N
maximize ∑ ok∑b
k=l n=l
subject to ≤ P
0 < Sk Vk,n
表示第 k个用户在第 n个子频带上所分配的功率;
°η表示第 k个用户线路在第 η个子频带上的传输系数;
G (j≠ k)表示第 j个用户对第 k个用户在第 η个子频带上的串扰系数; Pk表示第 k个用户的总功率门限;
°^表示第 k个用户的速率权系数;
σ2表示噪声功率;
Ν表示子频带总数;
Κ表示用户总数。
此问题为一非线性约束优化问题,目标函数及约束条件均为关于自变量的 非凸函数。因此不存在高效且完整的求解算法。在现有算法中,最为热门的为: 最优频谱平衡 ( OSB, Optimize Spectrum Balancing ), 迭代频谱平衡 ( ISB, Iterative Spectrum Balancing )、 迭代注水算法 ( IWF, Iterative Water - Filling )。 这三种算法都是通过引入拉格朗日乘子,将原始非线性约束优化问题转化为非 线性无约束优化问题。 因此, 上述公式转化为:
其中, λ¾为拉格朗曰乘子, 项表示为总功率约束(即等
同于上述公式中的 2 ≤^, \/k )0 Pk—般为一常数。 当 Pk为一常数时, 上述公式可以等价于: k=l
k=l n=l
故只需求解上述公式即可。
现有技术中的一种实现 DSM的方法为: 迭代注水算法。 这是一种會心算法, 只考虑 对第 k个用户线上速率的影响, 当用户线不 是当前用户 k时, 令其权重 Wj = 0,λ」 = 0,Vj≠k。 这样, 实际上的^可以简 化为 Jn≡j; = wkbk n - ksk n , 因为此时 = 0,入」= 0,^≠1^且 ,(>]≠10固定。 对于求解 = argmax J^,由于 J = wk¾ - λ 有唯一极值点,可以通过推导求 得其极值点的解析表达式。 整个算法不断迭代以求得各个用户在各个 tone上 的 max , 直到所有用户的功率分配都不再变化为止。
IWF的算法复杂度较低, 对于较大的 N和 K都可进行计算。 并且它是完 全自制的, 即各个用户只需优化自身的速率和满足自身的功率约束即可, 而不 需要不同用户之间进行数据信息的交互, 即不需要中心管理器, 易于在实际系 统中实现, 但是, 由于 IWF是會心算法, 所以无法保证最优解或近似最优解, 依赖于初始解。
具体实施方式
本发明实施例提供一种动态功率频谱管理方法、频语优化系统及客户端设 备, 用于提高频谱管理效果。
请参阅图 2,为本发明实施例动态功率频谱管理方法的流程图,包括步骤: 201、 预置税因子, 并初始化所有用户的功率谱;
其中,频谱管理中心随机设置一个满足预置条件的税因子, 并将该税因子 发送至 DSLAM。
其中, 预置条件可以根据实际情况进行确定, 本实施例中, 预置条件为非 负。
在 DSLAM中, 每个用户都有预先设置的原始功率谱, 在预置税因子时, 需要对这些原始功率谱进行初始化处理。
202、 判断是否所有用户的功率谱都已经收敛, 若是, 则执行步骤 207, 若否, 则执行步骤 203;
203、 确定需要处理的用户;
其中, 若该步骤为第一次执行, 则需要处理的用户为第一个用户, 若不是 第一次执行, 则需要处理的用户为当前用户的下一用户;
若当前用户是第 K个用户 (即最后一个用户), 则其下一个用户为第一个 用户。
204、 加载当前用户对应的注水参数;
其中, 注水参数包括税因子和串扰项,还包括计算过程中涉及到所需要的 其它参数。
串扰项可以在用户端直接测量(它会随着其他用户功率谱的改变而改变), 也可以根据后面的公式(5 )进行计算。
205、 计算当前用户的功率分配;
其中 ,在计算功率分配之前需要先进行公式的推导与演算,具体过程如下: 频语优化问题是在每一个用户发送的总功率不超过门限时,通过调整所有 用户分别在每一个子频带上的发送功率值, 来最大化所有用户的加权速率之 和。 其数学公式表示如下: maximize 公式( 1 )
subject to < Pj
0 < S Vk,n
表示第 k个用户在第 n个子频带上所分配的功率;
表示第 k个用户线路在第 η个子频带上的传输系数;
表示第 k个用户线路在第 η个子频带上所受到的串扰;
Pk表示第 k个用户的总功率门限;
表示第 k个用户的速率权系数;
2表示噪声功率;
Ν表示子频带总数;
Κ表示用户总数。
每个用户所受串扰的大小 不仅取决于其他用户在当前子频带上的功率: 而且有可能取决于其他用户在相邻子频带上的功率 (如果系统不同步的话)。 具体来说, i;=∑∑Y(n-m)G¾J
j≠k m=l (m)表示相邻频带串扰系数( Inter-Carrier-Interference )。
当系统是同步的, (m) = (m) = S(m)就是一个冲击函数,所以此时第 k 个用户线路在第 n , 也就是说没有相邻频
带串扰。 当系统不同步时, 最坏情况的相邻频带串扰系数可以表示为 m = 0 y(m) = YICI(m) 2
m≠ 0, < m <—
兀 2 2
N sin (― m
N 这里 , ICI(m) 具有对称性和循环性: 7ICI(m) = 7ICI(_m) , ICI(m) = 7ICI(N_m)。
通过拉格朗日乘子的方法把上面如公式(1)所示的有约束问题转化成无 约束问题, 如下式: 公式(2)
其中 和 分别为拉格朗日乘子。其中 l k Pk_∑ 项表示为总功率
Κ N
约束(即等同上述公式中的∑S ≤Pk, Vk)。 μ 项表示为发送功率
k=l n=l 非负约束(0≤S Vk,n)。 公式(2)所描述的问题等同于公式(1)的 问题, 依然是一个非凸问题, 没有高效的数字解决方法。 同时公式(2)所描 述的问题存在多个一阶极点 , 每一个极大值点就是这个问题的一个局部最优 解。 可以通过对公式(2)对 Si求偏导的方法得到下式: ω,
-λ,+μ^ =0 公式( 3 )
. k ^ N 7(n-m)GiiSj
其中 t = Y ω GJ Y ~ ^-^ 公式(4)
n 3 η^ί (GiiSJ +Γ +& )(Γ +& )
t叫作税因子。
如果能通过一定的方法从公式(3) 中找到一个 的解, 那么这个解是一 个局部最优解。从实际的角度来看如果这个方法的算法比较简单, 性能能有很 大的提升。 说明就通过一个简单的算法找到了公式(2)所描述的问题的一个 局部最优解。
当然这个除了要满足公式外, 同样还要满足以下约束条件:
pk-∑s;>o 公式( 6 )
实际上公式 (3)和公式 (6) —起构成了公式 (1 ) 所描述问题的 KKT 条件 ( Karush-Kuhn-Tucher conditions ), 从上面的描述可知这个条件对于满足 公式(1 )所描述的问题是一个必要条件。 意思是只要是最优解一定满足这个 条件, 但满足这个条件的解有可能是局部最优解, 就不一定是最优解。
可以进一步的通过公式演算把公式(3)和公式(6)构成的 KKT条件成 以下公式:
,Vn,k 公式( 7 )
, Pk-∑S = 0,Vk 公式( 8 )
+k ^ ( SINR SINRJ m
t =Y ro.GJ Y γ(η - m、) ~― ^ - 公式( 9 ) n 3 n^ GJJSJ 1 + SINR 其中公式( 7 ) 中 [x]+表示 max(0, x) , 其作用等效于拉格朗日乘子 , 所 以可以在(7) 中把 k省略掉。 公式(9)中, SINRj = . m m. , 其物理意义 n m Γ m +& m 为第 j个用户的第 m个子频带上的信号干扰噪声比,这个量可以在用户线上直 接测量。
针对单个用户的传统的优化注水算法, 其表达式如下:
∑ , 对这个公式就 sn 导得:
1 G
- λ = 0。 进一步化简得:
1 σ
S 公式( 10a )
λ G„ 其中 SNR = Gn :Sn , 其中丄表示注水(water level ) 水平, 对于同一个 σ λ 用户的所有子频带是一个常量。
传统的多用户的迭代注水是基于单用户的注水,但是因为有了多个用户之 间的串扰, 其频谱更新具体如下:
Vn,k 公式 ( 10b )
传统的多用户迭代注水的算法具体描述如下:
第一步: 对于用户 k,假设其他用户的频谱不变, 通过传统的排序算法的 、
方法找到一个 λ , 使得功率总和 Ρ
第二步:通过公式( 10b 十算出用户 k,各个子频带的 ,其中 n = 1,2,···,Ν。 第三步: 对所有的用户重复步骤一和步骤二, 直至收敛。 对比单用户的注水公式( 10a ), 传统的多用户的注水公式( 10b ), 以及 公式(7 ), 不难发现单用户注水(即当 tk = 0和 Ik = 0 )和传统的多用户注水
(即当 t =0 )都是公式(7)的特例。 为了解决公式(7)和(8)的问题, 则需要利用修正迭代注水的方法找 l 各个用户的发送功率谱。 其中步骤 205具体的步骤如下:
给定 t, k=l,2,...K, n=l,2,...N。 对当前用户 k执行以下步骤:
1、 定义功率总和公式
2、 如果 p(0)≤Pk, 设 =0, Sn k =0, 其中 η = 12... Ν
3、 如果 p(0)〉Pk, 使用二分法(bi-section)的方法寻找最优 λ¾和功率
3a )初始化入 k, 和
) < Pk
3b)重复以下步骤:
3bl)计算 2 k
3b2)如果 p k)〉Pk, mm = k
3b3)如果 P(^)<p k, λ
直至 ^收敛。
3c)通过公式(7)和收敛的^, 计算出当前用户各个子频带的 , 其中 n = l2 .. N 当 ^收敛后, 当前用户的总功率门限就会满足, 此时执行步骤 202 206、 更新税因子;
其中,根据测量的信号串扰噪声比与步骤 205的计算结果计算税因子,具 体有以下三种情况:
方案一
基于无记忆的一阶最优化条件的更新, 税因子的计算可根据如下表达式 k _ ^ ik N 、SINRj SINRJ ^ r 、 tk =Yro.GJky Y(n-m) ~—^- 公式(11 )
^ J - G S]m l + SINRJ m 此方案的优点在于更新简单,而且如果更新收敛,其必收敛于局部最优解。 但是, 此方案不保证收敛性。
方案二
基于有记忆的一阶最优化条件的更新, 也就是 ί τ
公式(12) 其中 ' 表示上一时刻的税因子, 是当前时刻的税因子; 是时间参数。 α(τ)表示一个和时间有关的记忆参数,其属于(0,1]。 当 α(τ) = 1时,公式(12) 就变化成公式(11 ), 也就是说无记忆更新是有记忆更新的一个特例。 此方案 的优点在于可以通过适当选取°^ )参数(比如反比于时间)使得税因子一定 收敛。 但缺点在于收敛的税因子不保证满足 ΚΚΤ条件。 也就是说此方案可能 收敛到一个次优解, 但不一定是局部最优解。
方案三
基于目标函数(所有用户的加权速率之和 )的更新。 实现原理是每次税因 子的更新都确保目标函数的值会增加。 具体方案如下:
(a)假设当前时刻的功率谱为 S 目标函数的值为 υτ, 未更新的税因子
为 可以通过公式(11)或者公式(12)计算出临时税因子 t'。
( b )根据临时税因子进一步推算出其对应的临时功率谱 S'以及临时目标 函数值 U'。
(c)如果 U'〉IT, 更新下一时刻的税因子为 τ _τ , 然后发送, 更新结
(d)如果 υ'≤ΙΓ, 更新临时税因子为 , 返回步骤(b)。
其中目标函数就是所有用户的加权速率之和, 即:
此方案的优点在于它能保证税因子收敛到一个局部最优解。缺点在于需要 更多的计算量。
可以通过方案一, 二, 三的任何一种对税因子进行更新。 具体使用那种方 法取决于很多因素, 比如实际的网络拓朴, 频语管理中心的计算能力, 以及用 户实际的串扰程度。
综上所述,可以通过结合税因子的更新和修正迭代注水来找到一个满足公 式(7)、 (8)和(9)的功率谱分布。 迭代的顺序可以有多种不同的组合, 例 如:
组合一: 首先所有用户同时更新税因子, 然后基于固定的税因子所有用户 进行修正迭代注水。 如此反复直到收敛。 组合二: 每个用户每次都更新自己的税因子, 拉格朗日乘子, 功率谱, 如 此反复直到收敛。
(ί ,λ^ )→ (ti,x2,s2 →■■■→ (tK a ,xK,sK a)→ (t )→ -
下面列举了组合一具体实现的算法。其中包括如何连续更新税因子, 以及 如何在给定税因子的情况下进行修正迭代注水直至收敛。 当税因子收敛后,用 修正迭代注水计算出来的各个用户的发送功率谱就是公式(1 )所描述问题的 一个局部最优解。 这个解的结果取决于税因子初始值的设置。
重复对所有用户执行以下步骤:
利用修正迭代注水的结果在线路上测试出或是计算出 SINR ,再通过公式 ( 11 )、 (12 )或者方案三计算出相应的 。
根据计算出的 进行多用户的修正迭代注水, 直至所有用户收敛。
207、 判断税因子是否收敛, 若收敛, 则执行步骤 208 , 若未收敛, 则执 行步骤 202;
208、 结束。
请参阅图 3 , 本发明频语优化系统实施例包括:
频语管理中心 301 , 数字用户线接入复用器 302以 ^户端设备 303;
所述频语管理中心 301 用于设置税因子并将其发送至数字用户线接入复 用器 302;
所述数字用户线接入复用器 302, 用于根据税因子计算当每个用户最大功 率约束小于或等于该用户总功率门限时对应的注水参数,并根据所述注水参数 或
所述频谱管理中心 301 用于设置税因子并将其发送至数字用户线接入复 用器 302;
所述数字用户线接入复用器 302用于将税因子转发至客户端设备 303; 所述客户端设备 303 用于根据税因子计算当每个用户最大功率约束小于 或等于该用户总功率门限时对应的注水参数,并根据所述注水参数计算出每个 用户在各子频带上所分配的功率直至所有用户收敛。
其中, 所述频语管理中心 301包括:
收集单元 3011 , 更新单元 3012以及传送单元 3013;
所述收集单元 3011用于获取不同用户的信号串扰噪声比, 功率谱以及相 关参数, 并将获取到的参数发送至更新单元 3012;
其中, 相关参数包括传输系数等, 或者其它所需参数, 或者不包括任何参 数。
所述更新单元 3012用于接收收集单元发送的参数, 并根据所述参数更新 税因子, 并将更新后的税因子发送至传送单元 3013;
所述传送单元 3013用于将接收到的税因子发送至数字用户线接入复用器
302。
其中, 由数字用户线接入复用器 302计算功率分配时,在所有用户计算完 功率分配之后数字用户线接入复用器 302将计算得到或检测到的信号串扰噪 声比发送至频语管理中心 301;
由客户端设备 303计算功率分配时,在所有用户计算完功率分配之后客户 端设备 303通过数字用户线接入复用器 302将计算得到或检测到的信号串扰噪 声比发送至频语管理中心 301;
所述频语管理中心 301根据接收到的串扰噪声比设置税因子。
其中, 所述数字用户线接入复用器 302包括:
频语管理中心感知器,用于当频谱管理中心 301失效时,将税因子设为 0。 具体流程为:
首先频谱管理中心设置一个满足预置条件的税因子 ,并通过税因子接口把 各个用户的税因子分发到各自的 DSLAM中 ,如果带税因子的迭代注水算法是 在客户端设备中完成,那么 DSLAM再通过 DSLAM和客户端设备间的内嵌通 道 兌因子分发到各用户所在的客户端设备中, 然后 DSLAM和 /或客户端设 备根据这个税因子进行计算, 当所有用户收敛后,各用户通过 DSLAM于频谱 管理中心的 SINR接口把各自的信号串扰噪声比上传到频谱管理中心, 由频谱 管理中心计算出新的税因子, 并重复上述过程, 直至收敛。
请参阅图 4, 本发明提供的客户端设备 303实施例包括:
税因子接收单元 401 , 第一注水单元 402以^送单元 403;
所述税因子接收单元 401 用于接收税因子并将接收到的税因子发送至注 水单元 402;
所述第一注水单元 402 用于根据接收到的税因子计算每个用户的总功率 小于或等于该用户总功率门限时对应的注水参数;并根据所述注水参数计算出 每个用户在各子频带上所分配的功率, 并将计算结果发送至发送单元;
所述发送单元用于向频谱管理中心发送计算结果。
优选的, 所述客户端设备还可以包括第二注水单元 404, 用于根据频谱管 理中心更新过税因子后再次计算每个用户在各子频带上所分配的功率直至收 敛, 并将计算结果发送至发送单元 403;
所述发送单元 403用于向频谱管理中心发送计算结果。
另外, 本发明实施例还提供一种频谱优化设备, 包括: 获取单元、 第一注 水参数计算单元和第一功率计算单元。其中,所述获取单元,用于获取税因子; 所述第一注水参数计算单元, 用于根据所述税因子计算每个用户的注水参数。 所述注水参数定义为每个用户的总功率小于或等于该用户总功率门限时对应 的值; 所述第一功率计算单元, 与第一注水参数计算单元相连, 用于根据所述 注水参数计算出每个用户在各子频带上所分配的功率。
优选的, 所述设备还包括: 税因子设置单元, 分别与获取单元和第一注水 参数计算单元相连, 用于设置每个用户满足预置条件的税因子, 所述获取单元 获取设置的税因子。
优选的, 所述设备还包括: 税因子更新单元, 分别与获取单元和第一功率 计算单元相连, 用于对每个用户的税因子进行更新; 所述获取单元获取更新后 的税因子。
优选的, 所述税因子更新单元至少包括下述一种: 串扰噪声比确定子单元 和第一税因子计算子单元; 串扰噪声比确定子单元、参数获取子单元和第二税 因子计算子单元; 目标函数获取子单元、临时目标函数计算子单元和税因子确 定子单元。所述第一税因子计算子单元和第二税因子计算子单元可以集成在一
起, 也可以单独存在。 其中,
所述串扰噪声比确定子单元,用于计算或测量每个用户在各子频带上的信 号串扰噪声比; 所述第一税因子计算子单元, 用于根据所述信号串扰噪声比以 及每个用户在各子频带上所分配的功率计算出对应的税因子。
所述串扰噪声比确定子单元,用于计算或测量每个用户在各子频带上的信 号串扰噪声比; 所述参数获取子单元, 用于获取记忆参数以及时间参数; 所述 第二税因子计算子单元,用于根据所述记忆参数、时间参数前一时刻的税因子、 信号串扰噪声比以及每个用户在各子频带上所分配的功率计算当前税因子。
所述目标函数获取子单元,用于获取所有用户的加权速率之和的数值并将 所述数值作为目标函数数值; 所述临时目标函数计算子单元, 用于计算临时税 因子并根据所述临时税因子计算临时目标函数数值; 所述税因子确定子单元, 用于将所述目标函数数值与临时目标函数数值进行比较,根据比较结果确定税 因子。
优选的, 在上述频谱优化设备的基础上, 还可以包括: 第二注水参数计算 单元和第二功率计算单元。 其中
所述第二注水参数计算单元, 用于根据所述第一税因子计算子单元、第二 税因子计算子单元或税因子确定子单元更新后的税因子计算每个用户的总功 率小于或等于该用户总功率门限时对应的注水参数; 所述第二功率计算单元, 用于根据所述注水参数计算出每个用户在各子频带上所分配的功率直至收敛。
对于所述频语优化设备中各个单元的功能和作用,详见上述方法、 系统及 客户端设备中对应部分的功能和作用, 在此不再赞述。
以上对本发明所提供的一种动态功率频谱管理方法及频语优化系统及客 户端设备进行了详细介绍,本文中应用了具体个例对本发明的原理及实施方式 进行了阐述; 对于本领域的一般技术人员, 依据本发明的思想, 在具体实施方 式及应用范围上均会有改变之处, 综上所述,本说明书内容不应理解为对本发 明的限制。