首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 46 毫秒
1.
In this paper we study the limit as \(\varepsilon \rightarrow 0\) of the singularly perturbed second order equation \(\varepsilon ^2 \ddot{u}_\varepsilon + \nabla _{\!x} V(t,u_\varepsilon (t))=0\), where V(tx) is a potential. We assume that \(u_0(t)\) is one of its equilibrium points such that \(\nabla _{\!x}V(t,u_0(t))=0\) and \(\nabla _{\!x}^2V(t,u_0(t))>0\). We find that, under suitable initial data, the solutions \(u_\varepsilon \) converge uniformly to \(u_0\), by imposing mild hypotheses on V. A counterexample shows that they cannot be weakened.  相似文献   

2.
A spatially two-dimensional sixth order PDE describing the evolution of a growing crystalline surface h(xyt) that undergoes faceting is considered with periodic boundary conditions, as well as its reduced one-dimensional version. These equations are expressed in terms of the slopes \(u_1=h_{x}\) and \(u_2=h_y\) to establish the existence of global, connected attractors for both equations. Since unique solutions are guaranteed for initial conditions in \(\dot{H}^2_{per}\), we consider the solution operator \(S(t): \dot{H}^2_{per} \rightarrow \dot{H}^2_{per}\), to gain our results. We prove the necessary continuity, dissipation and compactness properties.  相似文献   

3.
We consider a Riemann problem for the shallow water system \(u_{t} +\big (v+\textstyle \frac{1}{2}u^{2}\big )_{x}=0\), \(v_{t}+\big (u+uv\big )_{x}=0\) and evaluate all singular solutions of the form \(u(x,t)=l(t)+b(t)H\big (x-\gamma (t)\big )+a(t)\delta \big (x-\gamma (t)\big )\), \(v(x,t)=k(t)+c(t)H\big (x-\gamma (t)\big )\), where \(l,b,a,k,c,\gamma :\mathbb {R}\rightarrow \mathbb {R}\) are \(C^{1}\)-functions of time t, H is the Heaviside function, and \(\delta \) stands for the Dirac measure with support at the origin. A product of distributions, not constructed by approximation processes, is used to define a solution concept, that is a consistent extension of the classical solution concept. Results showing the advantage of this framework are briefly presented in the introduction.  相似文献   

4.
In this paper we focused our study on derived from Anosov diffeomorphisms (DA diffeomorphisms ) of the torus \(\mathbb {T}^3,\) it is, an absolute partially hyperbolic diffeomorphism on \(\mathbb {T}^3\) homotopic to a linear Anosov automorphism of the \(\mathbb {T}^3.\) We can prove that if \(f: \mathbb {T}^3 \rightarrow \mathbb {T}^3 \) is a volume preserving DA diffeomorphism homotopic to a linear Anosov A,  such that the center Lyapunov exponent satisfies \(\lambda ^c_f(x) > \lambda ^c_A > 0,\) with x belongs to a positive volume set, then the center foliation of f is non absolutely continuous. We construct a new open class U of non Anosov and volume preserving DA diffeomorphisms, satisfying the property \(\lambda ^c_f(x) > \lambda ^c_A > 0\) for \(m-\)almost everywhere \(x \in \mathbb {T}^3.\) Particularly for every \(f \in U,\) the center foliation of f is non absolutely continuous.  相似文献   

5.
We consider a family of linearly viscoelastic shells with thickness \(2\varepsilon\), clamped along their entire lateral face, all having the same middle surface \(S=\boldsymbol{\theta}(\bar{\omega})\subset \mathbb{R}^{3}\), where \(\omega\subset\mathbb{R}^{2}\) is a bounded and connected open set with a Lipschitz-continuous boundary \(\gamma\). We make an essential geometrical assumption on the middle surface \(S\), which is satisfied if \(\gamma\) and \(\boldsymbol{\theta}\) are smooth enough and \(S\) is uniformly elliptic. We show that, if the applied body force density is \(O(1)\) with respect to \(\varepsilon\) and surface tractions density is \(O(\varepsilon)\), the solution of the scaled variational problem in curvilinear coordinates, \(\boldsymbol{u}( \varepsilon)\), defined over the fixed domain \(\varOmega=\omega\times (-1,1)\) for each \(t\in[0,T]\), converges to a limit \(\boldsymbol{u}\) with \(u_{\alpha}(\varepsilon)\rightarrow u_{\alpha}\) in \(W^{1,2}(0,T,H ^{1}(\varOmega))\) and \(u_{3}(\varepsilon)\rightarrow u_{3}\) in \(W^{1,2}(0,T,L^{2}(\varOmega))\) as \(\varepsilon\to0\). Moreover, we prove that this limit is independent of the transverse variable. Furthermore, the average \(\bar{\boldsymbol{u}}= \frac{1}{2}\int_{-1}^{1} \boldsymbol{u}dx_{3}\), which belongs to the space \(W^{1,2}(0,T, V_{M}( \omega))\), where
$$V_{M}(\omega)=H^{1}_{0}(\omega)\times H^{1}_{0}(\omega)\times L ^{2}(\omega), $$
satisfies what we have identified as (scaled) two-dimensional equations of a viscoelastic membrane elliptic shell, which includes a long-term memory that takes into account previous deformations. We finally provide convergence results which justify those equations.
  相似文献   

6.
We investigate the influence of a shifting environment on the spreading of an invasive species through a model given by the diffusive logistic equation with a free boundary. When the environment is homogeneous and favourable, this model was first studied in Du and Lin (SIAM J Math Anal 42:377–405, 2010), where a spreading–vanishing dichotomy was established for the long-time dynamics of the species, and when spreading happens, it was shown that the species invades the new territory at some uniquely determined asymptotic speed \(c_0>0\). Here we consider the situation that part of such an environment becomes unfavourable, and the unfavourable range of the environment moves into the favourable part with speed \(c>0\). We prove that when \(c\ge c_0\), the species always dies out in the long-run, but when \(0<c<c_0\), the long-time behavior of the species is determined by a trichotomy described by (a) vanishing, (b) borderline spreading, or (c) spreading. If the initial population is written in the form \(u_0(x)=\sigma \phi (x)\) with \(\phi \) fixed and \(\sigma >0\) a parameter, then there exists \(\sigma _0>0\) such that vanishing happens when \(\sigma \in (0,\sigma _0)\), borderline spreading happens when \(\sigma =\sigma _0\), and spreading happens when \(\sigma >\sigma _0\).  相似文献   

7.
The presence of a finite tangential velocity on a hydrodynamically slipping surface is known to reduce vorticity production in bluff body flows substantially while at the same time enhancing its convection downstream and into the wake. Here, we investigate the effect of hydrodynamic slippage on the convective heat transfer (scalar transport) from a heated isothermal circular cylinder placed in a uniform cross-flow of an incompressible fluid through analytical and simulation techniques. At low Reynolds (\({\textit{Re}}\ll 1\)) and high Péclet (\({\textit{Pe}}\gg 1\)) numbers, our theoretical analysis based on Oseen and thermal boundary layer equations allows for an explicit determination of the dependence of the thermal transport on the non-dimensional slip length \(l_s\). In this case, the surface-averaged Nusselt number, Nu transitions gradually between the asymptotic limits of \(Nu \sim {\textit{Pe}}^{1/3}\) and \(Nu \sim {\textit{Pe}}^{1/2}\) for no-slip (\(l_s \rightarrow 0\)) and shear-free (\(l_s \rightarrow \infty \)) boundaries, respectively. Boundary layer analysis also shows that the scaling \(Nu \sim {\textit{Pe}}^{1/2}\) holds for a shear-free cylinder surface in the asymptotic limit of \({\textit{Re}}\gg 1\) so that the corresponding heat transfer rate becomes independent of the fluid viscosity. At finite \({\textit{Re}}\), results from our two-dimensional simulations confirm the scaling \(Nu \sim {\textit{Pe}}^{1/2}\) for a shear-free boundary over the range \(0.1 \le {\textit{Re}}\le 10^3\) and \(0.1\le {\textit{Pr}}\le 10\). A gradual transition from the lower asymptotic limit corresponding to a no-slip surface, to the upper limit for a shear-free boundary, with \(l_s\), is observed in both the maximum slip velocity and the Nu. The local time-averaged Nusselt number \(Nu_{\theta }\) for a shear-free surface exceeds the one for a no-slip surface all along the cylinder boundary except over the downstream portion where unsteady separation and flow reversal lead to an appreciable rise in the local heat transfer rates, especially at high \({\textit{Re}}\) and Pr. At a Reynolds number of \(10^3\), the formation of secondary recirculating eddy pairs results in appearance of additional local maxima in \(Nu_{\theta }\) at locations that are in close proximity to the mean secondary stagnation points. As a consequence, Nu exhibits a non-monotonic variation with \(l_s\) increasing initially from its lowermost value for a no-slip surface and then decreasing before rising gradually toward the upper asymptotic limit for a shear-free cylinder. A non-monotonic dependence of the spanwise-averaged Nu on \(l_s\) is observed in three dimensions as well with the three-dimensional wake instabilities that appear at sufficiently low \(l_s\), strongly influencing the convective thermal transport from the cylinder. The analogy between heat transfer and single-component mass transfer implies that our results can directly be applied to determine the dependency of convective mass transfer of a single solute on hydrodynamic slip length in similar configurations through straightforward replacement of Nu and \({\textit{Pr}}\) with Sherwood and Schmidt numbers, respectively.  相似文献   

8.
We consider the elliptic equation \(-\Delta u +u =0\) with nonlinear boundary condition \(\frac{\partial u}{\partial n}= \lambda u + g(\lambda ,x,u), \) where \(\frac{g(\lambda ,x,s)}{s} \rightarrow 0, \hbox { as }|s|\rightarrow \infty \) and g is oscillatory. We provide sufficient conditions on g for the existence of unbounded sequences of stable solutions, unstable solutions, and turning points, even in the absence of resonant solutions.  相似文献   

9.
We deal with one dimensional p-Laplace equation of the form
$$\begin{aligned} u_t = (|u_x|^{p-2} u_x )_x + f(x,u), \ x\in (0,l), \ t>0, \end{aligned}$$
under Dirichlet boundary condition, where \(p>2\) and \(f:[0,l]\times {\mathbb {R}}\rightarrow {\mathbb {R}}\) is a continuous function with \(f(x,0)=0\). We will prove that if there is at least one eigenvalue of the p-Laplace operator between \(\lim _{u\rightarrow 0} f(x,u)/|u|^{p-2}u\) and \(\lim _{|u|\rightarrow +\infty } f(x,u)/|u|^{p-2}u\), then there exists a nontrivial stationary solution. Moreover we show the existence of a connecting orbit between stationary solutions. The results are based on Conley index and detect stationary states even when those based on fixed point theory do not apply. In order to compute the Conley index for nonlinear semiflows deformation along p is used.
  相似文献   

10.
In the context of measure spaces equipped with a doubling non-trivial Borel measure supporting a Poincaré inequality, we derive local and global sup bounds of the nonnegative weak subsolutions of
$$\begin{aligned} (u^{q})_t-\nabla \cdot {(|\nabla u|^{p-2}\nabla u)}=0, \quad \mathrm {in} \ U_\tau = U \times (\tau _1, \tau _2] , \quad p>1,\quad q>1 \end{aligned}$$
and of its associated Dirichlet problem, respectively. For particular ranges of the exponents p and q, we show that any locally nonnegative weak subsolution, taken in \(Q (\subset \bar{Q}\subset U_\tau )\), is controlled from above by the \(L^\alpha (\bar{Q}) \)-norm, for \(\alpha = \max \{p, q+1\}\). As for the global setting, under suitable assumptions on the boundary datum g and on the initial datum, we obtain sup bounds for u, in \(U \times \{ t\}\), which depend on the \(\sup g\) and on the \(L^{q+1}(U \times (\tau _1, \tau _1+t])\)-norm of \((u-\sup g)_+\), for all \(t \in (0, \tau _2-\tau _1]\). On the critical ranges of p and q, a priori local and global \(L^\infty \) estimates require extra qualitative information on u.
  相似文献   

11.
This paper is concerned with time periodic traveling curved fronts for periodic Lotka–Volterra competition system with diffusion in two dimensional spatial space
$$\begin{aligned} {\left\{ \begin{array}{ll} \dfrac{\partial u_{1}}{\partial t}=\Delta u_{1} +u_{1}(x,y,t)\left( r_{1}(t)-a_{1}(t)u_{1}(x,y,t)-b_{1}(t)u_{2}(x,y,t)\right) ,\\ \dfrac{\partial u_{2}}{\partial t}=d\Delta u_{2} +u_{2}(x,y,t)\left( r_{2}(t)-a_{2}(t)u_{1}(x,y,t)-b_{2}(t)u_{2}(x,y,t)\right) , \end{array}\right. } \end{aligned}$$
where \(\Delta \) denotes \(\frac{\partial ^{2}}{\partial x^{2} }+ \frac{\partial ^{2}}{\partial y^{2} }\), \(x,y\in {\mathbb {R}}\) and \(d>0\) is a constant, the functions \(r_i(t),a_i(t)\) and \(b_i(t)\) are T-periodic and Hölder continuous. Under suitable assumptions that the corresponding kinetic system admits two stable periodic solutions (p(t), 0) and (0, q(t)), the existence, uniqueness and stability of one-dimensional traveling wave solution \(\left( \Phi _{1}(x+ct,t),\Phi _{2}(x+ct,t)\right) \) connecting two periodic solutions (p(t), 0) and (0, q(t)) have been established by Bao and Wang ( J Differ Equ 255:2402–2435, 2013) recently. In this paper we continue to investigate two-dimensional traveling wave solutions of the above system under the same assumptions. First, we establish the asymptotic behaviors of one-dimensional traveling wave solutions of the system at infinity. Using these asymptotic behaviors, we then construct appropriate super- and subsolutions and prove the existence and non-existence of two-dimensional time periodic traveling curved fronts. Finally, we show that the time periodic traveling curved front is asymptotically stable.
  相似文献   

12.
Motivated by some recent studies on the Allen–Cahn phase transition model with a periodic nonautonomous term, we prove the existence of complex dynamics for the second order equation
$$\begin{aligned} -\ddot{x} + \left( 1 + \varepsilon ^{-1} A(t)\right) G'(x) = 0, \end{aligned}$$
where A(t) is a nonnegative T-periodic function and \(\varepsilon > 0\) is sufficiently small. More precisely, we find a full symbolic dynamics made by solutions which oscillate between any two different strict local minima \(x_0\) and \(x_1\) of G(x). Such solutions stay close to \(x_0\) or \(x_1\) in some fixed intervals, according to any prescribed coin tossing sequence. For convenience in the exposition we consider (without loss of generality) the case \(x_0 =0\) and \(x_1 = 1\).
  相似文献   

13.
Let (XG) be a G-action topological dynamical system (t.d.s. for short), where G is a countably infinite discrete amenable group. In this paper, we study the topological pressure of the sets of generic points. We show that when the system satisfies the almost specification property, for any G-invariant measure \(\mu \) and any continuous map \(\varphi \),
$$\begin{aligned} P\left( X_{\mu },\varphi ,\{F_n\}\right) = h_{\mu }(X)+\int \varphi d\mu , \end{aligned}$$
where \(\{F_n\}\) is a Følner sequence, \(X_{\mu }\) is the set of generic points of \(\mu \) with respect to (w.r.t. for short) \(\{F_n\}\), \(P(X_{\mu },\varphi ,\{F_n\})\) is the topological pressure of \(X_{\mu }\) for \(\varphi \) w.r.t. \(\{F_n\}\) and \(h_{\mu }(X)\) is the measure-theoretic entropy.
  相似文献   

14.
The linear stability analysis of vertical throughflow of power law fluid for double-diffusive convection with Soret effect in a porous channel is investigated in this study. The upper and lower boundaries are assumed to be permeable, isothermal and isosolutal. The linear stability of vertical through flow is influenced by the interactions among the non-Newtonian Rayleigh number (Ra), Buoyancy ratio (N), Lewis number (Le), Péclet number (Pe), Soret parameter (Sr) and power law index (n). The results indicate that the Soret parameter has a significant influence on convective instability of power law fluid. It has also been noticed that buoyancy ratio has a dual effect on the instability of fluid flow. Further, it is noticed that the basic temperature and concentration profiles have singularities at \(Pe = 0\) and \(Le = 1\), the convective instability is looked into for the limiting case of \(Pe\rightarrow 0\) and \(Le \rightarrow 1\). For the case of pure thermal convection with no vertical throughflow, the present numerical results coincide with the solution of standard Horton–Rogers–Lapwood problem. The present results for critical Rayleigh number obtained using bvp4c and two-term Galerkin approximation are compared with those available in the literature and are tabulated.  相似文献   

15.
This study considers the quasilinear elliptic equation with a damping term,
$$\begin{aligned} \text {div}(D(u)\nabla u) + \frac{k(|{\mathbf {x}}|)}{|{\mathbf {x}}|}\,{\mathbf {x}}\cdot (D(u)\nabla u) + \omega ^2\big (|u|^{p-2}u + |u|^{q-2}u\big ) = 0, \end{aligned}$$
where \({\mathbf {x}}\) is an N-dimensional vector in \(\big \{{\mathbf {x}} \in \mathbb {R}^N: |{\mathbf {x}}| \ge \alpha \big \}\) for some \(\alpha > 0\) and \(N \in {\mathbb {N}}\setminus \{1\}\); \(D(u) = |\nabla u|^{p-2} + |\nabla u|^{q-2}\) with \(1 < q \le p\); k is a nonnegative and locally integrable function on \([\alpha ,\infty )\); and \(\omega \) is a positive constant. A necessary and sufficient condition is given for all radially symmetric solutions to converge to zero as \(|{\mathbf {x}}|\rightarrow \infty \). Our necessary and sufficient condition is expressed by an improper integral related to the damping coefficient k. The case that k is a power function is explained in detail.
  相似文献   

16.
In this paper we show a striking contrast in the symmetries of equilibria and extremisers of the total elastic energy of a hyperelastic incompressible annulus subject to pure displacement boundary conditions. Indeed upon considering the equilibrium equations, here, the nonlinear second order elliptic system formulated for the deformation \(u=(u_{1}, \ldots, u_{N})\):
$$ {\mathbb{E}} {\mathbb{L}}[u, {\mathbf {X}}] = \left \{ \textstyle\begin{array}{l@{\quad}l} \Delta u = \operatorname{div}(\mathscr{P} (x) \operatorname{cof} \nabla u) & \textrm{in }{\mathbf {X}},\\ \det\nabla u = 1 & \textrm{in }{\mathbf {X}},\\ u \equiv\varphi& \textrm{on }\partial{\mathbf {X}}, \end{array}\displaystyle \right . $$
where \({\mathbf {X}}\) is a finite, open, symmetric \(N\)-annulus (with \(N \ge2\)), \(\mathscr{P}=\mathscr{P}(x)\) is an unknown hydrostatic pressure field and \(\varphi\) is the identity mapping, we prove that, despite the inherent rotational symmetry in the system, when \(N=3\), the problem possesses no non-trivial symmetric equilibria whereas in sharp contrast, when \(N=2\), the problem possesses an infinite family of symmetric and topologically distinct equilibria. We extend and prove the counterparts of these results in higher dimensions by way of showing that a similar dichotomy persists between all odd vs. even dimensions \(N \ge4\) and discuss a number of closely related issues.
  相似文献   

17.
Consider a weakly nonlinear CGL equation on the torus \(\mathbb {T}^d\):
$$\begin{aligned} u_t+i\Delta u=\epsilon [\mu (-1)^{m-1}\Delta ^{m} u+b|u|^{2p}u+ ic|u|^{2q}u]. \end{aligned}$$
(*)
Here \(u=u(t,x)\), \(x\in \mathbb {T}^d\), \(0<\epsilon <<1\), \(\mu \geqslant 0\), \(b,c\in \mathbb {R}\) and \(m,p,q\in \mathbb {N}\). Define \(I(u)=(I_{\mathbf {k}},\mathbf {k}\in \mathbb {Z}^d)\), where \(I_{\mathbf {k}}=v_{\mathbf {k}}\bar{v}_{\mathbf {k}}/2\) and \(v_{\mathbf {k}}\), \(\mathbf {k}\in \mathbb {Z}^d\), are the Fourier coefficients of the function \(u\) we give. Assume that the equation \((*)\) is well posed on time intervals of order \(\epsilon ^{-1}\) and its solutions have there a-priori bounds, independent of the small parameter. Let \(u(t,x)\) solve the equation \((*)\). If \(\epsilon \) is small enough, then for \(t\lesssim {\epsilon ^{-1}}\), the quantity \(I(u(t,x))\) can be well described by solutions of an effective equation:
$$\begin{aligned} u_t=\epsilon [\mu (-1)^{m-1}\Delta ^m u+ F(u)], \end{aligned}$$
where the term \(F(u)\) can be constructed through a kind of resonant averaging of the nonlinearity \(b|u|^{2p}+ ic|u|^{2q}u\).
  相似文献   

18.
The first part of this paper is a general approach towards chaotic dynamics for a continuous map \(f:X\supset M\rightarrow X\) which employs the fixed point index and continuation. The second part deals with the differential equation
$$\begin{aligned} x'(t)=-\alpha \,x(t-d_{{\varDelta }}(x_t)). \end{aligned}$$
with state-dependent delay. For a suitable parameter \(\alpha \) close to \(5\pi /2\) we construct a delay functional \(d_{{\varDelta }}\), constant near the origin, so that the previous equation has a homoclinic solution, \(h(t)\rightarrow 0\) as \(t\rightarrow \pm \infty \), with certain regularity properties of the linearization of the semiflow along the flowline \(t\mapsto h_t\). The third part applies the method from the beginning to a return map which describes solution behaviour close to the homoclinic loop, and yields the existence of chaotic motion.
  相似文献   

19.
We consider positive classical solutions of
$$\begin{aligned} v_t=(v^{m-1}v_x)_x, \qquad x\in {\mathbb {R}}, \ t>0, \qquad (\star ) \end{aligned}$$
in the super-fast diffusion range \(m<-1\). Our main interest is in smooth positive initial data \(v_0=v(\cdot ,0)\) which decay as \(x\rightarrow +\infty \), but which are possibly unbounded as \(x\rightarrow -\infty \), having in mind monotonically decreasing data as prototypes. It is firstly proved that if \(v_0\) decays sufficiently fast only in one direction by satisfying
$$\begin{aligned} v_0(x) \le cx^{-\beta } \qquad \text{ for } \text{ all } ~x>0 \quad \hbox { with some }\quad \beta >\frac{2}{1-m} \end{aligned}$$
and some \(c>0\), then the so-called proper solution of (\(\star \)) vanishes identically in \({\mathbb {R}}\times (0,\infty )\), and accordingly no positive classical solution exists in any time interval in this case. Complemented by some sufficient criteria for solutions to remain positive either locally or globally in time, this condition for instantaneous extinction is shown to be optimal at least with respect to algebraic decay of the initial data. This partially extends some known nonexistence results for (\(\star \)) (Daskalopoulos and Del Pino in Arch Rat Mech Anal 137(4):363–380, 1997) in that it does not require any knowledge on the behavior of \(v_0(x)\) for \(x<0\). Next focusing on the phenomenon of extinction in finite time, we show that in this respect a mass influx from \(x=-\infty \) can interact with mass loss at \(x=+\infty \) in a nontrivial manner. Namely, we shall detect examples of monotone initial data, with critical decay as \(x\rightarrow +\infty \) and exponential growth as \(x\rightarrow -\infty \), that lead to solutions of (\(\star \)) which become extinct at a finite positive time, but which have empty extinction sets. This is in sharp contrast to known extinction mechanisms which are such that the corresponding extinction sets coincide with all of \({\mathbb {R}}\).
  相似文献   

20.
The flow induced by an impermeable flat surface executing orthogonal stretching and orthogonal shearing in a rotating fluid system is investigated. Both the stretching and shearing are linear in the coordinates. An exact similarity reduction of the Navier–Stokes equations gives rise to a pair of nonlinearly-coupled ordinary differential equations governed by three parameters. In this study we set one parameter and analyze the problem which leads to flow for an impermeable surface with shearing and stretching due to velocity u along the x-axis of equal strength a while the shearing and stretching due to velocity v along the y-axis of equal strength b. These solutions depend on two parameters—a Coriolis (rotation) parameter \(\sigma = \Omega /a\) and a stretching/shearing ratio \(\lambda =b/a\). A symmetry in solutions is found for \(\lambda = 1\). The exact solution for \(\sigma = 0\) and the asymptotic behavior of solutions for \(|\sigma | \rightarrow \infty\) are determined and compared with numerical results. Oscillatory solutions are found whose strength increases with increasing values of \(|\sigma |\). It is shown that these solutions tend to the well-known Ekman solution as \(|\sigma | \rightarrow \infty\).  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号