期刊文献+
共找到197篇文章
< 1 2 10 >
每页显示 20 50 100
A class of quasilinear equations with-1 powers
1
作者 ZHANG Heng SUN Yijing 《中国科学院大学学报(中英文)》 北大核心 2025年第1期13-19,共7页
This paper deals with quasilinear elliptic equations of singular growth like-Δu-uΔ(u^(2))=a(x)u^(-1).We establish the existence of positive solutions for general a(x)∈L^(p)(Ω),p>2,whereΩis a bounded domain inℝ... This paper deals with quasilinear elliptic equations of singular growth like-Δu-uΔ(u^(2))=a(x)u^(-1).We establish the existence of positive solutions for general a(x)∈L^(p)(Ω),p>2,whereΩis a bounded domain inℝ^(N)with N≥1. 展开更多
关键词 quasilinear singular equation -1 power elliptic equation
在线阅读 下载PDF
Normalized Solutions of Nonlinear Choquard Equations with Nonconstant Potential
2
作者 LI Nan XU Liping 《应用数学》 北大核心 2025年第1期14-29,共16页
In this paper,we mainly focus on a type of nonlinear Choquard equations with nonconstant potential.Under appropriate hypotheses on potential function and nonlinear terms,we prove that the above Choquard equation with ... In this paper,we mainly focus on a type of nonlinear Choquard equations with nonconstant potential.Under appropriate hypotheses on potential function and nonlinear terms,we prove that the above Choquard equation with prescribed 2-norm has some normalized solutions by introducing variational methods. 展开更多
关键词 Nonlinear Choquard equation Potential function Variational method Normalized solution
在线阅读 下载PDF
Two Second-Order Ecient Numerical Schemes for the Boussinesq Equations
3
作者 LIU Fang WANG Danxia ZHANG Jianwen 《应用数学》 北大核心 2025年第1期114-129,共16页
In this paper,we construct two fully decoupled,second-order semi-discrete numerical schemes for the Boussinesq equations based on the scalar auxiliary variable(SAV)approach.By introducing a scalar auxiliary variable,t... In this paper,we construct two fully decoupled,second-order semi-discrete numerical schemes for the Boussinesq equations based on the scalar auxiliary variable(SAV)approach.By introducing a scalar auxiliary variable,the original Boussinesq system is transformed into an equivalent one.Then we discretize it using the second-order backward di erentiation formula(BDF2)and Crank-Nicolson(CN)to obtain two second-order time-advanced schemes.In both numerical schemes,a pressure-correction method is employed to decouple the velocity and pressure.These two schemes possess the desired property that they can be fully decoupled with satisfying unconditional stability.We rigorously prove both the unconditional stability and unique solvability of the discrete schemes.Furthermore,we provide detailed implementations of the decoupling procedures.Finally,various 2D numerical simulations are performed to verify the accuracy and energy stability of the proposed schemes. 展开更多
关键词 Scalar auxiliary variable approach Pressure-correction method Fully decoupled Unconditional stability Boussinesq equations
在线阅读 下载PDF
Two Efficient Numerical Algorithms for the Thermally Coupled Incompressible MHD Equations
4
作者 LIAO Mingliang WANG Danxia 《应用数学》 北大核心 2025年第4期932-951,共20页
In this work,we construct two efficient fully decoupled,linear,unconditionally stable numerical algorithms for the thermally coupled incompressible magnetohydrodynamic equations.Firstly,in order to obtain the desired ... In this work,we construct two efficient fully decoupled,linear,unconditionally stable numerical algorithms for the thermally coupled incompressible magnetohydrodynamic equations.Firstly,in order to obtain the desired algorithm,we introduce a scalar auxiliary variable(SAV)to get a new equivalent system.Secondly,by combining the pressure-correction method and the explicit-implicit method,we perform semi-discrete numerical algorithms of first and second order,respectively.Then,we prove that the obtained algorithms follow an unconditionally stable law in energy,and we provide a detailed implementation process,which we only need to solve a series of linear differential equations with constant coefficients at each time step.More importantly,with some powerful analysis,we give the order of convergence of the errors.Finally,to illustrate theoretical results,some numerical experiments are given. 展开更多
关键词 Magnetohydrodynamic equation Fully decoupled Unconditionally stability Scalar auxiliary variable Error analysis
在线阅读 下载PDF
Investigating Solutions in Nonlinear Evolution Equations:A Focus on Local Existence in Mixed Types
5
作者 NAFFISA Toureche Trouba FAN Long ABDELGHANI Dahou 《应用数学》 北大核心 2025年第3期691-702,共12页
With the urgent need to resolve complex behaviors in nonlinear evolution equations,this study makes a contribution by establishing the local existence of solutions for Cauchy problems associated with equations of mixe... With the urgent need to resolve complex behaviors in nonlinear evolution equations,this study makes a contribution by establishing the local existence of solutions for Cauchy problems associated with equations of mixed types.Our primary contribution is the establishment of solution existence,illuminating the dynamics of these complex equations.To tackle this challenging problem,we construct an approximate solution sequence and apply the contraction mapping principle to rigorously prove local solution existence.Our results significantly advance the understanding of nonlinear evolution equations of mixed types.Furthermore,they provide a versatile,powerful approach for tackling analogous challenges across physics,engineering,and applied mathematics,making this work a valuable reference for researchers in these fields. 展开更多
关键词 Nonlinear evolution equation Contraction mapping principle Sobolev space Dissipative system
在线阅读 下载PDF
A New Class of Efficient Schemes for the Cahn-Hilliard-Navier-Stokes Equations
6
作者 WANG Lijing WANG Danxia ZHANG Jianwen 《应用数学》 北大核心 2025年第3期607-624,共18页
In this paper,we construct a new class of efficient and high-order schemes for the Cahn-Hilliard-Navier-Stokes equations with periodic boundary conditions.These schemes are based on two types of scalar auxiliary varia... In this paper,we construct a new class of efficient and high-order schemes for the Cahn-Hilliard-Navier-Stokes equations with periodic boundary conditions.These schemes are based on two types of scalar auxiliary variable approaches.By using a new pressure correction method,the accuracy of the pressure has been greatly improved.Furthermore,one only needs to solve a series of fully decoupled linear equations with constant coefficients at each time step.In addition,we prove the unconditional energy stability of the schemes,rigorously.Finally,plenty of numerical simulations are carried out to verify the convergence rates,stability,and effectiveness of the proposed schemes numerically. 展开更多
关键词 Cahn-Hilliard-Navier-Stokes equation Scalar auxiliary variable Pressurecorrection Unconditional energy stability
在线阅读 下载PDF
Normalized Positive Ground State Solutions for Nonhomogeneous Kirchhoff Equations
7
作者 ZHANG Xiaocang XU Liping 《应用数学》 北大核心 2025年第3期711-720,共10页
This article studies a class of nonlinear Kirchhoff equations with exponential critical growth,trapping potential,and perturbation.Under appropriate assumptions about f and h,the article obtained the existence of norm... This article studies a class of nonlinear Kirchhoff equations with exponential critical growth,trapping potential,and perturbation.Under appropriate assumptions about f and h,the article obtained the existence of normalized positive solutions for this equation via the Trudinger-Moser inequality and variational methods.Moreover,these solutions are also ground state solutions.Additionally,the article also characterized the asymptotic behavior of solutions.The results of this article expand the research of relevant literature. 展开更多
关键词 Normalized positive ground state solution Nonhomogeneous Kirchhoff equation Variational method Exponential critical growth Trapping potential
在线阅读 下载PDF
Smoluchowski-Kramers Approximation for Stochastic Differential Equations under Discretization
8
作者 Li Ge 《应用概率统计》 北大核心 2025年第4期622-635,共14页
This paper studies the Smoluchowski–Kramers approximation for a discrete-time dynamical system modeled as the motion of a particle in a force field.We show that the approximation holds for the drift-implicit Euler–M... This paper studies the Smoluchowski–Kramers approximation for a discrete-time dynamical system modeled as the motion of a particle in a force field.We show that the approximation holds for the drift-implicit Euler–Maruyama discretization and derive its convergence rate.In particular,the solution of the discretized system converges to the solution of the first-order limit equation in the mean-square sense,and this convergence is independent of the order in which the mass parameterμand the step size h tend to zero. 展开更多
关键词 stochastic differential equations Smoluchowski-Kramers approximation driftimplicit Euler-Maruyama scheme convergence rate
在线阅读 下载PDF
On some multiple solutions to a p(x)-Laplace equation with supercritical growth
9
作者 Lin Zhao 《中国科学技术大学学报》 北大核心 2025年第5期15-20,14,I0001,共8页
We consider the multiplicity of solutions to a p(x)-Laplacian problem involving supercritical Sobolev growth via Ricceri’s principle.By means of truncation combined with De Giorgi iteration,we can extend the results ... We consider the multiplicity of solutions to a p(x)-Laplacian problem involving supercritical Sobolev growth via Ricceri’s principle.By means of truncation combined with De Giorgi iteration,we can extend the results of subcritical and critical growth to supercritical growth and obtain at least three solutions to the p(x)-Laplacian problem. 展开更多
关键词 degenerate elliptic equations a priori estimates truncated technique supercritical Sobolev exponent
在线阅读 下载PDF
Evaluation of detonation performance of explosives ICM-101,ONC,and TNAZ based on improved VHL equation of state
10
作者 Yong Han Qin Liu +2 位作者 Yingliang Duan Yaqi Zhao Xinping Long 《Defence Technology(防务技术)》 2025年第2期83-97,共15页
Detonation performance is crucial for evaluating the power of high explosives(HEs),and the equation of state(EOS)that accurately describes the high-temperature,high-pressure,and high-temperature,medium-pressure states... Detonation performance is crucial for evaluating the power of high explosives(HEs),and the equation of state(EOS)that accurately describes the high-temperature,high-pressure,and high-temperature,medium-pressure states of detonation products is key to assessing the damage efficiency of these energetic materials.This article examines the limitations of the VLW EOS in representing the thermodynamic states of explosive detonation gas products under high-temperature and medium-to high-pressure conditions.A new gas EOS for detonation products,called VHL(Virial-Han-Long),is proposed.The accuracy of VHL in describing gas states under high-temperature and medium-to high-pressure conditions is verified,and its performance in evaluating explosive detonation and working capabilities is explored.The results demonstrate that VHL exhibits high precision in calculating detonation performance.Subsequently,the detonation performance of three new HEs(ICM-101,ONC,and TNAZ)was calculated and compared to traditional HEs(TATB,CL-20,and HMX).The results indicate that ONC has superior detonation performance compared to the other explosives,while ICM-101 shows a detonation velocity similar to CL-20 but with slightly lower detonation pressure.The detonation characteristics of TNAZ are comparable to those of the standard HE HMX.From the perspective of products,considering the comprehensive work performance(mechanical work and detonation heat),both ONC and ICM-101demonstrate relatively superior performance. 展开更多
关键词 equation of state Detonation performance Working capability THERMODYNAMICS High explosive
在线阅读 下载PDF
Development and experimental validation of the dynamic constitutive model and equation of state for Mo-10Cu alloy
11
作者 Bihui Hong Wenbin Li +3 位作者 Yiming Li Weihang Li Peng Chen Yalong Wang 《Defence Technology(防务技术)》 2025年第9期145-158,共14页
This study systematically investigates the mechanical response characteristics of Mo-10Cu pseudo-alloy under various conditions,including temperatures ranging from 298 K to 550 K,strain rates from1×10^(-2)s^(-1)t... This study systematically investigates the mechanical response characteristics of Mo-10Cu pseudo-alloy under various conditions,including temperatures ranging from 298 K to 550 K,strain rates from1×10^(-2)s^(-1)to 5.2×10^(3)s^(-1),and dynamic impact loads from 134 m/s to 837 m/s.The investigation is conducted using a combination of multi-method crossover experiment and numerical simulations,with accuracy validated through X-ray testing and static penetration test.Using a universal testing machine,Split-Hopkinson Pressure Bar(SHPB)system,and a light-gas gun,the dynamic constitutive behavior and shock adiabatic curves of the alloy under complex loading conditions are revealed.Experimental results demonstrate that the flow stress evolution of Mo-10Cu alloy exhibits significant strain hardening,and strain-rate strengthening.Based on these observations,a Johnson-Cook(J-C)constitutive model has been developed to describe the material's dynamic behavior.Through free-surface particle velocity measurements,the shock adiabatic relationship was obtained,and a Gruneisen equation of state was established.X-ray experimental results confirm that the Mo-10Cu liner can generate well-formed,cohesive jets.The penetration test results show that the maximum penetration depth can reach243.10 mm.The maximum error between the numerical simulation and the X-ray test is less than 7.70%,and the error with the penetration test is 4.73%,which confirms the accuracy of the constitutive parameters and the state equation.In conclusion,the proposed J-C model and Gruneisen equation effectively predict the dynamic response and jet formation characteristics of Mo-10Cu alloy under extreme loads.This work provides both theoretical support and experimental data for material design and performance optimization in shaped charge applications. 展开更多
关键词 Shaped charge Mo-10Cu alloy Johnson-cook Gruneisen equation of state X-RAY PENETRATION
在线阅读 下载PDF
Numerical Solution of Constrained Mechanical System Motions Equations and Inverse Problems of Dynamics 被引量:2
12
作者 R.G. Muharliamov (Russian Peoples’ Friendship University, 117198, Moscow, Mikluho Maklaya,6,Russia.) 《应用数学》 CSCD 北大核心 2001年第2期103-119,共17页
In this paper the method of design of kinematical and dynamical equations of mechanical systems, applied to numerical ealization, is proposed. The corresponding difference equations, which are obtained, give a guarant... In this paper the method of design of kinematical and dynamical equations of mechanical systems, applied to numerical ealization, is proposed. The corresponding difference equations, which are obtained, give a guarantee of computations with a given precision. The equations of programmed constraints and those of constraint perturbations are defined. The stability of the programmed manifold for numerical solutions of the kinematical and dynamical equations is obtained by corresponding construction of the constraint perturbation equations. The dynamical equations of system with programmed constraints are set up in the form of Lagrange’s equations in generalized coordinates. Certain inverse problems of rigid body dynamics are examined. 展开更多
关键词 Kinematies Dynamical equations CONSTRAINTS Lagrange’s equations Rigid body Numerical solution Differential algebraic equations
在线阅读 下载PDF
Contributions to Hom-Schunck optical flow equations-part I: Stability and rate of convergence of classical algorithm 被引量:3
13
作者 DONG Guo-hua AN Xiang-jing FANG Yu-qiang HU De-wen 《Journal of Central South University》 SCIE EI CAS 2013年第7期1909-1918,共10页
Globally exponential stability (which implies convergence and uniqueness) of their classical iterative algorithm is established using methods of heat equations and energy integral after embedding the discrete iterat... Globally exponential stability (which implies convergence and uniqueness) of their classical iterative algorithm is established using methods of heat equations and energy integral after embedding the discrete iteration into a continuous flow. The stability condition depends explicitly on smoothness of the image sequence, size of image domain, value of the regularization parameter, and finally discretization step. Specifically, as the discretization step approaches to zero, stability holds unconditionally. The analysis also clarifies relations among the iterative algorithm, the original variation formulation and the PDE system. The proper regularity of solution and natural images is briefly surveyed and discussed. Experimental results validate the theoretical claims both on convergence and exponential stability. 展开更多
关键词 optical flow Hom-Schunck equations globally exponential stability convergence convergence rate heat equations energy integral and estimate Gronwall inequality natural images REGULARITY
在线阅读 下载PDF
Boussinesq equations的新的精确解 被引量:1
14
作者 李伟 《重庆理工大学学报(自然科学)》 CAS 2015年第6期151-154,共4页
微分方程包含线性和非线性微分方程。微分方程研究的主体是非线性微分方程,特别是非线性偏微分方程。很多意义重大的自然科学和工程技术问题都可归结为非线性偏微分方程的研究。另外,随着研究的深入,有些原来可用线性偏微分方程近似处... 微分方程包含线性和非线性微分方程。微分方程研究的主体是非线性微分方程,特别是非线性偏微分方程。很多意义重大的自然科学和工程技术问题都可归结为非线性偏微分方程的研究。另外,随着研究的深入,有些原来可用线性偏微分方程近似处理的问题,也必须考虑非线性的影响。从传统的观点来看,求偏微分方程的精确解是十分困难的,但经过几十年的研究和探索,人们已经找到了一些构造精确解的方法。借助于行波变换法,改进的双曲函数法,齐次平衡法获法和拟解的方法,获得Boussinesq equations的新的精确解。这种方法可以解决一系列的偏微分方程。 展开更多
关键词 行波变换 精确解 BOUSSINESQ equationS
在线阅读 下载PDF
Solution of the HJI equations for nonlinear H_∞ control design by state-dependent Riccati equations approach 被引量:1
15
作者 Xueyan Zhao Feiqi Deng 《Journal of Systems Engineering and Electronics》 SCIE EI CSCD 2011年第4期654-660,共7页
The relationship between the technique by state- dependent Riccati equations (SDRE) and Hamilton-Jacobi-lsaacs (HJI) equations for nonlinear H∞ control design is investigated. By establishing the Lyapunov matrix ... The relationship between the technique by state- dependent Riccati equations (SDRE) and Hamilton-Jacobi-lsaacs (HJI) equations for nonlinear H∞ control design is investigated. By establishing the Lyapunov matrix equations for partial derivates of the solution of the SDREs and introducing symmetry measure for some related matrices, a method is proposed for examining whether the SDRE method admits a global optimal control equiva- lent to that solved by the HJI equation method. Two examples with simulation are given to illustrate the method is effective. 展开更多
关键词 nonlinear system robust control Hamilton-Jacobi-Isaacs (HJI) equation state-dependent Riccati equation (SDRE) global stabilization optimal control.
在线阅读 下载PDF
Improving Solution of Euler Equations by a Gas-Kinetic BGK Method
16
作者 Liu, Ya Gao, Chao Liu, F. 《西北工业大学学报》 EI CAS CSCD 北大核心 2009年第1期1-5,共5页
Aim.The well known JST(Jameson-Schmidt-Turkel) scheme requires the use of a dissipation term.We propose using gas-kinetic BGK(Bhatnagar-Gross-Krook) method,which is based on the more fundamental Boltzmann equation,in ... Aim.The well known JST(Jameson-Schmidt-Turkel) scheme requires the use of a dissipation term.We propose using gas-kinetic BGK(Bhatnagar-Gross-Krook) method,which is based on the more fundamental Boltzmann equation,in order to obviate the use of dissipation term and obtain,we believe,an improved solution.Section 1 deals essentially with three things:(1) as analytical solution of molecular probability density function at the cell interface has been obtained by the Boltzmann equation with BGK model,we can compute the flux term by integrating the density function in the phase space;eqs.(8) and(11) require careful attention;(2) the integrations can be expressed as the moments of Maxwellian distribution with different limits according to the analytical solution;eqs.(9) and(10) require careful attention;(3) the discrete equation by finite volume method can be solved using the time marching method.Computations are performed by the BGK method for the Sod′s shock tube problem and a two-dimensional shock reflection problem.The results are compared with those of the conventional JST scheme in Figs.1 and 2.The BGK method provides better resolution of shock waves and other features of the flow fields. 展开更多
关键词 Euler equations Boltzmann equation gas-kinetic BGK(Bhatnagar-Gross-Krook) scheme JST(Jameson-Schmidt-Turkel) scheme shock tube shock reflection
在线阅读 下载PDF
Generalization of integral inequalities and (c_1,c_1) stability of neutral differential equations with time-varying delays
17
作者 Shuli Guo Lina Han 《Journal of Systems Engineering and Electronics》 SCIE EI CSCD 2017年第2期347-360,共14页
A uniform stability analysis is developed for a type of neutral delays differential equations which depend on more general nonlinear integral inequalities. Many original investigations and results are obtained. Firstl... A uniform stability analysis is developed for a type of neutral delays differential equations which depend on more general nonlinear integral inequalities. Many original investigations and results are obtained. Firstly, generations of two integral nonlinear inequalities are presented, which are very effective in dealing with the complicated integro-differential inequalities whose variable exponents are greater than zero. Compared with existed integral inequalities, those proposed here can be applied to more complicated differential equations, such as time-varying delay neutral differential equations. Secondly, the notions of (ω, Ω) uniform stable and (ω, Ω) uniform asymptotically stable, especially (c1, c1) uniform stable and (c1, c1) uniform asymptotically stable, are presented. Sufficient conditions on about (c1, c1) uniform stable and (c1, c1) uniform asymptotically stable of time-varying delay neutral differential equations are established by the proposed integral inequalities. Finally, a complex numerical example is presented to illustrate the main results effectively. The above work allows to provide further applications on the proposed stability analysis and control system design for different nonlinear systems. © 2017 Beijing Institute of Aerospace Information. 展开更多
关键词 Control system stability Differential equations Nonlinear equations Time delay Time varying control systems
在线阅读 下载PDF
The Research of the Equation Model on the System-level Fault Diagnosis 被引量:5
18
作者 Hengnong Xuan1, Dafang Zhang2, Yue Wu3, Minwei He4 (1,4: Dept. of Computer, Wuyi University, Jiangmen, Guangdong, 529020, China (2: Dept. of Computer, Hunan University, Changsha, 410082, China ) (3: Institute of Computer Science and Engineering, Univers 《湖南大学学报(自然科学版)》 EI CAS CSCD 2000年第S2期148-157,共10页
In this paper we propose an equation model of system-level fault diagnoses, and construct corresponding theory and algorithms. People can turn any PMC model on ex-test into an equivalent equation (or a system of equat... In this paper we propose an equation model of system-level fault diagnoses, and construct corresponding theory and algorithms. People can turn any PMC model on ex-test into an equivalent equation (or a system of equations), and find all consistent fault patterns based on the equation model. We can also find all fault patterns, in which the fault node numbers are less than or equal to t without supposing t-diagnosable. It is not impossible for all graphic models. 展开更多
关键词 PMC MODEL equation MODEL ex-test system-level FAULT diagnosis DIAGNOSTIC algorithm
在线阅读 下载PDF
Vertical-to-horizontal response spectral ratio for offshore ground motions:Analysis and simplified design equation 被引量:15
19
作者 陈宝魁 王东升 +2 位作者 李宏男 孙治国 李超 《Journal of Central South University》 SCIE EI CAS CSCD 2017年第1期203-216,共14页
In order to study the differences in vertical component between onshore and offshore motions,the vertical-to-horizontal peak ground acceleration ratio(V/H PGA ratio) and vertical-to-horizontal response spectral ratio(... In order to study the differences in vertical component between onshore and offshore motions,the vertical-to-horizontal peak ground acceleration ratio(V/H PGA ratio) and vertical-to-horizontal response spectral ratio(V/H) were investigated using the ground motion recordings from the K-NET network and the seafloor earthquake measuring system(SEMS).The results indicate that the vertical component of offshore motions is lower than that of onshore motions.The V/H PGA ratio of acceleration time histories at offshore stations is about 50%of the ratio at onshore stations.The V/H for offshore ground motions is lower than that for onshore motions,especially for periods less than 0.8 s.Furthermore,based on the results in statistical analysis for offshore recordings in the K-NET,the simplified V/H design equations for offshore motions in minor and moderate earthquakes are proposed for seismic analysis of offshore structures. 展开更多
关键词 offshore ground motion vertical component simplified design equation vertical-to-horizontal response spectral ratio (V/H) vertical-to-horizontal peak ground acceleration ratio (V/H PGA ratio) K-NET network seafloor earthquake measuring system (SEMS)
在线阅读 下载PDF
Solving single-frequency phase ambiguity using parameter weights fitting and constrained equation ambiguity resolution methods 被引量:6
20
作者 阳仁贵 欧吉坤 +3 位作者 袁运斌 张克非 闻德保 Ron Grenfell 《Journal of Central South University of Technology》 EI 2006年第1期93-98,共6页
Based on the structural characteristics of the double-differenced normal equation, a new method was proposed to resolve the ambiguity float solution through a selection of parameter weights to construct an appropriate... Based on the structural characteristics of the double-differenced normal equation, a new method was proposed to resolve the ambiguity float solution through a selection of parameter weights to construct an appropriate regularized matrix, and a singular decomposition method was used to generate regularization parameters. Numerical test results suggest that the regularized ambiguity float solution is more stable and reliable than the least-squares float solution. The mean square error matrix of the new method possesses a lower correlation than the variancecovariance matrix of the least-squares estimation. The size of the ambiguity search space is reduced and the search efficiency is improved. The success rate of the integer ambiguity searching process is improved significantly when the ambiguity resolution by using constraint equation method is used to determine the correct ambiguity integervector. The ambiguity resolution by using constraint equation method requires an initial input of the ambiguity float solution candidates which are obtained from the LAMBDA method in the new method. In addition, the observation time required to fix reliable integer ambiguities can he significantly reduced. 展开更多
关键词 global position system ill-conditioned state parameter weight fitting method constraint equation integer ambiguity
在线阅读 下载PDF
上一页 1 2 10 下一页 到第
使用帮助 返回顶部