电波传播理论(徐立勤)db-1_第1页
电波传播理论(徐立勤)db-1_第2页
电波传播理论(徐立勤)db-1_第3页
电波传播理论(徐立勤)db-1_第4页
电波传播理论(徐立勤)db-1_第5页
已阅读5页,还剩28页未读 继续免费阅读

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

电波传播理论基于徐立勤《电磁场与电磁波理论》第三版·系统讲授电磁波辐射、传播与导行Contents课程目录从数学基础到工程应用,四段式构建电波传播理论的完整知识体系。01数学基础与电磁基本定律02静态场理论与边值问题03电磁波辐射与空间传播04导行电磁波与工程应用CHAPTER01数学基础与电磁基本定律矢量分析、场论基础与宏观电磁现象的基本规律Chapter02·Fundamentals矢量分析与坐标系基础矢量分析是电磁场理论的数学语言。标量场描述空间中的单值分布(如电位),矢量场描述有方向的物理量(如电场强度),而选择合适的正交坐标系是简化场方程求解的关键第一步。三维坐标系与矢量场教学模型标量场与梯度标量场在空间每点赋予一个实数值(如电位φ、电荷密度ρ),其空间变化率由梯度算子完整描述矢量场·散度与旋度矢量场在空间每点赋予大小和方向(如E、H),需通过散度和旋度两个独立算子刻画其源特性三种正交坐标系直角坐标系适合平面边界,柱坐标系适合同轴电缆等圆柱结构,球坐标系适合点源辐射等球对称场景坐标系选择原则使边界条件与坐标曲面吻合,将三维偏微分方程降维为可解的常微分方程VectorCalculus梯度、散度与旋度:场论三大算子梯度揭示标量场的最大变化方向,散度描述矢量场的源汇特性,旋度刻画矢量场的涡旋结构。三大算子与两大积分定理构成电磁场微分-积分对偶分析的完整框架。Gradient梯度(Gradient)方向导数极值:gradφ指向标量场增长最快的方向,其模等于该方向上的方向导数最大值电场关联:E=−gradφ,电位场的负梯度将标量电位与矢量电场直接关联E=−∇φDivergence散度(Divergence)源汇判定:divF表示单位体积内矢量场的净流出量,正值为源、负值为汇、零为无源场高斯定理:散度定理将体积分与闭合面积分等价,是积分与微分形式转换的核心桥梁∮F·dA=∫∇·FdVCurl旋度(Curl)涡旋源:curlF描述矢量场的微观旋转强度与旋转轴方向,非零旋度表示存在涡旋源斯托克斯定理:将曲面积分与边界线积分等价,奠定安培环路定律和法拉第定律基础∮F·dl=∫∇×F·dSELECTROMAGNETICTHEORY麦克斯韦方程组:电磁学的统一理论麦克斯韦方程组由四个偏微分方程构成,完整描述了电荷、电流与电磁场之间的因果关系。麦克斯韦引入位移电流假设使方程组自洽,由此推导出电磁波方程并预言光即为电磁波,实现了电、磁、光三大现象的历史性统一。JamesClerkMaxwell,1831–187901高斯电场定律(∇·D=ρ):电荷是电场的源,闭合曲面的电通量等于其包围的总自由电荷02高斯磁场定律(∇·B=0):不存在孤立磁荷,磁力线永远闭合,磁场是无散度场03法拉第定律(∇×E=-∂B/∂t):时变磁场产生涡旋电场,这是电磁感应与变压器工作的物理基础04安培-麦克斯韦定律(∇×H=J+∂D/∂t):传导电流和位移电流共同产生磁场,位移电流项是预言电磁波存在的关键ElectromagneticMediumParameters电磁场基本物理量与媒质参数介电常数ε、磁导率μ、电导率σ是描述媒质电磁特性的三大基本参数。它们共同决定了电磁波在媒质中的传播速度v=1/√(με)、波阻抗η=√(μ/ε)和衰减常数α,是连接电磁理论与工程应用的桥梁。介电常数ε描述材料极化储能能力,真空ε₀≈8.854×10⁻¹²F/m,常见介质εᵣ从2(聚乙烯)到80(水)Permittivityεᵣ:2–80磁导率μ描述材料导磁能力,真空μ₀=4π×10⁻⁷H/m,铁磁材料μᵣ可达数千至数万Permeabilityμᵣ≤10⁴电导率σ描述材料导电能力,良导体铜σ≈5.8×10⁷S/m,绝缘体σ<10⁻¹²S/m,跨越近20个数量级Conductivity×10²⁰量级本构关系D=εE、B=μH、J=σE将场量与媒质参数联系,是求解麦克斯韦方程组的必要补充条件ConstitutiveRelationsDBJCHAPTER02静态场理论与边值问题静电场、恒定电场与恒定磁场的分析方法及工程应用ELECTROSTATICS静电场:从库仑定律到电位理论静电场由静止电荷产生,遵循库仑定律的平方反比律。通过引入标量电位φ将矢量电场问题降维,利用泊松方程∇²φ=-ρ/ε统一求解,是静电场分析的核心方法论,也是后续所有边值问题解法的出发点。库仑定律F=kq₁q₂/r²描述点电荷间的作用力,是静电学的实验基础,其平方反比律与万有引力定律形式相同,奠定了静电场理论的基本框架F=kq₁q₂/r²电场强度E定义为E=F/q₀(试验电荷所受力除以电荷量),连续电荷分布通过体积分叠加计算各点场强,是描述电场力的空间分布特性E=F/q₀电位φ是标量辅助函数,满足E=-∇φ,将三维矢量场问题简化为标量场问题,大幅降低求解复杂度,引入电位是静电场分析的关键技巧E=-∇φ泊松方程∇²φ=-ρ/ε和拉普拉斯方程∇²φ=0是静电场的控制方程,配合边界条件可唯一确定空间电位分布,是边值问题求解的数学基础∇²φ=-ρ/εBoundaryValueProblems静电场边值问题:镜像法与分离变量法镜像法通过引入虚拟镜像电荷替代导体边界效应,分离变量法将偏微分方程分解为独立常微分方程求解。两种方法分别适用于特殊对称边界和规则几何边界,是求解拉普拉斯方程和泊松方程的两大经典工具。01·镜像法用假想镜像电荷替代导体边界,使原区域电位分布等价,导体平面处恰好满足零电位边界条件适用于点电荷与无限大接地导体平面、导体球等,可直接计算感应电荷分布和作用力仅适用于边界具有高度对称性的少数情况,复杂几何边界需借助数值方法求解02·分离变量法基本假设φ=X(x)Y(y)Z(z)代入拉普拉斯方程,利用边界条件确定分离常数和待定系数坐标系由边界几何决定:矩形用直角坐标,圆柱用柱坐标,球面用球坐标通解为无穷级数(傅里叶级数或贝塞尔函数展开),由边界条件确定各阶系数CHAPTER02·静态场恒定电场与恒定磁场恒定电场和恒定磁场构成电磁学的'静态双子星'。两者分别呈现'有散无旋'和'有旋无散'的对偶结构,为理解时变场的耦合奠定基础。01恒定电场由稳恒电流产生,满足J=σE和∇·J=0,在均匀导电媒质内部同样遵循拉普拉斯方程J=σE静电场分析方法(镜像法、分离变量法)可直接迁移,边界条件由电位连续和法向电流连续给出SeparationofVariables02恒定磁场毕奥-萨伐尔定律通过积分计算任意电流分布产生的磁场dB=μ₀Idl×r̂/(4πr²)安培环路定律在高对称情况下(直导线、螺线管、同轴电缆)可直接求解磁场分布∮H·dl=I磁场是有旋无散场:无磁荷且电流是磁场的旋度源,与静电场形成完美对偶∇·B=0∇×H=JChapter03电磁波辐射与空间传播从辐射源到远场传播,揭示电磁波的产生机制与空间行为RADIATIONTHEORY电磁波辐射原理与赫兹偶极子电磁波辐射源于加速电荷产生的时变电磁场向外传播。赫兹偶极子是最基本的辐射单元,其辐射场在近场区以储能为主(场量按1/r²和1/r³衰减),在远场区才形成真正的辐射场(场量按1/r衰减,功率按1/r²衰减),远场条件为r>>λ/(2π)。01波动方程∇²E−με∂²E/∂t²=0从麦克斯韦方程组推导而来,证明电磁扰动以v=1/√(με)的速度在空间中传播02赫兹偶极子是长度l<<λ的短导线辐射模型,其辐射电阻Rᵣ=80π²(l/λ)²Ω,表明短天线的辐射效率极低03近场区(r<<λ/2π)电场和磁场相位差90°,能量在源与场之间周期性交换而不向外传播,以储能为主04远场区(r>>λ/2π)电场E与磁场H互相垂直且同相,波阻抗η=E/H=√(μ/ε),能量沿径向单向传播偶极子天线辐射场实验室实物Radiation&Gain天线辐射方向性与增益天线增益是通信系统设计的核心参数,它综合反映了天线将辐射能量集中到特定方向的能力(方向性)和天线本身的损耗效率。增益每提高3dB,等效辐射功率翻倍,在相同条件下通信距离可增加约41%。01D方向性函数F(θ,φ)描述辐射功率密度的空间分布,方向性系数D=max[F(θ,φ)]衡量能量集中程度02ηₐ·D增益G=ηₐ·D将天线效率ηₐ与方向性结合,工程中用dBi表示:G(dBi)=10log₁₀(G),全向天线增益为0dBi03+3dB增益每提高3dB等效辐射功率翻倍,相同发射功率下通信距离增加√2≈1.41倍,这是链路预算的关键参数04权衡实际天线设计需在增益、波束宽度、旁瓣电平和物理尺寸之间权衡,不同应用场景(广播/雷达/卫星)需求各异ELECTROMAGNETICTHEORY均匀平面波在自由空间中的传播均匀平面波是远场电磁波的理想化模型,其等相位面为平面且面上各点场量相同。在自由空间中,平面波为TEM波(E、H、传播方向三者互相垂直),波阻抗η₀≈377Ω,传播速度为光速c,是最基础也是最重要的电磁波传播分析起点。01均匀平面波的等相位面是平面,每个面上电场和磁场的幅值与方向完全相同,是远场球面波的局部近似02TEM波特性:电场E⊥磁场H⊥传播方向k,三者构成右手螺旋关系,坡印廷矢量S=E×H指向能量传播方向03自由空间波阻抗η₀=√(μ₀/ε₀)≈120π≈377Ω,是电场与磁场幅值之比,也是电磁学最基本的常数之一04传播速度v=c=1/√(μ₀ε₀)≈3×10⁸m/s,波长λ=c/f,例如900MHz信号波长约33cm,2.4GHz信号波长约12.5cmTEM波传播示意·电场、磁场与传播方向三者正交关系WaveAttenuation平面波在有耗媒质中的衰减与趋肤效应有耗媒质使电磁波幅值随传播距离指数衰减。良导体中趋肤深度δ极小(如铜在1GHz仅2.1μm),高频电流被限制在表面薄层内流动;介质的损耗角正切tanδ决定电磁波在介质中的衰减速率,这直接影响通信链路设计和材料选择。趋肤深度公式δ=√(2/ωμσ),铜在1GHz时δ≈2.1μm,10GHz时仅≈0.66μm,高频电流集中在导体表面极薄层2.1μm电磁屏蔽效应良导体衰减常数α≈√(πfμσ),电磁波进入导体后迅速衰减,几个趋肤深度后场量几乎为零,形成电磁屏蔽α≈√(πfμσ)介质损耗机制由tanδ=σ/(ωε)描述,tanδ越大衰减越快;水在2.45GHz的tanδ较大,这是微波加热的物理基础2.45GHz传播常数γ=α+jβ综合描述衰减和相位变化,α为衰减常数(Np/m),β为相位常数(rad/m),vₚ=ω/β为相速度γ=α+jβPOLARIZATION电磁波极化:线极化、圆极化与椭圆极化极化描述电场矢量端点的时域轨迹,分为线极化、圆极化和椭圆极化三种形式。极化匹配是无线通信的关键——收发天线极化不一致会产生极化损耗,严重时可导致通信中断。线极化电场矢量始终沿一条固定直线振动,两个正交分量同相或反相,大多数地面广播天线采用线极化水平与垂直极化正交,理论隔离度无穷大,实际可达20–30dB,可用于频率复用20–30dB圆极化两个正交分量幅值相等且相位差±90°,电场矢量端点画出圆,旋向分左旋和右旋对天线姿态旋转不敏感,卫星通信和GPS广泛采用;反射后旋向反转,可抑制多径干扰AR=1椭圆极化最一般的极化形式,两正交分量幅值不等或相位差非90°,线极化和圆极化均为其特例轴比AR描述椭圆胖瘦程度,AR=1为圆极化,AR=∞为线极化,工程上AR<3dB视为合格AR<3dBWAVEPROPAGATION分界面上的反射与折射:菲涅尔公式电磁波在媒质分界面发生反射和折射,其行为由菲涅尔公式精确描述。布儒斯特角使平行极化波反射系数为零,全反射将电磁波约束在光密媒质内——这两个现象分别是偏光技术和光纤通信的物理基础。光纤中全反射光信号传输演示01斯涅尔定律n₁sinθ₁=n₂sinθ₂确定反射角和折射角方向,折射率n=√(εᵣμᵣ)决定波在不同媒质中的传播速度02菲涅尔反射系数取决于入射角和极化方式:TE波(垂直极化)和TM波(平行极化)的反射系数公式不同03布儒斯特角θB=arctan(n₂/n₁)使TM波反射系数为零,全部能量透射——偏光太阳镜利用此原理消除水面反射眩光04全反射条件:n₁>n₂且θ₁>θc=arcsin(n₂/n₁),光纤通信利用全反射将光信号约束在纤芯内实现低损耗长距离传输WAVEPROPAGATIONMECHANISMS实际电波传播的四种基本机制实际无线信道中电波传播由直射、反射、绕射和散射四种机制叠加构成,理解各机制的物理特征是无线信道建模的基础。城市环境中无线信号的多径传播场景直射直射波遵循自由空间路径损耗PL=32.44+20lgf+20lgd(dB),是视距传播的主要分量山区环境中信号绕射覆盖示意反射反射良导体全反射无衰减,绝缘体部分反射部分折射,非理想介质产生穿透衰落散射散射大量小于波长的小物体使信号能量分散到各方向,形成漫射分量Free-SpacePathLoss自由空间传播损耗与弗里斯传输公式弗里斯传输公式Pᵣ=PₜGₜGᵣ(λ/(4πd))²是无线通信链路预算的基石。每增加一倍距离损耗增加6dB,每增加一倍频率损耗也增加6dB。01弗里斯公式Pᵣ=PₜGₜGᵣ(λ/(4πd))²量化远场条件下接收功率与发射功率、天线增益、频率和距离的关系02分贝形式:Pᵣ=Pₜ+Gₜ+Gᵣ−32.44−20lgf(MHz)−20lgd(km),距离或频率每翻倍损耗增加6dB+6dB/倍频程03900MHz基站(43dBm/15dBi)覆盖1km处手机(0dBi),接收功率约−70.5dBm,链路余量约30dB−70.5dBm04自由空间损耗为理论下限,实际城市环境附加损耗20–40dB,需通过Okumura-Hata等实测模型修正20–40dB不同频率下自由空间路径损耗随距离变化相同距离下频率越高路径损耗越大,1800MHz比900MHz多约6dB损耗PropagationTheory高频与低频信号的室内传播差异高频信号穿透能力强但室内分布不均匀,低频信号穿透能力弱但反射后分布均匀——5G高频段需采用室分系统或小基站补充覆盖。低频800MHz穿透能力较弱,进入室内后在墙壁间多次反射,场强分布相对均匀,用户体感信号稳定,适合广域覆盖场景。800MHz高频1800MHz穿透能力强,部分信号穿透墙壁后又穿出室外,导致室内场强分布不均匀,用户体感信号波动大,覆盖质量不稳定。1800MHz室内信号是穿透分量与绕射分量的叠加,绕射分量占绝大部分;高频信号室内外电平差显著大于低频信号,衰减特性差异明显。DiffractionDominant5G高频段(3.5GHz/4.9GHz)室内覆盖挑战更大,需部署室分系统(DAS)或数字化小基站实现深度覆盖,保障用户体验。3.5/4.9GHzCHAPTER04导行电磁波与工程应用金属波导、传输线与史密斯圆图——从理论到微波工程实践WaveguideModes矩形波导中的TE/TM模式与截止频率矩形波导中电磁波只能以离散的TEₘₙ/TMₘₙ模式传播,每种模式有确定的截止频率fc。工作频率高于fc时模式可传播,低于fc时迅速衰减。工程设计通常让波导工作在主模TE₁₀的单模区域,以避免多模干扰并确保信号传输质量。矩形金属波导实物01TEₘₙ与TMₘₙ模—横电模电场无传播方向分量Ez=0;横磁模磁场无传播方向分量Hz=002截止频率公式—fc(m,n)=(c/2)√[(m/a)²+(n/b)²],a为宽边、b为窄边,m和n为半波数03主模TE₁₀—截止频率fc=c/(2a)最低;单模条件:fc(TE₁₀)<f<fc(下一高阶模)04波导波长与速度—λg=λ/√[1-(fc/f)²]>λ;相速度vp>c,群速度vg<cWaveguideStructures圆形波导、同轴线与微带线不同导波结构各有优势:矩形波导适合高功率微波系统,圆形波导TE₀₁模适合低损耗长距离传输,同轴线支持TEM模可传输从直流到高频的宽带信号,微带线则实现了微波电路的平面化集成,是现代射频PCB的基础。同轴电缆与射频连接器实物圆形波导与同轴线01圆形波导主模TE₁₁,TE₀₁模衰减随频率升高而单调下降,适合长距离低损耗传输和旋转关节应用02同轴线传播TEM模(fc=0),可传输直流到高频宽带信号;标准阻抗50Ω兼顾功率容量和损耗最优微带线与平面电路01微带线由介质基板上的导带和底部接地板构成,准TEM模传播,特性阻抗由导带宽度和介质厚度比决定02微带线易于与有源器件集成、可光刻批量生产,是现代射频PCB和微波集成电路(MMIC)的基础导波结构微波电路板PCB射频微带线实物ELECTROMAGNETICTHEORY传输线理论与电报方程当传输线长度与波长可比时,必须采用分布参数模型(R/L/G/C),特性阻抗Z₀是行波的电压电流比,传播常数γ决定衰减和相移。01·MODEL分布参数模型将传输线视为无穷多微小RLGC段的级联,适用于线长l与波长λ可比(l>λ/10)的高频场景。l>λ/1002·EQUATION电报方程∂V/∂z=−(R+jωL)I和∂I/∂z=−(G+jωC)V是传输线的基本微分方程,解为入射波与反射波的叠加。入射+反射03·IMPEDANCE特性阻抗Z₀=√[(R+jωL)/(G+jωC)],无耗线简化为Z₀=√(L/C),是纯实数且与频率无关。Z₀=√(L/C)04·PROPAGATION传播常数γ=α+jβ=√[(R+jωL)(G+jωC)],α为衰减常数(Np/m),β=ω√(LC)为相位常数(rad/m)。γ=α+jβELECTROMAGNETICTHEORY传输线反射、驻波与阻抗匹配负载与传输线阻抗不匹配时产生反射,反射系数Γ和驻波比VSWR是衡量匹配质量的核心指标。VSWR=1为完美匹配,工程中通常要求VSWR<2.0。01反射系数ΓΓ=(ZL−Z₀)/(ZL+Z₀),完全匹配时Γ=0,开路Γ=+1,短路Γ=−1,全反射|Γ|=102驻波比VSWRVSWR=(1+|Γ|)/(1−|Γ|),VSWR=1纯行波,VSWR=∞纯驻波,工程要求VSWR<1.5~2.003反射功率损耗Pr=|Γ|²Pi,VSWR=2.0时|Γ|=1/3,反射功率约11.1%,约1dB功率损耗04阻抗匹配方法串/并联电抗元件、λ/4变换器、渐变线匹配,史密斯圆图是匹配设计的核心工具不同VSWR对应的反射功率占比反射功率随VSWR升高急剧增加,VSWR=5.0时达44.4%MicrowaveEngineering史密斯圆图:微波工程的图形化利器史密斯圆图通过保角映射将复数阻抗变换到反射系数单位圆内,使传输线阻抗计算转化为直观的几何操作。01保角映射z=(1+Γ)/(1−Γ)将复阻抗平面映射到|Γ|≤1的单位圆内,每个点同时表示归一化阻抗和反射系数02等值圆族等电阻圆内切于右端点(r=0最左,r=∞最右),等电抗圆弧通过右端点(x>0上半圆,x<0下半圆)03沿线移动距离l对应旋转角度2βl:向源顺时针,向负载逆时针,半波长旋转一圈回到原点04阻抗匹配串联元件沿等电阻圆移动,并联元件沿等电导圆移动(导纳圆图),目标是将阻抗调至圆心(z=1)矢量网络分析仪上的史密斯圆图显示ImpedanceMatching阻抗匹配方法:λ/4变换器与枝节匹配阻抗匹配通过引入电抗元件或传输线段将负载阻抗变换为传输线特性阻抗。λ/4变换器利用四分之一波长线的阻抗倒置特性,枝节匹配通过并联导纳补偿实现匹配。方法选择取决于带宽需求、频率范围和实现工艺。CoreEquationsλ/4CHARACTERISTICIMPEDANCEZ₁=√(Z₀·ZL)IMPEDANCEINVERSIONZin=Z₁²/ZLMATCHINGCONDITIONYin=Y₀=1/Z₀Transformer在负载前插入长度λ/4、特性阻抗Z₁=√(Z₀·ZL)的传输线段,利用阻抗倒置Zin=Z₁²/ZL实现匹配Bandwidth仅适用于纯电阻负载,带宽有限;多级λ/4级联可扩展带宽,切比雪夫设计可实现等纹波最优宽带匹配StubMatching单枝节匹配:在距负载d处并联短路/开路枝节(长度l),调节d和l使总导纳Y=Y₀,可在史密斯圆图上图解设计LumpedNetworkL型/π型集总网络用LC元件实现匹配,适用于频率较低场景;微带线工艺中枝节匹配更易于实现和调试MULTIPATHPROPAGATION移动通信中的多径传播与衰落特性多径传播使接收信号产生小尺度衰落与时延扩展,均方根时延扩展决定相干带宽——宽带信号超过相干带宽时出现频率选择性衰落,这是4G/5G采用OFDM的根本动因。FADING多条路径分量相位随机叠加,无直射径时幅值服从瑞利分布,有直射径时服从莱斯分布DELAYSPREAD由不同路径传播距离差决定,城市典型值1-5μs,郊区0.1-1μs,室内10-50nsCOHERENTBWBc≈1/(5τrms):信号带宽小于Bc时平坦衰落,大于Bc时频率选择性衰落OFDM将宽带信号分为多个窄带子载波(每个<Bc),使每个子载波经历平坦衰落,简化均衡城市环境中的基站天线与多径传播场景EMPIRICALPROPAGATIONMODELS电波传播经验模型:Okumura-Hata与射线追踪经验传播模型基于海量实测数据拟合,弥补理论模型对复杂环境适应性的不足。5G时代射线追踪模型实现高精度预测,但经验模型仍是网络规划的首选快速工具。01Okumura-Hata模型:L=69.55+26.16lgf−13.82lghb−a(hm)+(44.9−6.55lghb)lgd,适用于150-1500MHz宏蜂窝场景02环境修正因子:大城市a(hm)=(1.1lgf−0.7)hm−(1.56lgf−0.8),中小城市a(hm)=(8.29lghm)²−1.103COST-231扩展:频率扩展至2GHz,增加大城市修正项CM=3dB,是GSM18

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论