针对旋转倒立摆自动起摆的控制问题,提出了基于BVP算法的自动起摆控制策略。该方法将倒立摆起摆控制问题转化成求解非线性方程的两点边值问题(Two-point BoundaryValue Problem BVP),构造了含参变量具有傅立叶级数形式的起摆力矩函数,...针对旋转倒立摆自动起摆的控制问题,提出了基于BVP算法的自动起摆控制策略。该方法将倒立摆起摆控制问题转化成求解非线性方程的两点边值问题(Two-point BoundaryValue Problem BVP),构造了含参变量具有傅立叶级数形式的起摆力矩函数,将力矩函数代入倒立摆系统,利用Matlab工具箱中的bvp4c函数求解两点边值条件,获得起摆过程的起摆控制的时间序列。基于BVP算法的起摆控制的求解,本质上属于开环前馈控制。为了抑制参数摄动,进行了平衡点附近的稳摆控制设计。稳摆设计是针对系统模型不稳定性和非最小相位特性分别进行的。对起摆、稳摆及其切换过程进行了仿真和实验研究,验证了所提出的自动起摆控制策略的有效性。展开更多
A conventional complex variable boundary integral equation (CVBIE) in plane elasticity is provided. After using the Somigliana identity between a particular fundamental stress field and a physical stress field, an a...A conventional complex variable boundary integral equation (CVBIE) in plane elasticity is provided. After using the Somigliana identity between a particular fundamental stress field and a physical stress field, an additional integral equality is obtained. By adding both sides of this integral equality to both sides of the conventional CVBIE, the amended boundary integral equation (BIE) is obtained. The method based on the discretization of the amended BIE is called the amended influence matrix method. With this method, for the Neumann boundary value problem (BVP) of an interior region, a unique solution for the displacement can be obtained. Several numerical examples are provided to prove the efficiency of the suggested method.展开更多
This paper is devoted to find the numerical solutions of one dimensional general nonlinear system of third-order boundary value problems (BVPs) for the pair of functions using Galerkin weighted residual method. We der...This paper is devoted to find the numerical solutions of one dimensional general nonlinear system of third-order boundary value problems (BVPs) for the pair of functions using Galerkin weighted residual method. We derive mathematical formulations in matrix form, in detail, by exploiting Bernstein polynomials as basis functions. A reasonable accuracy is found when the proposed method is used on few examples. At the end of the study, a comparison is made between the approximate and exact solutions, and also with the solutions of the existing methods. Our results converge monotonically to the exact solutions. In addition, we show that the derived formulations may be applicable by reducing higher order complicated BVP into a lower order system of BVPs, and the performance of the numerical solutions is satisfactory. .展开更多
The existence of positive solutions to second-order periodic BVPs-u'+Mu =j(t, u),t(0) = u(2π),u'(0) = '(2π) and u'+ Mu = I(t, u), u(0) = u(2π), u'(0) = u'(2π)is proved by a simple appliCati...The existence of positive solutions to second-order periodic BVPs-u'+Mu =j(t, u),t(0) = u(2π),u'(0) = '(2π) and u'+ Mu = I(t, u), u(0) = u(2π), u'(0) = u'(2π)is proved by a simple appliCation of a Fixed point Theorem in cones due to Krasnoselskii.展开更多
本文运用数值计算的方法研究了包含周期刺激项Acos(ω_0t)的BVP(Bonhoeffer-van der Pol)模型随控制参数ω_0变化出现的分岔、混沌和阵发混沌.结果表明,受环境周期性变化的影响,在适当的条件下,一类可兴奋细胞系统表现出的复杂的电生理...本文运用数值计算的方法研究了包含周期刺激项Acos(ω_0t)的BVP(Bonhoeffer-van der Pol)模型随控制参数ω_0变化出现的分岔、混沌和阵发混沌.结果表明,受环境周期性变化的影响,在适当的条件下,一类可兴奋细胞系统表现出的复杂的电生理活动是与混沌行为密切相关的.展开更多
The existence of positive radial solutions to the second order semilinear elliptic BVPΔu(X)+g(|X|)f(u(X))=0, R 1<|X|<R 2, u(X)=0, |X|=R 1 or |X|=R 2is considered. A general existence criterion and se...The existence of positive radial solutions to the second order semilinear elliptic BVPΔu(X)+g(|X|)f(u(X))=0, R 1<|X|<R 2, u(X)=0, |X|=R 1 or |X|=R 2is considered. A general existence criterion and several existence theorems of positive radial solution are established. Here it is not required that lim l→0f(l)/l and lim l→∞f(l)/l exist.展开更多
We explore the complicated bursting oscillations as well as the mechanism in a high-dimensional dynamical system.By introducing a periodically changed electrical power source in a coupled BVP oscillator, a fifth-order...We explore the complicated bursting oscillations as well as the mechanism in a high-dimensional dynamical system.By introducing a periodically changed electrical power source in a coupled BVP oscillator, a fifth-order vector field with two scales in frequency domain is established when an order gap exists between the natural frequency and the exciting frequency.Upon the analysis of the generalized autonomous system, bifurcation sets are derived, which divide the parameter space into several regions associated with different types of dynamical behaviors. Two typical cases are focused on as examples,in which different types of bursting oscillations such as sub Hopf/sub Hopf burster, sub Hopf/fold-cycle burster, and doublefold/fold burster can be observed. By employing the transformed phase portraits, the bifurcation mechanism of the bursting oscillations is presented, which reveals that different bifurcations occurring at the transition between the quiescent states(QSs) and the repetitive spiking states(SPs) may result in different forms of bursting oscillations. Furthermore, because of the inertia of the movement, delay may exist between the locations of the bifurcation points on the trajectory and the bifurcation points obtained theoretically.展开更多
In this paper, we discuss the uniform convergence of the simple upwind scheme on the Shishkin mesh and the Bakhvalov-Shishkin mesh for solving a singularly perturbed Robin boundary value problem, and investigate the m...In this paper, we discuss the uniform convergence of the simple upwind scheme on the Shishkin mesh and the Bakhvalov-Shishkin mesh for solving a singularly perturbed Robin boundary value problem, and investigate the midpoint upwind scheme on the Shishkin mesh and the Bakhvalov-Shishkin mesh to achieve better uniform convergence. The elaborate ε-uniform pointwise estimates are proved by using the comparison principle and barrier functions. The numerical experiments support the theoretical results for the schemes on the meshes.展开更多
This paper presents derivation of a priori error estimates and convergence rates of finite element processes for boundary value problems (BVPs) described by self adjoint, non-self adjoint, and nonlinear differential o...This paper presents derivation of a priori error estimates and convergence rates of finite element processes for boundary value problems (BVPs) described by self adjoint, non-self adjoint, and nonlinear differential operators. A posteriori error estimates are discussed in context with local approximations in higher order scalar product spaces. A posteriori error computational framework (without the knowledge of theoretical solution) is presented for all BVPs regardless of the method of approximation employed in constructing the integral form. This enables computations of local errors as well as the global errors in the computed finite element solutions. The two most significant and essential aspects of the research presented in this paper that enable all of the features described above are: 1) ensuring variational consistency of the integral form(s) resulting from the methods of approximation for self adjoint, non-self adjoint, and nonlinear differential operators and 2) choosing local approximations for the elements of a discretization in a subspace of a higher order scalar product space that is minimally conforming, hence ensuring desired global differentiability of the approximations over the discretizations. It is shown that when the theoretical solution of a BVP is analytic, the a priori error estimate (in the asymptotic range, discussed in a later section of the paper) is independent of the method of approximation or the nature of the differential operator provided the resulting integral form is variationally consistent. Thus, the finite element processes utilizing integral forms based on different methods of approximation but resulting in VC integral forms result in the same a priori error estimate and convergence rate. It is shown that a variationally consistent (VC) integral form has best approximation property in some norm, conversely an integral form with best approximation property in some norm is variationally consistent. That is best approximation property of the integral form and the VC of the integral form is equivalent, one cannot exist without the other, hence can be used interchangeably. Dimensional model problems consisting of diffusion equation, convection-diffusion equation, and Burgers equation described by self adjoint, non-self adjoint, and nonlinear differential operators are considered to present extensive numerical studies using Galerkin method with weak form (GM/WF) and least squares process (LSP) to determine computed convergence rates of various error norms and present comparisons with the theoretical convergence rates.展开更多
The boundary value problems of the third-order ordinary differential equation have many practical application backgrounds and their some special cases have been studied by many authors. However, few scholars have stud...The boundary value problems of the third-order ordinary differential equation have many practical application backgrounds and their some special cases have been studied by many authors. However, few scholars have studied the boundary value problems of the complete third-order differential equations u′′′(t) = f (t,u(t),u′(t),u′′(t)). In this paper, we discuss the existence and uniqueness of solutions and positive solutions of the fully third-order ordinary differential equation on [0,1] with the boundary condition u(0) = u′(1) = u′′(1) = 0. Under some inequality conditions on nonlinearity f some new existence and uniqueness results of solutions and positive solutions are obtained.展开更多
The present paper proposes a mathematical method to numerically treat a class of third-order linear Boundary Value Problems (BVPs). This method is based on the combination of the Adomian Decomposition Method (ADM) and...The present paper proposes a mathematical method to numerically treat a class of third-order linear Boundary Value Problems (BVPs). This method is based on the combination of the Adomian Decomposition Method (ADM) and, the modified shooting method. A complete derivation of the proposed method has been provided, in addition to its numerical implementation and, validation via the utilization of the Runge-Kutta method and, other existing methods. The method has been applied to diverse test problems and turned out to perform remarkably. Lastly, the simulated numerical results have been graphically illustrated and, also supported by some absolute error comparison tables.展开更多
文摘针对旋转倒立摆自动起摆的控制问题,提出了基于BVP算法的自动起摆控制策略。该方法将倒立摆起摆控制问题转化成求解非线性方程的两点边值问题(Two-point BoundaryValue Problem BVP),构造了含参变量具有傅立叶级数形式的起摆力矩函数,将力矩函数代入倒立摆系统,利用Matlab工具箱中的bvp4c函数求解两点边值条件,获得起摆过程的起摆控制的时间序列。基于BVP算法的起摆控制的求解,本质上属于开环前馈控制。为了抑制参数摄动,进行了平衡点附近的稳摆控制设计。稳摆设计是针对系统模型不稳定性和非最小相位特性分别进行的。对起摆、稳摆及其切换过程进行了仿真和实验研究,验证了所提出的自动起摆控制策略的有效性。
文摘A conventional complex variable boundary integral equation (CVBIE) in plane elasticity is provided. After using the Somigliana identity between a particular fundamental stress field and a physical stress field, an additional integral equality is obtained. By adding both sides of this integral equality to both sides of the conventional CVBIE, the amended boundary integral equation (BIE) is obtained. The method based on the discretization of the amended BIE is called the amended influence matrix method. With this method, for the Neumann boundary value problem (BVP) of an interior region, a unique solution for the displacement can be obtained. Several numerical examples are provided to prove the efficiency of the suggested method.
文摘This paper is devoted to find the numerical solutions of one dimensional general nonlinear system of third-order boundary value problems (BVPs) for the pair of functions using Galerkin weighted residual method. We derive mathematical formulations in matrix form, in detail, by exploiting Bernstein polynomials as basis functions. A reasonable accuracy is found when the proposed method is used on few examples. At the end of the study, a comparison is made between the approximate and exact solutions, and also with the solutions of the existing methods. Our results converge monotonically to the exact solutions. In addition, we show that the derived formulations may be applicable by reducing higher order complicated BVP into a lower order system of BVPs, and the performance of the numerical solutions is satisfactory. .
文摘The existence of positive solutions to second-order periodic BVPs-u'+Mu =j(t, u),t(0) = u(2π),u'(0) = '(2π) and u'+ Mu = I(t, u), u(0) = u(2π), u'(0) = u'(2π)is proved by a simple appliCation of a Fixed point Theorem in cones due to Krasnoselskii.
文摘The existence of positive radial solutions to the second order semilinear elliptic BVPΔu(X)+g(|X|)f(u(X))=0, R 1<|X|<R 2, u(X)=0, |X|=R 1 or |X|=R 2is considered. A general existence criterion and several existence theorems of positive radial solution are established. Here it is not required that lim l→0f(l)/l and lim l→∞f(l)/l exist.
基金Project supported by the National Natural Science Foundation of China(Grant No.21276115)
文摘We explore the complicated bursting oscillations as well as the mechanism in a high-dimensional dynamical system.By introducing a periodically changed electrical power source in a coupled BVP oscillator, a fifth-order vector field with two scales in frequency domain is established when an order gap exists between the natural frequency and the exciting frequency.Upon the analysis of the generalized autonomous system, bifurcation sets are derived, which divide the parameter space into several regions associated with different types of dynamical behaviors. Two typical cases are focused on as examples,in which different types of bursting oscillations such as sub Hopf/sub Hopf burster, sub Hopf/fold-cycle burster, and doublefold/fold burster can be observed. By employing the transformed phase portraits, the bifurcation mechanism of the bursting oscillations is presented, which reveals that different bifurcations occurring at the transition between the quiescent states(QSs) and the repetitive spiking states(SPs) may result in different forms of bursting oscillations. Furthermore, because of the inertia of the movement, delay may exist between the locations of the bifurcation points on the trajectory and the bifurcation points obtained theoretically.
文摘In this paper, we discuss the uniform convergence of the simple upwind scheme on the Shishkin mesh and the Bakhvalov-Shishkin mesh for solving a singularly perturbed Robin boundary value problem, and investigate the midpoint upwind scheme on the Shishkin mesh and the Bakhvalov-Shishkin mesh to achieve better uniform convergence. The elaborate ε-uniform pointwise estimates are proved by using the comparison principle and barrier functions. The numerical experiments support the theoretical results for the schemes on the meshes.
文摘This paper presents derivation of a priori error estimates and convergence rates of finite element processes for boundary value problems (BVPs) described by self adjoint, non-self adjoint, and nonlinear differential operators. A posteriori error estimates are discussed in context with local approximations in higher order scalar product spaces. A posteriori error computational framework (without the knowledge of theoretical solution) is presented for all BVPs regardless of the method of approximation employed in constructing the integral form. This enables computations of local errors as well as the global errors in the computed finite element solutions. The two most significant and essential aspects of the research presented in this paper that enable all of the features described above are: 1) ensuring variational consistency of the integral form(s) resulting from the methods of approximation for self adjoint, non-self adjoint, and nonlinear differential operators and 2) choosing local approximations for the elements of a discretization in a subspace of a higher order scalar product space that is minimally conforming, hence ensuring desired global differentiability of the approximations over the discretizations. It is shown that when the theoretical solution of a BVP is analytic, the a priori error estimate (in the asymptotic range, discussed in a later section of the paper) is independent of the method of approximation or the nature of the differential operator provided the resulting integral form is variationally consistent. Thus, the finite element processes utilizing integral forms based on different methods of approximation but resulting in VC integral forms result in the same a priori error estimate and convergence rate. It is shown that a variationally consistent (VC) integral form has best approximation property in some norm, conversely an integral form with best approximation property in some norm is variationally consistent. That is best approximation property of the integral form and the VC of the integral form is equivalent, one cannot exist without the other, hence can be used interchangeably. Dimensional model problems consisting of diffusion equation, convection-diffusion equation, and Burgers equation described by self adjoint, non-self adjoint, and nonlinear differential operators are considered to present extensive numerical studies using Galerkin method with weak form (GM/WF) and least squares process (LSP) to determine computed convergence rates of various error norms and present comparisons with the theoretical convergence rates.
文摘The boundary value problems of the third-order ordinary differential equation have many practical application backgrounds and their some special cases have been studied by many authors. However, few scholars have studied the boundary value problems of the complete third-order differential equations u′′′(t) = f (t,u(t),u′(t),u′′(t)). In this paper, we discuss the existence and uniqueness of solutions and positive solutions of the fully third-order ordinary differential equation on [0,1] with the boundary condition u(0) = u′(1) = u′′(1) = 0. Under some inequality conditions on nonlinearity f some new existence and uniqueness results of solutions and positive solutions are obtained.
文摘The present paper proposes a mathematical method to numerically treat a class of third-order linear Boundary Value Problems (BVPs). This method is based on the combination of the Adomian Decomposition Method (ADM) and, the modified shooting method. A complete derivation of the proposed method has been provided, in addition to its numerical implementation and, validation via the utilization of the Runge-Kutta method and, other existing methods. The method has been applied to diverse test problems and turned out to perform remarkably. Lastly, the simulated numerical results have been graphically illustrated and, also supported by some absolute error comparison tables.