期刊文献+
共找到34篇文章
< 1 2 >
每页显示 20 50 100
基于AR-Hankel矩阵的风力发电机早期故障诊断方法研究 被引量:6
1
作者 尹勋 张新燕 +2 位作者 常喜强 赵理威 王衡 《可再生能源》 CAS 北大核心 2016年第1期80-85,共6页
针对风力发电机故障早期阶段,故障特征不明显,继电保护装置检测不到相关电气量异常的问题,文章提出了基于AR模型、Hankel矩阵和奇异值分解的风力发电机早期故障诊断方法。该方法先对主轴径向振动信号进行总体平均经验模式分解,再按照互... 针对风力发电机故障早期阶段,故障特征不明显,继电保护装置检测不到相关电气量异常的问题,文章提出了基于AR模型、Hankel矩阵和奇异值分解的风力发电机早期故障诊断方法。该方法先对主轴径向振动信号进行总体平均经验模式分解,再按照互相关准则选择若干个固有模态函数建立AR模型,然后对自回归系数构建Hankel矩阵并作奇异值分解,将奇异值作为故障特征输入支持向量机判断发电机的运行状态。试验结果表明,该方法能对直驱风力发电机正常运行、定子线圈匝间短路、发电机主轴偏心、发电机轴承磨损4种状态进行准确诊断。 展开更多
关键词 总体平均经验模式分解 自回归模型 hankel矩阵 奇异值分解 早期故障诊断
在线阅读 下载PDF
基于最优Hankel范数近似的核电站电功率与冷却剂平均温度解耦控制方法 被引量:2
2
作者 钱虹 房振鲁 +1 位作者 杨祖魁 周蕾 《动力工程学报》 CAS CSCD 北大核心 2016年第5期383-388,共6页
采用最小二乘法获得了被控对象模型,针对核电站电功率和冷却剂平均温度控制耦合的特点,设计了对角阵完全解耦器,基于最优Hankel范数降阶模型近似的理论,对对角阵完全解耦器进行模型降阶,并在Matlab/Simulink中对所得到的模型进行仿真验... 采用最小二乘法获得了被控对象模型,针对核电站电功率和冷却剂平均温度控制耦合的特点,设计了对角阵完全解耦器,基于最优Hankel范数降阶模型近似的理论,对对角阵完全解耦器进行模型降阶,并在Matlab/Simulink中对所得到的模型进行仿真验证.结果表明:对对角阵完全解耦器降阶后的PID控制系统比单纯的PID控制系统的响应速度快,且进行功率控制时的超调量较小. 展开更多
关键词 核电站 最优hankel范数近似 解耦控制 模型降阶
在线阅读 下载PDF
高速列车万向轴动不平衡检测的EEMD-Hankel-SVD方法 被引量:9
3
作者 丁建明 林建辉 赵洁 《机械工程学报》 EI CAS CSCD 北大核心 2015年第10期143-151,159,共10页
针对聚合经验模式分解(Ensemble empirical model decomposition,EEMD)的等效滤波特性依然存在模式分量间频带重叠较大的根本缺陷,提出一种高速列车万向轴动不平衡动态检测的新方法。该方法的核心是对万向节安装机座的振动信号进行EEMD... 针对聚合经验模式分解(Ensemble empirical model decomposition,EEMD)的等效滤波特性依然存在模式分量间频带重叠较大的根本缺陷,提出一种高速列车万向轴动不平衡动态检测的新方法。该方法的核心是对万向节安装机座的振动信号进行EEMD分解得到基本模式分量,应用基本模式分量信号来构造Hankel矩阵,对该矩阵进行正交化奇异值(Singular value decomposition,SVD)分解,以奇异值关键叠层作为奇异值的选择准则对信号进行重构,应用重构信号的傅里叶谱来检测高速列车万向轴的动不平衡,消除EEMD分解模式频带重叠对故障特征的淹没和混淆效应,提高了谱的清晰度,凸显了故障特征。应用万向轴动不平衡试验数据对该方法进行试验验证,结果表明,该方法能够有效检测万向轴动不平衡引起的故障特征和万向轴的固有振动特征,与纯EEMD方法相比,该方法在谱的清晰度和故障表征力上得到了显著提高。 展开更多
关键词 高速列车 万向轴动不平衡 聚合经验模式分解(Ensemble empirical model decomposition EEMD) hankel矩阵 正交化奇异值(Singular value decomposition SVD) 动态检测
在线阅读 下载PDF
Hankel范数准则的应用:多变量模型匹配设计方法
4
作者 姚一新 王广雄 《宇航学报》 EI CAS CSCD 北大核心 1990年第4期15-22,共8页
本文在Hankel范数空间研究多变量控制系统模型匹配问题,其中包括完全和不完全模型匹配。应用本文提出的方法,可以获得一个Hankel范数意义下的最优模型匹配控制器。文中提出了一种多变量控制系统完全模型匹配的充要条件。如果这一条件满... 本文在Hankel范数空间研究多变量控制系统模型匹配问题,其中包括完全和不完全模型匹配。应用本文提出的方法,可以获得一个Hankel范数意义下的最优模型匹配控制器。文中提出了一种多变量控制系统完全模型匹配的充要条件。如果这一条件满足,可实现完全模型匹配。如果条件不满足,可以求得一个最优的近似解。文中还给出了设计例子。 展开更多
关键词 多变量控制 模型匹配 hankel范数
在线阅读 下载PDF
基于最优Hankel范数近似的线性相位ⅡR滤波器设计 被引量:4
5
作者 邓娜 邵世煌 顾国祥 《控制理论与应用》 EI CAS CSCD 北大核心 2007年第2期200-204,共5页
提出了一类基于最优Hankel范数近似的线性相位无限脉冲响应(IIR)滤波器设计方法.首先给出了Hankel范数的相关预备知识,然后给出了离散时间单输入单数出系统Hankel范数近似的定理及证明,最后给出了线性相位IIR滤波器的设计步骤.该方法... 提出了一类基于最优Hankel范数近似的线性相位无限脉冲响应(IIR)滤波器设计方法.首先给出了Hankel范数的相关预备知识,然后给出了离散时间单输入单数出系统Hankel范数近似的定理及证明,最后给出了线性相位IIR滤波器的设计步骤.该方法不但减小了逆矩阵求解过程中的计算量,同时给出了(?)2范数近似的误差边界,仿真结果验证了该方法的有效性. 展开更多
关键词 FIR滤波器 IIR滤波器 最优hankel范数近 模型降阶 线性相位 群时延
在线阅读 下载PDF
基于Hankel范数模型降价的控制对象的名义模型选择 被引量:2
6
作者 周彤 《自动化学报》 EI CSCD 北大核心 1997年第2期247-252,共6页
研究了控制对象具有多个模型时,求取其适合鲁棒控制器设计的名义模型的问题.提出了一种基于Hankel范数模型降阶的名义模型选择算法.
关键词 鲁棒控制 建模 模型误差 hankel范数
在线阅读 下载PDF
离散时间系统Hankel最优模型降阶及其应用
7
作者 邓娜 苏庆新 《控制工程》 CSCD 2006年第6期557-559,600,共4页
针对离散时间系统Hankel最优模型降阶问题,对标量情况由最优Hankel范数近似给出了状态空间解。将这一理论应用于无限脉冲响应(IIR)滤波器的设计,进一步提出了一种基于最优Hankel范数近似的线性相位IIR滤波器设计方法。通过采用线性相位... 针对离散时间系统Hankel最优模型降阶问题,对标量情况由最优Hankel范数近似给出了状态空间解。将这一理论应用于无限脉冲响应(IIR)滤波器的设计,进一步提出了一种基于最优Hankel范数近似的线性相位IIR滤波器设计方法。通过采用线性相位FIR滤波器和最优Hankel范数近似方法设计IIR滤波器,使其具有线性相位,减小了逆矩阵求解中的计算量。低通滤波器的仿真实例验证了该方法的有效性。 展开更多
关键词 离散时间系统 hankel最优模型降阶 线性相位 IIR滤波器 FIR滤波器
在线阅读 下载PDF
基于Hankel阵的荧光油膜灰度与厚度预测模型改进 被引量:3
8
作者 钱泓江 董秀成 张征宇 《航空动力学报》 EI CAS CSCD 北大核心 2021年第10期2061-2071,共11页
针对荧光油膜灰度与厚度关系数据采集方法较为繁琐、耗时耗力这一问题,研究了基于Hankel阵的系统辨识算法,并在此基础提出了Hankel阵误差修正模型和Hankel阵高阶迭代误差修正模型等两种改进方法,利用了极少数据量建立模型,实现了对其余... 针对荧光油膜灰度与厚度关系数据采集方法较为繁琐、耗时耗力这一问题,研究了基于Hankel阵的系统辨识算法,并在此基础提出了Hankel阵误差修正模型和Hankel阵高阶迭代误差修正模型等两种改进方法,利用了极少数据量建立模型,实现了对其余未知荧光油膜厚度值的预测,且保持了较高的精度。试验结果表明:基于极少数据量建模以预测该数据量外的数据点这一特殊背景下,插值法的外插能力显得并不适用。而传统Hankel阵预测模型的预测精度为76.69%,Hankel阵误差修正模型和Hankel阵高阶迭代误差修正模型的预测精度分别为85.69%和89.25%,较之传统方法预测精度分别提高了9%和12.56%,为荧光油流技术领域针对荧光油膜灰度与厚度建模问题提供了一种可行技术路线,具有一定的实际工程应用意义。 展开更多
关键词 荧光油膜 灰度 油膜厚度 hankel 预测模型
原文传递
一类截断Hankel算子的复对称性
9
作者 赖丽玲 梁金金 陈泳 《数学物理学报(A辑)》 CSCD 北大核心 2022年第4期961-968,共8页
截断Hankel算子是Hardy空间上的Hankel算子在模空间上的压缩.该文研究模空间上一类截断Hankel算子的复对称性,给出了完全刻画.所得结果表明,截断Hankel算子的复对称可能仅与模空间有关,也可能与模空间和算子的符号函数均相关.
关键词 hankel算子 模空间 复对称性 有限Blaschke乘积
在线阅读 下载PDF
基于Hankel范数与平衡截断法的联合降阶模型
10
作者 李久芹 陈佃浩 《科技视界》 2016年第3期187-188,共2页
为简化复杂系统控制器或观测器的设计,本文提出一种Hankel范数降阶与平衡截断降阶相结合的联合降阶法,得到了复杂系统的降阶模型;通过仿真分析,验证了该方法的有效性。
关键词 hankel范数降阶 平衡截断 电力系统
在线阅读 下载PDF
基于Hankel矩阵改进的TLS-ESPRIT多频带融合处理 被引量:3
11
作者 张硕 杨君 +2 位作者 葛鹏程 冯婷婷 杜钰 《电光与控制》 CSCD 北大核心 2022年第9期90-95,共6页
针对传统矩阵束算法在低信噪比条件下,估计全极点模型参数时存在计算数值不稳定、极点估计精度不高的缺点,提出先对回波信号引入Hankel矩阵以去信号相干性,然后利用改进的TLS-ESPRIT算法估计多频带信号极点,再利用基于Hankel矩阵改进的T... 针对传统矩阵束算法在低信噪比条件下,估计全极点模型参数时存在计算数值不稳定、极点估计精度不高的缺点,提出先对回波信号引入Hankel矩阵以去信号相干性,然后利用改进的TLS-ESPRIT算法估计多频带信号极点,再利用基于Hankel矩阵改进的TLS-ESPRIT替代传统的矩阵束求全频带全极点模型参数,以提高融合信号的精度。最终的仿真及实测数据结果显示,所提算法提高了多频带信号在低信噪比下的融合性能。 展开更多
关键词 多频带信号 全极点模型 改进的TLS-ESPRIT 非相干量估计 hankel矩阵
在线阅读 下载PDF
基于Koopman算子与卡尔曼滤波的同步发电机动态状态估计 被引量:6
12
作者 焦鹏悦 杨德友 蔡国伟 《电力系统保护与控制》 EI CSCD 北大核心 2024年第9期27-35,共9页
动态状态估计是监测同步发电机动态行为的重要手段,准确的动态状态估计结果对于指导电力系统安全运行与高效控制具有重要意义。从数据驱动的角度出发,提出了基于Koopman算子与卡尔曼滤波的同步发电机动态状态估计方法。该方法首先利用... 动态状态估计是监测同步发电机动态行为的重要手段,准确的动态状态估计结果对于指导电力系统安全运行与高效控制具有重要意义。从数据驱动的角度出发,提出了基于Koopman算子与卡尔曼滤波的同步发电机动态状态估计方法。该方法首先利用汉克尔动态模态分解算法从发电机动态响应数据中提取Koopman算子,进而以提取的Koopman算子为基础构建同步发电机状态空间模型,并利用卡尔曼滤波对同步发电机状态变量进行动态估计。该方法无须事先构建发电机模型及参数,实现了完全数据驱动的动态状态估计。仿真实验结果表明,在发电机模型及参数失配的情况下该方法估计精度明显高于传统以模型为基础的估计结果,具有较好的自适应性和鲁棒性。 展开更多
关键词 动态状态估计 模型 数据驱动 Koopman算子 卡尔曼滤波 汉克尔动态模态分解
在线阅读 下载PDF
循环荷载作用下分数阶饱和黏土地基的动力固结分析
13
作者 王立安 余云燕 +1 位作者 任新 陈辉 《岩土力学》 EI CAS CSCD 北大核心 2024年第8期2279-2289,共11页
在Biot固结理论框架下,利用分数阶Kelvin模型描述土骨架流变效应,考虑孔压消散和土骨架流变的耦合作用。通过建立周期性循环荷载的时空域解析函数,在柱坐标系下构建半空间饱和黏土地基的三维轴对称动力固结控制方程。采用Hankel-Laplac... 在Biot固结理论框架下,利用分数阶Kelvin模型描述土骨架流变效应,考虑孔压消散和土骨架流变的耦合作用。通过建立周期性循环荷载的时空域解析函数,在柱坐标系下构建半空间饱和黏土地基的三维轴对称动力固结控制方程。采用Hankel-Laplace联合变换及张量运算推导控制方程在变换域的解析解,然后利用数值反演得出物理场的时空域解。通过算例分析,对循环荷载作用下饱和黏土地基的动力固结特性开展研究。结果表明,饱和黏土在主固结阶段的沉降速率较慢,而次固结阶段的沉降速率加快;循环荷载作用下,土骨架流变性越强,土体的累积沉降发展越快,土体位移波动幅值随黏性阶数的增大而减小,而且黏性阶数越大位移滞后性越明显;土骨架流变性导致孔压响应滞后于有效应力,孔压与有效应力的螺旋曲线发生水平移动;在循环荷载的卸载阶段,正应力随外荷载的减小而减小,但切应力随之增大,土体因发生剪胀现象而产生负孔压。 展开更多
关键词 分数阶Kelvin模型 饱和黏土地基 动力固结 循环荷载 hankel-Laplace变换
原文传递
基于GA优化的LQR控制对并联机器人的控制研究 被引量:3
14
作者 王林军 史宝周 +2 位作者 张东 丁仕豪 刘建明 《振动与冲击》 EI CSCD 北大核心 2020年第20期82-90,96,共10页
并联机构由于非线性强、控制过程复杂,导致其控制困难,精度难以满足工作要求。建立了三自由度刚柔耦合并联机器人的空间结构简化模型和多体系统树模型,基于多体系统传递矩阵法(MSTMM)对其动力学和状态空间表示进行了分析,并对其状态空... 并联机构由于非线性强、控制过程复杂,导致其控制困难,精度难以满足工作要求。建立了三自由度刚柔耦合并联机器人的空间结构简化模型和多体系统树模型,基于多体系统传递矩阵法(MSTMM)对其动力学和状态空间表示进行了分析,并对其状态空间进行了Hankel模型约简;以此为理论基础,设计了LQR控制和基于遗传算法(GA)优化的LQR控制,分别对机器人工作情况进行了仿真对比分析。结果表明:基于遗传算法优化的LQR控制下机器人只需用时1 s便进入稳态。在稳态情况下,优化后与优化前相比,机器人动平台x、y和z方向最大位移分别下降了24.34%、13.51%和1.03%;各支腿上最大加速度分别下降了26.92%、33.96%和35.71%;各支腿上最大控制力分别下降了8.96%、7.37%和9.01%。由此可见在该优化控制方法下机器人响应快、精度高、稳定性好,充分证明了该方法的合理性和优越性,为进一步对并联机器人动力学性能和控制方法的研究提供了一定的理论依据。 展开更多
关键词 遗传算法(GA) 并联机器人 传递矩阵法 动力学 状态空间表示 hankel模型 LQR控制
在线阅读 下载PDF
考虑基岩作用的圆形粘弹性地基振动分析 被引量:2
15
作者 郑灶锋 蔡袁强 +1 位作者 徐长节 朱渠艳 《江南大学学报(自然科学版)》 CAS 2004年第5期512-515,543,共5页
针对Merchant粘弹性本构模型,从轴对称饱和地基动力方程出发,运用Hankel积分变换,研究了粘弹性地基的竖向振动问题,获得了位移、应力、孔隙水压力及渗流量的积分形式解.并由此结合数值算例,讨论了粘弹性地基受动力荷载作用的振动特性.... 针对Merchant粘弹性本构模型,从轴对称饱和地基动力方程出发,运用Hankel积分变换,研究了粘弹性地基的竖向振动问题,获得了位移、应力、孔隙水压力及渗流量的积分形式解.并由此结合数值算例,讨论了粘弹性地基受动力荷载作用的振动特性.结果表明:地基竖向振幅沿半径呈波动衰减,粘弹性土体竖向振幅随半径变化要比弹性土体的小. 展开更多
关键词 Merchant模型 hankel变换 饱和地基 振动分析
在线阅读 下载PDF
CSAMT一维正演可视化设计与激发极化效应 被引量:4
16
作者 钟华 唐新功 《物探与化探》 CAS 北大核心 2022年第2期459-466,共8页
针对可控源音频大地电磁法(CSAMT)开发了一套CSAMT一维数据正演可视化软件。其界面友好,操作简单;使用MATLAB进行核心算法的编写,正演计算过程引入线性滤波系数对汉克尔积分进行求解,并对不同线性滤波系数方案的精度进行了比较;基于Jav... 针对可控源音频大地电磁法(CSAMT)开发了一套CSAMT一维数据正演可视化软件。其界面友好,操作简单;使用MATLAB进行核心算法的编写,正演计算过程引入线性滤波系数对汉克尔积分进行求解,并对不同线性滤波系数方案的精度进行了比较;基于Java语言进行操作界面设计,提供多参数设置界面,自动化程度高;通过软件自带的绘图功能,将数据处理结果直观地提供给用户。在计算过程中引入Cole-Cole模型参数,以复电阻率代替不考虑地电体极化效应的直流电阻率,对CSAMT场源一维层状激化介质模型进行了正演模拟,讨论了极化层分别位于不同深度时一维层状介质的极化效应。研究发现极化层埋藏越浅,极化效应对于CSAMT的影响频带越宽。这对于认识含极化效应的CSAMT方法的电磁响应特征,提高数据处理便利度与工作效率方面都具有重要的作用。 展开更多
关键词 一维正演 汉克尔变换 可视化 激发极化 Java
在线阅读 下载PDF
CSAMT三维交错采样有限差分数值模拟 被引量:12
17
作者 邓居智 谭捍东 +2 位作者 陈辉 佟拓 林昌洪 《地球物理学进展》 CSCD 北大核心 2011年第6期2026-2032,共7页
系统分析大地电磁三维交错采样有限差分算法的基础上,根据可控源音频大地电磁场特征,采用将总场分解为一次场和二次场计算,一次场利用快速汉克尔变换,二次场利用数值模拟的思路.从CSAMT满足的麦克斯韦方程组积分形式出发,利用交错采样... 系统分析大地电磁三维交错采样有限差分算法的基础上,根据可控源音频大地电磁场特征,采用将总场分解为一次场和二次场计算,一次场利用快速汉克尔变换,二次场利用数值模拟的思路.从CSAMT满足的麦克斯韦方程组积分形式出发,利用交错采样有限差分算法推导了电场和磁场的离散关系式,提出了简洁的边界条件和合理的剖分方案,所实现的CSAMT三维正演算法具有迭代速度快、收敛稳定、计算精度高等特点.通过均匀半空间和三维低阻体两个理论模型的计算结果检验了算法的正确性和计算精度,为CSAMT三维反演奠定了基础. 展开更多
关键词 CSAMT 交错采样有限差分 三维正演 边界条件 汉克尔变换
在线阅读 下载PDF
瞬变电磁快速三维正演 被引量:5
18
作者 李瑞雪 王鹤 +2 位作者 席振铢 蒋欢 刘愿愿 《中南大学学报(自然科学版)》 EI CAS CSCD 北大核心 2016年第10期3477-3482,共6页
为了提高瞬变电磁法三维正演计算速度,一方面,采用矢量有限单元法与积分方程法相混合,使24维单元矩阵的阶数降低到12维,减少大型稀疏矩阵的计算时间,同时增强边界条件的连续性,避免出现伪解;另一方面,引进Flow-through Hankel transfor... 为了提高瞬变电磁法三维正演计算速度,一方面,采用矢量有限单元法与积分方程法相混合,使24维单元矩阵的阶数降低到12维,减少大型稀疏矩阵的计算时间,同时增强边界条件的连续性,避免出现伪解;另一方面,引进Flow-through Hankel transform快速计算技术,解决不同收发距的贝塞尔函数重复计算的问题,进一步节省计算时间。研究结果表明:引入矢量有限单元法的混合法和Flow-through Hankel transform快速计算技术不但能保证3维计算精度,而且能提高计算速度。 展开更多
关键词 瞬变电磁法 3维正演 FLOW-THROUGH hankel transform 矢量有限单元法
在线阅读 下载PDF
模块化多电平换流器的降阶小信号模型研究 被引量:17
19
作者 李正 郝全睿 +1 位作者 尹晓东 王淑颖 《中国电机工程学报》 EI CSCD 北大核心 2018年第12期3638-3650,共13页
建立模块化多电平换流器(modular multilevel converter,MMC)的小信号模型时,通常会考虑子模块电容电压波动、内部环流等状态变量,但却提高了小信号模型的阶数。为了减少多端MMC直流网络进行小信号分析时的模型阶数,该文首先建立MMC... 建立模块化多电平换流器(modular multilevel converter,MMC)的小信号模型时,通常会考虑子模块电容电压波动、内部环流等状态变量,但却提高了小信号模型的阶数。为了减少多端MMC直流网络进行小信号分析时的模型阶数,该文首先建立MMC全阶小信号模型;然后基于平衡理论和Hankel奇异值降阶思想,提出MMC全阶小信号模型的一般降阶方法;建立MMC整流站和逆变站的有效降阶小信号模型;并以三端MMC-HVDC系统为例对降阶前后系统的动态响应特性和不稳定特征根轨迹进行了对比。仿真结果表明:通过选取适当的阶数,MMC降阶模型与全阶模型具有一致性;直流系统在应用降阶小信号模型后,系统的降阶模型与原系统的全阶小信号模型具有相同的动态特性,并且同样能够对系统不稳定进行预估。 展开更多
关键词 模块化多电平换流器 降阶小信号模型 平衡理论 hankel奇异值 动态特性
原文传递
横观各向同性饱和土体三维粘弹性动力分析 被引量:4
20
作者 祝彦知 李冬霞 方志 《岩土力学》 EI CAS CSCD 北大核心 2005年第10期1557-1564,共8页
采用针对横观各向同性饱和土体u-w形式三维粘弹性动力方程,考虑土骨架的粘弹性性质且基于粘弹性理论,通过运用Fourier展开、Laplace和Hankel积分变换方法和引入中间变量,将含有粘弹性参数的六元二阶偏微分运动控制方程组,化为2组各含4... 采用针对横观各向同性饱和土体u-w形式三维粘弹性动力方程,考虑土骨架的粘弹性性质且基于粘弹性理论,通过运用Fourier展开、Laplace和Hankel积分变换方法和引入中间变量,将含有粘弹性参数的六元二阶偏微分运动控制方程组,化为2组各含4个未知变量的常微分方程组,从而给出了柱坐标系下粘弹性横观各向同性饱和土体在非轴对称动力荷载作用下的瞬态反应的土骨架位移分量、孔隙流体相对于土骨架的位移分量瞬态反应一般解。在此基础上,引入初始条件和边界条件,对垂直向和水平向动力荷载作用下半空间边值问题进行了求解。根据动力时域解答的一般解,利用Laplace和Hankel数值逆变换技术,编制了相应的数值计算程序。并进行了实例验证和弹性、粘弹性解的对比分析。结果表明,在进行横观各向同性饱和土体动力分析时,考虑土骨架的粘弹性是必要的。 展开更多
关键词 BIOT波动方程 横观各向同性 饱和土体 粘弹性模型 动力分析 LAPLACE-hankel变换 数值解
在线阅读 下载PDF
上一页 1 2 下一页 到第
使用帮助 返回顶部