基于演化博弈的制造企业服务创新扩散研究

制造企业服务化是我国实现制造强国的必经之路,本文以“中国制造2025”战略为背景,通过借鉴创新扩散和演化博弈理论,对社会系统中参与服务创新的相关主体进行分析,构建政府与企业之间演化博弈模型,结合并改进Bass经典创新扩散模型,探究政府规制行为对企业服务创新行为演化及扩散的影响,利用Matlab对复制动态方程和扩散模型进行数学推导和数值仿真分析。结果表明:当满足政府的社会总收益大于其调控成本和各项补贴总额且企业选择服务创新策略时收益的额外增加额大于其实际投入时,系统才会演化至政府调控且企业采纳服务创新的稳定状态;适当增加政府对企业的调控力度、投入补贴和税收补贴有助于推进服务创新在系统中的扩散,其中税收补贴的感知能力最强,投入补贴最弱;将抑制系数引入Bass创新扩散模型,构建互补型、竞争型和替代型三种服务创新的扩散模式,其中竞争型扩散符合中国现阶段市场规律,替代性扩散将成为未来制造企业的发展趋势。     

Persistence of centrality in random growing trees

We investigate properties of node centrality in random growing tree models. We focus on a measure of centrality that computes the maximum subtree size of the tree rooted at each node, with the most central node being the tree centroid. For random trees grown according to a preferential attachment model, a uniform attachment model, or a diffusion processes over a regular tree, we prove that a single node persists as the tree centroid after a finite number of steps, with probability 1. Furthermore, this persistence property generalizes to the top K ≥ 1 nodes with respect to the same centrality measure. We also establish necessary and sufficient conditions for the size of an initial seed graph required to ensure persistence of a particular node with probability , as a function of ϵ: In the case of preferential and uniform attachment models, we derive bounds for the size of an initial hub constructed around the special node. In the case of a diffusion process over a regular tree, we derive bounds for the radius of an initial ball centered around the special node. Our necessary and sufficient conditions match up to constant factors for preferential attachment and diffusion tree models.     

Respectively scaled HSS iteration methods for solving discretized spatial fractional diffusion equations

For the discrete linear systems resulted from the discretization of the one‐dimensional anisotropic spatial fractional diffusion equations of variable coefficients with the shifted finite‐difference formulas of the Grünwald–Letnikov type, we propose a class of respectively scaled Hermitian and skew‐Hermitian splitting iteration method and establish its asymptotic convergence theory. The corresponding induced matrix splitting preconditioner, through further replacements of the involved Toeplitz matrices with certain circulant matrices, leads to an economic variant that can be executed by fast Fourier transforms. Both theoretical analysis and numerical implementations show that this fast respectively scaled Hermitian and skew‐Hermitian splitting preconditioner can significantly improve the computational efficiency of the Krylov subspace iteration methods employed as effective linear solvers for the target discrete linear systems.     

Bayesian estimation of incompletely observed diffusions

We present a general framework for Bayesian estimation of incompletely observed multivariate diffusion processes. Observations are assumed to be discrete in time, noisy and incomplete. We assume the drift and diffusion coefficient depend on an unknown parameter. A data-augmentation algorithm for drawing from the posterior distribution is presented which is based on simulating diffusion bridges conditional on a noisy incomplete observation at an intermediate time. The dynamics of such filtered bridges are derived and it is shown how these can be simulated using a generalised version of the guided proposals introduced in Schauer, Van der Meulen and Van Zanten (2017, Bernoulli 23(4A)).     

Pullback exponential attractors for nonautonomous dynamical system in space of higher regularity

Under what condition, a process which exists a $(E,E)$-pullback exponential attractor implies the existence of $(E,V)$- pullback exponential attractor when $V$ embedded in $E$? We answer this question in this paper. As an application of this result, we prove the existence of pullback exponential attractor for a nonlinear reaction-diffusion equation with a polynomial growth nonlinearity in $L^q(\Omega)(\forall q\geq 2)$ and $H_0^1(\Omega)$.     

Attractors for the nonclassical diffusion equations with fading memory

We discuss long-time dynamical behavior of the nonclassical diffusion equation with fading memory when nonlinearity is critical. The existence and regularity of global attractors in weak topological space and strong topological space are obtained, while the forcing term only belongs to H⁻¹(Ω) and L²(Ω) respectively. The results in this part are new and appear to be optimal corresponding to the forcing term.     

New invariant sets to nonlinear diffusion equations with x-dependent convection and absorption

In this paper, we introduce a new invariant set

     

The use of Chebyshev cardinal functions for the solution of a partial differential equation with an unknown time-dependent coefficient subject to an extra measurement

A numerical technique is presented for the solution of a parabolic partial differential equation with a time-dependent coefficient subject to an extra measurement. The method is derived by expanding the required approximate solution as the elements of Chebyshev cardinal functions. Using the operational matrix of derivative, the problem can be reduced to a set of algebraic equations. From the computational point of view, the solution obtained by this method is in excellent agreement with those obtained by previous works and also it is efficient to use.     

Taylor‐Galerkin B‐spline finite element method for the one‐dimensional advection‐diffusion equation

The advection‐diffusion equation has a long history as a benchmark for numerical methods. Taylor‐Galerkin methods are used together with the type of splines known as B‐splines to construct the approximation functions over the finite elements for the solution of time‐dependent advection‐diffusion problems. If advection dominates over diffusion, the numerical solution is difficult especially if boundary layers are to be resolved. Known test problems have been studied to demonstrate the accuracy of the method. Numerical results show the behavior of the method with emphasis on treatment of boundary conditions. Taylor‐Galerkin methods have been constructed by using both linear and quadratic B‐spline shape functions. Results shown by the method are found to be in good agreement with the exact solution. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2010