求解非线性Black-Scholes模型的自适应小波精细积分法

针对非线性Black-Scholes方程,基于quasi-Shannon小波函数给出了一种求解非线性偏微分方程的自适应多尺度小波精细积分法.该方法首先利用插值小波理论构造了用于逼近连续函数的多尺度小波插值算子,利用该算子可以将非线性Black-Scholes方程自适应离散为非线性常微分方程组;然后将用于求解常微分方程组的精细积分法和小波变换的动态过程相结合,并利用非线性处理技术(如同伦分析技术)可有效求解非线性Black-Scholes方程.数值结果表明了该方法在数值精度和计算效率方面的优越性.     

基于拟Shannon小波浅水长波近似方程组的数值解

本文研究了浅水长波近似方程组初边值问题的数值解.利用小波多尺度分析和区间拟Shannon小波,对浅水长波近似方程组空间导数实施空间离散,用时间步长自适应精细积分法对其变换所的非线性常微分方程组进行求解,得到了浅水长波近似方程组的数值解,并将此方法计算的结果与其解析解进行比较和验证.     

一类非线性微分方程组的有理化Haar小波解法

利用有理化Haar小波性质和方法,建立了一类非线性微分方程组在任意区间[a,b）的求解算法．基于该算法,运用计算机代数系统Maple,给出了求解非线性微分方程组的程序．并运用此程序给出了一类微分方程组的计算实例,从数值模拟来看可以达到较高的精度,并对方程组的动力学行为给出较好的描述．     

基于Riccati-Bernoulli辅助常微分方程的Davey-Stewartson方程的行波解

Riccati-Bernoulli辅助常微分方程方法可以用来构造非线性偏微分方程的行波解.利用行波变换,将非线性偏微分方程化为非线性常微分方程, 再利用Riccati-Bernoulli方程将非线性常微分方程化为非线性代数方程组, 求解非线性代数方程组就能直接得到非线性偏微分方程的行波解.对Davey-Stewartson方程应用这种方法, 得到了该方程的精确行波解.同时也得到了该方程的一个Backlund变换.所得结果与首次积分法的结果作了比较.Riccati-Bernoulli辅助常微分方程方法是一种简单、有效地求解非线性偏微分方程精确解的方法.     

斜拉桥拉索-阻尼器系统非线性瞬态响应分析

考虑索的抗弯刚度、垂度及几何非线性的影响,得出了索-阻尼器系统的空间非线性振动偏微分方程,用中心差分法将偏微分方程在空间内离散,导出了系统的非线性振动常微分方程组．结合Newmark法及虚拟力法提出了一种用于求解非线性振动瞬态响应的杂交分析算法．并以典型的斜拉桥拉索为研究对象,给出了数值算例,并与Runge-Kutta直接积分法进行了比较,说明了杂交算法的准确性及有效性．为大跨斜拉桥拉索的振动控制研究提供了一种简便、有效、快速的时程分析方法．     

非线性动力系统的自适应显式Magnus数值方法

基于最近发展的矩阵李群上非线性微分方程的显式Magnus展式,给出了非线性动力系统的有效的数值算法,并且在数值求解过程中具有自适应的步长控制特点,可以显著地提高计算效率．最后,通过非线性动力系统典型问题Duffing方程和强刚性的Van derPol方程以及非线性振子的Hamilton方程的数值实验来说明方法的有效性．     

Navier—Stokes方程区域分解法的收敛性

０引言区域分解方法是近年来迅速发展的偏微分方程数值方法．区域分解方法及其收敛性的研究大多是在线性偏微分方程下得到的,对于非线性问题,经典的技巧在收敛性证明时遇到了困难．流体计算是一个较为复杂的非线性问题,数值模拟过程中因节点多．网格复杂,所以计算量很大．由于区域分解方法不但可以缩小求解规模,进行并行计算,而且可以在不同区域选取不同离散方法和模型,因此对Ｎ－Ｓ方程区域分解方法的研究会有较高的实用价值,也可以对其它非线性问题数值方法研究提供新的途径．本文首先给出了Ｎ－Ｓ方程的最优控制方法以及一些重要…     

正弦电磁场中铁磁材料数学模型

将电磁场理论与弹性力学理论相结合,建立了描述铁磁材料在正弦电磁场中的数学模型,并对该模型一类的4阶非线性偏微分方程的解进行了讨论．给出其一阶近似后得到的线性偏微分方程的解析表达式和数值计算方法．计算结果表明,本方法是有效的．     

解非线性椭圆型方程奇异摄动问题的一致收敛差分格式

给出了求解非线性椭圆型偏微分方程奇异摄动问题的广义OCI差分格式．证明了这种格式的解关于摄动参数一致收敛于连续问题的解．给出了数值例子．     

多步双变量Chebyshev谱配置方法求解初边值问题

提出了单步和多步双变量Chebyshev配置方法,用于求解非线性发展型偏微分方程的初边值问题.单步格式容易实施并且具有谱精度,并给出了多步方法的收敛性分析.数值实验表明:多步双变量Chebyshev谱配置方法在非线性发展型偏微分方程问题求解中是非常有效的,与理论分析一致,特别适合于长时间问题的数值模拟.     

Adaptive Wavelet Methods for Elliptic Operator Equations with Nonlinear Terms

We present an adaptive wavelet method for the numerical solution of elliptic operator equations with nonlinear terms. This method is developed based on tree approximations for the solution of the equations and adaptive fast reconstruction of nonlinear functionals of wavelet expansions. We introduce a constructive greedy scheme for the construction of such tree approximations. Adaptive strategies of both continuous and discrete versions are proposed. We prove that these adaptive methods generate approximate solutions with optimal order in both of convergence and computational complexity when the solutions have certain degree of Besov regularity.     

Numerical solution of nonlinear Volterra integro-differential equations of arbitrary order by CAS wavelets

A computational method for numerical solution of a nonlinear Volterra integro-differential equation of fractional (arbitrary) order which is based on CAS wavelets and BPFs is introduced. The CAS wavelet operational matrix of fractional integration is derived and used to transform the main equation to a system of algebraic equations. Some examples are included to demonstrate the validity and applicability of the technique.     

Computational inelasticity for loading conditions on multiple time scales by adaptive step size control

This contribution proposes an algorithm based on adaptive step size control for the simulation of inelastic solids and structures undergoing loading conditions at multiple time scales. Adaptivity in time integration of viscoelastic constitutive laws is directed by an refinement indicator which is constructed from integrators of different order, here a fourth-order Runge-Kutta (RK) method and linear Backward-Euler. The key novel aspect is that by virtue of an recently established consistency condition the higher order methods, p ≥ 2, can achieve their full nominal order without fulfilling the weak form of balance of linear momentum in the RK stages, but only at the end of the time interval. A representative numerical example illustrates the performance of the present adaptive method and underpins the computational savings compared with uniform time step sizes. (© 2017 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A CAS wavelet method for solving nonlinear Fredholm integro-differential equations of fractional order

In this paper we present a computational method for solving a class of nonlinear Fredholm integro-differential equations of fractional order which is based on CAS (Cosine And Sine) wavelets. The CAS wavelet operational matrix of fractional integration is derived and used to transform the equation to a system of algebraic equations. Some examples are included to demonstrate the validity and applicability of the technique.     

Solving fractional nonlinear Fredholm integro-differential equations by the second kind Chebyshev wavelet

In this paper, we first construct the second kind Chebyshev wavelet. Then we present a computational method based on the second kind Chebyshev wavelet for solving a class of nonlinear Fredholm integro-differential equations of fractional order. The second kind Chebyshev wavelet operational matrix of fractional integration is derived and used to transform the equation to a system of algebraic equations. The method is illustrated by applications and the results obtained are compared with the existing ones in open literature. Moreover, comparing the methodology with the known technique shows that the present approach is more efficient and more accurate.     

Asymptotics of generalized derangements

An adaptive wavelet-based method is proposed for solving TV(total variation)–Allen–Cahn type models for multi-phase image segmentation. The adaptive algorithm integrates (i) grid adaptation based on a threshold of the sparse wavelet representation of the locally-structured solution; and (ii) effective finite difference on irregular stencils. The compactly supported interpolating-type wavelets enjoy very fast wavelet transforms, and act as a piecewise constant function filter. These lead to fairly sparse computational grids, and relax the stiffness of the nonlinear PDEs. Equipped with this algorithm, the proposed sharp interface model becomes very effective for multi-phase image segmentation. This method is also applied to image restoration and similar advantages are observed.     

An optimal adaptive wavelet method without coarsening of the iterands

In this paper, an adaptive wavelet method for solving linear operator equations is constructed that is a modification of the method from [Math. Comp, 70 (2001), pp. 27-75] by Cohen, Dahmen and DeVore, in the sense that there is no recurrent coarsening of the iterands. Despite this, it will be shown that the method has optimal computational complexity. Numerical results for a simple model problem indicate that the new method is more efficient than an existing alternative adaptive wavelet method.

     

Uncertainty investigations in nonlinear aeroelastic systems

In this paper, the stochastic collocation method (SCM) is applied to investigate the nonlinear behavior of an aeroelastic system with uncertainties in the system parameter and the initial condition. Numerical case studies for problems with uncertainties are carried out. In particular, the performance of the SCM is compared with solutions based on other computational techniques such as Monte Carlo simulation, Wiener chaos expansion and wavelet chaos expansion. From the computational results, we conclude that the SCM is an effective tool to study a nonlinear aeroelastic system with random parameters.