Search

A local method is developed for estimating the effective conductivity for materials with varying deterministic or random conductivity on domains with mixed boundary conditions. The effective conductivity is the conductivity of a virtual homogeneous material which behaves globally as the original heterogeneous material. It is shown that the effective conductivity can be calculated from the values of the potential at a relatively small number of points in the specimen. A method is developed for calculating the potential at an arbitrary point in the specimen directly, rather than extracting its value from a field solution. The method is based on some technical concepts related to properties of diffusion processes and an extension of Itô’s formula for continuous semimartingales to the case of reflected diffusion processes. However, the application of the method is intuitive and uses elementary Monte Carlo simulation algorithms. The paper presents essential facts proving the validity of the proposed method and two numerical examples illustrating the calculation of the effective conductivity for specimens with deterministic and random conductivity fields. [ABSTRACT FROM AUTHOR]

This paper formulates a hybrid Monte Carlo implementation of the Fourier path integral (FPI-HMC) approach with partial averaging. Such a hybrid Monte Carlo approach allows one to generate collective moves through configuration space using molecular dynamics while retaining the computational advantages associated with the Fourier path integral Monte Carlo method. In comparison with the earlier Metropolis Monte Carlo implementations of the FPI algorithm, the present HMC method is shown to be significantly more efficient for quantum Lennard-Jones solids and suggests that such algorithms may prove useful for efficient simulations of a range of atomic and molecular systems. [ABSTRACT FROM AUTHOR]

This paper describes a new Monte Carlo method based on a novel stochastic potential switching algorithm. This algorithm enables the equilibrium properties of a system with potential V to be computed using a Monte Carlo simulation for a system with a possibly less complex stochastically altered potential V. By proper choices of the stochastic switching and transition probabilities, it is shown that detailed balance can be strictly maintained with respect to the original potential V. The validity of the method is illustrated with a simple one-dimensional example. The method is then generalized to multidimensional systems with any additive potential, providing a framework for the design of more efficient algorithms to simulate complex systems. A near-critical Lennard-Jones fluid with more than 20 000 particles is used to illustrate the method. The new algorithm produced a much smaller dynamic scaling exponent compared to the Metropolis method and improved sampling efficiency by over an order of magnitude. [ABSTRACT FROM AUTHOR]

The dynamical solution of a well-mixed, nonlinear stochastic chemical kinetic system, described by the Master equation, may be exactly computed using the stochastic simulation algorithm. However, because the computational cost scales with the number of reaction occurrences, systems with one or more “fast” reactions become costly to simulate. This paper describes a hybrid stochastic method that partitions the system into subsets of fast and slow reactions, approximates the fast reactions as a continuous Markov process, using a chemical Langevin equation, and accurately describes the slow dynamics using the integral form of the “Next Reaction” variant of the stochastic simulation algorithm. The key innovation of this method is its mechanism of efficiently monitoring the occurrences of slow, discrete events while simultaneously simulating the dynamics of a continuous, stochastic or deterministic process. In addition, by introducing an approximation in which multiple slow reactions may occur within a time step of the numerical integration of the chemical Langevin equation, the hybrid stochastic method performs much faster with only a marginal decrease in accuracy. Multiple examples, including a biological pulse generator and a large-scale system benchmark, are simulated using the exact and proposed hybrid methods as well as, for comparison, a previous hybrid stochastic method. Probability distributions of the solutions are compared and the weak errors of the first two moments are computed. In general, these hybrid methods may be applied to the simulation of the dynamics of a system described by stochastic differential, ordinary differential, and Master equations. © 2005 American Institute of Physics. [ABSTRACT FROM AUTHOR]

A framework, PhaseLift, was recently proposed to solve the phase retrieval problem. In this framework, the problem is solved by optimizing a cost function over the set of complex Hermitian positive semidefinite matrices. This approach to phase retrieval motivates a more general consideration of optimizing cost functions on semidefinite Hermitian matrices where the desired minimizers are known to have low rank. This paper considers an approach based on an alternative cost function defined on a union of appropriate manifolds. It is related to the original cost function in a manner that preserves the ability to find a global minimizer and is significantly more efficient computationally. A rank-based optimality condition for stationary points is given and optimization algorithms based on state-of-the-art Riemannian optimization and dynamically reducing rank are proposed. Empirical evaluations are performed using the PhaseLift problem. The new approach is shown to be an effective method of phase retrieval with computational efficiency increased substantially compared to a state-of-the-art algorithm, the Wirtinger flow algorithm. A preliminary version of this paper can be found in [W. Huang, K. A. Gallivan, and X. Zhang, Procedia Comput. Sci., 80 (2016), pp. 1125-1134]. [ABSTRACT FROM AUTHOR]

Modeling of a time-varying dynamical system provides insights into the functions of biological neural networks and contributes to the development of next-generation neural prostheses. In this paper, we have formulated a novel sparse multiwavelet-based generalized Laguerre-Volterra (sMGLV) modeling framework to identify the time-varying neural dynamics from multiple spike train data. First, the significant inputs are selected by using a group least absolute shrinkage and selection operator (LASSO) method, which can capture the sparsity within the neural system. Second, the multiwavelet-based basis function expansion scheme with an efficient forward orthogonal regression (FOR) algorithm aided by mutual information is utilized to rapidly capture the time-varying characteristics from the sparse model. Quantitative simulation results demonstrate that the proposed sMGLV model in this paper outperforms the initial full model and the state-of-the-art modeling methods in tracking performance for various time-varying kernels. Analyses of experimental data show that the proposed sMGLV model can capture the timing of transient changes accurately. The proposed framework will be useful to the study of how, when, and where information transmission processes across brain regions evolve in behavior. [ABSTRACT FROM AUTHOR]

An iterative QR-based soft feedback segment interference cancellation (QRSFSIC) detection and decoder algorithm for a Reed-Muller (RM) space-time turbo system is proposed in this paper. It forms the sufficient statistic for the minimum-mean-square error (MMSE) estimate according to QR decomposition-based soft feedback successive interference cancellation, stemmed from the a priori log-likelihood ratio (LLR) of encoded bits. Then, the signal originating from the symbols of the reliable segment, the symbol reliability metric, in terms of an a posteriori LLR of encoded bits which is larger than a certain threshold, is iteratively cancelled with the QRSFSIC in order to further obtain the residual signal for evaluating the symbols in the unreliable segment. This is done until the unreliable segment is empty, resulting in the extrinsic information for a RM turbo-coded bit with the greatest likelihood. Bridged by de-multiplexing and multiplexing, an iterative QRSFSIC detection is concatenated with an iterative trellis-based maximum a posteriori probability RM turbo decoder as if a principal Turbo detection and decoder is embedded with an iterative subordinate QRSFSIC detection and a RM turbo decoder, exchanging each other's detection and decoding soft-decision information iteratively. These three stages let the proposed algorithm approach the upper bound of the diversity. The simulation results also show that the proposed scheme outperforms the other suboptimum detectors considered in this paper. [ABSTRACT FROM AUTHOR]

In this paper, we propose a novel density-based fuzzy clustering algorithm based on Active Learning Method (ALM), which is a methodology of soft computing inspired by some hypotheses claiming that human brain interprets information in pattern-like images rather than numerical quantities. The proposed clustering algorithm, Fuzzy Unsupervised Active Learning Method (FUALM), is performed in two main phases. First, each data point spreads in the feature space just like an ink drop that spreads on a sheet of paper. As a result of this process, densely connected ink patterns are formed that represent clusters. In the second phase, a fuzzifying process is applied in order to summarize the effects of all members of each cluster. Finding arbitrary shaped clusters, noise robustness and proposing fuzzy clusters are some of the advantages of our proposed clustering algorithm. The algorithm is described in full details and its performance is evaluated and compared with well-known clustering algorithms on synthetic and real-world datasets. [ABSTRACT FROM AUTHOR]

Continuous top-$k$ query over streaming data is a fundamental problem in database. In this paper, we focus on the sliding window scenario, where a continuous top- $k$ query returns the top- $k$ objects within each query window on the data stream. Existing algorithms support this type of queries via incrementally maintaining a subset of objects in the window and try to retrieve the answer from this subset as much as possible whenever the window slides. However, since all the existing algorithms are sensitive to query parameters and data distribution, they all suffer from expensive incremental maintenance cost. In this paper, we propose a self-adaptive partition framework to support continuous top-$k$ query. It partitions the window into sub-windows and only maintains a small number of candidates with highest scores in each sub-window. Based on this framework, we have developed several partition algorithms to cater for different object distributions and query parameters. To our best knowledge, it is the first algorithm that achieves logarithmic complexity w.r.t. $k$ for incrementally maintaining the candidate set even in the worst case scenarios. [ABSTRACT FROM PUBLISHER]

This paper is devoted to the analysis of the sixth-order symplectic and symmetric explicit extended Runge-Kutta-Nyström (ERKN) schemes for solving multi-frequency oscillatory nonlinear Hamiltonian equations. Fourteen practical sixth-order symplectic and symmetric explicit ERKN schemes are constructed, and their phase properties are investigated. The paper is accompanied by five numerical experiments, including a nonlinear two-dimensional wave equation. The numerical results in comparison with the sixth-order symplectic and symmetric Runge-Kutta-Nyström methods and a Gautschi-type method demonstrate the efficiency and robustness of the new explicit schemes for solving multi-frequency oscillatory nonlinear Hamiltonian equations. [ABSTRACT FROM AUTHOR]

In those oceans where measured wave data are not available, numerical wave models are usually adopted to hindcast wave parameters in order to define design waves for marine structures. To utilize these hindcating data, it is very important to perform error corrections of model results for accurate estimation of the appropriate wave parameters. In this paper, a new non-parametric correction model is established to improve wave model accuracy through modifying a previous approach released by Caires and Sterl in 2005. The new correction model introduces a kernel algorithm to learn error information from both value magnitude and series trend through training datasets, and utilizes the information to correct potential errors in model outputs. It is shown that the two-dimensional learning method is more effective than the previous one-dimensional which only learns error information from the value magnitude. Furthermore, an error constraint parameter is initially adopted in the new correction model to decrease the possibility of overcorrection. The new correction model performs better than its predecessor, especially when modeling wave period and altimeter synchronized wave height. Though this paper evaluates the model correcting performance with WAVEWATCH III outputs, the modified model can be adopted to correct other kinds of time-series data. [ABSTRACT FROM AUTHOR]

Hydraulic fractures are tensile cracks that propagate in pre-stressed solid media due to the injection of a viscous fluid. Developing numerical schemes to model the propagation of these fractures is particularly challenging due to the degenerate, hypersingular nature of the coupled integropartial differential equations. These equations typically involve a singular free boundary whose velocity can only be determined by evaluating a distinguished limit. This review paper describes a class of numerical schemes that have been developed to use the multiscale asymptotic behaviour typically encountered near the fracture boundary as multiple physical processes compete to determine the evolution of the fracture. The fundamental concepts of locating the free boundary using the tip asymptotics and imposing the tip asymptotic behaviour in a weak form are illustrated in two quite different formulations of the governing equations. These formulations are the displacement discontinuity boundary integral method and the extended finite-element method. Practical issues are also discussed, including new models for proppant transport able to capture 'tip screen-out'; efficient numerical schemes to solve the coupled nonlinear equations; and fast methods to solve resulting linear systems. Numerical examples are provided to illustrate the performance of the numerical schemes. We conclude the paper with open questions for further research. [ABSTRACT FROM AUTHOR]

In computer aided geometric design a polynomial is usually represented in Bernstein form. This paper presents a family of compensated algorithms to accurately evaluate a polynomial in Bernstein form with floating point coefficients. The principle is to apply error-free transformations to improve the traditional de Casteljau algorithm. At each stage of computation, round-off error is passed on to first order errors, then to second order errors, and so on. After the computation has been "filtered" (K − 1) times via this process, the resulting output is as accurate as the de Casteljau algorithm performed in K times the working precision. Forward error analysis and numerical experiments illustrate the accuracy of this family of algorithms. [ABSTRACT FROM AUTHOR]

With the increasing complexity of engineering systems, reliability analysis and evaluation of systems with traditional methods can't meet practical engineering requirements. Based on limited experimental conditions, lack of data, complex structure models, insufficient cognitive abilities, and many other issues, people have to consider many uncertain factors in system reliability research. Besides, common cause failure (CCF) has become an important factor of system failure. In this paper, a discrete‐time Bayesian network (DTBN) associated with an eight‐rotor unmanned aerial vehicle (UAV) system is presented to discuss above problems. In this approach, the system is assumed as a two‐state system. After that, interval analysis theory is employed to deal with uncertainty. We consider the four sets of auxiliary propellers in the auxiliary power group as a 3/8 voting system, and β factor model is used to process CCF in the auxiliary power group. The proposed methods prove the validity of proposing interval analysis theory to solve uncertain problems and it is necessary to consider reducing or avoiding CCFs in system. [ABSTRACT FROM AUTHOR]

Tribological tests are adopted to estimate the kinetic friction coefficient of the material of interest according to the ASTM standard. Typically, for the measurement process, several replicas as well as a postprocessing data treatment are necessary to take into account the observed casual variability of the measurand. This paper describes a novel statistical approach aiming to highlight the most significant time intervals during each test replica. The experimental activity has been carried out with reference to the ball-on-flat configuration of the tribometer in order to verify the suggested procedure in terms of both accuracy and reliability. Numerical simulation and comparison of the measurement results with the output of a different outlier detection algorithm are also included to show the improvement allowed by the authors’ proposal about the dry friction coefficient measurement. [ABSTRACT FROM AUTHOR]

Highlights • The new algorithm is the second-order accurate, unconditionally stable (L-stable) and self-starting. • The new algorithm shares the identical effective stiffness matrices inside two sub-steps. • The new method does not involve any artificial parameters and additional variable, such as the Lagrange multipliers. • The new scheme achieves the same numerical properties as the Bathe algorithm, but requires less matrix-vector operations. Abstract This paper presents a new composite sub-steps algorithm for solving reliable numerical responses in structural dynamics. The newly developed algorithm is a two sub-steps, second-order accurate and unconditionally stable implicit algorithm with the same numerical properties as the Bathe algorithm. The detailed analysis of the stability and numerical accuracy is presented for the new algorithm, which shows that its numerical characteristics are identical to those of the Bathe algorithm. Hence, the new sub-steps scheme could be considered as an alternative to the Bathe algorithm. Meanwhile, the new algorithm possesses the following properties: (a) it produces the same accurate solutions as the Bathe algorithm for solving linear and nonlinear problems; (b) it does not involve any artificial parameters and additional variables, such as the Lagrange multipliers; (c) The identical effective stiffness matrices can be obtained inside two sub-steps; (d) it is a self-starting algorithm. Some numerical experiments are given to show the superiority of the new algorithm and the Bathe algorithm over the dissipative CH- α algorithm and the non-dissipative trapezoidal rule. [ABSTRACT FROM AUTHOR]

Fractional interpolation is used to provide sub-pixel level references for motion compensation in the interprediction of video coding, which attempts to remove temporal redundancy in video sequences. Traditional handcrafted fractional interpolation filters face the challenge of modeling discontinuous regions in videos, while existing deep learning-based methods are either designed for a single quantization parameter (QP), only generating half-pixel samples, or need to train a model for each sub-pixel position. In this paper, we present a one-for-all fractional interpolation method based on a grouped variation convolutional neural network (GVCNN). Our method can deal with video frames coded using different QPs and is capable of generating all sub-pixel positions at one sub-pixel level. Also, by predicting variations between integer-position pixels and sub-pixels, our network offers more expressive power. Moreover, we perform specific measurements in training data generation to simulate practical situations in video coding, including blurring the down-sampled sub-pixel samples to avoid aliasing effects and coding integer pixels to simulate reconstruction errors. In addition, we analyze the impact of the size of blur kernels theoretically. Experimental results verify the efficiency of GVCNN. Compared with HEVC, our method achieves 2.2% in bit saving on average and up to 5.2% under low-delay P configuration. [ABSTRACT FROM AUTHOR]

In this paper, we establish a framework for enabling ultra-reliable and low-latency communications in the control and non-payload communications (CNPC) links of the unmanned aerial vehicle (UAV) communication systems. We first derive the available range of the CNPC links between UAVs and a ground control station. The available range is defined as the maximal horizontal communication distance within which the round-trip delay and the overall packet loss probability can be ensured with a required probability. To exploit the macro-diversity gain of the distributed multi-antenna systems (DAS) and the array gain of the centralized multi-antenna systems (CAS), we consider a modified DAS (M-DAS), where the ground control station is equipped with the distributed access points (APs), and each AP can have multiple antennas. We then show that the available range can be maximized by judiciously optimizing the altitude of UAVs, the duration of the uplink and downlink phases, and the antenna configuration. To solve the non-convex problem, we propose an algorithm that can converge to the optimal solution in DAS and CAS, and then extend it into more general M-DAS. The simulation and numerical results validate our analysis and show that the available range of M-DAS can be significantly larger than those of the DAS and CAS. [ABSTRACT FROM AUTHOR]

2-D linear discriminant analysis (2DLDA) has been widely used in pattern recognition and image classification. 2DLDA selects discriminative features from the up and left corner of images. However, 2DLDA uses the Frobenius norm (F-norm), which is sensitive to noise or outliers in data, as a metric. In this paper, we propose a novel framework, called horizontal and vertical nuclear norm-based 2DLDA (HVNN-2DLDA) for image representation. In the proposed framework, HVNN-2DLDA methods (i.e., HNN-2DLDA and VNN-2DLDA) are proposed, and both use the nuclear norm as a criterion. The nuclear norm can provide more structure and global information for the reconstruction of noisy images. HNN-2DLDA and VNN-2DLDA represent images in the row and column directions, respectively. In addition, by combining the row and column directions, we propose a bilateral nuclear norm-based 2DLDA method called BNN-2DLDA. The advantage of BNN-2DLDA over HNN-2DLDA and VNN-2DLDA is that an image sample can be represented by both the row and the column directions instead of only the row or column direction. HVNN-2DLDA learns a set of local optimal projection vectors by maximizing the ratio of the nuclear norm of the between-class scatter matrix and the nuclear norm of the within-class scatter matrix. To verify the robustness and recognition performance in image classification of HVNN-2DLDA, six public image databases are used for experiments. The experimental results demonstrate the effectiveness and the feasibility of the proposed framework. [ABSTRACT FROM AUTHOR]

In this paper we consider the branched transportation problem in two dimensions associated with a cost per unit length of the form 1 + β ⁢ θ {1+\beta\,\theta} , where θ denotes the amount of transported mass and β > 0 {\beta>0} is a fixed parameter (notice that the limit case β = 0 {\beta=0} corresponds to the classical Steiner problem). Motivated by the numerical approximation of this problem, we introduce a family of functionals ( { ℱ ε } ε > 0 {\{\mathcal{F}_{\varepsilon}\}_{\varepsilon>0}}) which approximate the above branched transport energy. We justify rigorously the approximation by establishing the equicoercivity and the Γ-convergence of { ℱ ε } {\{\mathcal{F}_{\varepsilon}\}} as ε ↓ 0 {\varepsilon\downarrow 0}. Our functionals are modeled on the Ambrosio–Tortorelli functional and are easy to optimize in practice. We present numerical evidences of the efficiency of the method. [ABSTRACT FROM AUTHOR]

Effective monitoring and management applications on modern distribution networks (DNs) require a sound network model and the knowledge of line parameters. Network line impedances are used, among other things, for state estimation and protection relay setting. Phasor measurement units (PMUs) give synchronized voltage and current phasor measurements, referred to a common time reference (coordinated universal time). All synchrophasor measurements can thus be temporally aligned and coordinated across the network. This feature, along with high accuracy and reporting rates, could make PMUs useful for the evaluation of network parameters. However, instrument transformer behavior strongly affects the parameter estimation accuracy. In this paper, a new PMU-based iterative line parameter estimation algorithm for DNs, which includes in the estimation model systematic measurement errors, is presented. This method exploits the simultaneous measurements given by PMUs on different nodes and branches of the network. A complete analysis of uncertainty sources is also performed, allowing the evaluation of estimation uncertainty. Issues related to operating conditions, topology, and measurement uncertainty are thoroughly discussed and referenced to a realistic model of a DN to show how a full network estimator is possible. [ABSTRACT FROM AUTHOR]

Hyperspectral unmixing has attracted much attention in recent years. Single sparse unmixing assumes that a pixel in a hyperspectral image consists of a relatively small number of spectral signatures from large, ever-growing, and available spectral libraries. Joint-sparsity (or row-sparsity) model typically enforces all pixels in a neighborhood to share the same set of spectral signatures. The two sparse models are widely used in the literature. In this paper, we propose a joint-sparsity-blocks model for abundance estimation problem. Namely, the abundance matrix of size $m\times n$ is partitioned to have one row block and $s$ column blocks and each column block itself is joint-sparse. It generalizes both the single (i.e., $s=n$) and the joint (i.e., $s=1$) sparsities. Moreover, concatenating the proposed joint-sparsity-blocks structure and low rankness assumption on the abundance coefficients, we develop a new algorithm called joint-sparse-blocks and low-rank unmixing. In particular, for the joint-sparse-blocks regression problem, we develop a two-level reweighting strategy to enhance the sparsity along the rows within each block. Simulated and real-data experiments demonstrate the effectiveness of the proposed algorithm. [ABSTRACT FROM AUTHOR]

In this paper, we propose two groups of recursive codes, namely the Hadamard recursive carrier interferometry (HRCI) codes and diagonal recursive carrier interferometry (DRCI) codes, to simultaneously reduce the peak-to-average power ratio (PAPR) and suppress the phase noises and residual carrier frequency offset (RCFO) for high-speed orthogonal frequency division multiplexing (OFDM) systems. We exploit the recursion property of the Hadamard and diagonal matrices to constitute the proposed HRCI code and DRCI code using the carrier interferometry (CI) code, and present the system model in details to spread OFDM symbols in the frequency domain. Then, we prove that both HRCI and DRCI generation matrices keep their invertibility during the recursions. As a direct result of the new code design, the phase intervals are enlarged and the signals are shifted in the time domain, thus both the PAPR and the phase noises are reduced, and the RCFO can be mitigated. Moreover, the enlarged phase intervals enable the transceiver blocks including the channel estimator, which estimate the phase noise based on the received signals, to improve the performance. Furthermore, we present the expressions for the HRCI and DRCI embedded transmitted and received signals, and derive the analytical signal-to-interference-plus-noise ratio (SINR) expressions for bit error rate (BER) performance analysis. Then, we analyze the effects of presented codes on the phase differences and provide the computational complexity analysis. Numerical simulations verify the effectiveness of our theoretical analysis of SINR. Additionally, under different parameter settings, which consider the working conditions of practical systems, we demonstrate that the presented robust HRCI and DRCI coded OFDM systems can suppress the PAPR satisfactorily and better BER performances can be achieved by the recursive CI-aided spread OFDM systems with phase noise rejection and RCFO compensation when compared with counterpart systems. [ABSTRACT FROM AUTHOR]

Abstract Distributions of Monge type are a class of strongly regular bracket-generating distributions introduced by I. Anderson, Zh. Nie and P. Nurowski. Their symbol algebras prolong to simple graded Lie algebras, thus allowing one to associate a parabolic geometry to any given Monge distribution. This article is devoted to the classification problem for homogeneous models of Monge distributions of type C 3 in dimension eight, and is complementary to a paper by I. Anderson and P. Nurowski. The general classification algorithm, as well as most of its application to the particular problem, are joint work with Ian Anderson. [ABSTRACT FROM AUTHOR]

A new Interpolation based Reversible Data Hiding (IRDH) scheme is reported in this paper. For different applications of an IRDH scheme to the digital image, video, multimedia, big-data and biological data, the embedding capacity requirement usually varies. Disregarding this important consideration, existing IRDH schemes do not offer a better embedding rate-distortion performance for varying size payloads. To attain this varying capacity requirement with our proposed adaptive embedding, we formulate a capacity control parameter and propose to utilize it to determine a minimum set of embeddable bits in a pixel. Additionally, we use a logical (or bit-wise) correlation between the embeddable pixel and estimated versions of an embedded pixel. Thereby, while a higher range between an upper and lower limit of the embedding capacity is maintained, a given capacity requirement within that limit is also attained with a better-embedded image quality. Computational modeling of all new processes of the scheme is presented, and performance of the scheme is evaluated with a set of popular test-images. Experimental results of our proposed scheme compared to the prominent IRDH schemes have recorded a significantly better-embedding rate-distortion performance. [ABSTRACT FROM AUTHOR]

This paper proposes a novel algorithm to determine the path flow solution that satisfies the proportionality condition from a given origin-based user equilibrium (UE) link flow solution. The algorithm iteratively refines the origin-based bush representation of the UE solution through constructing a list of "node-based arrival bushes" and solving the entropy maximization subproblems defined on them. Thanks to the special structure of bushes, these subproblems can be solved efficiently. The proposed algorithm thus obviates the known hurdles to efficient implementation, such as enumerating all UE paths or collecting a set of paired alternative segments to cover them. We prove that the algorithm ensures convergence to a solution that perfectly satisfies the proportionality condition in general networks. Numerical experiments indicate that the proposed algorithm converges much faster than the known alternative, with a three- to eightfold speed-up on large networks. The online appendices are available at https://doi.org/10.1287/trsc.2018.0845. [ABSTRACT FROM AUTHOR]

Abstract Future Mars missions may require improved landed accuracy to facilitate the landing site selection and finally access a region of complex terrain with high scientific return. This paper is to develop a novel, robust, and precision entry guidance algorithm for Mars entry vehicles with low lift-to-drag ratios. In the presence of large uncertainties, the entry terminal point controller algorithm can encounter severe performance degradation due to: (1) the small perturbation assumption, and (2) theoretically ignoring deviations in the atmospheric-density model, aerodynamic-force model, etc. Based on numerical technologies and the classical variation method (VM), this work develops a semi-analytical (SA) algorithm, in which combined effects of several dynamic uncertainties now can be addressed. The terminal-downrange error is predicted by a numerical predictor such that the dependence on the reference trajectory can be reduced and then the issue caused by the small perturbation assumption can be addressed. Such a predicted terminal-downrange error is finally corrected by an analytical corrector, which is designed by the VM. It is indicative that there can be no numerical iterations in the SA algorithm. Numerical results demonstrate the effectiveness of the SA algorithm. [ABSTRACT FROM AUTHOR]

Highlights • Dynamical analysis has been provided to describe nonlinear behavior of supercavitating vehicles. • Robust H∞ loop-shaping synthesis with modified PID algorithm is proposed to control the dive plane maneuver of the HSSV. • Multi-objective control problems are solved using BMI optimization of an equivalent Schur formula. • Control scheme has the low order structure and provides robustness against uncertainties. • Integrated robust controller can deal with high parametric uncertainties and suppress exogenous disturbances and sensor noises. Abstract The high-speed supercavitating vehicle (HSSV) utilizes advanced technology that enables an underwater vehicle to reach its unprecedented high speed. The vertical motion control of the HSSV is challenging problem because of its complex dynamics with nonlinear planing force, parametric uncertainties, external disturbances, actuator saturation, and sensor noises. This paper deals with dynamical analysis and a robust H∞ loop-shaping synthesis with modified PID (proportional-integral-derivative) algorithm to control the dive plane maneuver of the HSSV. Typically, the control scheme has the low order structure and provides robustness against dynamic uncertainties, which can be implemented using the bilinear matrix inequality (BMI) optimization of an equivalent Schur formula. Simulation results show that the controlled vehicle system provides good performance and high robustness against uncertainties, ensuring no-overshoot and fast in time-domain responses. In addition, the control algorithm can decouple the dynamic interactions in the multi-input multi-output (MIMO) system, overcoming parametric uncertainty, planing force, and actuator saturation while minimizing the effect of the strong external disturbances and measurement noises. [ABSTRACT FROM AUTHOR]

In this paper, we propose a lattice-based robust distributed source coding system for two correlated sources and provide a detailed performance analysis under the high resolution assumption. It is shown, among other things, that, in the asymptotic regime where: 1) the side distortion approaches 0 and 2) the ratio between the central and side distortions approaches 0, our scheme is capable of achieving the information-theoretic limit of quadratic multiple description coding when the two sources are identical, whereas a variant of the random coding scheme by Chen and Berger with Gaussian codes has a performance loss of 0.5 bits relative to this limit. [ABSTRACT FROM AUTHOR]

Algorithms for secure encryption in a post-quantum world are currently receiving a lot of attention in the research community. One of the most promising such algorithms is the code-based scheme called QC-MDPC, which has excellent performance and a small public key size. In this paper, we present a very efficient key recovery attack on the QC-MDPC scheme using the fact that decryption uses an iterative decoding step, and this can fail with some small probability. We identify a dependence between the secret key and the failure in decoding. This can be used to build what we refer to as a distance spectrum for the secret key, which is the set of all distances between any two ones in the secret key. In a reconstruction step, we then determine the secret key from the distance spectrum. The attack has been implemented and tested on a proposed instance of QC-MDPC for 80-bit security. It successfully recovers the secret key in minutes. A slightly modified version of the attack can be applied on proposed versions of the QC-MDPC scheme that provides IND-CCA security. The attack is a bit more complex in this case, but still very much below the security level. The reason why we can break schemes with proved CCA security is that the model for these proofs typically does not include the decoding error possibility. At last, we present several algorithms for key reconstruction from an empirical distance spectrum. We first improve the naïve algorithm for key reconstruction by a factor of about 3 0000, when the parameters for 80-bit security are implemented. We further develop the algorithm to deal with errors in the distance spectrum. This ultimately reduces the requirement on the number of ciphertexts that need to be collected for a successful key recovery. [ABSTRACT FROM AUTHOR]

Dynamic robust principal component analysis (PCA) refers to the dynamic (time-varying) extension of robust PCA (RPCA). It assumes that the true (uncorrupted) data lie in a low-dimensional subspace that can change with time, albeit slowly. The goal is to track this changing subspace over time in the presence of sparse outliers. We develop and study a novel algorithm, which we call simple-ReProCS, based on the recently introduced Recursive Projected Compressive Sensing (ReProCS) framework. This paper provides the first guarantee for dynamic RPCA that holds under weakened versions of standard RPCA assumptions, slow subspace change, and a lower bound assumption on most outlier magnitudes. Our result is significant because: 1) it removes the strong assumptions needed by the two previous complete guarantees for ReProCS-based algorithms; 2) it shows that it is possible to achieve significantly improved outlier tolerance, compared with all existing RPCA or dynamic RPCA solutions by exploiting the above-mentioned two simple extra assumptions; and 3) it proves that simple-ReProCS is online (after initialization), fast, and, has near-optimal memory complexity. [ABSTRACT FROM AUTHOR]

Many statistical estimation procedures lead to nonconvex optimization problems. Algorithms to solve these problems are often guaranteed to output a stationary point of the optimization problem. Oracle inequalities are an important theoretical instrument to assess the statistical performance of an estimator. Oracle results have focused on the theoretical properties of the uncomputable (global) minimum or maximum. In this paper, a general framework used for convex optimization problems to derive oracle inequalities for stationary points is extended. A main new ingredient of these oracle inequalities is that they are sharp: they show closeness to the best approximation within the model plus a remainder term. We apply this framework to different estimation problems. [ABSTRACT FROM AUTHOR]

Intrusion detection systems are devoted to monitor a network with aims at finding and avoiding anomalous events. In particular, we focus on misuse detection systems, which are trained to identify several known types of attacks. These can be unauthorized accesses, or denial of service attacks, among others. Whenever it scans a trace of a suspicious event, it is programmed to trigger an alert and/or to block this dangerous access to the system. Depending on the security policies of the network, the administrator may seek different requirements that will have a strong dependency on the behavior of the intrusion detection system. For a given application, the cost of raising false alarms could be higher than carrying out a preventive access lock. In other scenarios, there could be a necessity of correctly identifying the exact type of cyber attack to proceed in a given way. In this paper, we propose a multi-objective evolutionary fuzzy system for the development of a system that can be trained using different metrics. By increasing the search space during the optimization of the model, more accurate solutions are expected to be obtained. Additionally, this scheme allows the final user to decide, among a broad set of solutions, which one is better suited for the current network characteristics. Our experimental results, using the well-known KDDCup'99 problem, supports the quality of this novel approach in contrast to the state-of-the-art for evolutionary fuzzy systems in intrusion detection, as well as the C4.5 decision tree. [ABSTRACT FROM AUTHOR]

Document representation is a key problem in document analysis and processing tasks, such as document classification, clustering and information retrieval. Especially for unstructured text data, the use of a suitable document representation method would affect the performance of the subsequent algorithms for applications and research. In this paper, we propose a novel document representation method called the conditional co-occurrence degree matrix document representation method (CCODM), which is based on word co-occurrence. CCODM not only considers the co-occurrence of terms but also considers the conditional dependencies of terms in a specific context, which leads to more available and useful structural and semantic information being retained from the original documents. Extensive experimental classification results with different supervised and unsupervised feature selection methods show that the proposed method, CCODM, achieves better performance than the vector space model, latent Dirichlet allocation, the general co-occurrence matrix representation method and the document embedding method. [ABSTRACT FROM AUTHOR]

Differential evolution (DE) is one of the most powerful and versatile evolutionary algorithms for efficiently solving complex real-world optimization problems in recent years. Since its introduction in 1995, the research focus in DE has mostly been on the variant side with so many new algorithms proposed based on the original DE algorithm. However, each new algorithm is only suitable for certain fitness landscapes, and, therefore, some types of optimization problems cannot be solved efficiently. To tackle this issue, this paper presents a new self-feedback DE algorithm, named the SFDE; its optimal variation strategy is selected by extracting the local fitness landscape characteristics in each generation population and combing the probability distributions of unimodality and multimodality in each local fitness landscape. The proposed algorithm is tested on a suite of 17 benchmark functions, and the experimental results demonstrated its advantages in a high search dimension in that it can ensure that the population moves to a better fitness landscape, then speeds up convergence to the global optimum, and avoids falling into local optima. [ABSTRACT FROM AUTHOR]

The use of the method of moments to solve surface integral equations is one of the most popular numerical techniques in electromagnetic modeling and analysis. This method requires the accurate and efficient numerical evaluation of iterated surface integrals over both source and testing domains. In this paper, we propose a scheme for evaluating these 4-D interaction integrals between pairs of arbitrarily positioned and oriented elements. The approach is based on applying the surface divergence theorem twice, once on the source and once on the test domain. When the integrations are reordered as two outer contour integrals plus two inner radial integrals, the initial radial integrations provide significant smoothing of the underlying singular integrands. The method is numerically validated for static and dynamic kernels arising in the electric field integral equation, i.e., for kernels with $1/R$ singularities, and linear basis functions. The proposed formula to evaluate 4-D reaction integrals can be extended to different kernels and to different elements, e.g., to curved or volumetric elements, and to basis functions of higher order. [ABSTRACT FROM AUTHOR]

Direct geolocation of radio emitters has traditionally been formulated using short measurement intervals over which delay and Doppler frequency shift are approximated as constant. In this paper, we consider maximum likelihood position estimation using long measurement intervals for which the Doppler is time varying. A series of assumptions and approximations are used to develop an efficient solution involving coherent combining of traditional, complex ambiguity function processing outputs. Simulation results are used to demonstrate potential gains. [ABSTRACT FROM AUTHOR]

This paper is devoted to the terrain aided navigation (TAN) systems based on the state estimation algorithms. In particular, the emphasis is laid on the design of the Bayesian Rao–Blackwellized point-mass filter for nonlinear state-space models of a specific structure typically used in the terrain aided navigation. The proposed filter preserves advantages of the point-mass filter, including the high estimation accuracy, robust initialization, deterministic nature of the algorithm, and predictable computational complexity, while the computational complexity is significantly reduced. The performance of the Rao–Blackwellized point-mass filter is illustrated using a simulated numerical example and a TAN system based on the real data. [ABSTRACT FROM AUTHOR]

The emerging of depth-camera technology is paving the way for a variety of new applications and it is believed that plane detection is one of them. In fact, planes are common in man-made living structures, thus their accurate detection can benefit many visual-based applications. The use of depth information allows detecting planes characterized by complex pattern and texture, where the texture-based plane detection algorithms usually fail. In this paper, we propose a robust depth-driven plane detection (DPD) algorithm which consists of two parts: the growing-based plane detection and a two-stage refinement. The proposed approach starts from the seed patch with the highest planarity and uses the estimated equation of the growing plane and a dynamic threshold function to steer the growing process. Aided with this mechanism, each seed patch can grow to its maximum extent, and then the next seed patch starts to grow. This process is iteratively repeated so as to detect all the planes. Moreover, the refinement is proposed to tackle two common problems suffered by growing-based approaches, the over-growing problem, and the under-growing problem. Validated by extensive experiments, the proposed DPD algorithm is able to accurately detect planes and robust to various testing conditions. In terms of applications, it can be used as the pre-processing step for a variety of applications, such as, planar object recognition, super-resolution of the time-of-flight depth images with intrinsically low resolution. [ABSTRACT FROM AUTHOR]

Frame rate up-conversion (FRUC) is a technology for generating a higher frame rate video from one with a lower frame rate. The multi-view plus depth map (MVD) system has been recently considered a promising solution for three-dimensional video systems. Accordingly, an FRUC efficient solution for the MVD system is greatly needed. In this paper, we thus propose an FRUC approach that considers the MVD system characteristics. It divides a frame into multiple layers using a depth map, and it reconstructs occluded regions using neighboring views and depth map information. With reconstructed information, it accurately performs motion estimation when a region occlusion occurs. Experimental results demonstrated that the proposed method significantly improved the quality of interpolated frames, and it can effectively complement other conventional FRUC approaches. [ABSTRACT FROM AUTHOR]

This paper formally constructs navigation functions with time-varying destinations on star worlds. The construction is based on appropriate diffeomorphic transformations and extends an earlier sphere-world formulation. A new obstacle modeling method is also introduced, reducing analytical complexity and offering unified expressions of common classes of $n$ -dimensional obstacles. The method allows for dynamic target tracking and is validated through simulations and experiments. [ABSTRACT FROM AUTHOR]

An alternative approach to standard computed torque with feedback linearization is proposed in this paper to control cable-driven parallel robots (CDPRs) with highly flexible cables. Exteroceptive feedback is used to measure the end-effector Cartesian position at a high sampling rate. Stability is demonstrated using singular perturbation theory. The proposed control scheme is experimentally validated on a planar 3-degree-of-freedom CDPR and its efficiency is assessed by comparison with a simple kinematic control law. [ABSTRACT FROM AUTHOR]

This paper is concerned with mathematical and numerical analysis of the system of radiative integral transfer equations. The existence and uniqueness of solution to the integral system is proved by establishing the boundedness of the radiative integral operators and proving the invertibility of the operator matrix associated with the system. A collocation-boundary element method is developed to discretize the differential-integral system. For the non-convex geometries, an element-subdivision algorithm is developed to handle the computation of the integrals containing the visibility factor. An efficient iterative algorithm is proposed to solve the nonlinear discrete system and its convergence is also established. Numerical experiment results are also presented to verify the effectiveness and accuracy of the proposed method and algorithm. [ABSTRACT FROM AUTHOR]

• An automatic algorithm for generating polyhedral meshes was systematically developed. • Hexahedral elements near the given STL surface were categorized into polyhedral volume topologies based on the MC algorithm. • Topology ambiguities of the MC algorithm were effectively resolved to provide an appropriate mesh for complex geometry. • The effectiveness was demonstrated through numerical analysis using the smoothed finite element method. Hexahedral meshes that exhibit the superiority in terms of solution accuracy and convergence rate are preferred to other types of meshes in the finite element analysis. However, the construction of the hexahedral meshes for complex geometries is still considered troublesome due to their poor geometric adaptability. This paper presents an efficient grid-based scheme to automatically generate polyhedral meshes including the hexahedral elements, and thus to provide hexahedral-dominant meshes for three-dimensional geometry with complex shapes. On the basis of the marching cube algorithm with a background grid composed of a regular arrangement of cubes, surface topologies for the background cubes are defined to represent the three-dimensional boundaries of a given domain. Then, in order to generate a three-dimensional finite element mesh, the surface topologies of the marching cube algorithm are systematically expanded to polyhedral volume topologies. Meanwhile, a topology ambiguity problem inherent in the marching cube algorithm is effectively resolved to generate an appropriate polyhedral mesh even for an arbitrary complex geometry. Several examples including biostructure modeling demonstrate that the proposed mesh generation scheme can easily discretize complex three-dimensional domains with hexahedral-dominant meshes, which are composed of the polyhedral elements near the domain boundaries and the hexahedral elements that come from the background cubes inside the domains. Furthermore, to show the applicability and effectiveness of polyhedral meshes in the finite element analysis, some structural analyses are performed using the smoothed finite element method that can be straightforwardly adapted to polyhedral elements of arbitrary shape. [ABSTRACT FROM AUTHOR]