无线/移动网络中基于动态阈值的呼叫准入控制方法

文档序号:7652277阅读:244来源:国知局
专利名称:无线/移动网络中基于动态阈值的呼叫准入控制方法
技术领域
无线/移动网络中基于动态阈值的呼叫准入控制方法属于无线/移动网络的服务质量控制算法领域。
背景技术
移动业务需求的高速增长促进了无线/移动网络的快速发展,现有的无线/移动网络如通用分组无线业务网络(GPRS)、无线局域网(WLAN)以及广带无线接入网(BWA)等,正朝着宽带化、智能化的方向发展,目的是为移动用户提供具有不同服务质量(QoS)保障的个性化、智能化的业务。然而在无线/移动网络中,由于存在拓扑和资源动态变化、带宽稀缺、多径衰落和阴影衰落等多种不利因素,使得在这样的网络环境下保证用户的QoS将比在有线网络中要复杂得多,尤其是在现有的组网技术中,为了获得更高的容量和更好的性能,无线/移动网络广泛采用微蜂窝/微微蜂窝(micro-cell/pico-cell)的结构,使得无线/移动子网的尺寸越变越小,从而导致更加频繁的切换和更加快速的网络负荷变化,最终使移动用户的QoS保障面临更为严峻的挑战。在这样的情况下,采用呼叫准入控制技术以提高网络资源利用率、保障已接纳呼叫的QoS,就显得尤为重要。
作为一种重要的QoS控制技术,呼叫准入控制的作用是对到来的呼叫连接请求作出是否接纳的判断,实现在保障已接纳呼叫的QoS的情况下,尽可能多地接纳新到达的呼叫、提高网络资源利用率的目的。在无线/移动网络中,存在两种类型的呼叫新呼叫和切换呼叫。其中,新呼叫是指由当前无线/移动子网的移动用户发起的呼叫,而切换呼叫是指由别的无线/移动子网的移动用户发起、并由于用户的移动而切换到当前无线/移动子网中的呼叫。无论是新呼叫还是切换呼叫在向网络发起连接请求时,它们是否能够享用网络资源、被网络所接纳,取决于网络所采用的呼叫准入控制技术,因此,新呼叫阻塞率和切换呼叫中断率成为衡量呼叫准入控制性能的两个重要参数,其中,新呼叫阻塞率是指在某一时段内,拒绝的新呼叫数量占该时段内到达的新呼叫数量的比值,而切换呼叫中断率则是指在某一时段内,拒绝的切换呼叫数量占该时段内到来的切换呼叫数量的比值。如何设计一种有效的呼叫准入控制方法,降低新呼叫阻塞率和切换呼叫中断率是当前呼叫准入控制研究中的关键问题。

发明内容
本发明的目的在于提供了一种能够减少新呼叫阻塞率和切换呼叫中断率的基于动态阈值的呼叫准入控制方法。
本发明的特征在于,它依次含有以下步骤步骤(1.)在系统即无线/移动网络中的基站、或接入路由器AR、或无线接入点AP中,记录业务类别总数K,K∈N以及剩余带宽Br,并初始化定时器,所述定时器用于将系统的时间分为等长的时段,时段长为20秒;步骤(2.)初始化计数器Ni,i=1,2,…,K为0,用于记录一个时段内系统接纳的第i,i=1,2,…,K类业务的新呼叫数,其中,新呼叫是指由当前无线/移动子网的移动用户发起的呼叫;初始化计数器Hi,i=1,2,…,K,用于记录一个时段内系统接纳的第i,i=1,2,…,K类业务的切换呼叫数,其中,切换呼叫是指由别的无线/移动子网的移动用户发起、并由于用户的移动而切换到当前无线/移动子网中的呼叫;步骤(3.)为了实施呼叫准入控制,系统在每个时段开始时为各类业务的新呼叫和切换呼叫设置呼叫准入控制阈值,其中第i类业务的新呼叫准入控制阈值rni和第i类业务的切换呼叫准入控制阈值rhi,i=1,2,…,K,分别是系统平均收益函数的参数,要在剩余带宽的约束下,使得系统平均收益函数最大化,所述的系统平均收益函数按照以下步骤获得步骤(3.1.)系统预测本时段内第i,i=1,2,…,K类业务的平均新呼叫到达率λni和平均切换呼叫到达率λhi,i=1,2,…,K,所述的平均新呼叫到达率和平均切换呼叫到达率分别指每秒达到系统的新呼叫数和切换呼叫数的平均值,这两个参数分别由系统根据呼叫到达的历史信息经统计分析得到,这里的呼叫到达的历史信息指本时段之前的各时段的新呼叫到达率和切换呼叫到达率的平均值;步骤(3.2.)系统根据本时段内的λni,i=1,2,…,K计算在本时段内到达系统的第i,i=1,2,…,K类业务的新呼叫个数Nni的概率Pn(Nni)=(λni×20)Nnie-(λni×20)/Nni!,i=1,2,···,K]]>Nni∈N,0<Nni<∞;系统根据本时段内的λhi,i=1,2,…,K计算在本时段内到达系统的第i类业务的切换呼叫个数的概率Ph(Nhi)=(λhi×20)Nnie-(λhi×20)/Nhi!,i=1,2,···,K,]]>Nhi∈N,0<Nhi<∞;步骤(3.3.)当本时段内到达系统的第i类业务的新呼叫个数Nni,0≤Nni≤rni,i=1,2,…,K且切换呼叫个数Nhi,0≤Nhi≤rhi,i=1,2,…,K时,系统的收益计算为pi×(Nni+Nhi), i=1,2,…,K,所述的pi是当一个第i类业务的呼叫被接纳时系统得到的奖赏值,且第i类业务的优先级越高,pi越大,0<pi<10,i=1,2,…,K;步骤(3.4.)当本时段内到达系统的第i类业务的新呼叫个数Nni,0≤Nni≤rni,i=1,2,…,K且切换呼叫个数Nhi,Nhi>rhi,i=1,2,…,K时,系统的收益计算为pi×(Nni+rhi)-chi×(Nhi-rhi)-αi×(rni-Nni),i=1,2,…,K,其中,Chi为拒绝一个第i类业务的切换呼叫时系统得到的惩罚值,且第i类业务的优先级越高,Chi越大,0<Chi<10,,i=1,2,…,K;而αi为系统为第i类业务的新呼叫分配过多的阈值,从而造成资源浪费时的惩罚值,第i类业务的优先级越高,αi越大,0<αi<10,i=1,2,…,K,;步骤(3.5.)当本时段内到达系统的第i类业务的新呼叫个数Nni,Nni>rni,i=1,2,…,K且切换呼叫个数Nhi,0≤Nhi≤rhi,i=1,2,…,K时,系统的收益计算为pi×(rni+Nhi)-cni×(Nni-rni)-βi×(rhi-Nhi),其中,Cni为拒绝一个第i类业务的新呼叫时系统得到的惩罚值且第i类业务的优先级越高,Cni越大,0<Cni<10,Cni<Chi,i=1,2,…,K;而βi为系统为第i类的切换呼叫分配过多的阈值,从而造成资源浪费时的惩罚值,且第i类业务的优先级越高,βi越大,0<βi<10,βi<αi,i=1,2,…,K;步骤(3.6.)当本时段内到达系统的第i类业务的新呼叫个数Nni,Nni>rni,i=1,2,…,K且切换呼叫个数Nhi,Nhi>rhi,i=1,2,…,K时,系统的收益计算为pi×(rni+rhi)-cni×(Nni-rni)-chi×(Nhi-rhi),i=1,2,…,K;步骤(3.7.)系统根据步骤(3.2.)得到的Pn(Nni)和Ph(Nhi),i=1,2,…,K以及步骤(3.3.)~步骤(3.6.)计算出的当Nni和Nhi,i=1,2,…,K处于不同范围时的系统的收益,得到系统的平均收益函数 其中,Gi(rni,rhi)为由于第i类业务的呼叫到达给系统带来的平均收益,它是阈值rni和rhi的函数,i=1,2,…,K,其计算为Gi(rni,rhi)=rni(pi+cni)+rhi(pi+chi)-20×(λnicni+λhichi)+ΣNni=0rniΣNhi=0rhi(αi(rni-Nni)+βi(rhi-Nhi))Pn(Nni)Ph(Nhi)+]]>ΣNni=0rni(pi+αi+cni)(Nni-rni)Pn(Nni)+ΣNhi=0rhi(pi+βi+chi)(Nhi-rhi)Ph(Nhi);]]>步骤(3.8.)系统在剩余带宽的约束下,计算使得平均收益函数最大化的rni和rhi,i=1,2,…,K;步骤(4.)系统启动定时器,并进入等待事件到达状态,所述的事件包括新呼叫到达事件、切换呼叫到达事件、呼叫离开事件和定时器超时事件;步骤(5.)在系统处于等待事件到达状态时,如果第i,i=1,2,…,K类业务的新呼叫到达,则系统判断接纳该呼叫是否会使该时段内接纳的第i,i=1,2,…,K类业务的新呼叫数目Ni,i=1,2,…,K超过阈值rni,i=1,2,…,K,如果是,则系统拒绝该呼叫,否则,系统接纳该呼叫,并执行如下操作步骤(5.1.)为该呼叫分配所需的带宽;步骤(5.2.)计数器Ni,i=1,2,…,K加1;步骤(5.3.)剩余带宽Br被更新为当前剩余带宽减去为该呼叫分配的带宽;步骤(5.4.)系统重新进入等待事件到达状态;步骤(6.)在系统处于等待事件到达状态时,如果第i,i=1,2,…,K类业务的切换呼叫到达,则系统判断接纳该呼叫是否会使该时段接纳的第i,i=1,2,…,K类业务的切换呼叫数目Hi,i=1,2,…,K超过阈值rhi,i=1,2,…,K,如果是,则系统拒绝该呼叫,否则,系统接纳该呼叫,并执行如下操作步骤(6.1.)为该呼叫分配所需的带宽;步骤(6.2.)计数器Hi,i=1,2,…,K加1;步骤(6.3.)剩余带宽Br被更新为当前剩余带宽减去为该呼叫分配的带宽;步骤(6.4)系统重新进入等待事件到达状态;步骤(7.)在系统处于等待事件到达状态时,如果有呼叫离开,则系统将剩余带宽Br更新为当前剩余带宽加上该呼叫释放的带宽,然后重新进入等待事件到达状态;
步骤(8.)在系统处于等待事件到达状态时,如果定时器超时,则从上述的步骤(2.)重新开始执行。
按照上述实施方式,我们在系统提供两类业务的情况下,以2000个仿真时段内的平均新呼叫阻塞率和平均切换呼叫中断率为衡量系统性能的参数,对比了本方法和其他两个呼叫准入控制方法——静态多阈值呼叫准入控制方法(SMCAC)和收益优化的呼叫准入控制方法(ROCAC)的性能,所述的平均新呼叫阻塞率(或平均切换呼叫中断率)是指2000个时段内的新呼叫阻塞率(切换呼叫中断率)之和除以2000,时段长即定时器长为20秒。图2和图3分别示出了当系统在初始时刻的剩余带宽Br=11兆比特/秒、第一类业务的新呼叫到达率和切换呼叫到达率,以及第二类业务的切换呼叫到达率均为0.01个呼叫/秒时,使用本方法、SMCAC和ROCAC时,各类业务的平均新呼叫阻塞率和平均切换呼叫中断率随第二类业务的平均新呼叫到达率的变化,在仿真中认为第一类业务的优先级要高于第二类业务的优先级,并设置如下奖赏值和惩罚值p1=8,p2=6,Cn1=4,Cn2=3,Ch1=8,Ch2=6,α1=6,α2=4,β1=3,β2=2。在每种方法中,系统的性能由第一类业务和第二类业务的平均新呼叫阻塞率和平均切换呼叫中断率共同决定,平均新呼叫阻塞率和平均切换呼叫中断率越低,系统的性能越好。
从两图可看出,随着第二类业务的平均新呼叫到达率的增加,SMCAC的第二类业务的平均新呼叫阻塞率急剧增加,虽然SMCAC的其它性能指标较低,但最高接近23%的第二类业务的平均新呼叫阻塞率使系统的整体性能急剧下降;此外,随着第二类业务的平均新呼叫到达率的增加,ROCAC的各类业务的平均新呼叫阻塞率均出现了较大幅度增长,最高接近6%,远远大于各类业务的平均切换中断率,虽然系统使用ROCAC的整体性能要比使用SMCAC的整体性能好,但ROCAC为了保障平均切换呼叫中断率而过分地牺牲平均新呼叫阻塞率的做法,是系统和用户所不希望的;最后,从两图可进一步看出,随着第二类业务的平均新呼叫到达率的增加,本方法的平均新呼叫阻塞率和平均切换中断率均较低,即使在第二类业务的平均新呼叫到达率大于等于0.025时,两类业务的平均切换呼叫中断率开始出现较大上升,但其最大值仍然被限制在0.53%和0.35%之内,而这些性能指标值是系统和用户可以接受的。
综上所述,本方法的性能比现有的两种呼叫准入控制方法——SMCAC和ROCAC的性能更优,因为它既不像SMCAC那样,由于不能动态地调整带宽,从而在平均呼叫到达率发生变化时,系统的整体性能变得很差;也不像ROCAC那样,虽能动态地调整带宽,但其动态调整的方式却是为了降低平均切换呼叫到达率而过分地牺牲平均新呼叫阻塞率;本方法不仅能动态地调整带宽,而且在平均新呼叫阻塞率和平均切换中断率这两个性能参数之间实现了平衡,达到了降低平均新呼叫阻塞率和平均切换呼叫中断率的目的。


图1.本发明的处理流程图;图2.平均新呼叫阻塞率曲线图;
本方法的第一类业务 本方法的第二类业务 ROCAC的第一类业务 ROCAC的第二类业务 SMCAC的第一类业务 SMCAC的第二类业务图3.平均切换呼叫中断率曲线图; 本方法的第一类业务 本方法的第二类业务 ROCAC的第一类业务 ROCAC的第二类业务 SMCAC的第一类业务 SMCAC的第二类业务图4.本发明的应用示例图。
具体实施例方式
无线/移动网络的最大特点是用户的移动性,用户的移动性带给网络的直接后果是拓扑和资源的动态变化,在具有动态特性的无线/移动网络中,降低新呼叫阻塞率和切换呼叫中断率具有重要的意义对网络来说,可以提高资源利用率;对用户来说,可以提升业务的QoS。根据动态变化的用户移动特性,设计一种能够降低新呼叫阻塞率和切换呼叫中断率的呼叫准入控制方法,是本发明的主要贡献。
该方法的基本思想是将系统时间划分为等长的时段(时段长为20秒),系统在每个时段初为各类业务设置新呼叫准入控制阈值和切换呼叫准入控制阈值,它们分别表示系统在一个时段内最多可接纳的各类业务的新呼叫数或切换呼叫数。系统统计每个时段内已接纳的各类业务的新呼叫数和切换呼叫数,当某类业务的新呼叫(或切换呼叫)向系统发出连接请求时,系统查看接纳该呼叫是否会导致该时段内已接纳的该类业务的新呼叫数(或切换呼叫数)超过该类业务的新呼叫准入控制阈值(或切换呼叫准入控制阈值),如果是,则拒绝该呼叫,否则接纳该呼叫,并为该呼叫分配所需的带宽。为了进行合理的带宽分配,系统根据剩余带宽以及各类业务的平均新呼叫到达率、平均切换呼叫到达率来计算最佳的呼叫准入控制阈值,从而实现减少新呼叫阻塞率和切换呼叫中断率的目的。为此,系统在接纳一个呼叫或当一个呼叫离开时,会更新系统的剩余带宽,同时还会对呼叫到达的历史信息进行统计分析,从而预测各类业务的平均新呼叫到达率和平均切换呼叫到达率。
本发明的总体流程图见图1.
本发明的应用示例图见图4。本发明所述的基于动态阈值的呼叫准入控制方法可应用于图中所示的接入点(AP)、接入路由器(AR)和基站中。
本发明将系统时间分为等长的时段,通过在各时段初设置合理的呼叫准入控制阈值,对到达系统的各类业务的新呼叫或切换呼叫进行是否接纳的判断,实现了在保障已有呼叫的QoS的同时,减少新呼叫阻塞率和切换呼叫中断率的目的,最终提高了系统的性能。
权利要求
1.无线/移动网络中基于动态阈值的呼叫准入控制方法,其特征在于,该方法依次含有以下步骤步骤(1.)在系统即无线/移动网络中的基站、或接入路由器AR、或无线接入点AP中,记录业务类别总数K,K∈N以及剩余带宽Br,并初始化定时器,所述定时器用于将系统的时间分为等长的时段,时段长为20秒;步骤(2.)初始化计数器Ni,i=1,2,…,K为0,用于记录一个时段内系统接纳的第i,i=1,2,…,K类业务的新呼叫数,其中,新呼叫是指由当前无线/移动子网的移动用户发起的呼叫;初始化计数器Hi,i=1,2,…,K,用于记录一个时段内系统接纳的第i,i=1,2,…,K类业务的切换呼叫数,其中,切换呼叫是指由别的无线/移动子网的移动用户发起、并由于用户的移动而切换到当前无线/移动子网中的呼叫;步骤(3.)为了实施呼叫准入控制,系统在每个时段开始时为各类业务的新呼叫和切换呼叫设置呼叫准入控制阈值,其中第i类业务的新呼叫准入控制阈值rni和第i类业务的切换呼叫准入控制阈值rhi,i=1,2,…,K,分别是系统平均收益函数的参数,要在剩余带宽的约束下,使得系统平均收益函数最大化,所述的系统平均收益函数按照以下步骤获得步骤(3.1.)系统预测本时段内第i,i=1,2,…,K类业务的平均新呼叫到达率λni和平均切换呼叫到达率λhi,i=1,2,…,K,所述的平均新呼叫到达率和平均切换呼叫到达率分别指每秒达到系统的新呼叫数和切换呼叫数的平均值,这两个参数分别由系统根据呼叫到达的历史信息经统计分析得到,这里的呼叫到达的历史信息指本时段之前的各时段的新呼叫到达率和切换呼叫到达率的平均值;步骤(3.2.)系统根据本时段内的λni,i=1,2,…,K计算在本时段内到达系统的第i,i=1,2,…,K类业务的新呼叫个数Nni的概率Pn(Nni)=(λni×20)Nnie-(λni×20)/Nni!]]>,i=1,2,…,KNni∈N,0<Nni<∞;系统根据本时段内的λhi,i=1,2,…,K计算在本时段内到达系统的第i类业务的切换呼叫个数的概率Ph(Nhi)=(λhi×20)Nnie-(λhi×20)/Nhi!]]>,i=1,2,…,K,Nhi∈N,0<Nhi<∞;步骤(3.3.)当本时段内到达系统的第i类业务的新呼叫个数Nni,0≤Nni≤rni,i=1,2,…,K且切换呼叫个数Nhi,0≤Nhi≤rhi,i=1,2,…,K时,系统的收益计算为pi×(Nni+Nhi),i=1,2,…,K,所述的pi是当一个第i类业务的呼叫被接纳时系统得到的奖赏值,且第i类业务的优先级越高,pi越大,0<pi<10,i=1,2,…,K;步骤(3.4.)当本时段内到达系统的第i类业务的新呼叫个数Nni,0≤Nni≤rni,i=1,2,…,K且切换呼叫个数Nhi,Nhi>rhi,i=1,2,…,K时,系统的收益计算为pi×(Nni+rhi)-chi×(Nhi-rhi)-αi×(rni-Nni),i=1,2,…,K,其中,chi为拒绝一个第i类业务的切换呼叫时系统得到的惩罚值,且第i类业务的优先级越高,chi越大,0<chi<10,,i=1,2,…,K;而αi为系统为第i类业务的新呼叫分配过多的阈值,从而造成资源浪费时的惩罚值,第i类业务的优先级越高,αi越大,0<αi<10,i=1,2,…,K,;步骤(3.5.)当本时段内到达系统的第i类业务的新呼叫个数Nni,Nni>rni,i=1,2,…,K且切换呼叫个数Nhi,0≤Nhi≤rhi,i=1,2,…,K时,系统的收益计算为pi×(rni+Nhi)-cni×(Nni-rni)-βi×(rhi-Nhi),其中,cni为拒绝一个第i类业务的新呼叫时系统得到的惩罚值且第i类业务的优先级越高,cni越大,0<cni<10,cni<chi,i=1,2,…,K;而βi为系统为第i类的切换呼叫分配过多的阈值,从而造成资源浪费时的惩罚值,且第i类业务的优先级越高,βi越大,0<βi<10,βi<αi,i=1,2,…,K;步骤(3.6.)当本时段内到达系统的第i类业务的新呼叫个数Nni,Nni>rni,i=1,2,…,K且切换呼叫个数Nhi,Nhi>rhi,i=1,2,…,K时,系统的收益计算为pi×(rni+rhi)-cni×(Nni-rni)-chi×(Nhi-rhi),i=1,2,…,K;步骤(3.7.)系统根据步骤(3.2.)得到的Pn(Nni)和Ph(Nhi),i=1,2,…,K以及步骤(3.3.)~步骤(3.6.)计算出的当Nni和Nhi,i=1,2,…,K处于不同范围时的系统的收益,得到系统的平均收益函数 其中,Gi(rni,rhi)为由于第i类业务的呼叫到达给系统带来的平均收益,它是阈值rni和rhi的函数,i=1,2,…,K,其计算为Gi(rni,rhi)=rni(pi+cni)+rhi(pi+chi)-20×(λnicni+λhichi)+ΣNni=0rniΣNhi=0rhi(αi(rni-Nni)+βi(rhi-Nhi))Pn(Nni)Ph(Nhi)+]]>ΣNni=0rni(pi+αi+cni)(Nni-rni)Pn(Nni)+ΣNhi=0rhi(pi+βi+chi)(Nhi-rhi)Ph(Nhi);]]>步骤(3.8.)系统在剩余带宽的约束下,计算使得平均收益函数最大化的rni和rhi,i=1,2,…,K;步骤(4.)系统启动定时器,并进入等待事件到达状态,所述的事件包括新呼叫到达事件、切换呼叫到达事件、呼叫离开事件和定时器超时事件;步骤(5.)在系统处于等待事件到达状态时,如果第i,i=1,2,…,K类业务的新呼叫到达,则系统判断接纳该呼叫是否会使该时段内接纳的第i,i=1,2,…,K类业务的新呼叫数目Ni,i=1,2,…,K超过阈值rni,i=1,2,…,K,如果是,则系统拒绝该呼叫,否则,系统接纳该呼叫,并执行如下操作步骤(5.1.)为该呼叫分配所需的带宽;步骤(5.2.)计数器Ni,i=1,2,…,K加1;步骤(5.3.)剩余带宽Br被更新为当前剩余带宽减去为该呼叫分配的带宽;步骤(5.4.)系统重新进入等待事件到达状态;步骤(6.)在系统处于等待事件到达状态时,如果第i,i=1,2,…,K类业务的切换呼叫到达,则系统判断接纳该呼叫是否会使该时段接纳的第i,i=1,2,…,K类业务的切换呼叫数目Hi,i=1,2,…,K超过阈值rhi,i=1,2,…,K,如果是,则系统拒绝该呼叫,否则,系统接纳该呼叫,并执行如下操作步骤(6.1.)为该呼叫分配所需的带宽;步骤(6.2.)计数器Hi,i=1,2,…,K加1;步骤(6.3.)剩余带宽Br被更新为当前剩余带宽减去为该呼叫分配的带宽;步骤(6.4)系统重新进入等待事件到达状态;步骤(7.)在系统处于等待事件到达状态时,如果有呼叫离开,则系统将剩余带宽Br更新为当前剩余带宽加上该呼叫释放的带宽,然后重新进入等待事件到达状态;步骤(8.)在系统处于等待事件到达状态时,如果定时器超时,则从上述的步骤(2.)重新开始执行。
全文摘要
本发明属于服务质量控制算法领域,其特征在于依次含有以下步骤将系统时间分为等长的时段,在每个时段初设置各类业务的新呼叫准入控制阈值和切换呼叫准入控制阈值。当某类业务的新呼叫或切换呼叫在某个时段内发起连接请求时,若接纳该呼叫会使该时段系统已接纳的该类业务的新呼叫数或切换呼叫数,超过该类业务的新呼叫准入控制阈值或切换呼叫准入控制阈值,则系统拒绝该呼叫,否则接纳该呼叫,并为其分配所需的带宽。为设置合理的呼叫准入控制阈值,使用了“奖惩策略”,即在系统接纳呼叫时给予奖赏,在系统拒绝呼叫或者由于阈值设置不当而造成资源浪费时给予惩罚。本发明实现了减小新呼叫阻塞率和切换呼叫中断率的目的、提高系统性能的目的。
文档编号H04W48/06GK101052232SQ20071009940
公开日2007年10月10日 申请日期2007年5月18日 优先权日2007年5月18日
发明者崔勇, 王胜灵 申请人:清华大学
网友询问留言 已有0条留言
  • 还没有人留言评论。精彩留言会获得点赞!
1