Two approximation laws of sliding mode for discrete-time variable structure control systems are proposed to overcome the limitations of the exponential approximation law and the variable rate approximation law. By app...Two approximation laws of sliding mode for discrete-time variable structure control systems are proposed to overcome the limitations of the exponential approximation law and the variable rate approximation law. By applying the proposed approximation laws of sliding mode to discrete-time variable structure control systems, the stability of origin can be guaranteed, and the chattering along the switching surface caused by discrete-time variable structure control can be restrained effectively. In designing of approximation laws, the problem that the system control input is restricted is also considered, which is very important in practical systems. Finally a simulation example shows the effectiveness of the two approximation laws proposed.展开更多
The numerical algorithms for finding the lines of branching and branching-off solutions of nonlinear problem on mean-square approximation of a real finite nonnegative function with respect to two variables by the modu...The numerical algorithms for finding the lines of branching and branching-off solutions of nonlinear problem on mean-square approximation of a real finite nonnegative function with respect to two variables by the modulus of double discrete Fourier transform dependent on two parameters, are constructed and justified.展开更多
We show that the best L_p-approximant to continuous functions by n-convex functions is the limit of discrete n-convex approximations.The techniques of the proof are then used to show the existence of near interpolants...We show that the best L_p-approximant to continuous functions by n-convex functions is the limit of discrete n-convex approximations.The techniques of the proof are then used to show the existence of near interpolants to discrete n-convex data by continuous n-convex functions if the data points are close.展开更多
A nonlinear problem of mean-square approximation of a real nonnegative continuous function with respect to two variables by the modulus of double Fourier integral dependent on two real parameters with use of the smoot...A nonlinear problem of mean-square approximation of a real nonnegative continuous function with respect to two variables by the modulus of double Fourier integral dependent on two real parameters with use of the smoothing functional is studied. Finding the optimal solutions of this problem is reduced to solution of the Hammerstein type two-dimensional nonlinear integral equation. The numerical algorithms to find the branching lines and branching-off solutions of this equation are constructed and justified. Numerical examples are presented.展开更多
In this paper, we apply homotopy analysis method to solve discrete mKdV equation and successfully obtain the bell-shaped solitary solution to mKdV equation. Comparison between our solution and the exact solution shows...In this paper, we apply homotopy analysis method to solve discrete mKdV equation and successfully obtain the bell-shaped solitary solution to mKdV equation. Comparison between our solution and the exact solution shows that homotopy analysis method is effective and validity in solving hybrid nonlinear problems, including solitary solution of difference-differential equation.展开更多
Nonlinearity has a crucial impact on the symmetry properties of dynamical systems. This paper studies a one-dimensional mixed Klein-Gordon/Fermi Pasta-Ulam diatomic chain using the expanded rotating plane-wave approxi...Nonlinearity has a crucial impact on the symmetry properties of dynamical systems. This paper studies a one-dimensional mixed Klein-Gordon/Fermi Pasta-Ulam diatomic chain using the expanded rotating plane-wave approximation and numerical calculations to determine the effect of cubic potentials on the symmetry properties of discrete breathers in this system. The results will be very useful to researchers in the field of numerical calculations on discrete breathers.展开更多
This paper deals with the inertial manifold and the approximate inertialmanifold concepts of the Navier-Stokes equations with nonhomogeneous boundary conditions and inertial algorithm. Furtheremore,we provide the erro...This paper deals with the inertial manifold and the approximate inertialmanifold concepts of the Navier-Stokes equations with nonhomogeneous boundary conditions and inertial algorithm. Furtheremore,we provide the error estimates of the approximate solutions of the Navier-Stokes Equations.展开更多
Consistency and asymptotic normality of the sieve estimator and an approximate maximum likelihood estimator of the drift coefficient of an interacting particles of diffusions are studied. For the sieve estimator, obse...Consistency and asymptotic normality of the sieve estimator and an approximate maximum likelihood estimator of the drift coefficient of an interacting particles of diffusions are studied. For the sieve estimator, observations are taken on a fixed time interval [0,T] and asymptotics are studied as the number of interacting particles increases with the dimension of the sieve. For the approximate maximum likelihood estimator, discrete observations are taken in a time interval [0,T] and asymptotics are studied as the number of interacting particles increases with the number of observation time points.展开更多
The aim of this study is to compare the Discrete wavelet decomposition and the modified Principal Analysis Component (PCA) decomposition to analyze the stabilogram for the purpose to provide a new insight about human ...The aim of this study is to compare the Discrete wavelet decomposition and the modified Principal Analysis Component (PCA) decomposition to analyze the stabilogram for the purpose to provide a new insight about human postural stability. Discrete wavelet analysis is used to decompose the stabilogram into several timescale components (i.e. detail wavelet coefficients and approximation wavelet coefficients). Whereas, the modified PCA decomposition is applied to decompose the stabilogram into three components, namely: trend, rambling and trembling. Based on the modified PCA analysis, the trace of analytic trembling and rambling in the complex plan highlights a unique rotation center. The same property is found when considering the detail wavelet coefficients. Based on this property, the area of the circle in which 95% of the trace’s data points are located, is extracted to provide important information about the postural equilibrium status of healthy subjects (average age 31 ± 11 years). Based on experimental results, this parameter seems to be a valuable parameter in order to highlight the effect of visual entries, stabilogram direction, gender and age on the postural stability. Obtained results show also that wavelets and the modified PCA decomposition can discriminate the subjects by gender which is particularly interesting in biometric applications and human stability simulation. Moreover, both techniques highlight the fact that male are less stable than female and the fact that there is no correlation between human stability and his age (under 60).展开更多
This paper proposes an extension of the algorithm in [1], as well as utilization of the wavelet transform in event detection, including High Impedance Fault (HIF). Techniques to analyze the abundant data of PMUs quick...This paper proposes an extension of the algorithm in [1], as well as utilization of the wavelet transform in event detection, including High Impedance Fault (HIF). Techniques to analyze the abundant data of PMUs quickly and effectively are paramount to increasing response time to events and unstable parameters. With the amount of data PMUs output, unstable parameters, tie line oscillations, and HIFs are often overlooked in the bulk of the data. This paper explores model-free techniques to attain stability information and determine events in real-time. When full system connectivity is unknown, many traditional methods requiring other bus measurements can be impossible or computationally extensive to apply. The traditional method of interest is analyzing the power flow Jacobian for singularities and system weak points, attained by applying singular value decomposition. This paper further develops upon the approach in [1] to expand the Discrete-Time Jacobian Eigenvalue Approximation (DDJEA), giving values to significant off-diagonal terms while establishing a generalized connectivity between correlated buses. Statistical linear models are applied over large data sets to prove significance to each term. Then the off diagonal terms are given time-varying weights to account for changes in topology or sensitivity to events using a reduced system model. The results of this novel method are compared to the present errors of the previous publication in order to quantify the degree of improvement that this novel method imposes. The effective bus eigenvalues are briefly compared to Prony analysis to check similarities. An additional application for biorthogonal wavelets is also introduced to detect event types, including the HIF, for PMU data.展开更多
In his series of three papers we study singularly perturbed (SP) boundary valueproblems for equations of elliptic and parabolic type. For small values of the pertur-bation parameter parabolic boundary and interior lay...In his series of three papers we study singularly perturbed (SP) boundary valueproblems for equations of elliptic and parabolic type. For small values of the pertur-bation parameter parabolic boundary and interior layers appear in these problems.If classical discretisation methods are used, the solution of the finite differencescheme and the approximation of the diffusive flux do not converge uniformly withrespect to this parameter. Using the method of special, adapted grids, we canconstruct difference schemes that allow approximation of the solution and the nor-malised diffusive flux uniformly with respect to the small parameter.We also consider sillgularly perturbed boundary value problems for convection-diffusion equations. Also for these problems we construct special finite differenceschemes, the solution of which converges ε-uniformly We study what problems ap-pear, when classical schemes are used for the approximation of the spatial deriva-tives. We compare the results with those obtained by the adapted approach. Re-sults of numerical experiments are discussed.In the three papers we first give an introduction on the general problem, andthen we consider respectively (i) Problems for SP parabolic equations, for whichthe solution and the normalised diffusive fluxes are required; (ii) Problems for SPelliptic equations with boundary conditions of Diriclilet, Neumann and RDbin type;(iii) Problems for SP parabolic equation with discontinuous boundary conditions-展开更多
We use the method of discrete dipole approximation with surface interaction to construct a model in which a plurality of nanoparticles is arranged on the surface of BK7 glass. Nanoparticles are in air medium illuminat...We use the method of discrete dipole approximation with surface interaction to construct a model in which a plurality of nanoparticles is arranged on the surface of BK7 glass. Nanoparticles are in air medium illuminated by evanescent wave generated from total internal reflection. The effects of the wavelength, the polarization of the incident wave, the number of nanoparticles and the spacing of multiple nanoparticles on the field enhancement and extinction efficiency are calculated by our model. Our work could pave the way to improve the field enhancement of multiple nanoparticles systems.展开更多
As a first step towards the numerical analysis of the stochastic primitive equations of the atmosphere and the oceans, the time discretization of these equations by an implicit Euler scheme is studied. From the determ...As a first step towards the numerical analysis of the stochastic primitive equations of the atmosphere and the oceans, the time discretization of these equations by an implicit Euler scheme is studied. From the deterministic point of view, the 3D primitive equations are studied in their full form on a general domain and with physically realistic boundary conditions. From the probabilistic viewpoint, this paper deals with a wide class of nonlinear, state dependent, white noise forcings which may be interpreted in either the Itor the Stratonovich sense. The proof of convergence of the Euler scheme,which is carried out within an abstract framework, covers the equations for the oceans, the atmosphere, the coupled oceanic-atmospheric system as well as other related geophysical equations. The authors obtain the existence of solutions which are weak in both the PDE and probabilistic sense, a result which is new by itself to the best of our knowledge.展开更多
This study tested a novel method designed to provide useful information for medical diagnosis and treatment. We measured electroencephalography (EEG) during a test of eye opening and closing, a common test in routine ...This study tested a novel method designed to provide useful information for medical diagnosis and treatment. We measured electroencephalography (EEG) during a test of eye opening and closing, a common test in routine EEG examination. This test is mainly used for measuring the degree of alpha blocking and sensitivity during eyes opening and closing. However, because these factors depend on the subject’s awareness, drowsiness can interfere with accurate diagnosis. We sought to determine the optimal EEG frequency band and optimal brain region for distinguishing healthy individuals from patients suffering from several neurophysiological diseases (including dementia, cerebrovascular disorder, schizophrenia, alcoholism, and epilepsy) while fully awake, and while in an early drowsy state. We tested four groups of subjects (awake healthy subjects, drowsy healthy subjects, awake patients and drowsy patients). The complexity of EEG band frequencies over five lobes in the human brain was analyzed using wavelet-based approximate entropy (ApEn). Two-way analysis of variance tested the effects of the two factors of interest (subjects’ health state, and subjects’ wakefulness state) on five different lobes of the brain during eyes opening and closing. The complexity of the theta and delta bands over frontal and central regions, respectively, was significantly greater in the healthy state during eyes opening. In contrast, patients exhibited increased complexity of gamma band activity over the temporal region only, during eyes-close. The early drowsy state and wakefulness state increased the complexity of theta band activity over the temporal region only during eyes-close and eyes-open states respectively, and this change was significantly greater in control subjects compared with patients. We propose that this method may be useful in routine EEG examination, to aid medical doctors and clinicians in distinguishing healthy individuals from patients, regardless of whether the subject is fully awake or in the early stages of drowsiness.展开更多
Every univariate Hermite interpolation problem can be written as a pointwise limit of Lagrange interpolants.However,this property is not preserved for the multivariate case.In this paper,the authors first generalize t...Every univariate Hermite interpolation problem can be written as a pointwise limit of Lagrange interpolants.However,this property is not preserved for the multivariate case.In this paper,the authors first generalize the result of P.Gniadek.As an application,the authors consider the discrete approximation problem for a special case when the interpolation condition contains all partial derivatives of order less than n and one nth order differential polynomial.In addition,for the case of n≥3,the authors use the concept of Cartesian tensors to give a sufficient condition to find a sequence of discrete points,such that the Lagrange interpolation problems at these points converge to the given Hermite-type interpolant.展开更多
We consider the scattering of light in participating media composed of sparsely and randomly distributed discrete particles.The particle size is expected to range from the scale of the wavelength to several orders of ...We consider the scattering of light in participating media composed of sparsely and randomly distributed discrete particles.The particle size is expected to range from the scale of the wavelength to several orders of magnitude greater,resulting in an appearance with distinct graininess as opposed to the smooth appearance of continuous media.One fundamental issue in the physically-based synthesis of such appearance is to determine the necessary optical properties in every local region.Since these properties vary spatially,we resort to geometrical optics approximation(GOA),a highly efficient alternative to rigorous Lorenz–Mie theory,to quantitatively represent the scattering of a single particle.This enables us to quickly compute bulk optical properties for any particle size distribution.We then use a practical Monte Carlo rendering solution to solve energy transfer in the discrete participating media.Our proposed framework is the first to simulate a wide range of discrete participating media with different levels of graininess,converging to the continuous media case as the particle concentration increases.展开更多
In this series of three papers we study singularly perturbed (SP) boundaryvalue problems for equations of eiliptic and parabolic type- For small values ofthe perturbation parameter parabolic boundary and interior laye...In this series of three papers we study singularly perturbed (SP) boundaryvalue problems for equations of eiliptic and parabolic type- For small values ofthe perturbation parameter parabolic boundary and interior layers appear in theseproblems. If classical discretisation methods are used, the solution of the finitedifference scheme and the approximation of the diffusive flux do not converge uniformly with respect to this parameter. Using the method of special, edapted grids,we can construct difference schemes that allow apprcximation of the solution andthe normalised diffusive flux uniformly with respect to the small parameter.We also consider singularly perturbed boundary value problems for convection-diffusion equations. Also for these problems we construct special finite differenceschemes, the solution of which converges ε-uniformly. We study what problems appear, when classical schemes are used for the approximation of the spatial derivatives. We compare the results with those obtained by the adapted approach. Re-sults of numerical experiments are discussed.In the three papers we first give an introduction on the general problem, andthen we consider respectively (i) Problems for SP parabolic equations, for whichthe solution and the normalised diffusive fluxes are required; (ii) Problems for SPelliptic equations with boundary conditions of Dirichlet, Neumann and Robin type;(iii) Problems for SP parabolic equation with discontinuous boundary conditions.展开更多
基金This work was supported by the National Natural Science Foundation of China (No.60274099) and the Foundation of Key Laboratory of Process Industry Automation, Ministry of Education
文摘Two approximation laws of sliding mode for discrete-time variable structure control systems are proposed to overcome the limitations of the exponential approximation law and the variable rate approximation law. By applying the proposed approximation laws of sliding mode to discrete-time variable structure control systems, the stability of origin can be guaranteed, and the chattering along the switching surface caused by discrete-time variable structure control can be restrained effectively. In designing of approximation laws, the problem that the system control input is restricted is also considered, which is very important in practical systems. Finally a simulation example shows the effectiveness of the two approximation laws proposed.
文摘The numerical algorithms for finding the lines of branching and branching-off solutions of nonlinear problem on mean-square approximation of a real finite nonnegative function with respect to two variables by the modulus of double discrete Fourier transform dependent on two parameters, are constructed and justified.
文摘We show that the best L_p-approximant to continuous functions by n-convex functions is the limit of discrete n-convex approximations.The techniques of the proof are then used to show the existence of near interpolants to discrete n-convex data by continuous n-convex functions if the data points are close.
文摘A nonlinear problem of mean-square approximation of a real nonnegative continuous function with respect to two variables by the modulus of double Fourier integral dependent on two real parameters with use of the smoothing functional is studied. Finding the optimal solutions of this problem is reduced to solution of the Hammerstein type two-dimensional nonlinear integral equation. The numerical algorithms to find the branching lines and branching-off solutions of this equation are constructed and justified. Numerical examples are presented.
基金the State Key Basic Research Program of China under Grant No.2004CB318000
文摘In this paper, we apply homotopy analysis method to solve discrete mKdV equation and successfully obtain the bell-shaped solitary solution to mKdV equation. Comparison between our solution and the exact solution shows that homotopy analysis method is effective and validity in solving hybrid nonlinear problems, including solitary solution of difference-differential equation.
基金Project supported by the National Natural Science Foundation of China (Grant No.10574011)the Foundation for Innovative Research Groups Foundation of Beijing Normal University
文摘Nonlinearity has a crucial impact on the symmetry properties of dynamical systems. This paper studies a one-dimensional mixed Klein-Gordon/Fermi Pasta-Ulam diatomic chain using the expanded rotating plane-wave approximation and numerical calculations to determine the effect of cubic potentials on the symmetry properties of discrete breathers in this system. The results will be very useful to researchers in the field of numerical calculations on discrete breathers.
文摘This paper deals with the inertial manifold and the approximate inertialmanifold concepts of the Navier-Stokes equations with nonhomogeneous boundary conditions and inertial algorithm. Furtheremore,we provide the error estimates of the approximate solutions of the Navier-Stokes Equations.
文摘Consistency and asymptotic normality of the sieve estimator and an approximate maximum likelihood estimator of the drift coefficient of an interacting particles of diffusions are studied. For the sieve estimator, observations are taken on a fixed time interval [0,T] and asymptotics are studied as the number of interacting particles increases with the dimension of the sieve. For the approximate maximum likelihood estimator, discrete observations are taken in a time interval [0,T] and asymptotics are studied as the number of interacting particles increases with the number of observation time points.
文摘The aim of this study is to compare the Discrete wavelet decomposition and the modified Principal Analysis Component (PCA) decomposition to analyze the stabilogram for the purpose to provide a new insight about human postural stability. Discrete wavelet analysis is used to decompose the stabilogram into several timescale components (i.e. detail wavelet coefficients and approximation wavelet coefficients). Whereas, the modified PCA decomposition is applied to decompose the stabilogram into three components, namely: trend, rambling and trembling. Based on the modified PCA analysis, the trace of analytic trembling and rambling in the complex plan highlights a unique rotation center. The same property is found when considering the detail wavelet coefficients. Based on this property, the area of the circle in which 95% of the trace’s data points are located, is extracted to provide important information about the postural equilibrium status of healthy subjects (average age 31 ± 11 years). Based on experimental results, this parameter seems to be a valuable parameter in order to highlight the effect of visual entries, stabilogram direction, gender and age on the postural stability. Obtained results show also that wavelets and the modified PCA decomposition can discriminate the subjects by gender which is particularly interesting in biometric applications and human stability simulation. Moreover, both techniques highlight the fact that male are less stable than female and the fact that there is no correlation between human stability and his age (under 60).
文摘This paper proposes an extension of the algorithm in [1], as well as utilization of the wavelet transform in event detection, including High Impedance Fault (HIF). Techniques to analyze the abundant data of PMUs quickly and effectively are paramount to increasing response time to events and unstable parameters. With the amount of data PMUs output, unstable parameters, tie line oscillations, and HIFs are often overlooked in the bulk of the data. This paper explores model-free techniques to attain stability information and determine events in real-time. When full system connectivity is unknown, many traditional methods requiring other bus measurements can be impossible or computationally extensive to apply. The traditional method of interest is analyzing the power flow Jacobian for singularities and system weak points, attained by applying singular value decomposition. This paper further develops upon the approach in [1] to expand the Discrete-Time Jacobian Eigenvalue Approximation (DDJEA), giving values to significant off-diagonal terms while establishing a generalized connectivity between correlated buses. Statistical linear models are applied over large data sets to prove significance to each term. Then the off diagonal terms are given time-varying weights to account for changes in topology or sensitivity to events using a reduced system model. The results of this novel method are compared to the present errors of the previous publication in order to quantify the degree of improvement that this novel method imposes. The effective bus eigenvalues are briefly compared to Prony analysis to check similarities. An additional application for biorthogonal wavelets is also introduced to detect event types, including the HIF, for PMU data.
文摘In his series of three papers we study singularly perturbed (SP) boundary valueproblems for equations of elliptic and parabolic type. For small values of the pertur-bation parameter parabolic boundary and interior layers appear in these problems.If classical discretisation methods are used, the solution of the finite differencescheme and the approximation of the diffusive flux do not converge uniformly withrespect to this parameter. Using the method of special, adapted grids, we canconstruct difference schemes that allow approximation of the solution and the nor-malised diffusive flux uniformly with respect to the small parameter.We also consider sillgularly perturbed boundary value problems for convection-diffusion equations. Also for these problems we construct special finite differenceschemes, the solution of which converges ε-uniformly We study what problems ap-pear, when classical schemes are used for the approximation of the spatial deriva-tives. We compare the results with those obtained by the adapted approach. Re-sults of numerical experiments are discussed.In the three papers we first give an introduction on the general problem, andthen we consider respectively (i) Problems for SP parabolic equations, for whichthe solution and the normalised diffusive fluxes are required; (ii) Problems for SPelliptic equations with boundary conditions of Diriclilet, Neumann and RDbin type;(iii) Problems for SP parabolic equation with discontinuous boundary conditions-
基金supported by the Zhejiang Provincial Natural Science Foundation of China(No.LGF20C050001)the National Nature Science Foundation of China(No.61805213)。
文摘We use the method of discrete dipole approximation with surface interaction to construct a model in which a plurality of nanoparticles is arranged on the surface of BK7 glass. Nanoparticles are in air medium illuminated by evanescent wave generated from total internal reflection. The effects of the wavelength, the polarization of the incident wave, the number of nanoparticles and the spacing of multiple nanoparticles on the field enhancement and extinction efficiency are calculated by our model. Our work could pave the way to improve the field enhancement of multiple nanoparticles systems.
基金supported by the National Science Foundation under the grants NSF-DMS-1206438 and NSF-DHS-1510249,the National Science Foundation under the grants NSF-DMS-1004638 and NSF-DMS-1313272the Research Fund of Indiana University
文摘As a first step towards the numerical analysis of the stochastic primitive equations of the atmosphere and the oceans, the time discretization of these equations by an implicit Euler scheme is studied. From the deterministic point of view, the 3D primitive equations are studied in their full form on a general domain and with physically realistic boundary conditions. From the probabilistic viewpoint, this paper deals with a wide class of nonlinear, state dependent, white noise forcings which may be interpreted in either the Itor the Stratonovich sense. The proof of convergence of the Euler scheme,which is carried out within an abstract framework, covers the equations for the oceans, the atmosphere, the coupled oceanic-atmospheric system as well as other related geophysical equations. The authors obtain the existence of solutions which are weak in both the PDE and probabilistic sense, a result which is new by itself to the best of our knowledge.
文摘This study tested a novel method designed to provide useful information for medical diagnosis and treatment. We measured electroencephalography (EEG) during a test of eye opening and closing, a common test in routine EEG examination. This test is mainly used for measuring the degree of alpha blocking and sensitivity during eyes opening and closing. However, because these factors depend on the subject’s awareness, drowsiness can interfere with accurate diagnosis. We sought to determine the optimal EEG frequency band and optimal brain region for distinguishing healthy individuals from patients suffering from several neurophysiological diseases (including dementia, cerebrovascular disorder, schizophrenia, alcoholism, and epilepsy) while fully awake, and while in an early drowsy state. We tested four groups of subjects (awake healthy subjects, drowsy healthy subjects, awake patients and drowsy patients). The complexity of EEG band frequencies over five lobes in the human brain was analyzed using wavelet-based approximate entropy (ApEn). Two-way analysis of variance tested the effects of the two factors of interest (subjects’ health state, and subjects’ wakefulness state) on five different lobes of the brain during eyes opening and closing. The complexity of the theta and delta bands over frontal and central regions, respectively, was significantly greater in the healthy state during eyes opening. In contrast, patients exhibited increased complexity of gamma band activity over the temporal region only, during eyes-close. The early drowsy state and wakefulness state increased the complexity of theta band activity over the temporal region only during eyes-close and eyes-open states respectively, and this change was significantly greater in control subjects compared with patients. We propose that this method may be useful in routine EEG examination, to aid medical doctors and clinicians in distinguishing healthy individuals from patients, regardless of whether the subject is fully awake or in the early stages of drowsiness.
基金supported by the National Natural Science Foundation of China under Grant Nos.11901402 and 11671169。
文摘Every univariate Hermite interpolation problem can be written as a pointwise limit of Lagrange interpolants.However,this property is not preserved for the multivariate case.In this paper,the authors first generalize the result of P.Gniadek.As an application,the authors consider the discrete approximation problem for a special case when the interpolation condition contains all partial derivatives of order less than n and one nth order differential polynomial.In addition,for the case of n≥3,the authors use the concept of Cartesian tensors to give a sufficient condition to find a sequence of discrete points,such that the Lagrange interpolation problems at these points converge to the given Hermite-type interpolant.
基金National Natural Science Foundation of China(Grant Nos.61972194 and 62032011)。
文摘We consider the scattering of light in participating media composed of sparsely and randomly distributed discrete particles.The particle size is expected to range from the scale of the wavelength to several orders of magnitude greater,resulting in an appearance with distinct graininess as opposed to the smooth appearance of continuous media.One fundamental issue in the physically-based synthesis of such appearance is to determine the necessary optical properties in every local region.Since these properties vary spatially,we resort to geometrical optics approximation(GOA),a highly efficient alternative to rigorous Lorenz–Mie theory,to quantitatively represent the scattering of a single particle.This enables us to quickly compute bulk optical properties for any particle size distribution.We then use a practical Monte Carlo rendering solution to solve energy transfer in the discrete participating media.Our proposed framework is the first to simulate a wide range of discrete participating media with different levels of graininess,converging to the continuous media case as the particle concentration increases.
文摘In this series of three papers we study singularly perturbed (SP) boundaryvalue problems for equations of eiliptic and parabolic type- For small values ofthe perturbation parameter parabolic boundary and interior layers appear in theseproblems. If classical discretisation methods are used, the solution of the finitedifference scheme and the approximation of the diffusive flux do not converge uniformly with respect to this parameter. Using the method of special, edapted grids,we can construct difference schemes that allow apprcximation of the solution andthe normalised diffusive flux uniformly with respect to the small parameter.We also consider singularly perturbed boundary value problems for convection-diffusion equations. Also for these problems we construct special finite differenceschemes, the solution of which converges ε-uniformly. We study what problems appear, when classical schemes are used for the approximation of the spatial derivatives. We compare the results with those obtained by the adapted approach. Re-sults of numerical experiments are discussed.In the three papers we first give an introduction on the general problem, andthen we consider respectively (i) Problems for SP parabolic equations, for whichthe solution and the normalised diffusive fluxes are required; (ii) Problems for SPelliptic equations with boundary conditions of Dirichlet, Neumann and Robin type;(iii) Problems for SP parabolic equation with discontinuous boundary conditions.