Advances in Applied Mathematics
Vol. 09  No. 05 ( 2020 ), Article ID: 35605 , 14 pages
10.12677/AAM.2020.95084

Numerical Simulation of the Brusselator Model with Spatial Spectral Interpolation Coordination Method

Tingting Ban, Yulan Wang

Inner Mongolia University of Technology, Hohhot Inner Mongolia

Received: Apr. 23rd, 2020; accepted: May 11th, 2020; published: May 18th, 2020

ABSTRACT

This paper studies the Brusselator model of a class of nonlinear reaction-diffusion equations, so it is very important to find a simple and effective numerical method for nonlinear reaction-diffusion systems. Based on this problem, some numerical examples are simulated by using this new spatial spectral interpolation collocation method.

Keywords:Brusselator Model, Turing Bifurcation Condition, Spatial Spectral Interpolation Coordination Method, The Approximate Solutions

Brusselator模型的空间谱插值配点方法 数值模拟

班亭亭,王玉兰

内蒙古工业大学,内蒙古 呼和浩特

收稿日期:2020年4月23日;录用日期:2020年5月11日;发布日期:2020年5月18日

摘 要

该文研究一类非线性反应扩散方程组Brusselator模型,因此寻找一种简单有效的非线性反应扩散系统的数值方法是非常重要的。基于这一问题,本文提出了采用这种新的空间谱插值配点方法模拟了一些数值算例,其结果和理论上的吻合较好,结果表明了该方法的有效性。

关键词 :Brusselator模型,图灵分叉条件,空间谱插值配点方法,近似解

Copyright © 2020 by author(s) and Hans Publishers Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).

http://creativecommons.org/licenses/by/4.0/

1. 引言

Brusselator模型是一个著名的模型 [1] [2],它通常是用来描述化学反应过程中化学元素变化的一类反应扩散方程组 [3]。其形式为:

{ u t = Δ ( u 11 u + u 12 v ) ( β + 1 ) u + u 2 v + α , ( x , t ) Ω × [ 0 , T ] v t = Δ ( u 21 u + u 22 v ) + β u u 2 v , ( x , t ) Ω × [ 0 , T ] u n = v n = 0 , ( x , t ) Ω × [ 0 , T ] u ( x , 0 ) = ψ 1 ( x ) , v ( x , 0 ) = ψ 2 ( x ) , x Ω (1)

其中: u ( x , y , t ) v ( x , y , t ) 是未知函数, u 11 , u 12 , u 21 , u 22 是扩散系数, Δ = 2 x 2 + 2 y 2 是拉普拉斯算子, ( x , y ) Ω , t > 0 在光滑的边界 Ω 上的齐次Neumann边界条件,即 u / n | n = v / n | n f( u,v ),g( u,v ) 是已知的光滑函数。

许多方法被用来解决Brusselator模型。这些方法包括有限差分法 [4],Galerkin有限元法 [5],最佳均匀(BURA)有理近似 [6]。本文基于一种新的空间谱插值配点方法,研究了Brusselator系统的动力学行为,分析了该系统的稳定性 [7] [8]。数值模拟结果和理论吻合较好,仿真结果表明了该方法的有效性。

2. 空间谱插值配点方法的描述

在本文中,我们使用空间谱插值配点方法来解决Brusselator模型(1)。

从参考文献 [9] [10] 中,我们可以得到离散序列 u 1 , u 2 , , u N 的插值函数 I N u ( x ) 可以写成:

u ( x ) I N u ( x ) = m = 1 N u m S N ( x x m ) (2)

其中 S N ( x ) = sin ( π x / h ) ( 2 π / h ) tan ( x / 2 ) I N 对于任何函数,都是这样的插值算子, u ( x ) 在区间 [ 0 , 2 π ] 定义为 u j = u ( j h ) , j = 1 , , N , x j x m = ( j m ) h ,插值空间为 s p a n { S N ( x j h ) , j = 1 , 2 , , N } 。由此推导出n阶导数的简单表达式并不难 I N u ( x ) x j = j h 处的n阶导数的简单表达式并不难,为:

I N u ( n ) ( x j ) = m = 1 N u m S N ( n ) ( x j x m ) (3)

其中 D N ( n ) = [ S N ( n ) ( x j x m ) ] i , j = 1 , 2 , , N 被称为第n个谱微分矩阵 [9]。

这里我们考虑有限空间域 Ω = [ 0 , 2 π ] × [ 0 , 2 π ] 。我们定义 N 2 个等间距的网格点在区域 Ω 上,因此

( x i , y j ) = ( i h , j h ) , i , j = 1 , 2 , 3 , , N .

其中 h = 2 π N 对于给定的有限自然数 N 。使用(2)则配点函数 I N u ( x , y , t ) I N v ( x , y , t ) 关于函数 u ( x , y , t ) v ( x , y , t ) 可以被写为:

u ( x , y , t ) I N u ( x , y , t ) = i = 1 N j = 1 N S N ( x x i ) S N ( y y i ) u ( x i , y j , t ) v ( x , y , t ) I N v ( x , y , t ) = i = 1 N j = 1 N S N ( x x i ) S N ( y y i ) v ( x i , y j , t ) (4)

其中 u i , j = u ( x i , y j , t ) , v i , j = v ( x i , y j , t ) , i , j = 1 , 2 , , N 。因此,下列关系在搭配点 ( x p , y q ) 处成立:

u ( x p , y q , t ) I N u ( x p , y q , t ) = i = 1 N j = 1 N S N ( x p x i ) S N ( y q y i ) u ( x i , y j , t ) v ( x p , y q , t ) I N v ( x p , y q , t ) = i = 1 N j = 1 N S N ( x p x i ) S N ( y q y i ) v ( x i , y j , t ) (5)

u ( 2 , 0 ) ( x p , y q , t ) I N u ( 2 , 0 ) ( x p , y q , t ) = 2 u ( x p , y q , t ) x 2 = i = 1 N j = 1 N S N ( 2 ) ( x p x i ) S N ( y q y i ) u ( x i , y j , t ) u ( 0 , 2 ) ( x p , y q , t ) I N u ( 0 , 2 ) ( x p , y q , t ) = 2 u ( x p , y q , t ) y 2 = i = 1 N j = 1 N S N ( x p x i ) S N ( 2 ) ( y q y i ) u ( x i , y j , t ) v ( 2 , 0 ) ( x p , y q , t ) I N v ( 2 , 0 ) ( x p , y q , t ) = 2 v ( x p , y q , t ) x 2 = i = 1 N j = 1 N S N ( 2 ) ( x p x i ) S N ( y q y i ) v ( x i , y j , t ) v ( 0 , 2 ) ( x p , y q , t ) I N v ( 0 , 2 ) ( x p , y q , t ) = 2 v ( x p , y q , t ) y 2 = i = 1 N j = 1 N S N ( x p x i ) S N ( 2 ) ( y q y i ) v ( x i , y j , t )

注意的是

u = [ u 11 , u 21 , , u N 1 , u 12 , u 22 , , u N 2 , u 1 N , , u N N ] T , v = [ v 11 , v 21 , , v N 1 , v 12 , v 22 , , v N 2 , v 1 N , , v N N ] T , (6)

因此,公式(5)可以被写成一下矩阵形式:

u ( 2 , 0 ) = D N ( 2 , 0 ) u , u ( 0 , 2 ) = D N ( 0 , 2 ) u , v ( 2 , 0 ) = D N ( 2 , 0 ) v , v ( 0 , 2 ) = D N ( 0 , 2 ) v . (7)

D N ( 2 , 0 ) u = D N ( 2 ) E N , D N ( 0 , 2 ) = E N D N ( 2 ) , D N ( 0 , 0 ) = E N E N ,

其中 E N 是N阶单位矩阵, 是矩阵的克罗内克积,相反地。二阶谱微分矩阵为:

D N ( 2 ) = ( π 2 3 h 2 1 6 1 2 csc 2 ( 1 h 2 ) 1 2 csc 2 ( 1 h 2 ) 1 2 csc 2 ( 2 h 2 ) 1 2 csc 2 ( 2 h 2 ) 1 2 csc 2 ( 2 h 2 ) 1 2 csc 2 ( 1 h 2 ) 1 2 csc 2 ( 3 h 2 ) π 2 3 h 2 1 6 1 2 csc 2 ( 3 h 2 ) 1 2 csc 2 ( 1 h 2 ) 1 2 csc 2 ( 2 h 2 ) 1 2 csc 2 ( 2 h 2 ) 1 2 csc 2 ( 1 h 2 ) 1 2 csc 2 ( 1 h 2 ) 1 2 csc 2 ( 1 h 2 ) π 2 3 h 2 1 6 ) (8)

结合等式(6)和等式(7),Equation (1)可以写成以下系统:

t [ u v ] = [ u 11 D u 12 D u 21 D u 22 D ] [ u v ] = [ f 1 ( u , v ) f 2 ( u , v ) ] (9)

这里

[ u , v ] = [ u 11 , , u N 1 , u 12 , , u N 2 , u 1 N , , u N N , v 11 , , v N 1 , v 12 , , v N 2 , v 1 N , , v N N ] ,

D = D N ( 2 , 0 ) + D N ( 0 , 2 ) = D N ( 2 ) E N + E N D N ( 2 ) ,

[ f 1 ( u , v ) , f 2 ( u , v ) ] = [ f 1 ( u 11 , v 11 ) , , f 1 ( u N N , v N N ) , f 2 ( u 11 , v 11 ) , , f 2 ( u N N , v N N ) , f 3 ( u 11 , v 11 ) , , f 3 ( u N N , v N N ) ]

利用Matlab中的ode45求解器求解系统(9),得到系统(1)的数值解。

3. 动力系统分析

在这一部分,我们给出了Brusselator模型的图灵分叉分析 [11]。如果不考虑扩散现象,即假设扩散系数为零时,则方程(1)就可以写为:

{ d u d t = α + u 2 v ( β + 1 ) u d v d t = β u u 2 v (10)

从动力学角度我们来分析系统(10)的解的性质 [12]。首先(10)只有唯一的平衡解 ( α , β / α ) 。其次,讨论这一平衡解的稳定性。为此考虑(10)的Jacobian矩阵:

D F ( u , v ) = [ 2 u v ( β + 1 ) u 2 β 2 u v u 2 ] ,

其在平衡解 ( α , β / α ) 的Jacobian矩阵为:

D F ( α , β / α ) = [ β 1 α 2 β α 2 ] (11)

显见矩阵(11)的迹 τ = t r a c D F ( α , β / α ) = α 2 + β 1 ,判别式 Δ = α 2 ( β 1 ) + α 2 β = α 2 。由于 Δ = α 2 > 0 ,所以(11)的特征根同号,故:

1) 若 τ < 0 ,即 β > α + 1 ,平衡解 ( α , β / α ) 是稳定的;

2) 若 τ > 0 ,即 β < α + 1 ,平衡解 ( α , β / α ) 是发散的。

下面讨论(10)的分岔。固定 α ,让 β 变动。由 τ = 0 得到 β H = α 2 + 1 。下面我们来证明 β H 是(10)的Hopf分岔点 [13] [14] [15],设 λ 为(11)的特征根,为此必须证明:(a) (11)在 β H 点只有一对纯虚根。(b) Re ( λ ) / β 0 。事实上,(11)的特征根为 λ 1 , 2 = ( τ ± τ 2 4 Δ ) / 2 ,由 τ = 0 Δ = α 2 λ 1 , 2 = ± α i 。所以(a)成立。其次,由 λ 1 , 2 = ( τ ± τ 2 4 Δ ) / 2 可得 Re ( λ ) = ( α 2 + β 1 ) / 2 ,所以 Re ( λ ) / B = 1 / 2 0 ,故(b)也成立。综上所述, β H 为系统(10)的Hopf分岔点,且由于 τ < 0 时系统的平衡不稳定,所以由 τ < 0 变化到 τ = 0 系统在Hopf分岔点处得到稳定性 [14]。

4. 数值实验

在这一节中,我们给出了一些数值例子,以便更好地解释使用不同的初始条件和参数的上述分析结果。

例1 考虑以下形式的Brusselator模型 [3]:

{ u t = Δ u 1 u ( β + 1 ) u + u 2 v + α , ( x , t ) Ω × [ 0 , T ] v t = Δ u 2 v + β u u 2 v , ( x , t ) Ω × [ 0 , T ] u n = v n = 0 , ( x , t ) Ω × [ 0 , T ] u ( x , 0 ) = ψ 1 ( x ) , v ( x , 0 ) = ψ 2 ( x ) , x Ω (12)

这里 α , β > 0 u 1 , u 2 为常数, t R + 为时间变量, x , y R 是空间变量, u , v 是t与 x , y 的函数。其中 u 1 = 0.0025 , u 2 = 0.0032 , α = 2 , β = 1 , v = o n e s ( N ) ,固定一个初始条件,当u变化时模拟的结果见图1~10,具体初始条件见表1

例2:考虑以下形式的Brusselator模型 [3]:

{ u t = Δ ( u 11 u + u 12 v ) ( β + 1 ) u + u 2 v + α , ( x , t ) Ω × [ 0 , T ] v t = Δ ( u 21 u + u 22 v ) + β u u 2 v , ( x , t ) Ω × [ 0 , T ] u n = v n = 0 , ( x , t ) Ω × [ 0 , T ] u ( x , 0 ) = ψ 1 ( x ) , v ( x , 0 ) = ψ 2 ( x ) , x Ω (13)

这里 α , β > 0 u 11 , u 12 , u 21 , u 22 为常数, t R + 为时间变量, x , y R 是空间变量, u , v 是t与 x , y 的函数。其中 u 11 = 0.4 , u 12 = 22.2665 , u 21 = 0.02 , u 22 = 2 , α = 6 , β = 1 ;本文方法得到的数值结果见图11~22,具体的初始条件见表2

Figure 1. Shows the numerical solution with the initial condition of u = sin ( 1000000 ( ( x 1 / 3 ) 2 ( y 1 / 2 ) 2 ) ) of example 1

图1. 例1初始条件为 u = sin ( 1000000 ( ( x 1 / 3 ) 2 ( y 1 / 2 ) 2 ) ) 的数值解

Figure 2. Shows the pattern with the initial condition of u = sin ( 1000000 ( ( x 1 / 3 ) 2 ( y 1 / 2 ) 2 ) ) example 1

图2. 例1初始条件为 u = sin ( 1000000 ( ( x 1 / 3 ) 2 ( y 1 / 2 ) 2 ) ) 的斑图

Figure 3. Shows the numerical solution with the initial condition of u = sin ( 1000 ( ( x 1 / 3 ) 2 ( y 1 / 2 ) 2 ) ) of example 1

图3. 例1初始条件是 u = sin ( 1000 ( ( x 1 / 3 ) 2 ( y 1 / 2 ) 2 ) ) 的数值解

Figure 4. Shows the pattern with the initial condition of u = sin ( 1000 ( ( x 1 / 3 ) 2 ( y 1 / 2 ) 2 ) ) example 1

图4. 例1初始条件是 u = sin ( 1000 ( ( x 1 / 3 ) 2 ( y 1 / 2 ) 2 ) ) 的斑图

Figure 5. Shows the numerical solution with the initial condition of u = sech ( 100 ( ( x 1 / 3 ) 2 ( y 1 / 2 ) 2 ) ) of example 1

图5. 例1初始条件是 u = sech ( 100 ( ( x 1 / 3 ) 2 ( y 1 / 2 ) 2 ) ) 的数值解

Figure 6. Shows the pattern with the initial condition of u = sech ( 100 ( ( x 1 / 3 ) 2 ( y 1 / 2 ) 2 ) ) example 1

图6. 例1初始条件是 u = sech ( 100 ( ( x 1 / 3 ) 2 ( y 1 / 2 ) 2 ) ) 的的斑图

Figure 7. Shows the numerical solution with the initial condition of u = sech ( x 2 / 0.02 2 π y 2 / 0.01 ) of example 1

图7. 例1初始条件是 u = sech ( x 2 / 0.02 2 π y 2 / 0.01 ) 的数值解

Figure 8. Shows the pattern with the initial condition of u = sech ( x 2 / 0.02 2 π y 2 / 0.01 ) example 1

图8. 例1初始条件是 u = sech ( x 2 / 0.02 2 π y 2 / 0.01 ) 的斑图

Figure 9. Shows the numerical solution with the initial condition of u = 12 / 25 + sin ( π x 2 510 ) + cos ( π y 220 ) of example 1

图9. 例1初始条件是 u = 12 / 25 + sin ( π x 2 510 ) + cos ( π y 220 ) 的数值解

Figure 10. Shows the pattern with the initial condition of u = 12 / 25 + sin ( π x 2 510 ) + cos ( π y 220 ) example 1

图10. 例1初始条件是 u = 12 / 25 + sin ( π x 2 510 ) + cos ( π y 220 ) 的斑图

Table 1. The different initial conditions corresponding to the numerical solution and pattern in example 1 are shown in Figures 1-10

表1. 例1中数值解和斑图所对应的不同的初始条件在图1~10

Figure 11. Shows the numerical solution with the initial condition of u = cos ( π ( ( x 0.4 ) 2 + ( y + 0.4 ) 2 ) ) + e sin ( ( x + 0.4 ) 2 + ( y 0.4 ) 2 ) of example 2

图11. 图为例2初始条件是 u = cos ( π ( ( x 0.4 ) 2 + ( y + 0.4 ) 2 ) ) + e sin ( ( x + 0.4 ) 2 + ( y 0.4 ) 2 ) 的数值解

Figure 12. Shows the pattern with the initial condition of u = cos ( π ( ( x 0.4 ) 2 + ( y + 0.4 ) 2 ) ) + e sin ( ( x + 0.4 ) 2 + ( y 0.4 ) 2 ) example 2

图12. 例2初始条件是 u = cos ( π ( ( x 0.4 ) 2 + ( y + 0.4 ) 2 ) ) + e sin ( ( x + 0.4 ) 2 + ( y 0.4 ) 2 ) 的斑图

Figure 13. Shows the numerical solution with the initial condition of u = sin ( π ( ( x 0.4 ) 2 + ( y + 0.4 ) 2 ) ) sin ( e ( x + 0.4 ) 2 + ( y 0.4 ) 2 ) of example 2

图13. 例2初始条件是 u = sin ( π ( ( x 0.4 ) 2 + ( y + 0.4 ) 2 ) ) sin ( e ( x + 0.4 ) 2 + ( y 0.4 ) 2 ) 的数值解

Figure 14. Shows the pattern with the initial condition of u = sin ( π ( ( x 0.4 ) 2 + ( y + 0.4 ) 2 ) ) sin ( e ( x + 0.4 ) 2 + ( y 0.4 ) 2 ) example 2

图14. 图为例2初始条件是 u = sin ( π ( ( x 0.4 ) 2 + ( y + 0.4 ) 2 ) ) sin ( e ( x + 0.4 ) 2 + ( y 0.4 ) 2 ) 的斑图

Figure 15. Shows the numerical solution with the initial condition of u = sech ( π ( x 2 + y 2 ) ) of example 2

图15. 例2初始条件是 u = sech ( π ( x 2 + y 2 ) ) 的数值解

Figure 16. Shows the pattern with the initial condition of u = sech ( π ( x 2 + y 2 ) ) example 2

图16. 图为例2初始条件是 u = sech ( π ( x 2 + y 2 ) ) 的斑图

Figure 17. Shows the numerical solution with the initial condition of u = sin ( π ( x 2 y 2 ) ) of example 2

图17. 例2初始条件是 u = sin ( π ( x 2 y 2 ) ) 的数值解

Figure 18. Shows the pattern with the initial condition of u = sin ( π ( x 2 y 2 ) ) example 2

图18. 图为例2初始条件是 u = sin ( π ( x 2 y 2 ) ) 的斑图

Figure 19. Shows the numerical solution with the initial condition of u = 12 / 25 sin ( 10 x 2 510 ) + sin ( π y 220 ) of example 2

图19. 例2初始条件是 u = 12 / 25 sin ( 10 x 2 510 ) + sin ( π y 220 ) 的数值解

Figure 20. Shows the pattern with the initial condition of u = 12 / 25 sin ( 10 x 2 510 ) + sin ( π y 220 ) example 2

图20. 图为例2初始条件是 u = 12 / 25 sin ( 10 x 2 510 ) + sin ( π y 220 ) 的斑图

Figure 21. Shows the numerical solution with the initial condition of u = 12 / 25 sin ( π x 2 510 ) + sin ( cos ( x 2 + y 2 22 ) ) of example 2

图21. 例2初始条件是 u = 12 / 25 sin ( π x 2 510 ) + sin ( cos ( x 2 + y 2 22 ) ) 的数值解

Figure 22. Shows the pattern with the initial condition of u = 12 / 25 sin ( π x 2 510 ) + sin ( cos ( x 2 + y 2 22 ) ) example 2

图22. 图为例2初始条件是 u = 12 / 25 sin ( π x 2 510 ) + sin ( cos ( x 2 + y 2 22 ) ) 的斑图

Table 2. The different initial conditions corresponding to the numerical solution and pattern in example 2 are shown in Figures 11-22

表2. 例2中数值解和斑图所对应的不同的初始条件在图 11~22

5. 结论

本文用谱插值配点方法求解Brusselator模型,数值结果表明了该方法与理论吻合较好。本文所有程序由matlab2017b算得。

致谢

感谢王玉兰老师的支持与帮助。

文章引用

班亭亭,王玉兰. Brusselator模型的空间谱插值配点方法数值模拟
Numerical Simulation of the Brusselator Model with Spatial Spectral Interpolation Coordination Method[J]. 应用数学进展, 2020, 09(05): 708-721. https://doi.org/10.12677/AAM.2020.95084

参考文献

  1. 1. Kang, H. and Pesin, Y. (2005) Dynamics of a Discrete Brusselator Model: Escape to Infinity and Julia Set. Milan Journal of Mathematics, 73, 1-17. https://doi.org/10.1007/s00032-005-0036-y

  2. 2. Mittal, R.C. and Jiwari, R. (2011) Numerical Solution of Two-Dimensional Reaction-C Diffusion Brusselator System. Applied Mathematics & Computation, 217, 5404-5415. https://doi.org/10.1016/j.amc.2010.12.010

  3. 3. Tyson, J.J. (1976) The Belousov-Zhabotinskii Reaction. In: Lecture Notes in Biomathematics 10, Springer-Verlag, Heidelberg. https://doi.org/10.1007/978-3-642-93046-1

  4. 4. Wang, Y.M. and Zhang, H.B. (2009) Higher-Order Compact Finite Difference Method for Systems of Reaction Diffusion Equations. Journal of Computational and Applied Mathematics, 233, 502-518. https://doi.org/10.1016/j.cam.2009.07.052

  5. 5. Macha, J., Benes, M. and Strachota, P. (2017) Nonlinear Galerkin Finite Element Method Applied to the System of Reaction-C Diffusion Equations in One Space Dimension. Computers and Mathematics with Applications, 73, 2053- 2065. https://doi.org/10.1016/j.camwa.2017.02.032

  6. 6. Harizanov, S., Lazarov, R., Margenov, S., et al. (2019) Numerical Solution of Fractional Diffusion Creaction Problems Based on BURA. Computers and Mathematics with Applications. https://doi.org/10.1016/j.camwa.2019.07.002

  7. 7. 唐晓栋. 多反馈反应扩散系统斑图动力学研究[D]: [博士学位论文]. 徐州: 中国矿业大学, 2014.

  8. 8. Siraj-ul-Islam, A.A. and Haq, S. (2010) A Computational Modeling of the Behavior of the Two-Dimensional Reaction-Diffusion Brusselator System. Applied Mathematical Modelling, 34, 3896-3909. https://doi.org/10.1016/j.apm.2010.03.028

  9. 9. Zhang, X. (2015) Efficient Solution of Differential Equation Based on MATLAB: Spectral Method Principle and Implementation. Mechanical Industry Press, Beijing.

  10. 10. 李庆阳, 王能超, 易大义. 数值分析[M]. 第5版. 北京: 清华大学出版社, 2008.

  11. 11. Peng, Y.H. and Ling, H. (2018) Pattern formation in a Ratio-Dependent Predator-Prey Model with Cross-Diffusion. Applied Mathematics and Computation, 331, 307-318. https://doi.org/10.1016/j.amc.2018.03.033

  12. 12. 熊晓华, 周天寿. Brusselator模型的定性分析及数值模拟[J]. 江西师范大学学报(自然科学版), 2007, 31(2): 164-168.

  13. 13. Guckenheimer, J. and Holmes, P. (1999) Nonlinear Oscillations, Dynamical Systems, and Bifurcations of Vector Fields. Springer-Verlag, New York.

  14. 14. Sey del, R. (1999) Practical Bifurcation and Stability Analysis. From Equilibrium to Chaos, Second Edition, Springer-Verlag, New York.

  15. 15. 周天寿. 线性耦合振子中的回声波及其稳定性[J]. 江西师范大学学报: 自然科学版, 2006, 30(1): 1-6.

期刊菜单