期刊文献+
共找到99篇文章
< 1 2 5 >
每页显示 20 50 100
An efficient cubic trigonometric B-spline collocation scheme for the time-fractional telegraph equation 被引量:1
1
作者 Muhammad Yaseen Muhammad Abbas 《Applied Mathematics(A Journal of Chinese Universities)》 SCIE CSCD 2020年第3期359-378,共20页
In this paper,a proficient numerical technique for the time-fractional telegraph equation(TFTE)is proposed.The chief aim of this paper is to utilize a relatively new type of B-spline called the cubic trigonometric B-s... In this paper,a proficient numerical technique for the time-fractional telegraph equation(TFTE)is proposed.The chief aim of this paper is to utilize a relatively new type of B-spline called the cubic trigonometric B-spline for the proposed scheme.This technique is based on finite difference formulation for the Caputo time-fractional derivative and cubic trigonometric B-splines based technique for the derivatives in space.A stability analysis of the scheme is presented to confirm that the errors do not amplify.A convergence analysis is also presented.Computational experiments are carried out in addition to verify the theoretical analysis.Numerical results are contrasted with a few present techniques and it is concluded that the presented scheme is progressively right and more compelling. 展开更多
关键词 Time-fractional telegraph equation finite difference method cubic trigonometric b-splines collocation method Stability CONVERGENCE
下载PDF
Performance of global look-up table strategy in digital image correlation with cubic B-spline interpolation and bicubic interpolation 被引量:4
2
作者 Zhiwei Pan Wei Chen +3 位作者 Zhenyu Jiang Liqun Tang Yiping Liu Zejia Liu 《Theoretical & Applied Mechanics Letters》 CAS CSCD 2016年第3期126-130,共5页
Global look-up table strategy proposed recently has been proven to be an efficient method to accelerate the interpolation, which is the most time-consuming part in the iterative sub-pixel digital image correlation (... Global look-up table strategy proposed recently has been proven to be an efficient method to accelerate the interpolation, which is the most time-consuming part in the iterative sub-pixel digital image correlation (DIC) algorithms. In this paper, a global look-up table strategy with cubic B-spline interpolation is developed for the DIC method based on the inverse compositional Gauss-Newton (IC-GN) algorithm. The performance of this strategy, including accuracy, precision, and computation efficiency, is evaluated through a theoretical and experimental study, using the one with widely employed bicubic interpolation as a benchmark. The global look-up table strategy with cubic B-spline interpolation improves significantly the accuracy of the IC-GN algorithm-based DIC method compared with the one using the bicubic interpolation, at a trivial price of computation efficiency. 展开更多
关键词 Digital image correlation Inverse compositional Gauss-Newton algorithm cubic b-spline interpolation Bicubic interpolation Global look-up table
下载PDF
Redefined Extended Cubic B-Spline Functions for Numerical Solution of Time-Fractional Telegraph Equation
3
作者 Muhammad Amin Muhammad Abbas +2 位作者 Dumitru Baleanu Muhammad Kashif Iqbal Muhammad Bilal Riaz 《Computer Modeling in Engineering & Sciences》 SCIE EI 2021年第4期361-384,共24页
This work is concerned with the application of a redefined set of extended uniform cubic B-spline(RECBS)functions for the numerical treatment of time-fractional Telegraph equation.The presented technique engages finit... This work is concerned with the application of a redefined set of extended uniform cubic B-spline(RECBS)functions for the numerical treatment of time-fractional Telegraph equation.The presented technique engages finite difference formulation for discretizing the Caputo time-fractional derivatives and RECBS functions to interpolate the solution curve along the spatial grid.Stability analysis of the scheme is provided to ensure that the errors do not amplify during the execution of the numerical procedure.The derivation of uniform convergence has also been presented.Some computational experiments are executed to verify the theoretical considerations.Numerical results are compared with the existing schemes and it is concluded that the present scheme returns superior outcomes on the topic. 展开更多
关键词 Extended cubic b-spline redefined extended cubic b-spline time fractional telegraph equation caputo fractional derivative finite difference method CONVERGENCE
下载PDF
Estimation of the water–oil–gas relative permeability curve from immiscible WAG coreflood experiments using the cubic B-spline model
4
作者 Dai-Gang Wang Yong-Le Hu +1 位作者 Jing-Jing Sun Yong Li 《Petroleum Science》 SCIE CAS CSCD 2016年第3期507-516,共10页
Immiscible water-alternating-gas(WAG) flooding is an EOR technique that has proven successful for water drive reservoirs due to its ability to improve displacement and sweep efficiency.Nevertheless,considering the c... Immiscible water-alternating-gas(WAG) flooding is an EOR technique that has proven successful for water drive reservoirs due to its ability to improve displacement and sweep efficiency.Nevertheless,considering the complicated phase behavior and various multiphase flow characteristics,gas tends to break through early in production wells in heterogeneous formations because of overriding,fingering,and channeling,which may result in unfavorable recovery performance.On the basis of phase behavior studies,minimum miscibility pressure measurements,and immiscible WAG coreflood experiments,the cubic B-spline model(CBM) was employed to describe the three-phase relative permeability curve.Using the Levenberg-Marquardt algorithm to adjust the vector of unknown model parameters of the CBM sequentially,optimization of production performance including pressure drop,water cut,and the cumulative gas-oil ratio was performed.A novel numerical inversion method was established for estimation of the water-oil-gas relative permeability curve during the immiscible WAG process.Based on the quantitative characterization of major recovery mechanisms,the proposed method was validated by interpreting coreflood data of the immiscible WAG experiment.The proposed method is reliable and can meet engineering requirements.It provides a basic calculation theory for implicit estimation of oil-water-gas relative permeability curve. 展开更多
关键词 cubic b-spline model Immiscible WAG flooding Relative permeability Numerical inversion
下载PDF
Shape Control and Modification of Rational Cubic B-Spline Curves
5
作者 Zhao Hongsheng Zhang Mandong(Department of Mechanical EngineeringTaiyuan University of Technology, Taiyuan 030024, P. R. China) 《Computer Aided Drafting,Design and Manufacturing》 1999年第1期17-23,共7页
This paper considers the construction of a rational cubic B-spline curve that willinterpolate a sequence of data points x'+ith specified tangent directions at those points. It is emphasisedthat the constraints are... This paper considers the construction of a rational cubic B-spline curve that willinterpolate a sequence of data points x'+ith specified tangent directions at those points. It is emphasisedthat the constraints are purely geometrical and that the pararnetric tangent magnitudes are notassigned as in many' curl'e manipulation methods. The knot vector is fixed and the unknowns are thecontrol points and x'eightsf in this respect the technique is fundamentally different from otherswhere knot insertion is allowed.First. the theoretical result3 for the uniform rational cubic B-spline are presented. Then. in theplanar case. the effect of changes to the tangent at a single point and the acceptable bounds for thechange are established so that all the weights and tangent magnitUdes remain positive. Finally, aninteractive procedure for controlling the shape of a planar rational cubic B-spline curve is presented. 展开更多
关键词 rational cubic b-spline shape control tangent modification
全文增补中
Numerical Solution of Obstacle Problems by B-Spline Functions
6
作者 G.B Loghmani F Mahdifar S.R Alavizadeh 《American Journal of Computational Mathematics》 2011年第2期55-62,共8页
In this study, we use B-spline functions to solve the linear and nonlinear special systems of differential equations associated with the category of obstacle, unilateral, and contact problems. The problem can easily c... In this study, we use B-spline functions to solve the linear and nonlinear special systems of differential equations associated with the category of obstacle, unilateral, and contact problems. The problem can easily convert to an optimal control problem. Then a convergent approximate solution is constructed such that the exact boundary conditions are satisfied. The numerical examples and computational results illustrate and guarantee a higher accuracy for this technique. 展开更多
关键词 Least SQUARE Method uniform b-splines Boundary Value PROBLEMS OBSTACLE PROBLEMS
下载PDF
星载共口径多波束赋形反射面天线设计
7
作者 杨小凤 杨承坤 郭庆功 《强激光与粒子束》 CAS CSCD 北大核心 2024年第9期52-59,共8页
为了满足星载通信中多任务并行需求,提出了一种由格里高利型反射面和三个馈源喇叭天线组成的共口径多波束反射面天线,该天线能够产生两个固定轮廓波束和一个点波束。天线设计是采用射线追踪法确定点波束馈源位置以建立天线基础结构,并通... 为了满足星载通信中多任务并行需求,提出了一种由格里高利型反射面和三个馈源喇叭天线组成的共口径多波束反射面天线,该天线能够产生两个固定轮廓波束和一个点波束。天线设计是采用射线追踪法确定点波束馈源位置以建立天线基础结构,并通过Zernike多项式和Cubic B-splines函数共同对主、副反射面进行赋形优化来完成。为了验证该方法的有效性,对口径为1.1 m的天线进行仿真设计,结果表明两个赋形轮廓波束在Ku收、发频段边缘增益(EOC gain)分别为27.7、28.0、28.0、28.2 dBi,固定点波束在服务区EOC gain不低于34 dBi并且在0~6.5°范围内的扫描波束增益不低于35 dBi。 展开更多
关键词 星载天线 赋形波束 ZERNIKE多项式 cubic b-splines函数 射线追踪法
下载PDF
基于改进人工势场法的RRT^(*)无人船路径规划算法 被引量:1
8
作者 周卫祥 许继强 《中北大学学报(自然科学版)》 CAS 2024年第2期123-131,共9页
为了使RRT^(*)能够更好地适应不同复杂程度的环境,并快速生成一条平滑的较优路径,本文在RRT^(*)算法的基础上引入人工势场法,设计了基于改进人工势场法的RRT^(*)算法。首先,对全局地图进行划分,并进行分区偏置采样;然后,改进节点拓展方... 为了使RRT^(*)能够更好地适应不同复杂程度的环境,并快速生成一条平滑的较优路径,本文在RRT^(*)算法的基础上引入人工势场法,设计了基于改进人工势场法的RRT^(*)算法。首先,对全局地图进行划分,并进行分区偏置采样;然后,改进节点拓展方式,引入障碍物大小因子来改进斥力势场函数,引导新节点的生成;同时,引入自适应变步长策略,根据距离障碍物的远近,以不同的步长拓展路径点。为了使规划路径更符合无人船的航行特性,采用三次非均匀B样条对改进算法生成的路径进行了平滑处理。为了验证本文改进算法的优势,通过设计特殊障碍物环境、简单障碍物环境以及复杂障碍物环境,对比分析了RRT、RRT^(*)、人工势场法和本文算法,发现了本文改进算法生成的路径平均长度短于RRT、RRT^(*)和人工势场法所规划的路径长度,路径规划效率更高,面对不同障碍物环境有更好的适用性。 展开更多
关键词 无人船 RRT^(*) 人工势场法 三次非均匀B样条 路径规划
下载PDF
A New Isogeometric Finite Element Method for Analyzing Structures
9
作者 Pan Su Jiaxing Chen +1 位作者 Ronggang Yang Jiawei Xiang 《Computer Modeling in Engineering & Sciences》 SCIE EI 2024年第11期1883-1905,共23页
High-performance finite element research has always been a major focus of finite element method studies.This article introduces isogeometric analysis into the finite element method and proposes a new isogeometric fini... High-performance finite element research has always been a major focus of finite element method studies.This article introduces isogeometric analysis into the finite element method and proposes a new isogeometric finite element method.Firstly,the physical field is approximated by uniform B-spline interpolation,while geometry is represented by non-uniform rational B-spline interpolation.By introducing a transformation matrix,elements of types C^(0)and C^(1)are constructed in the isogeometric finite element method.Subsequently,the corresponding calculation formats for one-dimensional bars,beams,and two-dimensional linear elasticity in the isogeometric finite element method are derived through variational principles and parameter mapping.The proposed method combines element construction techniques of the finite element method with geometric construction techniques of isogeometric analysis,eliminating the need for mesh generation and maintaining flexibility in element construc-tion.Two elements with interpolation characteristics are constructed in the method so that boundary conditions and connections between elements can be processed like the finite element method.Finally,the test results of several examples show that:(1)Under the same degree and element node numbers,the constructed elements are almost consistent with the results obtained by traditional finite element method;(2)For bar problems with large local field variations and beam problems with variable cross-sections,high-degree and multi-nodes elements constructed can achieve high computational accuracy with fewer degrees of freedom than finite element method;(3)The computational efficiency of isogeometric finite element method is higher than finite element method under similar degrees of freedom,while as degrees of freedom increase,the computational efficiency between the two is similar. 展开更多
关键词 Finite element method isogeometric analysis uniform b-spline non-uniform rational b-spline beam and bar
下载PDF
A Certified Cubic B-Spline Interpolation Method with Tangential Direction Constraints
10
作者 HE Shitao SHEN Liyong +1 位作者 WU Qin YUAN Chunming 《Journal of Systems Science & Complexity》 SCIE EI CSCD 2024年第3期1271-1294,共24页
Curve interpolation with B-spline is widely used in various areas. This problem is classic and recently raised in application scenario with new requirements such as path planning following the tangential vector field ... Curve interpolation with B-spline is widely used in various areas. This problem is classic and recently raised in application scenario with new requirements such as path planning following the tangential vector field under certified error in CNC machining. This paper proposes an algorithm framework to solve Hausdorff distance certified cubic B-spline interpolation problem with or without tangential direction constraints. The algorithm has two stages: The first stage is to find the initial cubic B-spine fitting curve which satisfies the Hausdorff distance constraint;the second stage is to set up and solve the optimization models with certain constraints. Especially, the sufficient conditions of the global Hausdorff distance control for any error bound are discussed, which can be expressed as a series of linear and quadratic constraints. A simple numerical algorithm to compute the Hausdorff distance between a polyline and its B-spline interpolation curve is proposed to reduce our computation.Experimental results are presented to show the advantages of the proposed algorithms. 展开更多
关键词 cubic b-spline Hausdorff distance INTERPOLATION tangential direction
原文传递
带多局部形状参数的三次扩展均匀B样条曲线 被引量:9
11
作者 胡钢 刘哲 +1 位作者 秦新强 戴芳 《西安交通大学学报》 EI CAS CSCD 北大核心 2008年第10期1245-1249,共5页
为了构造带局部形状控制参数的B样条曲线,给出了一组含有λi、μi2个形状参数的四次多项式调配函数,它是三次均匀B样条基函数的新扩展.同时,分析了这组调配函数的性质,并基于调配函数定义了一种新的带有λi、μi2个局部形状控制参数的... 为了构造带局部形状控制参数的B样条曲线,给出了一组含有λi、μi2个形状参数的四次多项式调配函数,它是三次均匀B样条基函数的新扩展.同时,分析了这组调配函数的性质,并基于调配函数定义了一种新的带有λi、μi2个局部形状控制参数的分段多项式样条曲线,其以三次均匀B样条曲线为特殊情形.最后,讨论了新曲线在曲线造型中的应用,并给出了相应扩展曲面的定义.造型实例表明,新曲线不仅具有灵活的局部形状可调性和更强的描述能力,而且可以在不改变曲线G1连续性和不影响曲线其他各段形状的同时,通过改变局部形状参数对曲线每段的形状进行多种方式的局部调整,为曲线和曲面的设计提供了一种有效的新方法. 展开更多
关键词 三次均匀B样条 调配函数 局部形状参数 曲线设计
下载PDF
三次均匀B样条曲线的新扩展及应用 被引量:15
12
作者 胡钢 刘哲 徐华楠 《计算机工程与应用》 CSCD 北大核心 2008年第32期161-164,共4页
给出了一组含有2个形状参数λ_i,μ_i的三次多项式调配函数,它是三次均匀B样条基函数的扩展;分析了这组调配函数的性质,基于此组调配函数定义了一种带2个局部形状控制参数λ_i,μ_i的分段多项式样条曲线,它以三次均匀B样条曲线为特殊情... 给出了一组含有2个形状参数λ_i,μ_i的三次多项式调配函数,它是三次均匀B样条基函数的扩展;分析了这组调配函数的性质,基于此组调配函数定义了一种带2个局部形状控制参数λ_i,μ_i的分段多项式样条曲线,它以三次均匀B样条曲线为特殊情形。新曲线不仅具有灵活的局部形状可调性和更强的描述能力,而且可以在不改变曲线G^1连续性和不影响曲线其他各段形状的同时,通过改变局部形状参数对曲线每段的形状进行多种方式的局部调整。最后讨论了新曲线在曲线造型中的应用,并给出了一个扩展曲面的定义。实例表明,新扩展曲线为曲线/曲面的设计提供了一种有效的新方法。 展开更多
关键词 三次均匀B样条 调配函数 局部形状参数 扩展 曲线设计
下载PDF
带多个形状参数的三次均匀B样条曲线的扩展 被引量:17
13
作者 夏成林 邬弘毅 +1 位作者 郑兴国 彭凯军 《工程图学学报》 CSCD 北大核心 2011年第2期73-79,共7页
通过构造两类带多个形状参数的调配函数,生成三次均匀B样条基函数的扩展。基于给出的调配函数定义了两类带多个形状参数的分段多项式曲线。这些曲线具有三次均匀B样条曲线的绝大多数重要性质,能达到GC1或GC2连续。改变形状参数的值可以... 通过构造两类带多个形状参数的调配函数,生成三次均匀B样条基函数的扩展。基于给出的调配函数定义了两类带多个形状参数的分段多项式曲线。这些曲线具有三次均匀B样条曲线的绝大多数重要性质,能达到GC1或GC2连续。改变形状参数的值可以独立地调控各子段的端点的位置及其切矢的长度,对曲线进行整体或局部调整,甚至直接插值任何所需的控制点。 展开更多
关键词 三次均匀B样条 调配函数 端点位置参数 切矢 插值
下载PDF
基于离散曲率的三次均匀B样条的局部光顺算法 被引量:10
14
作者 喻德生 程程 《浙江大学学报(理学版)》 CAS CSCD 北大核心 2011年第5期511-517,共7页
在研究B样条曲线节点的曲率和对应控制点的离散曲率之间关系的基础上,引入了一种新的离散曲率——第二离散曲率的概念,得出了三次均匀B样条曲线节点的曲率和对应控制点的第二离散曲率成正比的结论,并给出了基于第二离散曲率三次均匀B样... 在研究B样条曲线节点的曲率和对应控制点的离散曲率之间关系的基础上,引入了一种新的离散曲率——第二离散曲率的概念,得出了三次均匀B样条曲线节点的曲率和对应控制点的第二离散曲率成正比的结论,并给出了基于第二离散曲率三次均匀B样条曲线的光顺算法.该算法通过直接调整控制点的第二离散曲率进行曲线的光顺,从而使光顺过程更为简洁、更具几何直观性.算例表明,该算法具有较好的光顺效果. 展开更多
关键词 离散曲率 差分 光顺 三次均匀B样条
下载PDF
一种三次均匀B样条曲线的轨迹规划方法 被引量:20
15
作者 毛征宇 刘中坚 《中国机械工程》 EI CAS CSCD 北大核心 2010年第21期2569-2572,2577,共5页
针对复杂曲线的数控加工,提出一种新的插补方法。首先采用一种三次B样条曲线的重叠拼接算法实时地对复杂曲线进行拟合,进而综合弓高误差、速度和加速度等因素,给出具有自适应调整能力的插补步长确定算法,该算法在提高轮廓加工精度的同时... 针对复杂曲线的数控加工,提出一种新的插补方法。首先采用一种三次B样条曲线的重叠拼接算法实时地对复杂曲线进行拟合,进而综合弓高误差、速度和加速度等因素,给出具有自适应调整能力的插补步长确定算法,该算法在提高轮廓加工精度的同时,可减小加工过程中的冲击。最后,对插补周期内节点的轨迹进行了规划,以确保运动轨迹满足速度、加速度以及加加速度的平滑约束条件。仿真实验结果表明,该方法在实时插补过程中,可以保证复杂曲线插补加工的高速与高精度,且具有很好的速度、加速度以及加加速度平滑性。 展开更多
关键词 数控加工 复杂曲线 三次均匀B样条 插补 轨迹规划
下载PDF
带局部形状参数的三次均匀B样条曲线的扩展 被引量:16
16
作者 徐岗 汪国昭 《计算机研究与发展》 EI CSCD 北大核心 2007年第6期1032-1037,共6页
带形状参数的B样条曲线的构造已成为计算机辅助几何设计中的热点问题.为了使形状参数具有局部修改功能,给出了两类带局部形状参数的调配函数,它们都是三次均匀B样条基函数的扩展.基于给出的调配函数,定义了两种带局部形状参数的分段多... 带形状参数的B样条曲线的构造已成为计算机辅助几何设计中的热点问题.为了使形状参数具有局部修改功能,给出了两类带局部形状参数的调配函数,它们都是三次均匀B样条基函数的扩展.基于给出的调配函数,定义了两种带局部形状参数的分段多项式曲线.可以通过改变局部形状参数的取值对曲线进行局部调整.调整形状参数可使三次多项式曲线在三次均匀B样条曲线远离控制多边形的一侧摆动,而四次多项式曲线在三次均匀B样条曲线的两侧摆动.最后讨论了它们在曲线设计及曲线插值中的应用.造型实例表明,该类曲线在计算机辅助几何设计中具有重要的应用价值. 展开更多
关键词 三次均匀B样条 局部形状参数 曲线设计 插值
下载PDF
超光谱成像仪图像均匀性校正 被引量:7
17
作者 张春雷 向阳 《中国光学》 EI CAS 2013年第4期584-590,共7页
为获得高质量新型超光谱成像仪图像,根据新型超光谱成像仪图像特点,对其图像均匀性校正进行研究。利用转台旋转模拟超光谱成像仪对地面推扫,对景物成像获得超光谱成像仪三维数据;采用太阳作为新型超光谱成像仪图像均匀性校正光源,应用... 为获得高质量新型超光谱成像仪图像,根据新型超光谱成像仪图像特点,对其图像均匀性校正进行研究。利用转台旋转模拟超光谱成像仪对地面推扫,对景物成像获得超光谱成像仪三维数据;采用太阳作为新型超光谱成像仪图像均匀性校正光源,应用相应校正系数对超光谱成像仪的数据进行处理。结果表明:所采用校正方法对新型超光谱成像仪具有很好校正效果,同时也验证了新型成像光谱仪具有较高成像质量;单像元校正时间为2×10-6s,大数据量图像校正需较长时间;利用太阳照射漫反射板作为扩展源对亮度均匀物体图像进行校正,校正均匀性小于4%。 展开更多
关键词 超光谱 成像光谱仪 推扫 数据立方 非均匀性 太阳
下载PDF
基于均匀三次B样条曲面的斜齿轮副精确建模 被引量:4
18
作者 李学艺 李超超 +1 位作者 李三帅 江守波 《机械传动》 CSCD 北大核心 2012年第8期44-47,共4页
基于均匀三次B样条曲面的基本理论,在ANSYS系统中得到了精确的斜齿轮轮齿曲面,建立了完整的斜齿轮副装配模型,为斜齿轮有限元接触分析和疲劳分析的可靠性和精度提供了保证。基于B样条方法的曲面造型方法研究,对于在计算机上建立一些复... 基于均匀三次B样条曲面的基本理论,在ANSYS系统中得到了精确的斜齿轮轮齿曲面,建立了完整的斜齿轮副装配模型,为斜齿轮有限元接触分析和疲劳分析的可靠性和精度提供了保证。基于B样条方法的曲面造型方法研究,对于在计算机上建立一些复杂自由型曲面的零件模型具有重要的参考意义。 展开更多
关键词 斜齿轮 均匀三次B样条 ANSYS 造型方法
下载PDF
立方抛物线形渠道水力最优断面的计算 被引量:34
19
作者 魏文礼 杨国丽 《武汉大学学报(工学版)》 CAS CSCD 北大核心 2006年第3期49-51,共3页
根据明渠均匀流理论,推求了立方抛物线形渠道正常水深与流量关系的计算公式及其水力最优断面所满足的条件.与其他形式断面渠道的比较结果表明,立方抛物线形断面渠道也是一种水力较优的输水渠道.
关键词 立方抛物线形渠道 明渠均匀流 水力最优断面
下载PDF
基于小波的不均匀光照的校正 被引量:14
20
作者 王宇 龙兴武 《光学技术》 EI CAS CSCD 北大核心 2005年第5期726-728,共3页
给出了一种减小图像中的不均匀光照的方法。利用基于小波的非线性回归来估计降质图像的光照分布,再通过从原图中减去光照分布图达到减小不均匀光照的目的。对此方法和常用的双立方插值法的计算量进行了分析比较,证明此方法计算量小。最... 给出了一种减小图像中的不均匀光照的方法。利用基于小波的非线性回归来估计降质图像的光照分布,再通过从原图中减去光照分布图达到减小不均匀光照的目的。对此方法和常用的双立方插值法的计算量进行了分析比较,证明此方法计算量小。最后取三幅图像进行了实验验证。 展开更多
关键词 不均匀光照 小波非线性回归 双立方插值
下载PDF
上一页 1 2 5 下一页 到第
使用帮助 返回顶部