期刊文献+
共找到581篇文章
< 1 2 30 >
每页显示 20 50 100
三维抛物方程模型海底边界处理与并行算法研究
1
作者 瞿宇轩 程兴华 +1 位作者 周晓兰 邹劭芬 《声学技术》 北大核心 2025年第2期181-188,共8页
构建一种迅速而准确的三维声场模型以满足实际海洋条件,始终是当前科学研究的重点目标。作为一种步进求解的手段,抛物方程近似的方法在已知初始地的声场后以间距为单位步进,最终能得到整个波导中的声场,这对于处理声场三维变化和声道水... 构建一种迅速而准确的三维声场模型以满足实际海洋条件,始终是当前科学研究的重点目标。作为一种步进求解的手段,抛物方程近似的方法在已知初始地的声场后以间距为单位步进,最终能得到整个波导中的声场,这对于处理声场三维变化和声道水平变化相关的场景是十分方便的。但相对地,三维空间中的计算复杂度和计算范围要远大于二维空间,这就严重制约了模型的计算精度和速度。为了提高现有抛物方程模型的计算精度和运算速度,文章对同时考虑距离、深度和水平方位角的水下声传播模型,即三维抛物方程模型进行了推导并分析其数学特性。同时通过在海底边界插入完全匹配层的方法,对模型进行了改进,提出一种计算精度更高的抛物方程模型。此外,为了减小三维声场计算的时间,文章研究改进了模型的可并行性,在FOR3D串行算法的基础上基于多线程的共享内存技术OpenMP开发了三维抛物方程模型的并行算法,并与普通吸收层在不同线程数下的运行时间进行对比。结果表明,基于完全匹配层模型的并行算法与原模型相比,无论是在加速比还是运算速度上都有很大的提高,具有重要的工程意义。 展开更多
关键词 三维抛物方程 人工吸收层 完全匹配层 并行算法
在线阅读 下载PDF
基于deluxeFETI-DP并行算法的矩形太阳翼展开动力学高效仿真
2
作者 兰靖杰 张世雄 +1 位作者 荣吉利 刘铖 《宇航学报》 CSCD 北大核心 2024年第12期1897-1909,共13页
随着柔性太阳翼等可展开空间结构向轻量化和大型化方向发展,其展开过程常面临大转动和大变形等非线性动力学问题。通常采用柔性多体系统动力学方法构建高维非线性模型,以精确预测此类复杂动力学行为。然而,计算效率低下已成为限制该数... 随着柔性太阳翼等可展开空间结构向轻量化和大型化方向发展,其展开过程常面临大转动和大变形等非线性动力学问题。通常采用柔性多体系统动力学方法构建高维非线性模型,以精确预测此类复杂动力学行为。然而,计算效率低下已成为限制该数值模拟方法的主要瓶颈。为解决这一问题,基于有限元撕裂与连接-对偶原始(FETI-DP)并行算法,引入deluxe缩放方法用于柔性多体系统动力学仿真。该方法通过耦合子域界面的刚度矩阵,利用矩阵组合构造加权因子,平衡材料系数差异导致的刚度矩阵量级不均,从而优化界面方程的条件数,减少迭代次数。最后,基于deluxe FETI-DP并行算法对矩形太阳翼展开动力学进行高效仿真,并与串行算法及传统FETI-DP算法进行对比研究。数值结果表明,deluxe FETI-DP算法的仿真计算效率较串行算法提高了70.98%,较传统FETI-DP算法提高了64.16%。上述并行算法研究为大型矩形太阳翼展开动力学的高效仿真提供了理论参考。 展开更多
关键词 deluxe FETI-DP 并行算法 矩形太阳翼 展开动力学
在线阅读 下载PDF
混凝土细观力学分析程序中的快速算法与并行算法设计 被引量:6
3
作者 吴建平 王正华 +2 位作者 朱星明 马怀发 李晓梅 《计算力学学报》 EI CAS CSCD 北大核心 2008年第3期352-358,共7页
针对一套混凝土细观力学分析程序,在分析其计算方法与计算效率的不足之后,提出了采用稀疏矩阵与稀疏向量技术来高效实现有限元刚度矩阵装配过程的算法,并采用双门槛不完全Cholesky分解预条件技术与CG法相结合来高效地求解稀疏线性方程... 针对一套混凝土细观力学分析程序,在分析其计算方法与计算效率的不足之后,提出了采用稀疏矩阵与稀疏向量技术来高效实现有限元刚度矩阵装配过程的算法,并采用双门槛不完全Cholesky分解预条件技术与CG法相结合来高效地求解稀疏线性方程组。之后,从整体上提出了一个将有限单元分布与未知量分布有机结合的并行算法设计方案,并分别针对刚度矩阵装配、双门槛不完全Cholesky分解、稀疏矩阵与稠密向量相乘、稀疏向量相加等核心算法,进行了相应的并行算法设计。最后,在由每节点2 CPU的8个Intel Xeon节点采用千兆以太网连成的机群上,针对两个混凝土数值试样进行了数值实验,第一个试样含44117个网格点与53200个有限单元,第二个试样含71013个网格点与78800个有限单元;对第一个试样,原串行程序进行全程567次加载计算需要984.83小时约41天,采用文中串行算法后,模拟时间减少到22531秒约6.26小时,采用并行算法在16个CPU上的模拟时间进一步降为3860秒约1.07小时。对第二个试样,原串行程序进行全程94次加载计算需要467.19小时约19.5天,采用文中串行算法后,模拟时间减少到11453秒约3.18小时,采用并行算法在16个CPU上的模拟时间进一步降为1704秒约28.4分钟。串行算法的改进与并行算法的设计大大缩短了计算时间,对加快混凝土力学性能的分析研究具有重要意义。 展开更多
关键词 混凝土细观力学 有限元 稀疏线性方程组 预条件 并行算法
在线阅读 下载PDF
五点差分格式求解泊松方程并行算法的研究 被引量:15
4
作者 廖臣 祝大军 刘盛纲 《电子科技大学学报》 EI CAS CSCD 北大核心 2008年第1期81-83,127,共4页
以二维静电场泊松方程数值求解的串行算法(雅可比迭代、超松弛迭代)为基础,提出了五点差分格式超松弛迭代(SOR)求解二维静电场泊松方程的并行算法,通过与雅可比迭代(Jacobi)并行算法的时间复杂度、加速比和空间复杂度进行对比,得出超松... 以二维静电场泊松方程数值求解的串行算法(雅可比迭代、超松弛迭代)为基础,提出了五点差分格式超松弛迭代(SOR)求解二维静电场泊松方程的并行算法,通过与雅可比迭代(Jacobi)并行算法的时间复杂度、加速比和空间复杂度进行对比,得出超松弛迭代的并行算法具有更低的时间复杂度、空间复杂度和更高的加速比与效率。通过实验验证,CHIPIC软件的泊松模块宜采用超松弛迭代并行算法。 展开更多
关键词 雅可比迭代 并行算法 泊松 超松弛迭代
在线阅读 下载PDF
三维非结构网格DSMC并行算法及应用研究 被引量:7
5
作者 王学德 伍贻兆 +1 位作者 夏健 谭俊杰 《宇航学报》 EI CAS CSCD 北大核心 2007年第6期1500-1505,共6页
基于PC-CLUSTER机群并行体系结构与消息传递库MPI并行环境,研究了三维非结构网格DSMC并行算法。提出一类基于结构背景网格上的非结构网格动态分区策略,保证各子区域的分子数量大致相等,实现计算进程间的动态负载平衡。利用MPI库函数构... 基于PC-CLUSTER机群并行体系结构与消息传递库MPI并行环境,研究了三维非结构网格DSMC并行算法。提出一类基于结构背景网格上的非结构网格动态分区策略,保证各子区域的分子数量大致相等,实现计算进程间的动态负载平衡。利用MPI库函数构造了两类符合DSMC并行原理的通讯法:单步通讯法与多步通讯法。采用单控制多数据流(SPMD)以及Master/Slave并行模式,设计了三维非结构网格DSMC整体并行算法。给出了跟踪模拟分子在四面体网格间迁移运动的详细计算过程。最后对全尺寸航天飞机高超声速绕流进行了并行模拟,验证并行算法的有效性。 展开更多
关键词 非结构网格 直接物理模拟 DSMC 并行算法 动态负载平衡 MPI
在线阅读 下载PDF
非线性动力方程精细积分级数解的并行算法 被引量:9
6
作者 李渊印 金先龙 +1 位作者 张晓云 郭毅之 《上海交通大学学报》 EI CAS CSCD 北大核心 2006年第10期1809-1812,共4页
非线性动力方程通过变量变换可以转化为一阶微分方程,该方程的解由表示初值影响的齐次方程解和反映荷载作用的积分之和组成.其中:第一项用指数矩阵计算;第二项在文中采用级数解计算(设计了3种相应的并行算法),算法1对级数解的每一项先... 非线性动力方程通过变量变换可以转化为一阶微分方程,该方程的解由表示初值影响的齐次方程解和反映荷载作用的积分之和组成.其中:第一项用指数矩阵计算;第二项在文中采用级数解计算(设计了3种相应的并行算法),算法1对级数解的每一项先做若干个向量的线性组合,再做矩阵向量乘1次;算法2与算法1原理相同,只是将矩阵的幂运算转换成乘积;算法3先做若干个矩阵向量乘,再做若干个向量的线性组合.算法1的并行效率最好,但存储空间需求大,不利于大型结构的求解.算法2、3利用动力方程的稀疏变换改善了算法1的不足,算法3中级数解每一项计算均在其前一项基础上进行,一般能比算法2节省时间.最后,给出了算例验证,三种算法都获得了较好的加速比. 展开更多
关键词 非线性动力系统 精细积分法 显式级数解 并行算法
在线阅读 下载PDF
一类虚拟边界预条件多重网格并行算法 被引量:7
7
作者 卫加宁 王伟沧 +1 位作者 皮新明 章社生 《武汉理工大学学报(交通科学与工程版)》 北大核心 2001年第1期4-7,共4页
对网格方程组作多重网格区域分裂并行计算 ,以消去法为预条件 ,讨论了 Gauss消去法的迭代性质 ,指出该迭代法有效改善通常迭代在子域信息传递过程中的障碍性问题 ,提高多重网格区域分裂并行计算效率 .
关键词 多重网格 区域分裂 预条件 并行算法 Gauss消去法 预处理迭代法 边值问题
在线阅读 下载PDF
电各向异性介质FDTD并行算法的研究 被引量:9
8
作者 杨利霞 葛德彪 +1 位作者 郑奎松 魏兵 《电波科学学报》 EI CSCD 北大核心 2006年第1期43-48,共6页
基于消息传递模式的网络并行计算系统和区域分割技术成功地实现了电各向异性介质FDTD并行算法,并用此程序代码计算了电大目标的RCS。经测试,该程序代码并行效率达到87%。根据电各向异性介质FDTD迭代式,电场某一节点的计算,需牵涉到其周... 基于消息传递模式的网络并行计算系统和区域分割技术成功地实现了电各向异性介质FDTD并行算法,并用此程序代码计算了电大目标的RCS。经测试,该程序代码并行效率达到87%。根据电各向异性介质FDTD迭代式,电场某一节点的计算,需牵涉到其周围28个节点的电磁场值,详细分析了该情况下并行算法中的数据通讯规律,并实现了局域网内各节点机的协同并行计算。数值结果表明了该算法和程序的正确性及优越性。 展开更多
关键词 各向异性介质 各向同性介质 并行算法 FDTD方法
在线阅读 下载PDF
自适应多层快速多极子算法及其并行算法 被引量:7
9
作者 袁军 邱扬 +2 位作者 刘其中 郭景丽 谢拥军 《电波科学学报》 EI CSCD 北大核心 2008年第3期455-459,共5页
在多层快速多极子算法的基础上提出了一种改进的电大问题电磁散射快速算法及其并行算法。算法中将基函数和权函数分别用不同空间位置上的点源函数展开,使各部分的积分计算得到大大简化,所有转移过程可由快速傅里叶变换计算完成,同时还... 在多层快速多极子算法的基础上提出了一种改进的电大问题电磁散射快速算法及其并行算法。算法中将基函数和权函数分别用不同空间位置上的点源函数展开,使各部分的积分计算得到大大简化,所有转移过程可由快速傅里叶变换计算完成,同时还给出了该算法的并行化算法。数值结果说明了算法的有效性。 展开更多
关键词 多层快速多极子 自适应 电磁散射 并行算法
在线阅读 下载PDF
高效并行算法的设计与实现 被引量:6
10
作者 刘兴平 莫则尧 +2 位作者 雷光耀 张宝琳 张景琳 《高技术通讯》 EI CAS CSCD 1998年第8期15-18,共4页
在三类不同体系结构的并行多处理机系统上,通过求解大型稀疏线性代数方程组的高效并行迭代算法的设计与实现,指出了高效并行算法的设计必须与某类高性能并行计算机相结合,才能使得同一解算问题在不同类型的高性能并行机上获得最佳性... 在三类不同体系结构的并行多处理机系统上,通过求解大型稀疏线性代数方程组的高效并行迭代算法的设计与实现,指出了高效并行算法的设计必须与某类高性能并行计算机相结合,才能使得同一解算问题在不同类型的高性能并行机上获得最佳性能价格比。 展开更多
关键词 并行计算机 并行算法 并行效率 设计
在线阅读 下载PDF
多层嵌套动态规划并行算法在梯级水库优化调度中的应用 被引量:10
11
作者 蒋志强 纪昌明 +2 位作者 孙平 王丽萍 张验科 《中国农村水利水电》 北大核心 2014年第9期70-75,共6页
为增强动态规划算法对梯级水库群联合优化调度问题的求解能力,基于多层嵌套的思想,提出了多层嵌套多维动态规划算法。通过与传统多维动态规划算法的对比分析,发现该方法能在内存占用量以及程序复杂度2方面有效缓解多维动态规划算法应用... 为增强动态规划算法对梯级水库群联合优化调度问题的求解能力,基于多层嵌套的思想,提出了多层嵌套多维动态规划算法。通过与传统多维动态规划算法的对比分析,发现该方法能在内存占用量以及程序复杂度2方面有效缓解多维动态规划算法应用于梯级水库群联合优化调度时的维数灾问题。针对该方法在运行时间上的缺陷,引入并行处理技术,利用状态离散点间的计算独立性构建了多层嵌套多维动态规划并行算法。以李仙江流域三库梯级为背景进行了实例研究,并从运行时间、并行加速比、并行效率以及CPU核数等方面对该并行算法进行了详细的性能分析。结果表明在现有的计算条件下该并行算法能有效提高计算效率和缓解多维动态规划维数灾问题。 展开更多
关键词 梯级水库 多维动态规划 并行算法 优化调度 李向江流域
在线阅读 下载PDF
二维叠前深度偏移的并行算法 被引量:6
12
作者 王真理 杨长春 +4 位作者 张洪宙 张得顺 符基昕 匡兵 王建立 《石油地球物理勘探》 EI CSCD 北大核心 1997年第6期784-788,817,共6页
王真理,杨长春,张洪审,张得顺,符基昕,匡兵,王建立.二维叠前深度偏移的并行算法.石油地球物理勘探,1997,32(6):784~788,817鉴于PVM软件是一种几乎能够在所有流行的计算机体系结构上运行的并行环境,本文提出了单程序多... 王真理,杨长春,张洪审,张得顺,符基昕,匡兵,王建立.二维叠前深度偏移的并行算法.石油地球物理勘探,1997,32(6):784~788,817鉴于PVM软件是一种几乎能够在所有流行的计算机体系结构上运行的并行环境,本文提出了单程序多数据流(SPMD)方式的二维叠前深度偏移的并行其法,即走时表计算采用炮点循环并行方式,成像过程采用炮检距并行算法。采用PVM并行环境提供的消息传递方式进行进程间通讯,在SGIPowerChallenge和IBMSP2两种并行机上实现了其二维叠前深度偏移的并行其法。测试结果表明:本文的二维叠前深度偏移并行化算法效率是比较高的,在SGIPowerChallenge机上,当CPU数为4时,测得的加速比达3.91,并行效率达到0.98;在IBMSP2机上,当书点数为3时,测得的加速比为2.94,效率达到0.98。 展开更多
关键词 叠前深度偏移 并行算法 地质体 区域构造
在线阅读 下载PDF
USB数据传输中CRC校验码的并行算法实现 被引量:11
13
作者 程军 陈贵灿 姜飞 《微电子学与计算机》 CSCD 北大核心 2003年第3期77-80,共4页
文章介绍了用于USB总线数据传输的CRC校验的原理和算法,并且采用并行电路实现USB2.0中的CRC产生和CRC校验,与传统的串行电路实现相比,并行电路实现方法虽然在芯片面积上大于串行电路实现,但由于降低了时钟频率,电路更容易综合实现,并且... 文章介绍了用于USB总线数据传输的CRC校验的原理和算法,并且采用并行电路实现USB2.0中的CRC产生和CRC校验,与传统的串行电路实现相比,并行电路实现方法虽然在芯片面积上大于串行电路实现,但由于降低了时钟频率,电路更容易综合实现,并且大大降低了功耗,有利于低功耗电路设计。 展开更多
关键词 USB 数据传输 CRC校验码 并行算法 通用串行总线 循环冗余校验码
在线阅读 下载PDF
二维自适应非结构网格DSMC并行算法研究 被引量:5
14
作者 王学德 伍贻兆 +1 位作者 夏健 林晓宏 《计算力学学报》 CAS CSCD 北大核心 2009年第2期276-281,共6页
研究了二维自适应非结构网格DSMC并行算法实现的过程。首先提出了一类非结构网格自适应策略,有效降低了网格尺度对计算结果的影响,提高了流场的分辨率;然后基于PC-CLUSTER群机并行体系结构与消息传递库MPI并行环境,利用分区并行思想,设... 研究了二维自适应非结构网格DSMC并行算法实现的过程。首先提出了一类非结构网格自适应策略,有效降低了网格尺度对计算结果的影响,提高了流场的分辨率;然后基于PC-CLUSTER群机并行体系结构与消息传递库MPI并行环境,利用分区并行思想,设计了非结构网格DSMC并行算法,节约了计算时间。利用For-tran90的动态分配内存技术编制了通用计算程序;最后对过渡流域高超声绕流进行了数值模拟,计算结果初步验证了算法的可行性与有效性。 展开更多
关键词 自适应非结构网格 DSMC 并行算法 MPI
在线阅读 下载PDF
整体最小二乘参数估计的并行算法 被引量:9
15
作者 康正九 石毅明 胡保生 《西安交通大学学报》 EI CAS CSCD 北大核心 1998年第3期1-4,共4页
基于递推最小二乘的逆QR分解方法,给出一种整体最小二乘参数估计的递推算法,并利用已获得的并行实现方法进行并行实现,最后给出了仿真结果.文中给出的并行计算方法使得在实际中应用整体最小二乘进行参数估计变得可行.
关键词 整体最小二乘 参数估计 并行算法 控制论
在线阅读 下载PDF
从图像中快速检测直线的并行算法 被引量:5
16
作者 虞凡 吴惠思 +1 位作者 覃征 杨博 《西安交通大学学报》 EI CAS CSCD 北大核心 2006年第12期1370-1373,1387,共5页
提出了一种在具有可扩展机群体系结构的通用超级并行机环境下的快速直线检测算法.采用最优域划分法将原图像进行N×N等面积划分,并映射到并行系统各节点的局部坐标下进行直线检测.通过先后两次在Hough变换的扫描过程中提前引入阈值... 提出了一种在具有可扩展机群体系结构的通用超级并行机环境下的快速直线检测算法.采用最优域划分法将原图像进行N×N等面积划分,并映射到并行系统各节点的局部坐标下进行直线检测.通过先后两次在Hough变换的扫描过程中提前引入阈值作用,可有效减小计算的复杂度.同时证明了,划分后的子图像采用原阈值的1/(2N-1)进行检测,可满足原图像中的直线目标不丢失的要求.实验结果表明,所提算法具备高鲁棒性和强抗噪能力,能有效提升加速比,该比值的最大值可为25.320. 展开更多
关键词 并行算法 图像处理 直线检测
在线阅读 下载PDF
遥感图像PCA融合的并行算法研究与实现 被引量:8
17
作者 胡冰 周海芳 +1 位作者 王攀峰 刘衡竹 《微电子学与计算机》 CSCD 北大核心 2006年第10期153-155,161,共4页
文章针对已有的遥感图像PCA融合串行算法,提出了一种基于数据并行的新的PCA融合并行算法P-PCA,并对该并行算法进行了通信优化,在机群系统上进行了实现。针对IKONOS图像进行实验的结果表明该算法可获得良好的并行加速比,并行效率较高,为... 文章针对已有的遥感图像PCA融合串行算法,提出了一种基于数据并行的新的PCA融合并行算法P-PCA,并对该并行算法进行了通信优化,在机群系统上进行了实现。针对IKONOS图像进行实验的结果表明该算法可获得良好的并行加速比,并行效率较高,为遥感图像的实际应用提供了有益的指导和借鉴。 展开更多
关键词 遥感图像融合 PCA 并行算法
在线阅读 下载PDF
海量数据上挖掘关联规则的并行算法 被引量:5
18
作者 张兆功 李建中 张艳秋 《哈尔滨工业大学学报》 EI CAS CSCD 北大核心 2004年第5期561-565,共5页
针对目前关联规则挖掘算法中数据库规模很大时算法执行时间太长的问题.指出了并行计算是解决该问题的一个有效方法.利用新提出的可以忽略仅仅在少于1/4的结点机上的局部频繁项集,给出了一种新的并行随机抽样方法,并利用机群并行计算机... 针对目前关联规则挖掘算法中数据库规模很大时算法执行时间太长的问题.指出了并行计算是解决该问题的一个有效方法.利用新提出的可以忽略仅仅在少于1/4的结点机上的局部频繁项集,给出了一种新的并行随机抽样方法,并利用机群并行计算机的自治能力和I/O高度并行的特点,提高了抽样算法对海量数据的处理能力和效率.理论分析和实验数据显示,该算法的加速比接近于处理机的个数p,通信复杂性为处理机的个数p的对数,具有良好的扩展性和海量处理能力,且精确度较高. 展开更多
关键词 海量数据 关联规则 并行算法 数据挖掘 数据库
在线阅读 下载PDF
基于投影寻踪的高光谱图像异常检测并行算法 被引量:5
19
作者 王维 赵慧洁 董超 《北京航空航天大学学报》 EI CAS CSCD 北大核心 2009年第3期342-346,共5页
投影寻踪方法能有效提取数据中的非高斯结构凸显异常信息,但在求解最优投影方向时存在计算量大、运行时间长的问题,为提高处理效率,提出一种机群环境下的并行算法.选用偏度和峰度组合作为投影指标,将所有像素光谱作为特定投影方向集依... 投影寻踪方法能有效提取数据中的非高斯结构凸显异常信息,但在求解最优投影方向时存在计算量大、运行时间长的问题,为提高处理效率,提出一种机群环境下的并行算法.选用偏度和峰度组合作为投影指标,将所有像素光谱作为特定投影方向集依次搜索,求解最优投影方向.在并行计算各候选方向投影指标时,分割图像数据分布存储于各机群结点,数据子块朝候选方向并行投影后,将指标计算式变形分解,使各结点在指标计算过程中所需数据均为本地数据,解决数据局部性问题,并采用一种"轮流作主"的机制提高算法负载均衡程度.利用实用型模块化高光谱仪数据在机群系统上进行测试,达到了较好的加速效果,表明该并行算法具有良好的并行性能. 展开更多
关键词 遥感 图像异常检测 并行算法 投影寻踪 特定投影方向集
在线阅读 下载PDF
FDTD并行算法研究:电和磁本构参数均为各向异性情形 被引量:3
20
作者 杨利霞 葛德彪 +2 位作者 魏兵 郑奎松 葛宁 《电子学报》 EI CAS CSCD 北大核心 2006年第9期1703-1707,共5页
研究并实现了三维电与磁本构参数均为各向异性介质时的时域有限差分(ANI-FDTD)并行算法.根据电磁各向异性介质时域有限差分(ANI-FDTD)迭代式,电场和磁场某一节点场值的计算均涉及到其周围邻近28个场分量节点值.详细分析了该情况下各向... 研究并实现了三维电与磁本构参数均为各向异性介质时的时域有限差分(ANI-FDTD)并行算法.根据电磁各向异性介质时域有限差分(ANI-FDTD)迭代式,电场和磁场某一节点场值的计算均涉及到其周围邻近28个场分量节点值.详细分析了该情况下各向异性介质时域有限差分(ANI-FDTD)并行算法中的数据通讯规律.通过采用相邻子域重叠半个网格的区域划分方式,传递分界面两边相关节点的场值,实现了局域网内各节点机的协同并行计算.数值计算结果表明了该算法和程序的正确性. 展开更多
关键词 各向异性介质 各向同性介质 并行算法 时域有限差分(FDTD)方法
在线阅读 下载PDF
上一页 1 2 30 下一页 到第
使用帮助 返回顶部