逼近已知函数微商的广义Lanczos算法

给出逼近已知函数微商的广义Lanczos 算法, 构造了一列逼近算子$D_{h}^{n}$以提高稳定近似解的收敛速率. 当$n=2$时, 逼近精度达到$O(\delta^{6 \over 7})$, 而对一般的自然数$n$逼近精度为$O(\delta^{\frac{2n+2}{2n+3}})$, 这里$\delta$是近似函数的误差界.     

一类$L^2(\mathbb{R}^n)$的紧支撑不可分正交小波

本篇文章给出一类$L^{2}(\mathbb{R}^{n})$, $n\geq2$的紧支撑不可分正交小波基的具体构造算法,其中正交小波的伸缩矩阵为$\alpha I_{n}~(\alpha\geq2,\ \alpha \in \mathbb{Z})$, $I_{n}$是$n$阶单位矩阵.最后给出两个不可分正交小波基的构造算例.     

在Freud正交多项式零点处的Gr\"{u}nwald插值算子的收敛性

设$W_{\beta}(x)=\exp(-\frac{1}{2}|x|^{\beta})~(\beta > 7/6)$ 为Freud权, Freud正交多项式定义为满足下式$\int_{- \infty}^{\infty}p_{n}(x)p_{m}(x)W_{\beta}^{2}(x)\rd x=\left \{ \begin{array}{ll} 0 & \hspace{3mm} n \neq m , \\ 1 & \hspace{3mm}n = m \end{array} \right.$的     

一类具有奇异灵敏度的logistic趋化系统解的渐近行为

本文在无边界流的光滑有界区域$\Omega\subset\mathbb{R}^n~(n>2)$上研究了具有奇异灵敏度及logistic源的抛物-椭圆趋化系统$$\left\{\begin{array}{ll}u_t=\Delta u-\chi\nabla\cdot(\frac{u}{v}\nabla v)+r u-\mu u^k,&x\in\Omega,\,t>0,\\ 0=\Delta v-v+u,&x\in\Omega,\,t>0\end{array}\right.$$ 其中$\chi$, $r$, $\mu>0$, $k\geq2$. 证明了若当$r$适当大, 则当$t\rightarrow\infty$时该趋化系统全局有界解呈指数收敛于$((\frac{r}{\mu})^{\frac{1}{k-1}}, (\frac{r}{\mu})^{\frac{1}{k-1}})$.     

严格上三角矩阵李代数上的积零导子

设$F$ 为域, $n\geq 3$, $\bf{N}$$(n,\mathbb{F})$ 为域$\mathbb{F}$ 上所有$n\times n$ 阶严格上三角矩阵构成的严格上三角矩阵李代数, 其李运算为$[x,y]=xy-yx$. $\bf{N}$$(n, \mathbb{F})$ 上一线性映射$\varphi$ 称为积零导子,如果由$[x,y]=0, x,y\in \bf{N}$$(n,\mathbb{F})$,总可推出 $[\varphi(x), y]+[x,\varphi(y)]=0$. 本文证明 $\bf{N}$$(n,\mathbb{F})$上一线性映射 $\varphi$ 为积零导子当且仅当 $\varphi$ 为$\bf{N}$$(n,\mathbb{F})$ 上内导子, 对角线导子, 极端导子, 中心导子和标量乘法的和.     

无线网络中全调度问题的一种随机分布式算法

无线网络中的全调度,要确保网络中每个节点所可能的链路信息和广播信息都能无冲突地进行传输.通过简单的构造方法,证明了多项式时间内,能找到一个长度为$O(\bigtriangleup_{\rm out}^2\bigtriangleup_{\rm in})$的全调度;并且给出了全调度问题的一种随机分布式算法,证明了这种随机分布式算法,对任意的常数$h$,~$0     

关于奇次矩形元恢复导数的强超收敛性的进一步研究

通过推广林群等的超收敛结果及Green函数估计, 对矩形元利用SPR技巧给出了一种强超收敛方法, 证明了在局部对称点上导数具有$O(h^{k+3})$~($k\geq 3$ 为奇数)的强超收敛阶及位移具有$O(h^{k+4})$~($k\geq 4$ 为偶数)的强超收敛性.     

变系数混合效应模型

本文研究了下列变系数混合效应模型: $y_{ij}=z_{ij}^{\tau}b_i+x_{ij}^{\tau}\beta(w_{ij}) +\xe_{ij},\;i=1,\cdots,m;\;j=1,\cdots,n_i$, 其中$b_i$为i.i.d.期望为$\xt$, 协方差阵为$\xs^2_bI_q$的随机效应向量, $\xe_{ij}$是i.i.d.期望为零, 具有有限方差的随机误差. 文中我们不仅给出了函数系数向量$\xb(\cdot)$的局部多项式估计, 同时给出了随机效应期望、方差和随机误差方差的估计, 并给出了这些估计量的渐进正态性和相合性, 研究结果表明了这些估计量的可靠性.     

在Banach空间中推广的 Roper-Suffridge算子(III)

假定 $X$ 是具有范数$\|\cdot\|$的复 Banach 空间, $n$ 是一个满足 $\dim X\geq n\geq2$的正整数. 本文考虑由下式定义的推广的Roper-Suffridge算子 $\Phi_{n,\beta_2, \gamma_2, \ldots , \beta_{n+1}, \gamma_{n+1}}(f)$: \begin{equation} \begin{array}{lll} \Phi _{n, \beta_2, \gamma_2, \ldots, \beta_{n+1},\gamma_{n+1}}(f)(x) &;\hspace{-3mm}=&;\hspace{-3mm}\dl\he{j=1}{n}\bigg(\frac{f(x^*_1(x))}{x^*_1(x)})\bigg)^{\beta_j}(f''(x^*_1(x))^{\gamma_j}x^*_j(x) x_j\\ &;&;+\bigg(\dl\frac{f(x^*_1(x))}{x^*_1(x)}\bigg)^{\beta_{n+1}}(f''(x^*_1(x)))^{\gamma_{n+1}}\bigg(x-\dl\he{j=1}{n}x^*_j(x) x_j\bigg),\nonumber \end{array} \end{equation} 其中 $x\in\Omega_{p_1, p_2, \ldots, p_{n+1}}$, $\beta_1=1, \gamma_1=0$ 和 \begin{equation} \begin{array}{lll} \Omega_{p_1, p_2, \ldots, p_{n+1}}=\bigg\{x\in X: \dl\he{j=1}{n}| x^*_j(x)|^{p_j}+\bigg\|x-\dl\he{j=1}{n}x^*_j(x)x_j\bigg\|^{p_{n+1}}<1\bigg\},\nonumber \end{array} \end{equation} 这里 $p_j>1 \,( j=1, 2,\ldots, n+1$), 线性无关族 $\{x_1, x_2, \ldots, x_n \}\subset X $ 与 $\{x^*_1, x^*_2, \ldots, x^*_n \}\subset X^* $ 满足 $x^*_j(x_j)=\|x_j\|=1 (j=1, 2, \ldots, n)$ 和 $x^*_j(x_k)=0 \, (j\neq k)$, 我们选取幂函数的单值分支满足 $(\frac{f(\xi)}{\xi})^{\beta_j}|_{\xi=0}= 1$ 和 $(f''(\xi))^{\gamma_j}|_{\xi=0}=1, \, j=2, \ldots , n+1$. 本文将证明: 对某些合适的常数$\beta_j, \gamma_j$, 算子$\Phi_{n,\beta_2, \gamma_2, \ldots, \beta_{n+1}, \gamma_{n+1}}(f)$ 在$\Omega_{p_1, p_2, \ldots , p_{n+1}}$上保持$\alpha$阶的殆$\beta$型螺形映照和 $\alpha$阶的$\beta$型螺形映照.     

强连通双圈有向图的无符号拉普拉斯谱刻画

设$\overrightarrow{G}$ 是一个强连通双圈有向图, $A(\overrightarrow{G})$是其邻接矩阵.设$D(\overrightarrow{G})$ 是$\overrightarrow{G}$的顶点出度的对角矩阵, $Q(\overrightarrow{G})=D(\overrightarrow{G})+A(\overrightarrow{G})$是$\overrightarrow{G}$ 的无符号拉普拉斯矩阵. $Q(\overrightarrow{G})$的谱半径称为$\overrightarrow{G}$的无符号拉普拉斯谱半径.在这篇文章中, 确定了在所有强连通双圈有向图中达到最大或最小无符号拉普拉斯谱半径的唯一有向图. 此外,还证明了任意一个强连通双圈有向图是由它的无符号拉普拉斯谱所确定的.     

A spatially second-order accurate implicit numerical method for the space and time fractional Bloch-Torrey equation

In recent years, it has been found that many phenomena in engineering, physics, chemistry and other sciences can be described very successfully by models using mathematical tools from Fractional Calculus. Recently, a new space and time fractional Bloch-Torrey equation (ST-FBTE) has been proposed (Magin et al., J. Magn. Reson. 190(2), 255–270, 2008), and successfully applied to analyse diffusion images of human brain tissues to provide new insights for further investigations of tissue structures. In this paper, we consider the ST-FBTE with a nonlinear source term on a finite domain in three-dimensions. The time and space derivatives in the ST-FBTE are replaced by the Caputo and the sequential Riesz fractional derivatives, respectively. Firstly, we propose a spatially second-order accurate implicit numerical method (INM) for the ST-FBTE whereby we discretize the Riesz fractional derivative using a fractional centered difference. Secondly, we prove that the implicit numerical method for the ST-FBTE is uniquely solvable, unconditionally stable and convergent, and the order of convergence of the implicit numerical method is \(O\left (\tau ^{2-\alpha }+\tau +h_{x}^{2}+h_{y}^{2}+h_{z}^{2}\right )\) . Finally, some numerical results are presented to support our theoretical analysis.     

LIMIT CYCLES OF THE GENERALIZED POLYNOMIAL LINARD DIFFERENTIAL SYSTEMS

Using the averaging theory of first and second order we study the maximum number of limit cycles of generalized Linard differential systems{x = y + εh_l~1(x) + ε~2h_l~2(x),y=-x- ε(f_n~1(x)y~(2p+1) + g_m~1(x)) + ∈~2(f_n~2(x)y~(2p+1) + g_m~2(x)),which bifurcate from the periodic orbits of the linear center x = y,y=-x,where ε is a small parameter.The polynomials h_l~1 and h_l~2 have degree l;f_n~1and f_n~2 have degree n;and g_m~1,g_m~2 have degree m.p ∈ N and[·]denotes the integer part function.     

Compact finite difference schemes of the time fractional Black-Scholes model

In this paper, three compact difference schemes for the time-fractional Black-Scholes model governing European option pricing are presented. Firstly, in order to obtain the fourth-order accuracy in space by applying the Pad\''{e} approximation, we eliminate the convection term of the B-S equation by an exponential transformation. Then the time fractional derivative is approximated by $L1$ formula, $L2 - 1_\sigma$ formula and $L1 - 2$ formula respectively, and three compact difference schemes with oders $O(\Delta t^{2-\alpha}+h ^4)$, $O(\Delta t^{2}+h ^4)$ and $O(\Delta t^{3-\alpha}+h ^4)$ are constructed. Finally, numerical example is carried out to verify the accuracy and effectiveness of proposed methods, and the comparisons of various schemes are given. The paper also provides numerical studies including the effect of fractional orders and the effect of different parameters on option price in time-fractional B-S model.     

A Nonhomogeneous Boundary Value Problem for the Boussinesq Equation on a Bounded Domain

In this paper, we study the well-posedness of an initial-boundary-value problem (IBVP) for the Boussinesq equation on a bounded domain,\begin{cases} &u_{tt}-u_{xx}+(u^2)_{xx}+u_{xxxx}=0,\quad x\in (0,1), \;\;t>0,\\ &u(x,0)=\varphi(x),\;\;\; u_t(x,0)=ψ(x),\\ &u(0,t)=h_1(t),\;\;\;u(1,t)=h_2(t),\;\;\;u_{xx}(0,t)=h_3(t),\;\;\;u_{xx}(1,t)=h_4(t).\\ \end{cases} It is shown that the IBVP is locally well-posed in the space $H^s (0,1)$ for any $s\geq 0$ with the initial data $\varphi,$ $\psi$ lie in $H^s(0,1)$ and $ H^{s-2}(0,1)$, respectively, and the naturally compatible boundary data $h_1,$ $h_2$ in the space $H_{loc}^{(s+1)/2}(\mathbb{R}^+)$, and $h_3 $, $h_4$ in the the space of $H_{loc}^{(s-1)/2}(\mathbb{R}^+)$ with optimal regularity.     

Boundedness Characterization of Maximal Commutators on Orlicz Spaces in the Dunkl Setting

On the real line, the Dunkl operators$$D_{\nu}(f)(x):=\frac{d f(x)}{dx} + (2\nu+1) \frac{f(x) - f(-x)}{2x}, ~~ \quad\forall \, x \in \mathbb{R}, ~ \forall \, \nu \ge -\tfrac{1}{2}$$are differential-difference operators associated with the reflection group $\mathbb{Z}_2$ on $\mathbb{R}$, and on the $\mathbb{R}^d$ the Dunkl operators $\big\{D_{k,j}\big\}_{j=1}^{d}$ are the differential-difference operators associated with the reflection group $\mathbb{Z}_2^d$ on $\mathbb{R}^{d}$.In this paper, in the setting $\mathbb{R}$ we show that $b \in BMO(\mathbb{R},dm_{\nu})$ if and only if the maximal commutator $M_{b,\nu}$ is bounded on Orlicz spaces $L_{\Phi}(\mathbb{R},dm_{\nu})$. Also in the setting $\mathbb{R}^{d}$ we show that $b \in BMO(\mathbb{R}^{d},h_{k}^{2}(x) dx)$ if and only if the maximal commutator $M_{b,k}$ is bounded on Orlicz spaces $L_{\Phi}(\mathbb{R}^{d},h_{k}^{2}(x) dx)$.     

Uniformly Accurate Multiscale Time Integrators for Highly Oscillatory Second Order Differential Equations

In this paper, two multiscale time integrators (MTIs), motivated from two types of multiscale decomposition by either frequency or frequency and amplitude, are proposed and analyzed for solving highly oscillatory second order differential equations with a dimensionless parameter $0 < \varepsilon≤ 1.$ In fact, the solution to this equation propagates waves with wavelength at $O(\varepsilon^2)$ when $0<\varepsilon≪1,$ which brings significantly numerical burdens in practical computation. We rigorously establish two independent error bounds for the two MTIs at $O(\tau^2/\varepsilon^2)$ and $O(\varepsilon^2)$ for $\varepsilon ∈ (0,1]$ with $\tau > 0$ as step size, which imply that the two MTIs converge uniformly with linear convergence rate at $O(\tau)$ for $ε ∈ (0,1]$ and optimally with quadratic convergence rate at $O(\tau^2)$ in the regimes when either $ε=O(1)$ or $0<ε≤\tau.$ Thus the meshing strategy requirement (or $ε$-scalability) of the two MTIs is $\tau =O(1)$ for $0<ε≪1,$ which is significantly improved from $\tau =O(ε^3)$ and $\tau =O(ε^2)$ requested by finite difference methods and exponential wave integrators to the equation, respectively. Extensive numerical tests and comparisons with those classical numerical integrators are reported, which gear towards better understanding on the convergence and resolution properties of the two MTIs. In addition, numerical results support the two error bounds very well.     

A Trilayer Difference Scheme for One-Dimensional Parabolic Systems

In order to obtain the numerical solution for a one-dimensional parabolic system, an unconditionally stable difference method is investigated in [1]. If the number of unknown functions is M, for each time step only M times of calculation are needed. The rate of convergence is $O(\tau+h^2)$. On the basis of [1], an alternating calculation difference scheme is presented in [2]; the rate of the convergence is $O(\tau^2+h^2)$. The difference schemes in [1] and [2] are economic ones. For the $\alpha$-$th$ equation, only $U_{\alpha}$ is an unknown function; the others $U_{\beta}$ are given evaluated either in the last step or in the present step. So the practical calculation is quite convenient. The purpose of this paper is to derive a trilayer difference scheme for one-dimensional parabolic systems. It is known that the scheme is also unconditionally stable and the rate of convergence is $O(\tau^2+h^2)$.     

EXPANSIONS OF STEP-TRANSITION OPERATORS OF MULTI-STEP METHODS AND ORDER BARRIERS FOR DAHLQUIST PAIRS

Using least parameters, we expand the step-transition operator of any linear multi-step method （LMSM） up to O（τ^s＋5） with order s = 1 and rewrite the expansion of the steptransition operator for s = 2 （obtained by the second author in a former paper）. We prove that in the conjugate relation G3^λτ o G1^τ =G2^τ o G3^λτ with G1 being an LMSM,（1） theorder of G2 can not be higher than that of G1; （2） if G3 is also an LMSM and G2 is a symplectic B-series, then the orders of G1, G2 and G3 must be 2, 2 and 1 respectively.     

A Multiplicity Result for a Singular and Nonhomogeneous Elliptic Problem in Rn

We establish sufficient conditions under which the quasilinear equation $$-div(|∇u|^{n-2}∇u)+V(x)|u|^{n-2}u=\frac{f(x,u)}{|x|^β}+εh(x) in \mathbb{R}^n,$$ has at least two nontrivial weak solutions in $W^{1,n} (\mathbb{R}^n)$ when ε > 0 is small enough, 0≤β < n, V is a continuous potential, f(x,u) behaves like $exp{γ|u|^{n/(n-1)}}$ as $|u|→∞$ for some γ > 0 and h≢ 0 belongs to the dual space of $W^{1,n} (\mathbb{R}^n)$.     

A Fourth Order Finite Difference Approximation to the Eigenvalues Approximation to the Eigenvalues of a Clamped Plate

In a 21-point finite difference scheme, assign suitable interpolation values to the fictitious node points. The numerical eigenvalues are then of $O(h^2)$ precision. But the corrected value $\hat{λ}_h=λ_h+\frac{h^2}{6}λ_h^{\frac{3}{2}}$ and extrapolation $\hatλ_h=\frac{4}{3}λ_{\frac{λ}{2}}-\frac{1}{3}λ_h$can be proved to have $O(h^4)$ precision.