期刊文献+
共找到19篇文章
< 1 >
每页显示 20 50 100
微光数字观瞄镜的强杂光抑制技术研究
1
作者 梁国龙 张明超 +3 位作者 黄剑波 丁浩 白晶 张尧禹 《吉林大学学报(信息科学版)》 CAS 2024年第1期81-85,共5页
针对观瞄镜强杂光干扰问题,提出一套完整的强杂光抑制技术方案。首先,对物镜内表面采用绒布消光;然后,在软件图像处理上结合使用相邻图像累加积分、直方图统计和宽动态灰度增强等算法完成强杂光抑制;在天空照度为1×10^(-3)lx以下... 针对观瞄镜强杂光干扰问题,提出一套完整的强杂光抑制技术方案。首先,对物镜内表面采用绒布消光;然后,在软件图像处理上结合使用相邻图像累加积分、直方图统计和宽动态灰度增强等算法完成强杂光抑制;在天空照度为1×10^(-3)lx以下的野外暗夜环境下,加入强杂光干扰后进行实验测试。实验结果表明,该技术方案能有效完成强杂光抑制和图像细节增强,提升图像质量。软件算法在FPGA(Field Programmable Gate Array)平台上运行,2 ms内完成图像处理工作,满足实时性要求。 展开更多
关键词 数字观瞄镜 强杂光 干扰 抑制 消光绒布 直方图 宽动态灰度增强
下载PDF
无源声呐水下多目标融合跟踪方法 被引量:1
2
作者 梁国龙 张博宇 +3 位作者 齐滨 郝宇 杜致尧 李想 《声学学报》 EI CAS CSCD 北大核心 2024年第3期501-512,共12页
针对海洋环境噪声导致弱目标在不同子频带检测结果差异较大,致使以全频带探测结果为输入的跟踪算法出现性能退化的问题,提出一种子带融合跟踪方法。该方法利用改进的高斯混合概率假设密度滤波器对各频率子带输出的方位估计结果进行跟踪... 针对海洋环境噪声导致弱目标在不同子频带检测结果差异较大,致使以全频带探测结果为输入的跟踪算法出现性能退化的问题,提出一种子带融合跟踪方法。该方法利用改进的高斯混合概率假设密度滤波器对各频率子带输出的方位估计结果进行跟踪,并采用广义协方差交集准则对子带跟踪结果进行融合,以获得综合各子带信息的跟踪结果。仿真结果表明,所提方法可以提高弱目标在各子带信噪比不均衡情况下的跟踪能力,且运算时间与对比方法较为接近。海试数据处理结果进一步验证了所提方法的有效性。 展开更多
关键词 无源声呐 广义协方差交集 高斯混合概率假设密度滤波器 子带融合跟踪
原文传递
联合多传感器的水下多目标无源声学定位
3
作者 李想 王燕 +3 位作者 齐滨 郝宇 梁国龙 张涵 《声学学报》 EI CAS CSCD 北大核心 2024年第1期16-27,共12页
在多传感器无源声学定位问题中,不同传感器接收的来自同一目标的信号均对应于此目标位置,依据这一物理基础,提出了一种基于粒子滤波的无源声学定位方法,以有效融合多传感器数据,进而提高定位性能。该方法将粒子滤波中的似然函数定义为... 在多传感器无源声学定位问题中,不同传感器接收的来自同一目标的信号均对应于此目标位置,依据这一物理基础,提出了一种基于粒子滤波的无源声学定位方法,以有效融合多传感器数据,进而提高定位性能。该方法将粒子滤波中的似然函数定义为粒子状态所对应不同传感器信号之间互相关输出的乘积。该似然函数的设计确保所提方法可以充分获取多传感器的处理增益。此外,所提方法摆脱了传统定位范式,因此可以规避传统定位范式必须面对的测量−跟踪关联问题。湖上试验表明,在强多途干扰的条件下,传统定位方法的平均定位误差为7.2 m,而所提方法的平均定位误差为1.2 m,具有更好的性能。 展开更多
关键词 粒子滤波 检测前跟踪 无源声学定位 互相关
原文传递
多站联合观测非合作目标的定位批号关联算法
4
作者 梁国龙 李硕 +2 位作者 邹男 李想 李宇昂 《声学学报》 EI CAS CSCD 北大核心 2024年第5期1005-1018,共14页
在多站联合对非合作多目标进行定位场景中,提出了一种批号关联算法。该方法利用模糊数学理论描述已关联集和待关联量测中方位信息、连续谱特征信息和线谱特征信息之间的模糊隶属度,构建综合模糊关系矩阵进行二维关联,通过设置隶属度门... 在多站联合对非合作多目标进行定位场景中,提出了一种批号关联算法。该方法利用模糊数学理论描述已关联集和待关联量测中方位信息、连续谱特征信息和线谱特征信息之间的模糊隶属度,构建综合模糊关系矩阵进行二维关联,通过设置隶属度门限消除低隶属度关联的影响。仿真结果表明,所提算法克服了传统算法在目标增多关联成功率下降的问题,在一定的漏报虚警时关联结果召回率可达到0.8以上,关联结果精度也保持在较高水平。仿真结果及试验数据均验证了该方法在复杂水下环境中较好的目标关联定位性能。 展开更多
关键词 目标关联 特征信息 多站联合观测 非合作目标 模糊隶属度函数
原文传递
一种适用于半互质阵的高精度波达方向估计方法
5
作者 梁国龙 滕远鑫 +1 位作者 王晋晋 付进 《电子与信息学报》 EI CAS CSCD 北大核心 2024年第8期3228-3237,共10页
在半互质阵列(SCA)下,经典波达方向估计(DoA)算法在面对邻近相干信源时估计性能退化。为了解决该问题,该文提出一种适用于半互质阵列的高精度波达方向估计方法。首先,将半互质阵划分为3个均匀线阵并利用常规波束形成算法对子阵1、子阵2... 在半互质阵列(SCA)下,经典波达方向估计(DoA)算法在面对邻近相干信源时估计性能退化。为了解决该问题,该文提出一种适用于半互质阵列的高精度波达方向估计方法。首先,将半互质阵划分为3个均匀线阵并利用常规波束形成算法对子阵1、子阵2和子阵3的阵列输出信号进行处理;然后,对子阵3的输出信号进行加权后与子阵1,2的输出信号加和构建和波束,利用子阵1与子阵2的输出信号做差构建差波束;最后,将上述和差波束做差得到最终输出信号,计算最终信号的功率得到方位谱。该方法基于半互质阵列特点构建和差波束,充分利用了3个子阵的重叠阵元实现估计精度的提高。通过仿真和湖上实测数据验证表明,该方法可以适用于半互质阵列实现波达方向估计,并且在面对邻近相干信源时其波达方向估计性能优于最小方差无失真估计算法(MVDR)和最小处理算法(MP)。 展开更多
关键词 波达方向估计 半互质阵列 邻近相干信源 波束形成
下载PDF
一种基于组合子阵协方差矩阵的阵列扩展方法
6
作者 梁国龙 罗钧戈 +1 位作者 郝宇 付进 《兵工学报》 EI CAS CSCD 北大核心 2024年第5期1717-1724,共8页
针对阵列在低信噪比条件下信号检测概率低且方位估计性能下降的问题,提出一种基于组合子阵协方差矩阵的线列阵扩展方法。该方法将阵列划分为奇偶子阵,根据子阵的互协方差与自协方差矩阵构建扩展接收数据,对扩展接收数据进行组合作为最... 针对阵列在低信噪比条件下信号检测概率低且方位估计性能下降的问题,提出一种基于组合子阵协方差矩阵的线列阵扩展方法。该方法将阵列划分为奇偶子阵,根据子阵的互协方差与自协方差矩阵构建扩展接收数据,对扩展接收数据进行组合作为最终的阵列接收数据进行处理,由此在突破阵列物理孔径的基础上避免仅使用奇偶阵元互协方差重构接收数据时产生的栅瓣问题。仿真结果表明:新方法可以降低波束输出的旁瓣,具有较好的弱目标检测能力和方位分辨力;相比于基于奇偶阵元互协方差的阵列扩展方法,新方法可有效提升方位估计精度,且不存在栅瓣问题。 展开更多
关键词 阵列扩展 弱目标检测 波达方向估计 栅瓣 协方差矩阵
下载PDF
针对周期偏移声信标的二阶高精度水声定位模型 被引量:2
7
作者 孙思博 明瑞和 +2 位作者 杨卓 梁国龙 付进 《声学学报》 EI CAS CSCD 北大核心 2023年第1期41-49,共9页
针对传统长基线水声定位模型在以飞机黑匣子等周期偏移声信标为目标时定位精度显著下降的问题,提出了一种不依赖于信号周期信息的新型长基线定位模型。该模型将传统时延差定位模型拓展至二阶,通过求取到达时间的二阶导数,消去未知声信... 针对传统长基线水声定位模型在以飞机黑匣子等周期偏移声信标为目标时定位精度显著下降的问题,提出了一种不依赖于信号周期信息的新型长基线定位模型。该模型将传统时延差定位模型拓展至二阶,通过求取到达时间的二阶导数,消去未知声信标周期。为避免由阶数提高引起的解算非收敛问题,引入了改进牛顿迭代算法,有效提高了新模型的收敛概率。所提出新模型的仿真试验及湖上试验的定位误差分别为3.14 m和1.19 m,试验结果表明,该模型实现了周期偏移声信标的高精度定位解算。 展开更多
关键词 水声定位 牛顿迭代算法 信号周期 收敛问题 二阶导数 黑匣子 时延差 定位误差
原文传递
双基地有源探测柔性水平阵阵形实时校准方法 被引量:1
8
作者 梁国龙 周奕辰 邱龙皓 《声学学报》 EI CAS CSCD 北大核心 2023年第4期831-842,共12页
双基地有源模式下接收基地搭载的柔性水平阵在运动过程中,阵形会发生实时的畸变。针对阵形误差的存在会严重影响水听器阵列探测性能的问题,本文提出一种基于发射信号为线性调频(LFM)信号,利用直达波实现阵形参数估计的方法。该方法首先... 双基地有源模式下接收基地搭载的柔性水平阵在运动过程中,阵形会发生实时的畸变。针对阵形误差的存在会严重影响水听器阵列探测性能的问题,本文提出一种基于发射信号为线性调频(LFM)信号,利用直达波实现阵形参数估计的方法。该方法首先利用分数阶傅里叶变换(FRFT)对LFM信号在特定阶数的分数阶傅里叶变换域上的能量聚集性,实现直达波与目标回波的分离,避免了目标回波对于阵形参数估计的影响。而后利用各阵元直达波FRFT后对应峰值点之间的相位差实现了阵形参数估计,并提出了一种基于拖曳船辐射噪声的阵形解模糊方法。计算机仿真表明该方法对于噪声及干扰影响、入射角度误差以及多普勒频移均有一定的宽容性,湖上实验验证了阵形估计方法的有效性,是一种可实时校准、高精度、稳健的柔性水平阵阵形估计算法。 展开更多
关键词 双基地 有源探测 阵形校准 分数阶傅里叶变换
原文传递
水下单信标导航航路规划和最优航路接近方法
9
作者 张新宇 王燕 +2 位作者 邹男 梁国龙 付进 《导航定位与授时》 CSCD 2023年第5期18-30,共13页
相比于传统的长基线和超短基线等导航方式,水下单信标导航具有布放简单的优点,但其导航精度有待进一步提高。为此,提出了单信标导航的航路规划方案,通过泰勒级数展开推导了水平位置精度因子的表达式,并分析了导航点和声信标的相对几何... 相比于传统的长基线和超短基线等导航方式,水下单信标导航具有布放简单的优点,但其导航精度有待进一步提高。为此,提出了单信标导航的航路规划方案,通过泰勒级数展开推导了水平位置精度因子的表达式,并分析了导航点和声信标的相对几何位置关系对导航精度的影响,最终提出了航路规划方案。在此基础上,还提出了自主水下航行器(autonomous underwater vehicle,AUV)从其他位置接近最优航路的方法,包括两部分:一是建立以时延为观测量的滤波模型,利用滤波算法实时获取AUV的位置估计值;二是基于可观测度分析结果对最优航路接近过程的轨迹进行了设计。二者的结合使得AUV高效率且高精度地逼近最优航路。仿真证明了采用所提出的航路规划方案和最优航路接近方法可以提高导航精度,航路规划方案使得AUV位置估计的均方根误差近似为2 m。 展开更多
关键词 航路规划 最优航路接近 精度分析 可观测度 单信标导航
下载PDF
信标信号的随机共振稳态跃变检测方法
10
作者 万光明 梁国龙 +1 位作者 王晋晋 董文峰 《哈尔滨工程大学学报》 EI CAS CSCD 北大核心 2023年第9期1605-1612,共8页
针对飞行数据记录器水下搜探中的微弱信标信号检测问题,本文提出了一种基于稳态跃变分布的信标信号随机共振检测方法。通过选取合理的混频器载波频率和非线性系统参数,实现了水下声信标信号这类高频窄脉宽信号的随机共振。引入势垒高度... 针对飞行数据记录器水下搜探中的微弱信标信号检测问题,本文提出了一种基于稳态跃变分布的信标信号随机共振检测方法。通过选取合理的混频器载波频率和非线性系统参数,实现了水下声信标信号这类高频窄脉宽信号的随机共振。引入势垒高度变换系数,通过研究非线性系统输出波形的稳态跃变分布特性,构建了检测统计量即稳态跃变适配次数。该方法可实现微弱信标信号的随机共振检测,调整势垒高度变换系数可使输出时域波形的稳态跃变分布产生变化,减少虚警。通过仿真及湖上试验表明:该方法检测性能优于频域能量检测方法。 展开更多
关键词 飞行数据记录器 水下定位信标 黑匣子 随机共振 微弱信号检测 势垒高度 非线性系统 检测统计量
下载PDF
高速垂向运动声源声场干涉特性及其可观测性
11
作者 赵春鹏 梁国龙 +2 位作者 邱龙皓 付进 柳国龙 《声学学报》 EI CAS CSCD 北大核心 2023年第6期1081-1097,共17页
以高速垂向运动声源宽带连续谱噪声为研究对象,基于射线理论分析了声场干涉现象的产生机理,推导了干涉条纹的轨迹方程,讨论了影响干涉结构特性的因素,揭示了声场干涉结构的特性和规律。理论和仿真结果表明:干涉周期与多途信号时延差成... 以高速垂向运动声源宽带连续谱噪声为研究对象,基于射线理论分析了声场干涉现象的产生机理,推导了干涉条纹的轨迹方程,讨论了影响干涉结构特性的因素,揭示了声场干涉结构的特性和规律。理论和仿真结果表明:干涉周期与多途信号时延差成反比关系;减小接收深度或增大接收距离可以提高干涉条纹的可观测性;较大的海面起伏易使干涉条纹形成断点,可观测性降低。湖试数据分析结果表明,干涉条纹主要由直达声-海面反射声干涉条纹构成,干涉条纹宽度和干涉周期随水平距离的增大而增大。该结果可为高速垂向运动声源运动状态分析、环境参数反演等提供参考。 展开更多
关键词 声场干涉特性 干涉条纹 水下动目标观测 浅海波导 随机海面起伏
原文传递
光电瞄准镜电池监测系统研究
12
作者 梁国龙 白晶 +3 位作者 黄剑波 张明超 丁浩 张尧禹 《电测与仪表》 北大核心 2021年第4期194-200,共7页
为了实现便携式光电瞄准镜电池信息监测,考虑到光电瞄准镜低功耗、小尺寸的要求,提出一种仅使用单片FPGA作为处理器的电池监测系统结构。系统对2个串联的18650电池进行实时状态监测,FPGA对电池信息和图像信息进行同步融合叠加,在OLED屏... 为了实现便携式光电瞄准镜电池信息监测,考虑到光电瞄准镜低功耗、小尺寸的要求,提出一种仅使用单片FPGA作为处理器的电池监测系统结构。系统对2个串联的18650电池进行实时状态监测,FPGA对电池信息和图像信息进行同步融合叠加,在OLED屏上显示。显示屏上可以直观看到剩余电量、电池温度、电压、电流信息,当电池电量不足或者电流、电压异常时,通过字体颜色变红给出相应警示。经过多次反复系统实验,光电瞄准镜的电池监测系统具有精度高、稳定性高等特点,具有工程实用性。 展开更多
关键词 光电瞄准镜 FPGA DS2438Z OLED
下载PDF
利用夹角几何关系的超短基线定位方法 被引量:12
13
作者 梁国龙 张毅锋 付进 《哈尔滨工程大学学报》 EI CAS CSCD 北大核心 2019年第8期1474-1479,共6页
针对传统超短基线的定位方法定位精度有限的问题,本文提出了一种利用夹角几何关系的超短基线高精度定位方法。该方法构建超短基线定位误差分析模型,依据超短基线中存在的夹角几何关系,研究了深度信息对超短基线定位精度的影响。通过有... 针对传统超短基线的定位方法定位精度有限的问题,本文提出了一种利用夹角几何关系的超短基线高精度定位方法。该方法构建超短基线定位误差分析模型,依据超短基线中存在的夹角几何关系,研究了深度信息对超短基线定位精度的影响。通过有效利用深度信息,减小主要误差源对定位精度的影响,从而达到提高定位精度的目的。研究表明:该方法能够显著提高超短基线定位精度,尤其是对以相位差估计误差为主的情况效果更为明显,且在绝大多数区域均能提高定位精度。 展开更多
关键词 超短基线 定位方法 误差分析 深度信息 夹角几何关系 高精度 相位差估计误差 误差源
下载PDF
基于均匀圆阵中心对称性的相干源方位估计 被引量:4
14
作者 梁国龙 田蕴琦 +1 位作者 付进 邹男 《哈尔滨工程大学学报》 EI CAS CSCD 北大核心 2019年第12期1958-1964,共7页
均匀圆阵的相干源方位估计一直是阵列信号处理中备受关注的问题,而现行算法的分辨率仍有进一步提高的空间。针对这一问题,本文提出了一种基于均匀圆阵中心对称性的解相干算法。首先利用模式空间变换将均匀圆阵等效为虚拟均匀线阵,然后... 均匀圆阵的相干源方位估计一直是阵列信号处理中备受关注的问题,而现行算法的分辨率仍有进一步提高的空间。针对这一问题,本文提出了一种基于均匀圆阵中心对称性的解相干算法。首先利用模式空间变换将均匀圆阵等效为虚拟均匀线阵,然后利用均匀圆阵的中心对称性对模式空间数据进行共轭平均,最后计算数据的相关函数并据此构造Hermitian Toeplitz矩阵,实现相干源的方位估计。本算法实现简单,与传统算法相比具有更高的分辨率和方位估计精度。仿真实验证实了所提算法的有效性和性能优势。 展开更多
关键词 均匀圆阵 相干源 中心对称性 模式空间 HERMITIAN TOEPLITZ矩阵 方位估计
下载PDF
基于Chan算法的水下测控设备组网数据融合技术研究 被引量:6
15
作者 梁国龙 赵天白 +1 位作者 邹男 张博轩 《电子与信息学报》 EI CSCD 北大核心 2018年第5期1181-1186,共6页
依据水下测控设备组网系统的工作特点,为了有效挖掘现有测控设备的使用效能,构建高精度的综合性水下测控网络,该文提出基于Chan算法的水下测控设备组网集中式数据融合定位算法。该算法首先利用基于波达时间的加权最小二乘算法粗测目标位... 依据水下测控设备组网系统的工作特点,为了有效挖掘现有测控设备的使用效能,构建高精度的综合性水下测控网络,该文提出基于Chan算法的水下测控设备组网集中式数据融合定位算法。该算法首先利用基于波达时间的加权最小二乘算法粗测目标位置,然后依据该目标位置和测量时延信息的关系,构造新的误差矢量,利用该误差矢量再次加权最小二乘估计解算目标位置。研究结果表明,该算法可实现多套水下测控设备的数据融合,可有效提高全域范围的定位精度,且精度高于纯基于波达时间的融合定位算法。 展开更多
关键词 水下定位 CHAN算法 数据融合 定位精度
下载PDF
Frequency-invariant robust adaptive beamforming based on interference covariance matrix reconstruction 被引量:6
16
作者 FAN Zhan liang guolong 《Chinese Journal of Acoustics》 2014年第4期359-368,共10页
Consider the problems of frequency-invariant beampattern optimization and robustness in broadband beamforming.Firstly,a global optimization algorithm,which is based on phase compensation of the array manifolds,is used... Consider the problems of frequency-invariant beampattern optimization and robustness in broadband beamforming.Firstly,a global optimization algorithm,which is based on phase compensation of the array manifolds,is used to construct the frequency-invariant beampattern.Compared with some methods presented recently,the proposed algorithm is not only available to get the global optimal solution,but also simple for physical realization.Meanwhile,a robust adaptive broadband beamforming algorithm is also derived by reconstructing the covariance matrix.The essence of the proposed algorithm is to estimate the space-frequency spectrum using Capon estimator firstly,then integrate over a region separated from the desired signal direction to reconstruct the interference-plus-noise covariance matrix,and finally caleulate the adaptive beamformer weights with the reconstructed matrix.The design of beamformer is formulated as a convex optimization problem to be solved.Simulation results show that the performance of the proposed algorithm is almost always close to the optimal value across a wide range of signal to noise ratios. 展开更多
关键词 Frequency-invariant robust adaptive beamforming based on interference covariance matrix reconstruction
原文传递
Variance-of-instantaneous-frequency detector and its performance evaluation 被引量:2
17
作者 liang guolong HUI Junying(Harbin Engineering University Harbin 150001) 《Chinese Journal of Acoustics》 1999年第1期19-28,共10页
The optimum detector for a determined signal in the stationary white noise background In ideal channels is a matched filter.But its performance degenerates seriously in a practical underwater acoustic system for the ... The optimum detector for a determined signal in the stationary white noise background In ideal channels is a matched filter.But its performance degenerates seriously in a practical underwater acoustic system for the complexityof underwater acoustic channels.Based on the estimation of instantaneous frequency sequence and its statistical characteristics,a new type ofrobust detector,Varlance-of-Instantaneous-Frequency Detector(VIFD for short),is presented for the detection of CW pulse according to the features of underwater acousticchannels and operating environment of active sonars. It is proved through theoretical analysis,simulation and experiments that VIFD is promising for its robustness in nonstatlonary interference background and channels and that it is alsorobust for the fluctuations and decny of signal or mlsmatch caused by Doppler-frequency shift. 展开更多
关键词 Variance-of-instantaneous-frequency detector and its performance evaluation KHZ
原文传递
A subspace fitting algorithm of acoustic vector sensor array and corresponding matrix pre-filter design 被引量:1
18
作者 WANG Yan WU Wenfeng +1 位作者 FAN Zhan liang guolong 《Chinese Journal of Acoustics》 2014年第3期267-278,共12页
In order to ease the pass-band response distortion of the matrix pre-filter,a simple approach for designing matrix spatial filter is proposed,which minimizes the sum of the k maximal distortion norm(k is the number o... In order to ease the pass-band response distortion of the matrix pre-filter,a simple approach for designing matrix spatial filter is proposed,which minimizes the sum of the k maximal distortion norm(k is the number of the constraint points)within the pass-band,while constraining the filter response within the stop-band.Considering the costly amount of calculation of the high-resolution methods,an algorithm with small amount of calculation based on matrix pre-filtering and subspace fitting using acoustic vector array(MF-VSSF)is proposed.Through joint processing of signal subspace of both pressure and particle velocity,the pre-filtering matrix and the signal subspace is decreased to M-dimensional(M is the number of array-element),hence reduces the time-consumption of the matrix pre-filter design and DOA searching.Simulation results show that,the method offers the same performance as MUSIC with pre-filtering,but has much lesser amount of calculation.Moreover,the designed prefilter can efficiently suppress the interference in the stop-band and improve the estimation and resolution performance of successive DOA estimators. 展开更多
关键词 MUSIC A subspace fitting algorithm of acoustic vector sensor array and corresponding matrix pre-filter design
原文传递
Pseudo bearing sequence analysis 被引量:1
19
作者 liang guolong HUI Junying CAI Ping (Harbin Engineering University Harbin 150001) 《Chinese Journal of Acoustics》 2000年第1期22-34,共13页
A signal analyzing method, named pseudo bearing sequence analysis (PBSA), is proposed. The theoretical foundation and basic ideas of the method, as well as its practical applications in the area of underwater acoustic... A signal analyzing method, named pseudo bearing sequence analysis (PBSA), is proposed. The theoretical foundation and basic ideas of the method, as well as its practical applications in the area of underwater acoustic signal processing, are introduced. A signal detector named variance-of-pseudo-bearing detector (VPBD) is presented as well, which is superior to variance-of-instantaneous frequency detector (VIFD ) in t he reverberation background. The validity of the proposed method bas been confirmed by means of both simulations in computer and the practical processing of trial data. 展开更多
关键词 PBSA Pseudo bearing sequence analysis
原文传递
上一页 1 下一页 到第
使用帮助 返回顶部