/ 笔记  

【无线通信学习笔记(一)】路径损耗与阴影衰落

Wireless Communications (Andrea Goldsmith)学习笔记。

Path Loss and Shadowing

  1. 主要讲述 path loss(路径损耗) and shadowing(阴影衰弱)相对于距离对接收信号功率的影响。并考虑这两方面的影响建立信道模型。

    • Path loss is caused by dissipation of the power radiated by the transmitter as well as effects of the propagation channel. (be the same at a given transmit-receive distance)
    • Shadowing is caused by obstacles between the transmitter and receiver that attenuate signal power through absorption, reflection, scattering, and diffraction. (When the attenuation is very strong, the signal is blocked)
  2. 大尺度传播效应和小尺度传播效应

    • large-scale propagation effects: occur over relatively large distances. e.g. variations due to path loss and shadowing
    • small-scale propagation effects: occurs over very short distances, on the order of the signal wavelength
  3. 开发了一些近似方法来描述信号传播的特点,而不需要借助麦克斯韦方程。最常见的近似方法是使用 ray-tracing techniques(射线跟踪技术)。

Transmit and Receive Signal Models

  • transmitted signal:

    s(t)={u(t)e2πfct}=x(t)cos(2πfct)y(t)sin(2πfct) s(t)=\Re\left\{u(t)e^{2\pi f_c t}\right\}=x(t)\cos(2\pi f_ct)-y(t)\sin(2\pi f_c t)

    其中,信号 u(t)=x(t)+jy(t)u(t)=x(t)+jy(t) 是复包络,带宽为 BuB_u,,功率为 PuP_u. 发射机发射功率为 P(t)=Pu2P(t)=\frac{ P_u}{2}

  • received signal:

    r(t)={v(t)ej2πfct}={u(t)c(t)ej2πfct} r(t)=\Re{\left\{v(t)e^{j2\pi f_c t}\right\}}=\Re{\left\{u(t)*c(t)\cdot e^{j2\pi f_c t}\right\}}

    其中,c(t)c(t)是等效低通信道的脉冲响应。接收信号可能会存在多普勒频移(Doppler shift)fD=vcosθλf_D=\frac{v\cos\theta}{\lambda}.

Path Loss Models

由free-space path loss出发建立最基础的模型。由于处理RCS过于复杂,在本章中都是从ray tracing的角度来建立复杂模型分析问题的,包括two-ray model, 10-ray model进而推广到一般的general ray tracing。

路径损耗的定义

  • 路径损耗(Path Loss)

PL=PtPrP_L=\frac{P_t}{P_r}

  • 路径增益(Path gain)

    PG dB=PL dBP_G \text{ dB}=-P_L \text{ dB}

自由空间路径损耗

  • 视距信道(LOS,line-of-sight)

    • 没有障碍物
    • 信号沿直线传播
  • 接收信号

    r(t)={λGlej2πd/λ4πdu(t)ej2πfct}r(t)=\Re\left\{\frac{\lambda \sqrt{G_{l}} e^{-j 2 \pi d / \lambda}}{4 \pi d} u(t) e^{j 2 \pi f_{c} t}\right\}

    其中,Gl\sqrt{G_l}是发射天线和接收天线增益的乘积。ej2πd/λe^{-j 2 \pi d / \lambda}是距离dd带来的相移。功率是信号的平方。

    注:这里Gl\sqrt{G_l}是发射天线和接收天线的增益乘积,而非原书的field radiation(直译为场强方向图F(θ,ϕ)F(\theta,\phi)),也并非方向系数DD(在一定意义上可以认为是方向系数)。具体可以看到下面路径增益的表达式实际上就是《天线原理》中的Friis传输方程

    PR=Pi(λ4πr)2DiDrP_R=P_i\left(\frac{\lambda}{4\pi r}\right)^2D_iD_r

    上式是在满足极化匹配、共轭匹配和收发最大辐射方向对准时的表达式。在天线原理书后习题中也被写作

    PR=Pi(λ4πr)2GiGrP_R=P_i\left(\frac{\lambda}{4\pi r}\right)^2G_iG_r

    显然这是在考虑天线效率情况下的Friis传输方程。注意到下面路径增益中将收发天线的增益看作是一样的也就是GlG_l

  • 路径增益

    PrPt=[Glλ4πd]2\frac {P_r}{P_t}=\left[\frac{\sqrt{G_l}\lambda}{4\pi d}\right]^2

    • 接收机接收功率随着收发天线之间距离dd平方成反比。
    • 当载波频率增加,接收功率下降。
  • 自由空间路径损耗

    PL dB=10lgPtPr=10lg[Glλ4πd]2P_L\text{ dB}=10 \lg\frac{P_t}{P_r}=-10\lg\left[\frac{\sqrt{G_l}\lambda}{4\pi d}\right]^2

射线跟踪(Ray tracing)

  • multipath signal components:additional copies of the transmitted signal
  • multipath signal components相对于LOS信道的接收信号有功率衰减,延时,相移和/或频移。
  • 射线追踪技术通过将波面表示为简单的粒子来逼近电磁波的传播。
  • smallest when the receiver is many wavelengths from the nearest scatterer, and all the scatterers are large relative to a wavelength and fairly smooth.
  • 射线跟踪不能很好的反映多径时延扩展(delay spread)

Two-Ray Model

当单一地面反射主导多径效应时,采用Two-Ray Model。
Two-Ray Model

  • 由LOS路径的接受的直射信号和经过x,xx,x'反射信号的叠加,忽略反射面衰减,有

    r2ray (t)={λ4π[Glu(t)ej2πl/λl+RGru(tτ)ej2π(x+x)/λx+x]ej2πfct} r_{2 \text {ray }}(t)=\Re\left\{\frac{\lambda}{4 \pi}\left[\frac{\sqrt{G_{l}} u(t) e^{-j 2 \pi l / \lambda}}{l}+\frac{R \sqrt{G_{r}} u(t-\tau) e^{-j 2 \pi\left(x+x^{\prime}\right) / \lambda}}{x+x^{\prime}}\right] e^{j 2 \pi f_{c} t}\right\}

  • 进一步假设接收信号相对于delay spread(τ\tau)是窄带的,即τBu1\tau\ll B_u^{-1},有u(t)u(tτ)u(t)\approx u(t-\tau),有

    Pr=Pt[λ4π]2Gll+RGrejΔϕx+x2 P_{r}=P_{t}\left[\frac{\lambda}{4 \pi}\right]^{2}\left|\frac{\sqrt{G_{l}}}{l}+\frac{R \sqrt{G_{r}} e^{-j \Delta \phi}}{x+x^{\prime}}\right|^{2}

    其中,Δϕ=2π(x+xl)/λ\Delta\phi=2\pi (x+x'-l)/\lambda是接收到的两个信号的相位差。

  • dd相当大时,x+xld,θ0,GlGr,R=1x+x^{\prime} \approx l \approx d, \theta \approx 0, G_{l} \approx G_{r},R=-1,有

    Pr[λGl4πd]2[4πhthrλd]2Pt=[Glhthrd2]2Pt P_{r} \approx\left[\frac{\lambda \sqrt{G_{l}}}{4 \pi d}\right]^{2}\left[\frac{4 \pi h_{t} h_{r}}{\lambda d}\right]^{2} P_{t}=\left[\frac{\sqrt{G_{l}} h_{t} h_{r}}{d^{2}}\right]^{2} P_{t}

  • 一些说明:

    • 接收功率与λ\lambda无关。

    • 近距离(d<htd<h_t):两路径分量正向相加,路径损耗粗略的是平坦的,接收功率近似与1d2+ht2\frac {1}{d^2+h_t^2}成正比。

    • 远距离(ht<d<dch_t<d<d_c):两路径分量相加此消彼长,造成一系列极大极小值,接收功率近似于d4d^{-4}成正比。

    • d>dcd>d_c:critical distance(dcd_c)时到达了最后一个极大值。超过dcd_c后,接收功率急剧下降。dcd_c可由Δϕ=π\Delta \phi=\pi估计,即dc=4hthrλd_c=\frac{4h_th_r}{\lambda}

    • 蜂窝系统(cellular systems)中传播很少符合two-ray model。

      Received Power versus Distance for Two-Ray Model

Ten-Ray Model (Dielectric Canyon)

当单一路径反射超过三次可以被忽略。后向反射也可以被忽略。

Overhead View of the Ten-Ray Model
  • Ten-Ray Model下的接收信号为

    r10ray (t)={λ4π[Glu(t)ej2πl/λl+i=19RiGxiu(tτi)ej2πxi/λxi]ej2πfct}r_{10 \text {ray }}(t)=\Re\left\{\frac{\lambda}{4 \pi}\left[\frac{\sqrt{G_{l}} u(t) e^{-j 2 \pi l / \lambda}}{l}+\sum_{i=1}^{9} \frac{R_{i} \sqrt{G_{x_{i}}} u\left(t-\tau_{i}\right) e^{-j 2 \pi x_{i} / \lambda}}{x_{i}}\right] e^{j 2 \pi f_{c} t}\right\}

  • 进一步假设接收信号相对于delay spread(τ\tau)是窄带的,有u(t)u(tτ)u(t)\approx u(t-\tau),有Ten-Ray Model的接收功率

    Pr=Pt[λ4π]2Gll+i=19RiGxiejΔϕixi2 P_{r}=P_{t}\left[\frac{\lambda}{4 \pi}\right]^{2}\left|\frac{\sqrt{G_{l}}}{l}+\sum_{i=1}^{9} \frac{R_{i} \sqrt{G_{x_{i}}} e^{-j \Delta \phi_{i}}}{x_{i}}\right|^{2}

  • Ten-Ray Model的接收功率随d2d^{2}下降

General Ray Tracing(GRT)

  • 对于GRT模型,建筑物参数(高度、位置和介电性质)以及发射机和接收机相对于建筑物的位置必须准确确定。

  • 一般来说,由于绕射(diffraction)和散射(scattering)的损耗很大,所以在接收信号中起主导作用的十LOS信号和反射信号。

  • diffraction(绕射)分量可由Fresnel knife edge diffraction model求得,接收信号为

    r(t)={L(v)Gdu(tτ)ej2π(d+d)/λej2πfct}r(t)=\Re\left\{L(v) \sqrt{G_{d}} u(t-\tau) e^{-j 2 \pi\left(d+d^{\prime}\right) / \lambda} e^{j 2 \pi f_{c} t}\right\}

    其中,Gd\sqrt{G_d}是天线增益,τ=Δdc\tau=\frac{\Delta d}{c}绕射分量相对于直射分量的时延

  • scattered ray(散射)分量可有双基地雷达方程(bistatic radar equation)获得,接收信号为

  • 结合上述各种分量的影响,假设有一个LOS路径,NrN_r个反射路径,NdN_d个绕射路径,NsN_s个散射路径,则总的接收信号为

    rtotal (t)={[λ4π][Glu(t)ej2πl/λl+i=1NrRxiGxiu(tτi)ej2πxi/λxi+j=1NdLj(v)Gdju(tτj)ej2π(dj+dj)/λ+k=1NsGskσku(tτk)ej2π(sk+sk)/λsksk]ej2πfct}\begin{aligned} r_{\text {total }}(t) &=\Re\left\{\left[ \frac { \lambda } { 4 \pi } \right] \left[\frac{\sqrt{G_{l}} u(t) e^{j 2 \pi l / \lambda}}{l}+\sum_{i=1}^{N_{r}} \frac{R_{x_{i}} \sqrt{G_{x_{i}}} u\left(t-\tau_{i}\right) e^{-j 2 \pi x_{i} / \lambda}}{x_{i}}\right.\right.\\ &+\sum_{j=1}^{N_{d}} L_{j}(v) \sqrt{G_{d_{j}}} u\left(t-\tau_{j}\right) e^{-j 2 \pi\left(d_{j}+d_{j}^{\prime}\right) / \lambda} \\ &\left.\left.+\sum_{k=1}^{N_{s}} \frac{\sqrt{G_{s_{k}} \sigma_{k}} u\left(t-\tau_{k}\right) e^{j 2 \pi\left(s_{k}+s_{k}^{\prime}\right) / \lambda}}{s_{k} s_{k}^{\prime}}\right] e^{j 2 \pi f_{c} t}\right\} \end{aligned}

    如果其传播路径被建筑物或其他物体阻挡,任何一个多径分量都可能有一个额外的衰减系数。

  • local mean received power(Pˉr\bar P_r,本地平均接收功率):射线跟踪中所有路径的幅度平方和,可以较好地反映link quality(链路质量)。

Empirical Path Loss Model

为了消除多径效应,路径损耗的经验测量通常将其接收功率测量值和相应的路径损耗在给定距离上的几个波长上求平均值。这个平均路径损耗称为距离dd处的本地平均衰减(LMA),由于自由空间的路径损耗和信号障碍物的影响,一般会随着dd的增大而减小。

  • The Okumura Model: One of the most common models for signal prediction in large urban macrocells is the Okumura model . This model is applicable over distances of 1-100 km and frequency ranges of 150-1500 MHz. The base station heights
    for these measurements were 30-100 m, the upper end of which is higher than typical base stations today

    PL(d)dB=L(fc,d)+Amu(fc,d)G(ht)G(hr)GAREAP_{L}(d) \mathrm{dB}=L\left(f_{c}, d\right)+A_{m u}\left(f_{c}, d\right)-G\left(h_{t}\right)-G\left(h_{r}\right)-G_{A R E A}

  • Hata Model: This empirical model simplifies calculation of path loss since it is a closed-form formula and is not based on empirical curves for the different parameters. the same range of frequencies, 150-1500 MHz.

    PL,urban(d)dB=69.55+26.16log10(fc)13.82log10(ht)a(hr)+(44.96.55log10(ht))log10(d)P_{L, u r b a n}(d) \mathrm{dB}=69.55+26.16 \log _{10}\left(f_{c}\right)-13.82 \log _{10}\left(h_{t}\right)-a\left(h_{r}\right)+\left(44.9-6.55 \log _{10}\left(h_{t}\right)\right) \log _{10}(d)

  • COST 231 Extension to Hata Model: The Hata model was extended to 2 GHz .

    PL,urban(d)dB=46.3+33.9log10(fc)13.82log10(ht)a(hr)+(44.96.55log10(ht))log10(d)+CMP_{L, u r b a n}(d) \mathrm{dB}=46.3+33.9 \log _{10}\left(f_{c}\right)-13.82 \log _{10}\left(h_{t}\right)-a\left(h_{r}\right)+\left(44.9-6.55 \log _{10}\left(h_{t}\right)\right) \log _{10}(d)+C_{M}

    This model is referred to as the COST 231 extension to the Hata model, and is restricted
    to the following range of parameters: 1.5GHz < fcf_c < 2 GHz, 30m < hth_t < 200 m, 1m < hrh_r < 10 m, and 1km < dd < 20 km.

  • Piecewise Linear (Multi-Slope) Model

  • Indoor Attenuation Factors: it is difficult to find generic models that can be accurately applied to determine empirical path loss in a specific indoor setting.

Simplified Path Loss Model

  • 对于一般的衰减分析,可以使用路径损耗的简化模型:

    Pr=PtK[d0d]γP_r=P_tK\left[\frac{d_0}{d}\right]^\gamma

    用分贝表示

    PrdBm=PtdBm+K dB10γlog10[dd0]P_{r} \mathrm{dBm}=P_{t} \mathrm{dBm}+K \mathrm{~dB}-10 \gamma \log _{10}\left[\frac{d}{d_{0}}\right]

    其中,KK是一个无量纲常数,由天线特性和平均信道特性决定,d0d_0是天线远场的参考距离,γ\gamma是路径衰落指数。

  • KK的取值:为了简化计算,将K<1K<1取为全向天线自由空间内在d0d_0处的路径增益:

    K dB=20log10λ4πd0K \mathrm{~dB}=20 \log _{10} \frac{\lambda}{4 \pi d_{0}}

    这里假设天线效率ηA\eta_A为1,则天线增益G=ηADG=\eta_A D与天线方向系数DD就是一致的。

    对于全向辐射的假想元:

    D=4πo2π0πsinθdθdϕ=1=0 dBiD=\frac{4\pi}{\int_o^{2\pi}\int_{0}^{\pi}\sin \theta \mathrm{d}\theta\mathrm{d}\phi}=1=0\text{ dBi}

    代入路径增益的表达式即为上式。

  • γ\gamma的取值:γ\gamma在自用空间一般取2,在Two-Ray model中取4,对于更复杂的环境可通过对实测数据进行MMSE(minimum mean square error)拟合。

    F(γ)=i=1n[MmeasuredMmodel]2F(\gamma)=\sum_{i=1}^{n}\left[M_\text{measured}-M_\text{model}\right]^2

    其中MmeasuredM_\text{measured}是实测的10lgPrPt dB10\lg\frac{P_r}{P_t} \text{ dB},而MmodelM_\text{model}是上述路径损耗的简化模型计算得到的10lgPrPt dB=K dB10γlog10[dd0]10\lg\frac{P_r}{P_t} \text{ dB}=K \mathrm{~dB}-10 \gamma \log _{10}\left[\frac{d}{d_{0}}\right]。对F(γ)F(\gamma)求极值点:

    F(γ)γ=0γ\frac{\partial F(\gamma)}{\gamma}=0\to\gamma

    也可以通过经验数据获得:

Environment γ range
Urban macrocells 3.7-6.5
Urban microcells 2.7-3.5
Office Building (same floor) 1.6-3.5
Office Building (multiple floors) 2-6
Store 1.8-2.2
Factory 1.6-3.3
Home 3

Shadow Fading

通过无线信道传输的信号通常会因为信号路径中物体的阻挡而发生随机变化,从而引起给定距离内接收功率的随机变化。由于参数是未知的,所以只能通过统计模型来表征这种随机变化。

log-normal shadowing model

  • 最常用来描述这种加性衰减的模型是log-normal shadowing。

  • 在log-normal shadowing model中,发射接收功率比ψ=PtPr\psi=\frac{P_t}{P_r}是服从对数正态分布的

    p(ψ)=ξ2πσψdBψexp[(10log10ψμψdB)22σψdB2],ψ>0p(\psi)=\frac{\xi}{\sqrt{2 \pi} \sigma_{\psi_\text{dB}} \psi} \exp \left[-\frac{\left(10 \log _{10} \psi-\mu_{\psi_\text{dB}}\right)^{2}}{2 \sigma_{\psi_\text{dB}}^{2}}\right], \psi>0

    其中,ξ=10ln10\xi=\frac{10}{\ln 10}ψdB=10lgψ\psi_\text{dB}=10\lg\psi。同时,μψdB,σψdB\mu_{\psi_\text{dB}},\sigma_{\psi_\text{dB}}都是ψdB\psi_\text{dB}按分贝计算的均值和方差。

  • ψ=PtPr\psi=\frac{P_t}{P_r}服从对数正态分布,则接收功率和接收信噪比都服从对数正态分布。

  • 注意:在经验路径损耗模型中,由于无法拆分,μψdB\mu_{\psi_\text{dB}}已经包含Shadowing的影响。而对于理论分析,必须单独考虑Shadowing的影响。

  • 通常我们对测量值的分贝值取平均来确定平均的路径损耗和方差。

correlation of shadowing over distance

  • 假设Shadowing ψ(d)\psi(d)是一个一阶自回归过程,相距δ\delta两点Shadowing fading(阴影衰落)的协方差为

    A(δ)=E[(ψdB(d)μψdB)(ψdB(d+δ)μψdB)]=σψdB2ρDδ/DA(\delta)=\mathrm{E}\left[\left(\psi_{d B}(d)-\mu_{\psi_{d B}}\right)\left(\psi_{d B}(d+\delta)-\mu_{\psi_{d B}}\right)\right]=\sigma_{\psi_{d B}}^{2} \rho_{D}^{\delta / D}

    其中,ρD\rho_D是相距为DD时的归一化协方差,且必须通过测量得到。

  • 假设固定D=XcD=X_c时,ρD=1e\rho_D=\frac 1e,则模型可简化为

    A(δ)=σψdB2eδ/XcA(\delta)=\sigma^2_{\psi_\text{dB}}e^{-\delta/X_c}

    其中,XcX_c是decorrelation distance(去相关距离),指信号自相关等于其最大值的1e\frac 1e时的距离,并且与障碍物或障碍物群的大小成正比。

  • 当用户以速度vv运动,δ=vτ\delta=v\tau代入上式可求得在τ\tau时刻的阴影-去相关关系。

Combined Path Loss and Shadowing

the ratio of received to transmitted power

综合上述路径损耗的简化模型和对数正态分布的阴影衰减模型,可得接收发射功率比为

PrPt(dB)=10log10K10γlog10dd0ψdB\frac{P_{r}}{P_{t}}\text{(dB)}=10 \log _{10} K-10 \gamma \log _{10} \frac{d}{d_{0}}-\psi_{d B}

  • ψdBN(0,σψdB2)\psi_\text{dB}\sim N(0,\sigma^2_{\psi_\text{dB}})
  • 路径损耗随着lgd\lg d10γ dB/decade10\gamma\text{ dB/decade}的斜率线性递减
  • 阴影衰减导致的变化非常快,和去相关距离XcX_c在同一数量级

Outage Probability

无线通信中有目标允许最小接收功率PminP_{\min}

  • 中断概率(Outage Probability):在距离dd上,接收功率Pr(d)P_r(d)衰减到小于一个服务质量不可接受的临界值PminP_{\min}的概率:

    pout(Pmin,d)=p(Pr(d)<Pmin)p_\text{out}(P_{\min},d)=p\left(P_r(d)<P_{\min}\right)

  • 对于路径损耗的简化模型和对数正态分布的阴影衰减模型,上式可写为

    p(Pr(d)Pmin)=1Q(Pmin(Pt+10log10K10γlog10(d/d0))σψdB)p\left(P_{r}(d) \leq P_{\min }\right)=1-Q\left(\frac{P_{\min }-\left(P_{t}+10 \log _{10} K-10 \gamma \log _{10}\left(d / d_{0}\right)\right)}{\sigma_{\psi_\text{dB}}}\right)

    其中,

    Q(z)p(x>z)=z12πey2/2dy=12erfc(z2)Q(z) \triangleq p(x>z)=\int_{z}^{\infty} \frac{1}{\sqrt{2 \pi}} e^{-y^{2} / 2} d y=\frac 12\mathrm{erfc}\left(\frac{z}{\sqrt{2}}\right)

Cell Coverage Area

Contours of Constant Received Power
  • 考虑路径损耗和平均的阴影衰减,等功率线是以基站(Base station, BS)为圆心的圆。而考虑路径损耗和随机的阴影损耗,等功率线就变成一个不规则图形。

  • 不能保证在cell的边界上的人接收相同的功率量级,两种解决方式

    • transmit extra power to insure users affected by shadowing receive their minimum required powerPminP_{\min}(会影响邻近的cell)
    • some users within the cell will not meet their minimum received power requirement.
  • outage probability of the cell(小区中的中断概率):小区内未达到最低功率要求的面积百分比

    poutcell=1Cp_\text{out}^\text{cell}=1-C

    C=1πR2cell area PAdA=1πR202π0RPArdrdθC=\frac{1}{\pi R^{2}} \int_{\text {cell area }} P_{A} d A=\frac{1}{\pi R^{2}} \int_{0}^{2 \pi} \int_{0}^{R} P_{A} r d r d \theta

    其中,PAP_A表示在元dA\mathrm{d}A内接收功率达到最小接收功率的概率,即PA=E[1[Pr(r)>Pmin in dA]]P_{A}=\mathrm{E}\left[1\left[P_{r}(r)>P_{\min } \text { in } d A\right]\right]1[]1[\bullet]是示性函数。

  • 对于路径损耗的简化模型和对数正态分布的阴影衰减模型,有

    p(Pr(r)>Pmin)=1pout(Pmin,r)p(P_{r}(r)>P_{\min })=1-p_\text{out}(P_{\min},r)

    小区中接收功率低于PminP_{\min}的位置被称为outage location(中断位置)。

  • 若目标最小接收功率PminP_{\min}等于小区边界处的接收功率Pˉr\bar P_r,再次假设下小区边界上的中断概率为0.5.

在整理的时候,真的发现好多和《天线原理》和《雷达原理与系统》交叉的内容。

参考文献

  1. ANDREA GOLDSMITH.WIRELESS COMMUNICATIONS.2005.
  2. (美)Andrea Goldsmith著;杨鸿文,李卫东,郭文彬等译. 无线通信. 北京:人民邮电出版社, 2007.06.
  3. 魏文元等编. 天线原理. 北京:国防工业出版社, 1985.06.
  4. 无线通信学习笔记(二) 大尺度路径损耗.https://www.jianshu.com/p/5e9141e6cd76

本文标题:【无线通信学习笔记(一)】路径损耗与阴影衰落

文章作者:Levitate_

发布时间:2021年04月15日 - 22:05:31

原始链接:https://levitate-qian.github.io/2021/04/16/Wireless-Communications-Ch-2/

许可协议: 署名-非商业性使用-禁止演绎 4.0 国际 转载请保留原文链接及作者。