期刊文献+
共找到3,007篇文章
< 1 2 151 >
每页显示 20 50 100
Linear numerical calculation method for obtaining critical point,pore fluid,and framework parameters of gas-bearing media 被引量:3
1
作者 牛滨华 孙春岩 +2 位作者 闫国英 杨维 刘畅 《Applied Geophysics》 SCIE CSCD 2009年第4期319-326,393,共9页
Up to now, the primary method for studying critical porosity and porous media are experimental measurements and data analysis. There are few references on how to numerically calculate porosity at the critical point, p... Up to now, the primary method for studying critical porosity and porous media are experimental measurements and data analysis. There are few references on how to numerically calculate porosity at the critical point, pore fluid-related parameters, or framework-related parameters. So in this article, we provide a method for calculating these elastic parameters and use this method to analyze gas-bearing samples. We first derive three linear equations for numerical calculations. They are the equation of density p versus porosity Ф, density times the square of compressional wave velocity p Vp^2 versus porosity, and density times the square of shear wave velocity pVs^2 versus porosity. Here porosity is viewed as an independent variable and the other parameters are dependent variables. We elaborate on the calculation steps and provide some notes. Then we use our method to analyze gas-bearing sandstone samples. In the calculations, density and P- and S-velocities are input data and we calculate eleven relative parameters for porous fluid, framework, and critical point. In the end, by comparing our results with experiment measurements, we prove the viability of the method. 展开更多
关键词 linear equation numerical calculation gas-bearing media critical point pore fluid and framework elastic parameters
下载PDF
The Order of Approximation by Complex Rational Type Interpolating Operators at Fejer's Points 被引量:1
2
作者 王子玉 朱长青 《Chinese Quarterly Journal of Mathematics》 CSCD 1992年第3期66-70,共5页
In this paper, supose Γ be a boundary of a Jordan domain D and Γ satisfied Альпер condition, the order that rational type interpolating operators at Fejer's points of f(z)∈C(Γ) converge to f(z) in the se... In this paper, supose Γ be a boundary of a Jordan domain D and Γ satisfied Альпер condition, the order that rational type interpolating operators at Fejer's points of f(z)∈C(Γ) converge to f(z) in the sense of uniformly convergence is obtained. 展开更多
关键词 rational type interpolation Fejer's points order of approximation Jordan domain
下载PDF
The Order of Approximation by(0,1,…,q) Hermite-Fejer Interpolation Polynomials at Nearly Fejer's Points
3
作者 沈燮昌 王鸿飞 朱长青 《Chinese Quarterly Journal of Mathematics》 CSCD 1992年第1期81-88,共8页
Suppose that the outer mapping function of domain D has its second continuous derivatives. In this paper, the order proximation by (0,1,…,q) Hermite-Fejer interpolating polynomials at nearly Fejer's points of fun... Suppose that the outer mapping function of domain D has its second continuous derivatives. In this paper, the order proximation by (0,1,…,q) Hermite-Fejer interpolating polynomials at nearly Fejer's points of function of class A(D) are presented. Moreover in general the order of approximation is sharp. 展开更多
关键词 Hermite-Fejer interpolating Fejer's point
下载PDF
An interpolating reproducing kernel particle method for two-dimensional scatter points 被引量:2
4
作者 秦义校 刘营营 +1 位作者 李中华 杨明 《Chinese Physics B》 SCIE EI CAS CSCD 2014年第7期238-241,共4页
An interpolating reproducing kernel particle method for two-dimensional (2D) scatter points is introduced. It elim- inates the dependency of gridding in numerical calculations. The interpolating shape function in th... An interpolating reproducing kernel particle method for two-dimensional (2D) scatter points is introduced. It elim- inates the dependency of gridding in numerical calculations. The interpolating shape function in the interpolating repro- ducing kernel particle method satisfies the property of the Kronecker delta function. This method offers a mathematics basis for recognition technology and simulation analysis, which can be expressed as simultaneous differential equations in science or project problems. Mathematical examples are given to show the validity of the interpolating reproducing kernel particle method. 展开更多
关键词 interpolating reproducing kernel particle method point interpolating characteristic scatter points
下载PDF
COMPLEX RATIONAL TYPE INTERPOLATION AT DISTURBED FEJR POINTS ON A CURVE 被引量:1
5
作者 涂天亮 莫烔 《Acta Mathematica Scientia》 SCIE CSCD 2012年第2期441-454,共14页
The results of accurate order of uniform approximation and simultaneous approximation in the early work "Jackson Type Theorems on Complex Curves" are improved from Fejer points to disturbed Fejer points in this arti... The results of accurate order of uniform approximation and simultaneous approximation in the early work "Jackson Type Theorems on Complex Curves" are improved from Fejer points to disturbed Fejer points in this article. Furthermore, the stability of convergence of Tn,∈(f,z) with disturbed sample values f(z^*) + Sk are also proved in this article. 展开更多
关键词 Disturbed Fejer points complex interpolation order of approximation closed Jordan curve uniform convergence
下载PDF
Effective point kinetic parameters calculation in Tehran research reactor using deterministic and probabilistic methods 被引量:1
6
作者 M.Kheradmand Saadi A.Abbaspour 《Nuclear Science and Techniques》 SCIE CAS CSCD 2017年第12期182-192,共11页
The exact calculation of point kinetic parameters is very important in nuclear reactor safety assessment, and most sophisticated safety codes such as RELAP5, PARCS,DYN3D, and PARET are using these parameters in their ... The exact calculation of point kinetic parameters is very important in nuclear reactor safety assessment, and most sophisticated safety codes such as RELAP5, PARCS,DYN3D, and PARET are using these parameters in their dynamic models. These parameters include effective delayed neutron fractions as well as mean generation time.These parameters are adjoint-weighted, and adjoint flux is employed as a weighting function in their evaluation.Adjoint flux calculation is an easy task for most of deterministic codes, but its evaluation is cumbersome for Monte Carlo codes. However, in recent years, some sophisticated techniques have been proposed for Monte Carlo-based point kinetic parameters calculation without any need of adjoint flux. The most straightforward scheme is known as the ‘‘prompt method'' and has been used widely in literature. The main objective of this article is dedicated to point kinetic parameters calculation in Tehran research reactor(TRR) using deterministic as well as probabilistic techniques. WIMS-D5B and CITATION codes have been used in deterministic calculation of forward and adjoint fluxes in the TRR core. On the other hand, the MCNP Monte Carlo code has been employed in the ‘‘prompt method''scheme for effective delayed neutron fraction evaluation.Deterministic results have been cross-checked with probabilistic ones and validated with SAR and experimental data. In comparison with experimental results, the relativedifferences of deterministic as well as probabilistic methods are 7.6 and 3.2%, respectively. These quantities are10.7 and 6.4%, respectively, in comparison with SAR report. 展开更多
关键词 point kinetic parameters TEHRAN research reactor ADJOINT flux Prompt METHOD DETERMINISTIC METHOD Probabilistic METHOD
下载PDF
Interpolating Isogeometric Boundary Node Method and Isogeometric Boundary Element Method Based on Parameter Space 被引量:1
7
作者 Hongyin Yang Jiwei Zhong +2 位作者 Ying Wang Xingquan Chen Xiaoya Bian 《Computer Modeling in Engineering & Sciences》 SCIE EI 2020年第9期807-824,共18页
In this paper,general interpolating isogeometric boundary node method(IIBNM)and isogeometric boundary element method(IBEM)based on parameter space are proposed for 2D elasticity problems.In both methods,the integral c... In this paper,general interpolating isogeometric boundary node method(IIBNM)and isogeometric boundary element method(IBEM)based on parameter space are proposed for 2D elasticity problems.In both methods,the integral cells and elements are defined in parameter space,which can reproduce the geometry exactly at all the stages.In IIBNM,the improved interpolating moving leastsquare method(IIMLS)is applied for field approximation and the shape functions have the delta function property.The Lagrangian basis functions are used for field approximation in IBEM.Thus,the boundary conditions can be imposed directly in both methods.The shape functions are defined in 1D parameter space and no curve length needs to be computed.Besides,most methods for the treatment of the singular integrals in the boundary element method can be applied in IIBNM and IBEM directly.Numerical examples have demonstrated the accuracy of the proposed methods. 展开更多
关键词 interpolating isogeometric boundary node method isogeometric boundary element method parameter space improved interpolating moving least-square method Lagrangian basis functions
下载PDF
PARAMETER IDENTIFICATION PROBLEM OF THE FRACTAL INTERPOLATION FUNCTIONS 被引量:4
8
作者 阮火军 沙震 苏维宜 《Numerical Mathematics A Journal of Chinese Universities(English Series)》 SCIE 2003年第2期205-213,共9页
Parameter identification problem is one of essential problem in order to model effectively experimental data by fractal interpolation function.In this paper,we first present an example to explain a relationship betwee... Parameter identification problem is one of essential problem in order to model effectively experimental data by fractal interpolation function.In this paper,we first present an example to explain a relationship between iteration procedure and fractal function.Then we discuss conditions that vertical scaling factors must obey in one typical case. 展开更多
关键词 分形插值函数 参数鉴定 吸引子 垂直定标因数
下载PDF
Three-dimensional(3D)parametric measurements of individual gravels in the Gobi region using point cloud technique
9
作者 JING Xiangyu HUANG Weiyi KAN Jiangming 《Journal of Arid Land》 SCIE CSCD 2024年第4期500-517,共18页
Gobi spans a large area of China,surpassing the combined expanse of mobile dunes and semi-fixed dunes.Its presence significantly influences the movement of sand and dust.However,the complex origins and diverse materia... Gobi spans a large area of China,surpassing the combined expanse of mobile dunes and semi-fixed dunes.Its presence significantly influences the movement of sand and dust.However,the complex origins and diverse materials constituting the Gobi result in notable differences in saltation processes across various Gobi surfaces.It is challenging to describe these processes according to a uniform morphology.Therefore,it becomes imperative to articulate surface characteristics through parameters such as the three-dimensional(3D)size and shape of gravel.Collecting morphology information for Gobi gravels is essential for studying its genesis and sand saltation.To enhance the efficiency and information yield of gravel parameter measurements,this study conducted field experiments in the Gobi region across Dunhuang City,Guazhou County,and Yumen City(administrated by Jiuquan City),Gansu Province,China in March 2023.A research framework and methodology for measuring 3D parameters of gravel using point cloud were developed,alongside improved calculation formulas for 3D parameters including gravel grain size,volume,flatness,roundness,sphericity,and equivalent grain size.Leveraging multi-view geometry technology for 3D reconstruction allowed for establishing an optimal data acquisition scheme characterized by high point cloud reconstruction efficiency and clear quality.Additionally,the proposed methodology incorporated point cloud clustering,segmentation,and filtering techniques to isolate individual gravel point clouds.Advanced point cloud algorithms,including the Oriented Bounding Box(OBB),point cloud slicing method,and point cloud triangulation,were then deployed to calculate the 3D parameters of individual gravels.These systematic processes allow precise and detailed characterization of individual gravels.For gravel grain size and volume,the correlation coefficients between point cloud and manual measurements all exceeded 0.9000,confirming the feasibility of the proposed methodology for measuring 3D parameters of individual gravels.The proposed workflow yields accurate calculations of relevant parameters for Gobi gravels,providing essential data support for subsequent studies on Gobi environments. 展开更多
关键词 Gobi gravels three-dimensional(3D)parameters point cloud 3D reconstruction Random Sample Consensus(RANSAC)algorithm Density-Based Spatial Clustering of Applications with Noise(DBSCAN)
下载PDF
Saddle Point Solution for System with Parameter Uncertainties
10
作者 Abiola Bankole T. C. Obiwuru 《Advances in Pure Mathematics》 2013年第8期685-688,共4页
In this paper, we consider dynamical system, in the presence of parameter uncertainties. We apply max-min principles to determine the saddle point solution for the class of differential game arising from the associate... In this paper, we consider dynamical system, in the presence of parameter uncertainties. We apply max-min principles to determine the saddle point solution for the class of differential game arising from the associated dynamical system. We also provide sufficient condition for the existence of this saddle point. 展开更多
关键词 parameter Uncertainty MIN-MAX PRINCIPLES SADDLE point DIFFERENTIAL GAMES
下载PDF
C^1 C^2INTERPOLATION OF SCATTERED DATA POINTS 
11
作者 WANG JIAYE AND ZHANG CAIMING(Department of Computer Science,Shandong University Jinan 250100) 《Applied Mathematics(A Journal of Chinese Universities)》 SCIE CSCD 1994年第1期1-9,共9页
In this paper an error in[4]is pointed out and a method for constructingsurface interpolating scattered data points is presented.The main feature of the methodin this paper is that the surface so constructed is polyno... In this paper an error in[4]is pointed out and a method for constructingsurface interpolating scattered data points is presented.The main feature of the methodin this paper is that the surface so constructed is polynomial,which makes the construction simple and the calculation easy. 展开更多
关键词 interpolATION Scattered Data points TRIANGLE POLYNOMIAL BarycentricCoordinate.
下载PDF
Convexity-preserving interpolation of trigonometric polynomial curves with a shape parameter
12
作者 PAN Yong-juan WANG Guo-jin 《Journal of Zhejiang University-Science A(Applied Physics & Engineering)》 SCIE EI CAS CSCD 2007年第8期1199-1209,共11页
In computer aided geometric design(CAGD) ,it is often needed to produce a convexity-preserving interpolating curve according to the given planar data points. However,most existing pertinent methods cannot generate con... In computer aided geometric design(CAGD) ,it is often needed to produce a convexity-preserving interpolating curve according to the given planar data points. However,most existing pertinent methods cannot generate convexity-preserving in-terpolating transcendental curves;even constructing convexity-preserving interpolating polynomial curves,it is required to solve a system of equations or recur to a complicated iterative process. The method developed in this paper overcomes the above draw-backs. The basic idea is:first to construct a kind of trigonometric polynomial curves with a shape parameter,and interpolating trigonometric polynomial parametric curves with C2(or G1) continuity can be automatically generated without having to solve any system of equations or do any iterative computation. Then,the convexity of the constructed curves can be guaranteed by the appropriate value of the shape parameter. Performing the method is easy and fast,and the curvature distribution of the resulting interpolating curves is always well-proportioned. Several numerical examples are shown to substantiate that our algorithm is not only correct but also usable. 展开更多
关键词 Computer aided geometric design (CAGD) α-trigonometric polynomial curves interpolATION Convexity-preserving Shape parameter
下载PDF
Along Road Height Interpolation Based on Discrete Elevation Points
13
作者 Tran Thong Nhat 《Journal of Earth Science and Engineering》 2012年第11期691-695,共5页
The simulating exactly compared with realty of ground surface to run a model is more and more highly required. In the real, terrain of the earth surface is always complicated by the natural and human made ground objec... The simulating exactly compared with realty of ground surface to run a model is more and more highly required. In the real, terrain of the earth surface is always complicated by the natural and human made ground objects. Because of limitation of collecting and storing technologies in the past time, data are usually not detailed so that the data can not be full for the simulation. Besides computing power and simulation increase more day by day, the increasing requirements more detailed of topography surface simulation is a demand. In simulated flooding phenomenon or phenomena related to energy and momentum of water flow, the linear objects of ground surface such as roads, dikes, dams, etc. need to have their vertical dimension along continuously. However, these datas have often no height information alternately, there are only discrete elevation points that are extracted from topographic maps. Consequently, the demand of a suitable method for linear objects height interpolation is necessary. This paper aims to provide a method and evaluate its accuracy to meet this requirement. 展开更多
关键词 Height interpolation GIS road interpolation discrete elevation points.
下载PDF
A MAXIMUM PRINCIPLE FOR DISTRIBUTED PARAMETER SYSTEMS WITH MIXED PHASE-CONTROL CONSTRAINTS ANDENDPOINT CONSTRAINTS
14
作者 张平健 《Acta Mathematica Scientia》 SCIE CSCD 1997年第2期151-158,共8页
This paper considered the optimal control problem for distributed parameter systems with mixed phase-control constraints and end-point constraints. Pontryagin's maximum principle for optimal control are derived vi... This paper considered the optimal control problem for distributed parameter systems with mixed phase-control constraints and end-point constraints. Pontryagin's maximum principle for optimal control are derived via Duboviskij-Milujin theorem. 展开更多
关键词 maximum principle distributed parameter systems end-point constraints Duboviskij-Milujin theorem additive vector measure
下载PDF
A Method Determining Vertical Scaling Parameters of Fractal Interpolation
15
作者 Dai Fang Xu Xiaoge 1. Xi’an University of Technology, Xi’an 710048, China 2. Beijing Information Technology Institute, Beijing 100101, China 《Computer Aided Drafting,Design and Manufacturing》 2002年第1期37-41,共5页
A method determining vertical scaling parameters of fractal interpolation is given in this paper. By computer experiments, it is clear that this method is very effective.
关键词 FRACTAL fractal interpolation scaling parameter
全文增补中
Ground-based/UAV-LiDAR data fusion for quantitative structure modeling and tree parameter retrieval in subtropical planted forest 被引量:5
16
作者 Reda Fekry Wei Yao +1 位作者 Lin Cao Xin Shen 《Forest Ecosystems》 SCIE CSCD 2022年第5期674-691,共18页
Light detection and ranging(LiDAR)has contributed immensely to forest mapping and 3D tree modelling.From the perspective of data acquisition,the integration of LiDAR data from different platforms would enrich forest i... Light detection and ranging(LiDAR)has contributed immensely to forest mapping and 3D tree modelling.From the perspective of data acquisition,the integration of LiDAR data from different platforms would enrich forest information at the tree and plot levels.This research develops a general framework to integrate ground-based and UAV-LiDAR(ULS)data to better estimate tree parameters based on quantitative structure modelling(QSM).This is accomplished in three sequential steps.First,the ground-based/ULS LiDAR data were co-registered based on the local density peaks of the clustered canopy.Next,redundancy and noise were removed for the ground-based/ULS LiDAR data fusion.Finally,tree modeling and biophysical parameter retrieval were based on QSM.Experiments were performed for Backpack/Handheld/UAV-based multi-platform mobile LiDAR data of a subtropical forest,including poplar and dawn redwood species.Generally,ground-based/ULS LiDAR data fusion outperforms ground-based LiDAR with respect to tree parameter estimation compared to field data.The fusion-derived tree height,tree volume,and crown volume significantly improved by up to 9.01%,5.28%,and 18.61%,respectively,in terms of rRMSE.By contrast,the diameter at breast height(DBH)is the parameter that has the least benefits from fusion,and rRMSE remains approximately the same,because stems are already well sampled from ground data.Additionally,particularly for dense forests,the fusion-derived tree parameters were improved compared to those derived from ground-based LiDAR.Ground-based LiDAR can potentially be used to estimate tree parameters in low-stand-density forests,whereby the improvement owing to fusion is not significant. 展开更多
关键词 Ground/aerial view mobile LiDAR point cloud CO-REGISTRATION FUSION QSM Tree parameter retrieval
下载PDF
MODIFIABLE QUARTIC AND QUINTIC CURVES WITH SHAPE-PARAMETERS 被引量:2
17
作者 邬弘毅 朱功勤 《Numerical Mathematics A Journal of Chinese Universities(English Series)》 SCIE 2001年第1期9-19,共11页
This paper presents a method for creating modificable quartic and quintic curves with shape parameters. The curves can achieve C 2 even C 3 continuity and unify both interpolation and approximation to the control poin... This paper presents a method for creating modificable quartic and quintic curves with shape parameters. The curves can achieve C 2 even C 3 continuity and unify both interpolation and approximation to the control points without solving a system of equations or inserting additional control points. They have the local properties like the cubic B spline. Besides, the quintic curve would be able globally to tend the control polygon. 展开更多
关键词 modifiable quartic and quintic curves shape parameters C 2 or C 3 continuity interpolation and approximation.
下载PDF
Non-parametric machine learning methods for interpolation of spatially varying non-stationary and non-Gaussian geotechnical properties 被引量:3
18
作者 Chao Shi Yu Wang 《Geoscience Frontiers》 SCIE CAS CSCD 2021年第1期339-350,共12页
Spatial interpolation has been frequently encountered in earth sciences and engineering.A reasonable appraisal of subsurface heterogeneity plays a significant role in planning,risk assessment and decision making for g... Spatial interpolation has been frequently encountered in earth sciences and engineering.A reasonable appraisal of subsurface heterogeneity plays a significant role in planning,risk assessment and decision making for geotechnical practice.Geostatistics is commonly used to interpolate spatially varying properties at un-sampled locations from scatter measurements.However,successful application of classic geostatistical models requires prior characterization of spatial auto-correlation structures,which poses a great challenge for unexperienced engineers,particularly when only limited measurements are available.Data-driven machine learning methods,such as radial basis function network(RBFN),require minimal human intervention and provide effective alternatives for spatial interpolation of non-stationary and non-Gaussian data,particularly when measurements are sparse.Conventional RBFN,however,is direction independent(i.e.isotropic)and cannot quantify prediction uncertainty in spatial interpolation.In this study,an ensemble RBFN method is proposed that not only allows geotechnical anisotropy to be properly incorporated,but also quantifies uncertainty in spatial interpolation.The proposed method is illustrated using numerical examples of cone penetration test(CPT)data,which involve interpolation of a 2D CPT cross-section from limited continuous 1D CPT soundings in the vertical direction.In addition,a comparative study is performed to benchmark the proposed ensemble RBFN with two other non-parametric data-driven approaches,namely,Multiple Point Statistics(MPS)and Bayesian Compressive Sensing(BCS).The results reveal that the proposed ensemble RBFN provides a better estimation of spatial patterns and associated prediction uncertainty at un-sampled locations when a reasonable amount of data is available as input.Moreover,the prediction accuracy of all the three methods improves as the number of measurements increases,and vice versa.It is also found that BCS prediction is less sensitive to the number of measurement data and outperforms RBFN and MPS when only limited point observations are available. 展开更多
关键词 Spatial interpolation Multiple point statistics Bayesian compressive sampling Compressive sensing Sparse measurement
下载PDF
Factors affecting accuracy of radial point interpolation meshfree method for 3-D solid mechanics
19
作者 彭翀 袁会娜 +1 位作者 张丙印 张琰 《Journal of Central South University》 SCIE EI CAS 2013年第11期3229-3246,共18页
Recently,the radial point interpolation meshfree method has gained popularity owing to its advantages in large deformation and discontinuity problems,however,the accuracy of this method depends on many factors and the... Recently,the radial point interpolation meshfree method has gained popularity owing to its advantages in large deformation and discontinuity problems,however,the accuracy of this method depends on many factors and their influences are not fully investigated yet.In this work,three main factors,i.e.,the shape parameters,the influence domain size,and the nodal distribution,on the accuracy of the radial point interpolation method(RPIM)are systematically studied and conclusive results are obtained.First,the effect of shape parameters(R,q)of the multi-quadric basis function on the accuracy of RPIM is examined via global search.A new interpolation error index,closely related to the accuracy of RPIM,is proposed.The distribution of various error indexes on the R q plane shows that shape parameters q[1.2,1.8]and R[0,1.5]can give good results for general 3-D analysis.This recommended range of shape parameters is examined by multiple benchmark examples in 3D solid mechanics.Second,through numerical experiments,an average of 30 40 nodes in the influence domain of a Gauss point is recommended for 3-D solid mechanics.Third,it is observed that the distribution of nodes has significant effect on the accuracy of RPIM although it has little effect on the accuracy of interpolation.Nodal distributions with better uniformity give better results.Furthermore,how the influence domain size and nodal distribution affect the selection of shape parameters and how the nodal distribution affects the choice of influence domain size are also discussed. 展开更多
关键词 meshfree method radial point interpolation method shape parameter influence domain size nodal distribution
下载PDF
An improved local radial point interpolation method for transient heat conduction analysis
20
作者 王峰 林皋 +1 位作者 郑保敬 胡志强 《Chinese Physics B》 SCIE EI CAS CSCD 2013年第6期127-134,共8页
The smoothing thin plate spline (STPS) interpolation using the penalty function method according to the optimization theory is presented to deal with transient heat conduction problems. The smooth conditions of the ... The smoothing thin plate spline (STPS) interpolation using the penalty function method according to the optimization theory is presented to deal with transient heat conduction problems. The smooth conditions of the shape functions and derivatives can be satisfied so that the distortions hardly occur. Local weak forms are developed using the weighted residual method locally from the partial differential equations of the transient heat conduction. Here the Heaviside step function is used as the test function in each sub-domain to avoid the need for a domain integral. Essential boundary conditions can be implemented like the finite element method (FEM) as the shape functions possess the Kronecker delta property. The traditional two-point difference method is selected for the time discretization scheme. Three selected numerical examples are presented in this paper to demonstrate the availability and accuracy of the present approach comparing with the traditional thin plate spline (TPS) radial basis functions. 展开更多
关键词 thin plate splines transient heat conduction penalty function method local radial point interpolation method
下载PDF
上一页 1 2 151 下一页 到第
使用帮助 返回顶部