In this article we use high-throughput epigenomics, transcriptomics, and proteomics data to construct fine-graded models of the "protein-coding units"gathering all transcript isoforms and chromatin accessibility peaks associated with more than 4000 genes in humans. Each protein-coding unit has the structure of a directed acyclic graph (DAG) and can be represented as a Bayesian network. The factorization of the joint probability distribution induced by the DAGs imposes a number of conditional independence relationships among the variables forming a protein-coding unit, corresponding to the missing edges in the DAGs. We show that a large fraction of these conditional independencies are indeed verified by the data. Factors driving this verification appear to be the structural and functional annotation of the transcript isoforms, as well as a notion of structural balance (or frustration-free) of the corresponding sample correlation graph, which naturally leads to reduction of correlation (and hence to independence) upon conditioning.

GNSS receivers are vulnerable to spoofing attacks in which false satellite signals deceive receivers to compute false position and/or time estimates. This work derives and evaluates algorithms that perform spoofing mitigation by utilizing double differences of pseudorange or carrier phase measurements from multiple receivers. The algorithms identify pseudorange and carrier-phase measurements originating from spoofing signals, and omit these from the position and time computation. The algorithms are evaluated with simulated and live-sky meaconing attacks. The simulated spoofing attacks show that mitigation using pseudoranges is possible in these tests when the receivers are separated by five meters or more. At 20 meters, the pseudorange algorithm correctly authenticates six out of seven pseudoranges within 30 seconds in the same simulator tests. Using carrier phase allows mitigation with shorter distances between receivers, but requires better time synchronization between the receivers. Evaluations with live-sky meaconing attacks show the validity of the proposed mitigation algorithms.

We propose a linear programming method that is based on active-set changes and proximal-point iterations. The method solves a sequence of least-distance problems using a warm-started quadratic programming solver that can reuse internal matrix factorizations from the previously solved least-distance problem. We show that the proposed method terminates in a finite number of iterations and that it outperforms state-of-the-art LP solvers in scenarios where an extensive number of small/medium scale LPs need to be solved rapidly, occurring in, for example, multi-parametric programming algorithms. In particular, we show how the proposed method can accelerate operations such as redundancy removal, computation of Chebyshev centers and solving linear feasibility problems.

Motivation: The simultaneous availability of ATAC-seq and RNA-seq experiments allows to obtain a more in-depth knowledge on the regulatory mechanisms occurring in gene regulatory networks. In this article, we highlight and analyze two novel aspects that leverage on the possibility of pairing RNA-seq and ATAC-seq data. Namely we investigate the causality of the relationships between transcription factors, chromatin and target genes and the internal consistency between the two omics, here measured in terms of structural balance in the sample correlations along elementary length-3 cycles. Results: We propose a framework that uses the a priori knowledge on the data to infer elementary causal regulatory motifs (namely chains and forks) in the network. It is based on the notions of conditional independence and partial correlation, and can be applied to both longitudinal and non-longitudinal data. Our analysis highlights a strong connection between the causal regulatory motifs that are selected by the data and the structural balance of the underlying sample correlation graphs: strikingly, >97% of the selected regulatory motifs belong to a balanced subgraph. This result shows that internal consistency, as measured by structural balance, is close to a necessary condition for 3-node regulatory motifs to satisfy causality rules.

Introduction: By means of adding more sensor technology, modern hearing aids (HAs) strive to become better, more personalized, and self-adaptive devices that can handle environmental changes and cope with the day-to-day fitness of the users. The latest HA technology available in the market already combines sound analysis with motion activity classification based on accelerometers to adjust settings. While there is a lot of research in activity tracking using accelerometers in sports applications and consumer electronics, there is not yet much in hearing research. Objective: This study investigates the feasibility of activity tracking with ear-level accelerometers and how it compares to waist-mounted accelerometers, which is a more common measurement location. Method: The activity classification methods in this study are based on supervised learning. The experimental set up consisted of 21 subjects, equipped with two XSens MTw Awinda at ear-level and one at waist-level, performing nine different activities. Results: The highest accuracy on our experimental data as obtained with the combination of Bagging and Classification tree techniques. The total accuracy over all activities and users was 84% (ear-level), 90% (waist-level), and 91% (ear-level + waist-level). Most prominently, the classes, namely, standing, jogging, laying (on one side), laying (face-down), and walking all have an accuracy of above 90%. Furthermore, estimated ear-level step-detection accuracy was 95% in walking and 90% in jogging. Conclusion: It is demonstrated that several activities can be classified, using ear-level accelerometers, with an accuracy that is on par with waist-level. It is indicated that step-detection accuracy is comparable to a high-performance wrist device. These findings are encouraging for the development of activity applications in hearing healthcare.

We demonstrate that finite impulse response (FIR) models can be applied to analyze the time evolution of an epidemic with its impact on deaths and healthcare strain. Using time series data for COVID-19-related cases, ICU admissions and deaths from Sweden, the FIR model gives a consistent epidemiological trajectory for a simple delta filter function. This results in a consistent scaling between the time series if appropriate time delays are applied and allows the reconstruction of cases for times before July 2020, when RT-PCR testing was not widely available. Combined with randomized RT-PCR study results, we utilize this approach to estimate the total number of infections in Sweden, and the corresponding infection-to-fatality ratio (IFR), infection-to-case ratio (ICR), and infection-to-ICU admission ratio (IIAR). Our values for IFR, ICR and IIAR are essentially constant over large parts of 2020 in contrast with claims of healthcare adaptation or mutated virus variants importantly affecting these ratios. We observe a diminished IFR in late summer 2020 as well as a strong decline during 2021, following the launch of a nation-wide vaccination program. The total number of infections during 2020 is estimated to 1.3 million, indicating that Sweden was far from herd immunity.

Objectives: The investigation of auditory cognitive processes recently moved from strictly controlled, trial-based paradigms toward the presentation of continuous speech. This also allows the investigation of listening effort on larger time scales (i.e., sustained listening effort). Here, we investigated the modulation of sustained listening effort by a noise reduction algorithm as applied in hearing aids in a listening scenario with noisy continuous speech. The investigated directional noise reduction algorithm mainly suppresses noise from the background. Design: We recorded the pupil size and the EEG in 22 participants with hearing loss who listened to audio news clips in the presence of background multi-talker babble noise. We estimated how noise reduction (off, on) and signal-to-noise ratio (SNR; +3 dB, +8 dB) affect pupil size and the power in the parietal EEG alpha band (i.e., parietal alpha power) as well as the behavioral performance. Results: Our results show that noise reduction reduces pupil size, while there was no significant effect of the SNR. It is important to note that we found interactions of SNR and noise reduction, which suggested that noise reduction reduces pupil size predominantly under the lower SNR. Parietal alpha power showed a similar yet nonsignificant pattern, with increased power under easier conditions. In line with the participants reports that one of the two presented talkers was more intelligible, we found a reduced pupil size, increased parietal alpha power, and better performance when people listened to the more intelligible talker. Conclusions: We show that the modulation of sustained listening effort (e.g., by hearing aid noise reduction) as indicated by pupil size and parietal alpha power can be studied under more ecologically valid conditions. Mainly concluded from pupil size, we demonstrate that hearing aid noise reduction lowers sustained listening effort. Our study approximates to real-world listening scenarios and evaluates the benefit of the signal processing as can be found in a modern hearing aid.

In this work we consider a collective decision-making process in a network of agents described by a nonlinear interconnected dynamical model with sigmoidal nonlinearities and signed interaction graph. The decisions are encoded in the equilibria of the system. The aim is to investigate this multiagent system when the signed graph representing the community is not structurally balanced and in particular as we vary its frustration, i.e., its distance to structural balance. The model exhibits bifurcations, and a "social effort" parameter, added to the model to represent the strength of the interactions between the agents, plays the role of bifurcation parameter in our analysis. We show that, as the social effort increases, the decision-making dynamics exhibits a pitchfork bifurcation behavior where, from a deadlock situation of "no decision" (i.e., the origin is the only globally stable equilibrium point), two possible (alternative) decision states for the community are achieved (corresponding to two nonzero locally stable equilibria). The value of social effort for which the bifurcation is crossed (and a decision is reached) increases with the frustration of the signed network. 

This article deals with the codesign of an output-dependent switching function and a full-order affine filter for discrete-time switched affine systems. More specifically, from the measured output, the switched filter has the role of providing essential information for the switching function, which must assure global practical stability of a desired equilibrium point. The design conditions are based on a general quadratic Lyapunov function and are expressed in terms of linear matrix inequalities. Moreover, whenever the system is quadratically detectable, the solution to the output feedback problem coincides with the one for the state feedback case and the associated filter admits the observer form. To the best of the authors knowledge, this is the first time that the dynamic output feedback control problem is treated in the context of discrete-time switched affine systems. The results can be used to cope with sampled-data control and have the property of assuring global asymptotic stability when the sampling period tends to zero. A practical application concerning the velocity control of a dc motor driven by a buck-boost converter illustrates the theoretical results.

This tutorial reviews a series of reinforcement learning (RL) methods implemented in a reproducing kernel Hilbert space (RKHS) developed to address the challenges imposed on decoder design. RL-based decoders enable the user to learn the prosthesis control through interactions without desired signals and better represent the subjects goal to complete the task. The numerous actions in complex tasks and nonstationary neural states form a vast and dynamic state-action space, imposing a computational challenge in the decoder to detect the emerging neural patterns as well as quickly establish and adjust the globally optimal policy.

Recently we showed that higher reward results in increased pupil dilation during listening (listening effort). Remarkably, this effect was not accompanied with improved speech reception. Still, increased listening effort may reflect more in-depth processing, potentially resulting in a better memory representation of speech. Here, we investigated this hypothesis by also testing the effect of monetary reward on recognition memory performance. Twenty-four young adults performed speech reception threshold (SRT) tests, either hard or easy, in which they repeated sentences uttered by a female talker masked by a male talker. We recorded the pupil dilation response during listening. Participants could earn a high or low reward and the four conditions were presented in a blocked fashion. After each SRT block, participants performed a visual sentence recognition task. In this task, the sentences that were presented in the preceding SRT task were visually presented in random order and intermixed with unfamiliar sentences. Participants had to indicate whether they had previously heard the sentence or not. The SRT and sentence recognition were affected by task difficulty but not by reward. Contrary to our previous results, peak pupil dilation did not reflect effects of reward. However, post-hoc time course analysis (GAMMs) revealed that in the hard SRT task, the pupil response was larger for high than low reward. We did not observe an effect of reward on visual sentence recognition. Hence, the current results provide no conclusive evidence that the effect of monetary reward on the pupil response relates to the memory encoding of speech. (c) 2021 The Author(s). Published by Elsevier B.V. This is an open access article under the CC BY license ( http://creativecommons.org/licenses/by/4.0/ )

A novel method for accurate speed estimation of a vehicle using a deep learning convolutional neural network (CNN), with accelerometer and gyroscope measurements as input, is presented. It does not suffer from the fundamental drift problem present in all dead reckoning methods, and yet yields about 2 m/s in accuracy. Efficient drift-free vehicle speed estimates are essential in many automotive applications, where internal wheel speed sensors or GPS are unavailable. Using extensive experimental data, the proposed CNN method is compared to an existing frequency analysis method. The proposed method is shown to perform significantly better, particularly during low speed and rapid speed changes where the frequency method struggles.

Amaximumlikelihood estimator is presented for self-calibrating both accelerometers and gyroscopes in an inertial sensor array, including scale factors, misalignments, biases, and sensor positions. By simultaneous estimation of the calibration parameters and the motion dynamics of the array, external equipment is not required for the method. A computational efficient iterative optimizationmethod is proposed where the calibration problem is divided into smaller subproblems. Further, an identifiability analysis of the calibration problem is presented. The analysis shows that it is sufficient to know the magnitude of the local gravity vector and the average scale factor gain of the gyroscopes, and that the array is exposed to two types ofmotions for the calibration problemto bewell defined. The proposedestimator is evaluatedby real-worldexperimentsand byMonteCarlo simulations. The results show that the parameters can be consistently estimated and that the calibration significantly improves the accuracy of the motion estimation. This enables on-the-fly calibration of small inertial sensors arrays by simply twisting them by hand.

Performing multiple experiments is common when learning internal mechanisms of complex systems. These experiments can include perturbations of parameters or external disturbances. A challenging problem is to efficiently incorporate all collected data simultaneously to infer the underlying dynamic network. This paper addresses the reconstruction of dynamic networks from heterogeneous datasets under the assumption that the underlying networks share the same Boolean structure across all experiments. Parametric models are derived for dynamical structure functions, which describe causal interactions between measured variables. Multiple datasets are integrated into one regression problem with additional demands on group sparsity to assure network sparsity and structure consistency. To acquire structured group sparsity, we propose a sampling-based method, together with extended versions of l1-methods and sparse Bayesian learning. The performance of the proposed methods is benchmarked in numerical simulation. In summary, this paper presents efficient methods on network reconstruction from multiple experiments, and reveals practical experience that could guide applications. (c) 2020 Elsevier Ltd. All rights reserved.

 This paper presents a unified optimization-based path planning approach to efficiently compute locally optimal solutions to optimal path planning problems in unstructured environments. The approach is motivated by showing that a lattice-based planner can be cast and analyzed as a bilevel optimization problem. This insight is used to integrate a lattice-based planner and an optimal control-based method in a novel way. The lattice-based planner is applied to the problem in a first step using a discretized search space. In a second step, an optimal control-based method is applied using the lattice-based solution as an initial iterate. In contrast to prior work, the system dynamics and objective function used in the first step are chosen to coincide with those used in the second step. As an important consequence, the lattice planner provides a solution which is highly suitable as a warm-start to the optimal control step. This proposed combination makes, in a structured way, benefit of sampling-based methods ability to solve combinatorial parts of the problem and optimal control-based methods ability to obtain locally optimal solutions. Compared to previous work, the proposed approach is shown in simulations to provide significant improvements in terms of computation time, numerical reliability and objective function value.

When solving a quadratic program (QP), one can improve the numerical stability of any QP solver by performing proximal-point outer iterations, resulting in solving a sequence of better conditioned QPs. In this letter we present a method which, for a given multi-parametric quadratic program (mpQP) and any polyhedral set of parameters, determines which sequences of QPs will have to be solved when using outer proximal-point iterations. By knowing this sequence, bounds on the worst-case complexity of the method can be obtained, which is of importance in, for example, real-time model predictive control (MPC) applications. Moreover, we combine the proposed method with previous work on complexity certification for active-set methods to obtain a more detailed certification of the proximal-point methods complexity, namely the total number of inner iterations.

A new method for nonintrusive elevator fault detection is presented. A computationally efficient algorithm for implementing the method is also proposed. The method is employed to detect when the elevator has been stationary for an unusually long period of time compared to historical traffic load patterns. This information can be used for fault detection but also indirectly to monitor the condition of the doors. The traffic load on the elevator is modeled as a nonhomogeneous Poisson process, and a generalized linear model is used to describe how the intensity of the process varies over time. A statistical hypothesis test is then used to determine if the elevator has been stationary for an unusually long time. The application of the proposed method is illustrated by an example where the detected faults are compared with the elevator service log. All faults were detected long before the service company was notified by the facility owner. Furthermore, based on the evaluation of 30 weeks of data, the method achieves a precision of 0.82 at a recall probability of 0.80.

In this letter we propose a method to exactly certify the complexity of an active-set method which is based on reformulating strictly convex quadratic programs to nonnegative least-squares problems. The exact complexity of the method is determined by proving the correspondence between the method and a standard primal active-set method for quadratic programming applied to the dual of the quadratic program to be solved. Once this correspondence has been established, a complexity certification method which has already been established for the primal active-set method is used to also certify the complexity of the nonnegative least-squares method. The usefulness of the proposed method is illustrated on a multi-parametric quadratic program originating from model predictive control of an inverted pendulum.

In multiagent dynamical systems, privacy protection corresponds to avoid disclosing the initial states of the agents while accomplishing a distributed task. The system-theoretic framework described in this paper for this scope, denoted dynamical privacy, relies on introducing output maps which act as masks, rendering the internal states of an agent indiscernible by the other agents. Our output masks are local (i.e., decided independently by each agent), time-varying functions asymptotically converging to the true states. The resulting masked system is also time-varying, and has the original unmasked system as its limit system. It is shown that dynamical privacy is not compatible with the existence of equilibria. Nevertheless the masked system retains the same convergence properties of the original system: the equilibria of the original systems become attractors for the masked system but lose the stability property. Application of dynamical privacy to popular examples of multiagent dynamics, such as models of social opinions, average consensus and synchronization, is investigated in detail. (c) 2020 Elsevier Ltd. All rights reserved.

Objectives Selectively attending to a target talker while ignoring multiple interferers (competing talkers and background noise) is more difficult for hearing-impaired (HI) individuals compared to normal-hearing (NH) listeners. Such tasks also become more difficult as background noise levels increase. To overcome these difficulties, hearing aids (HAs) offer noise reduction (NR) schemes. The objective of this study was to investigate the effect of NR processing (inactive, where the NR feature was switched off,vs.active, where the NR feature was switched on) on the neural representation of speech envelopes across two different background noise levels [+3 dB signal-to-noise ratio (SNR) and +8 dB SNR] by using a stimulus reconstruction (SR) method. Design To explore how NR processing supports the listeners selective auditory attention, we recruited 22 HI participants fitted with HAs. To investigate the interplay between NR schemes, background noise, and neural representation of the speech envelopes, we used electroencephalography (EEG). The participants were instructed to listen to a target talker in front while ignoring a competing talker in front in the presence of multi-talker background babble noise. Results The results show that the neural representation of the attended speech envelope was enhanced by the active NR scheme for both background noise levels. The neural representation of the attended speech envelope at lower (+3 dB) SNR was shifted, approximately by 5 dB, toward the higher (+8 dB) SNR when the NR scheme was turned on. The neural representation of the ignored speech envelope was modulated by the NR scheme and was mostly enhanced in the conditions with more background noise. The neural representation of the background noise was modulated (i.e., reduced) by the NR scheme and was significantly reduced in the conditions with more background noise. The neural representation of the net sum of the ignored acoustic scene (ignored talker and background babble) was not modulated by the NR scheme but was significantly reduced in the conditions with a reduced level of background noise. Taken together, we showed that the active NR scheme enhanced the neural representation of both the attended and the ignored speakers and reduced the neural representation of background noise, while the net sum of the ignored acoustic scene was not enhanced. Conclusion Altogether our results support the hypothesis that the NR schemes in HAs serve to enhance the neural representation of speech and reduce the neural representation of background noise during a selective attention task. We contend that these results provide a neural index that could be useful for assessing the effects of HAs on auditory and cognitive processing in HI populations.

Regularized system identification of linear time invariant systems in the presence of outliers is investigated. The finite impulse response (FIR) model and the Gaussian scale mixture are chosen to be the system model and the noise model, respectively. Two special cases of the noise model are considered: the well-known Students t distribution and a proposed G-confluent distribution. Both the FIR model parameter and the latent variables in the noise model are treated as parameters of our statistical model and moreover, the scale of the noise variance is treated as a hyper-parameter besides the hyper-parameters used to parameterize the priors of the impulse response and the latent variables. Then a variational expectation-maximization algorithm is proposed for inference of the parameters and hyper-parameters of the statistical model, and the algorithm is guaranteed to converge to a stationary point. Monte Carlo numerical simulations show that when the relative size of outliers is small, the proposed approach performs comparably to a state-of-the-art method and when the relative size of outliers and/or the occurrence probability of outliers is large, the proposed approach outperforms the state-of-the-art method. (C) 2020 Elsevier Ltd. All rights reserved.

The angular wheel speed of a vehicle is estimated by tracking the frequency of chassis vibrations measured with an accelerometer. A Bayesian filtering framework is proposed, allowing for straightforward incorporation of supporting information. The framework is evaluated on a large number of experimental test drives, showing comparable performance to the standard periodogram method. We then demonstrate the flexibility of the framework using accelerometer information in two ways, combining the high-frequency vibrations with low-frequency information about the vehicle acceleration. This is shown to improve robustness and resolve many cases where stand-alone frequency tracking fails.

This article provides an external stakeholder perspective on the influence of higher education in Sweden, exploring their views on curriculum development and qualitywork at the programme level. Semi-structured interviews with a selected number of representatives of external stakeholders involved in various educational areas were conducted at seven higher education institutions. The participants argued that changes intheir business sectors, and subsequent changes in the knowledge and skills in the labour needed, should encourage higher education institutions to adjust and develop their programmes. They did not anticipate or demand immediate changes in response to their comments, nor did they see themselves as a part of any quality assurance scheme. Uncertainties about the internal decision-making process and organisation in higher education institutions apparently do not facilitate external stakeholders’ understanding of their role in the larger scheme. However, all informants had comments on quality in higher education, perceiving it predominantly as something connected to the world of work. The practical implication of this study is that curriculum development at higher education institutions would benefit from communicating the internal decision-making processes to external stakeholders and agreeing on the expectations with them, in collaboration.

Individuals with hearing loss allocate cognitive resources to comprehend noisy speech in everyday life scenarios. Such a scenario could be when they are exposed to ongoing speech and need to sustain their attention for a rather long period of time, which requires listening effort. Two well-established physiological methods that have been found to be sensitive to identify changes in listening effort are pupillometry and electroencephalography (EEG). However, these measurements have been used mainly for momentary, evoked or episodic effort. The aim of this study was to investigate how sustained effort manifests in pupillometry and EEG, using continuous speech with varying signal-to-noise ratio (SNR). Eight hearing-aid users participated in this exploratory study and performed a continuous speech-in-noise task. The speech material consisted of 30-second continuous streams that were presented from loudspeakers to the right and left side of the listener (+/- 30 degrees azimuth) in the presence of 4-talker background noise (+180 degrees azimuth). The participants were instructed to attend either to the right or left speaker and ignore the other in a randomized order with two different SNR conditions: 0 dB and -5 dB (the difference between the target and the competing talker). The effects of SNR on listening effort were explored objectively using pupillometry and EEG. The results showed larger mean pupil dilation and decreased EEG alpha power in the parietal lobe during the more effortful condition. This study demonstrates that both measures are sensitive to changes in SNR during continuous speech.

Novel features for joint classification of gait and device modes are proposed and multiple machine learning methods are adopted to jointly classify the modes. The classification accuracy as well as the F₁ score of two standard classification algorithms, K-nearest neighbor (KNN) and Gaussian process (GP), are evaluated and compared against a proposed neural network (NN)-based classifier. The proposed features are the correlation scores of a detected gait cycle relative to a set of unique gait signatures as well as the gait cycle time, all extracted from handheld inertial measurement units (IMUs). Each gait signature is defined such that it contains one full cycle of the human gait. In order to take the temporal correlation between classes into account, the initial classifiers’ estimates are fed into a hidden Markov model (HMM) unit to obtain the final class estimates. The performance of the proposed method is evaluated on a large data set, including two classes of gait modes (walking and running) and four classes of device modes (fixed and face-up in the hand, swinging in the hand, in the pocket, and in the backpack). The experimental results validate the reliability of the considered features and effectiveness of the HMM unit. The initial classification accuracy of the NN-based approach is 91%, which is further improved to 99% after the smoothing stage on the validation set and 98% on the test set.

The constrained consensus problem considered in this paper, denoted interval consensus, is characterized by the fact that each agent can impose a lower and upper bound on the achievable consensus value. Such constraints can be encoded in the consensus dynamics by saturating the values that an agent transmits to its neighboring nodes. We show in the paper that when the intersection of the intervals imposed by the agents is nonempty, the resulting constrained consensus problem must converge to a common value inside that intersection. In our algorithm, convergence happens in a fully distributed manner, and without need of sharing any information on the individual constraining intervals. When the intersection of the intervals is an empty set, the intrinsic nonlinearity of the network dynamics raises new challenges in understanding the node state evolution. Using Brouwer fixed-point theorem we prove that in that case there exists at least one equilibrium, and in fact the possible equilibria are locally stable if the constraints are satisfied or dissatisfied at the same time among all nodes. For graphs with sufficient sparsity it is further proven that there is a unique equilibrium that is globally attractive if the constraint intervals are pairwise disjoint.

We propose a Bayesian framework for the received-signal-strength-based cooperative localization problem with unknown path loss exponent. Our purpose is to infer the marginal posterior of each unknown parameter: the position or the path loss exponent. This probabilistic inference problem is solved using message passing algorithms that update messages and beliefs iteratively. For numerical tractability, we combine the variable discretization and Monte-Carlo-based numerical approximation schemes. To further improve computational efficiency, we develop an auxiliary importance sampler that updates the beliefs with the help of an auxiliary variable. An important ingredient of the proposed auxiliary importance sampler is the ability to sample from a normalized likelihood function. To this end, we develop a stochastic sampling strategy that mathematically interprets and corrects an existing heuristic strategy. The proposed message passing algorithms are analyzed systematically in terms of computational complexity, demonstrating the computational efficiency of the proposed auxiliary importance sampler. Various simulations are conducted to validate the overall good performance of the proposed algorithms.

An Informed Path Planning algorithm for multiple agents is presented. It can be used to efficiently utilize available agents when surveying large areas, when total coverage is unattainable. Internally the algorithm has a Probability Hypothesis Density (PHD) representation, inspired by modern multi-target tracking methods, to represent unseen objects. Using the PHD, the expected number of observed objects is optimized. In a sequential manner, each agent maximizes the number of observed new targets, taking into account the probability of undetected objects due to previous agents actions and the probability of detection, which yields a scalable algorithm. Algorithm properties are evaluated in simulations, and shown to outperform a greedy base line method. The algorithm is also evaluated by applying it to a sea ice tracking problem, using two datasets collected in the Arctic, with reasonable results. An implementation is provided under an Open Source license.

System identification is a mature research area with well established paradigms, mostly based on classical statistical methods. Recently, there has been considerable interest in so called kernel-based regularisation methods applied to system identification problem. The recent literature on this is extensive and at times difficult to digest. The purpose of this contribution is to provide an accessible account of the main ideas and results of kernel-based regularisation methods for system identification. The focus is to assess the impact of these new techniques on the field and traditional paradigms.

A Bayesian zero-velocity detector for foot-mounted inertial navigation systems is presented. The detector extends existing zero-velocity detectors based on the likelihood-ratio test and allows, possibly time-dependent, prior information about the two hypotheses-the sensors being stationary or in motion-to be incorporated into the test. It is also possible to incorporate information about the cost of a missed detection or a false alarm. Specifically, we consider a hypothesis prior based on the velocity estimates provided by the navigation system and an exponential model for how the cost of a missed detection increases with the time since the last zero-velocity update. Thereby, we obtain a detection threshold that adapts to the motion characteristics of the user. Thus, the proposed detection framework efficiently solves one of the key challenges in current zero-velocity-aided inertial navigation systems: the tuning of the zero-velocity detection threshold. A performance evaluation on data with normal and fast gait demonstrates that the proposed detection framework outperforms any detector that chooses two separate fixed thresholds for the two gait speeds.

Nonlinear system identification is an extremely broad topic, since every system that is not linear is nonlinear. That makes it impossible to give a full overview of all aspects of the fi eld. For this reason, the selection of topics and the organization of the discussion are strongly colored by the personal journey of the authors in this nonlinear universe.

Posterior Cramér-Rao bounds (CRBs) are derived for the estimation performance of three Gaussian process-based state-space models. The parametric CRB is derived for the case with a parametric state transition and a Gaussian process-based measurement model. We illustrate the theory with a target tracking example and derive both parametric and posterior filtering CRBs for this specific application. Finally, the theory is illustrated with a positioning problem, with experimental data from an office environment where the obtained estimation performance is compared to the derived CRBs.

Maneuvering a general 2-trailer with a car-like tractor in backward motion is a task that requires a significant skill to master and is unarguably one of the most complicated tasks a truck driver has to perform. This paper presents a path planning and path-following control solution that can be used to automatically plan and execute difficult parking and obstacle avoidance maneuvers by combining backward and forward motion. A lattice-based path planning framework is developed in order to generate kinematically feasible and collision-free paths and a path-following controller is designed to stabilize the lateral and angular path-following error states during path execution. To estimate the vehicle state needed for control, a nonlinear observer is developed, which only utilizes information from sensors that are mounted on the car-like tractor, making the system independent of additional trailer sensors. The proposed path-planning and path-following control framework is implemented on a full-scale test vehicle and results from simulations and real-world experiments are presented.

INTRODUCTION: The HeartMate 3^™ has shown lower rates of adverse events compared to previous devices due to the design and absence of mechanical bearings. For previous devices, sound analysis emerged as a way to assess pump function. The aims of this study were to determine if sound analysis can be applied to the HeartMate 3 in vivo and in vitro and to evaluate an electronic stethoscope.

METHOD: Sound recordings were performed with microphones and clinical accessible electronic stethoscope. The recordings were studied in both the time and the frequency domains. Recordings from four patients were performed to determine if in vivo and in vitro recordings are comparable.

RESULTS: The results show that it is possible to detect sound from HeartMate 3 and the sound spectrum is clear. Pump frequency and frequency of the pulsatile mode are easily determined. Frequency spectra from in vitro and in vivo recordings have the same pattern, and the major proportion (96.7%) of signal power is located at the pump speed frequency ±40 Hz. The recordings from the patients show low inter-individual differences except from location of peaks originating from pump speed and harmonics. Electronic stethoscopes could be used for sound recordings, but the dedicated equipment showed a clearer sound spectrum.

DISCUSSION: The results show that acoustic analysis can also be performed with the HeartMate 3 and that in vivo and in vitro sound spectrum is similar. The frequency spectra are different from previous devices, and methods for assessing pump function or thrombosis need further evaluation.

Objectives: The aim of this study was to compare listening effort, as estimated via pupillary response, during a speech-in-noise test in bone-anchored hearing system (BAHS) users wearing three different sound processors. The three processors, Ponto Pro (PP), Ponto 3 (P3), and Ponto 3 SuperPower (P3SP), differ in terms of maximum force output (MFO) and MFO algorithm. The hypothesis was that listeners would allocate lower listening effort with the P3SP than with the PP, as a consequence of a higher MFO and, hence, fewer saturation artifacts in the signal. Design: Pupil dilations were recorded in 21 BAHS users with a conductive or mixed hearing loss, during a speech-in-noise test performed at positive signal-to-noise ratios (SNRs), where the speech and noise levels were individually adjusted to lead to 95% correct intelligibility with the PP. The listeners had to listen to a sentence in noise, retain it for 3 seconds and then repeat it, while an eye-tracking camera recorded their pupil dilation. The three sound processors were tested in random order with a single-blinded experimental design. Two conditions were performed at the same SNR: Condition 1, where the speech level was designed to saturate the PP but not the P3SP, and condition 2, where the overall sound level was decreased relative to condition 1 to reduce saturation artifacts. Results: The P3SP led to higher speech intelligibility than the PP in both conditions, while the performance with the P3 did not differ from the performance with the PP and the P3SP. Pupil dilations were analyzed in terms of both peak pupil dilation (PPD) and overall pupil dilation via growth curve analysis (GCA). In condition 1, a significantly lower PPD, indicating a decrease in listening effort, was obtained with the P3SP relative to the PP. The PPD obtained with the P3 did not differ from the PPD obtained with the other two sound processors. In condition 2, no difference in PPD was observed across the three processors. The GCA revealed that the overall pupil dilation was significantly lower, in both conditions, with both the P3SP and the P3 relative to the PP, and, in condition 1, also with the P3SP relative to the P3. Conclusions: The overall effort to process a moderate to loud speech signal was significantly reduced by using a sound processor with a higher MFO (P3SP and P3), as a consequence of fewer saturation artifacts. These findings suggest that sound processors with a higher MFO may help BAHS users in their everyday listening scenarios, in particular in noisy environments, by improving sound quality and, thus, decreasing the amount of cognitive resources utilized to process incoming speech sounds.

Modeling and failure prediction are important tasks in many engineering systems. For these tasks, the machine learning literature presents a large variety of models such as classification trees, random forest, artificial neural networks, among others. Standard statistical models such as the logistic regression, linear discriminant analysis, k-nearest neighbors, among others, can be applied. This work evaluates advantages and limitations of statistical and machine learning methods to predict failures in industrial robots. The work is based on data from more than five thousand robots in industrial use. Furthermore, a new approach combining standard statistical and machine learning models, named hybrid gradient boosting, is proposed. Results show that the hybrid gradient boosting achieves significant improvement as compared to statistical and machine learning methods. Furthermore, local joint information has been identified as the main driver for failure detection, whereas failure classification can be improved using additional information from different joints and hybrid models. (C) 2019 Elsevier Ltd. All rights reserved.

Differential graphical games have been introduced in the literature to solve state synchronization problem for linear homogeneous agents. When the agents are heterogeneous, the previous notion of graphical games cannot be used anymore and a new definition is required. In this paper, we define a novel concept of differential graphical games for linear heterogeneous agents subject to external unmodeled disturbances, which contain the previously introduced graphical game for homogeneous agents as a special case. Using our new formulation, we can solve both the output regulation and H-infinity output regulation problems. Our graphical game framework yields coupled Hamilton-Jacobi-Bellman equations, which are, in general, impossible to solve analytically. Therefore, we propose a new actor-critic algorithm to solve these coupled equations numerically in real time. Moreover, we find an explicit upper bound for the overall L2-gain of the output synchronization error with respect to disturbance. We demonstrate our developments by a simulation example.

A signed network is a network in which each link is associated with a positive or negative sign. Models for nodes interacting over such signed networks arise from various biological, social, political, and economic systems. As modifications to the conventional DeGroot dynamics for positive links, two basic types of negative interactions along negative links, namely, the opposing rule and the repelling rule, have been proposed and studied in the literature. This paper reviews a few fundamental convergence results for such dynamics over deterministic or random signed networks under a unified algebraic-graphical method. We show that a systematic tool for studying node state evolution over signed networks can be obtained utilizing generalized Perron-Frobenius theory, graph theory, and elementary algebraic recursions.

We consider a linear state estimation problem, where, in addition to the usual timestamped measurements, observations with uncertain timestamps are available. Such observations could, e.g., come from traces left by a target in a tracking scenario or from witnesses of an event and have the potential to improve the estimation accuracy significantly. We derive the posterior distribution and point estimators for a linear Gaussian smoothing formulation of this problem and illustrate with two numerical examples.

This tutorial provides a gentle introduction to the particle Metropolis-Hastings (PMH) algorithm for parameter inference in nonlinear state-space models together with a software implementation in the statistical programming language R. We employ a step-by-step approach to develop an implementation of the PMH algorithm (and the particle filter within) together with the reader. This final implementation is also available as the package pmhtutorial from the Comprehensive R Archive Network (CRAN) repository. Throughout the tutorial, we provide some intuition as to how the algorithm operates and discuss some solutions to problems that might occur in practice. To illustrate the use of PMH, we consider parameter inference in a linear Gaussian state-space model with synthetic data and a nonlinear stochastic volatility model with real-world data.

Migraine is a neurological disorder characterized by persisting attacks, underlined by the sensitivity to light. One of the leading reasons that make migraine a bigger issue is that it cannot be diagnosed easily by physicians because of the numerous overlapping symptoms with other diseases, such as epilepsy and tension-headache. Consequently, studies have been growing on how to make a computerized decision support system for diagnosis of migraine. In most laboratory studies, flash stimulation is used during the recording of electroencephalogram (EEG) signals with different frequencies and variable (seconds) time windows. The main contribution of this study is the investigation of the effects of flash stimulation on the classification accuracy, and how to find the effective window length for EEG signal classification. To achieve this, we tested different machine learning algorithms on the EEG signals features extracted by using discrete wavelet transform. Our tests on the real-world dataset, recorded in the laboratory, show that the flash stimulation can improve the classification accuracy for more than 10%. Not surprisingly, it is seen that the same holds for the selection of time window length, i.e. the selection of the proper window length is crucial for the accurate migraine identification. (C) 2018 Elsevier Ltd. All rights reserved.

This paper proposes a new method to identify nonlinear state-space systems from heterogeneous datasets. The method is described in the context of identifying biochemical/gene networks (i.e., identifying both reaction dynamics and kinetic parameters) from experimental data. Simultaneous integration of various datasets has the potential to yield better performance for system identification. Data collected experimentally typically vary depending on the specific experimental setup and conditions. Typically, heterogeneous data are obtained experimentally through 1) replicate measurements from the same biological system or 2) application of different experimental conditions such as changes/perturbations in biological inductions, temperature, gene knock-out, gene over-expression, etc. We formulate here the identification problem using a Bayesian learning framework that makes use of “sparse group” priors to allow inference of the sparsest model that can explain the whole set of observed heterogeneous data. To enable scale up to large number of features, the resulting nonconvex optimization problem is relaxed to a reweighted Group Lasso problem using a convex-concave procedure. As an illustrative example of the effectiveness of our method, we use it to identify a genetic oscillator (generalized eight species repressilator). Through this example we show that our algorithm outperforms Group Lasso when the number of experiments is increased, even when each single time-series dataset is short. We additionally assess the robustness of our algorithm against noise by varying the intensity of process noise and measurement noise.

The behavior of a person during a conversation typically involves both auditory and visual attention. Visual attention implies that the person directs his or her eye gaze toward the sound target of interest, and hence, detection of the gaze may provide a steering signal for future hearing aids. The steering could utilize a beamformer or the selection of a specific audio stream from a set of remote microphones. Previous studies have shown that eye gaze can be measured through electrooculography (EOG). To explore the precision and real-time feasibility of the methodology, seven hearing-impaired persons were tested, seated with their head fixed in front of three targets positioned at -30 degrees, 0 degrees, and +30 degrees azimuth. Each target presented speech from the Danish DAT material, which was available for direct input to the hearing aid using head-related transfer functions. Speech intelligibility was measured in three conditions: a reference condition without any steering, a condition where eye gaze was estimated from EOG measures to select the desired audio stream, and an ideal condition with steering based on an eye-tracking camera. The "EOG-steering" improved the sentence correct score compared with the "no-steering" condition, although the performance was still significantly lower than the ideal condition with the eye-tracking camera. In conclusion, eye-gaze steering increases speech intelligibility, although real-time EOG-steering still requires improvements of the signal processing before it is feasible for implementation in a hearing aid.

The stable spline (SS) kernel and the diagonal correlated (DC) kernel are two kernels that have been applied and studied extensively for kernel-based regularized LTI system identification. In this note, we show that similar to the derivation of the SS kernel, the continuous-time DC kernel can be derived by applying the same "stable" coordinate change to a "generalized" first-order spline kernel, and thus, can be interpreted as a stable generalized first-order spline kernel. This interpretation provides new facets to understand the properties of the DC kernel. In particular, we derive a new orthonormal basis expansion of the DC kernel and the explicit expression of the norm of the reproducing kernel Hilbert space associated with the DC kernel. Moreover, for the nonuniformly sampled DC kernel, we derive its maximum entropy property and show that its kernel matrix has tridiagonal inverse.

Filtering and smoothing algorithms for linear discrete-time state-space models with skew-t-distributed measurement noise are proposed. The algorithms use a variational Bayes based posterior approximation with coupled location and skewness variables to reduce the error caused by the variational approximation. Although the variational update is done suboptimally using an expectation propagation algorithm, our simulations show that the proposed method gives a more accurate approximation of the posterior covariance matrix than an earlier proposed variational algorithm. Consequently, the novel filter and smoother outperform the earlier proposed robust filter and smoother and other existing low-complexity alternatives in accuracy and speed. We present both simulations and tests based on real-world navigation data, in particular the global positioning system data in an urban area, to demonstrate the performance of the novel methods. Moreover, the extension of the proposed algorithms to cover the case where the distribution of the measurement noise is multivariate skew-t is outlined. Finally, this paper presents a study of theoretical performance bounds for the proposed algorithms.

The models of collective decision-making considered in this paper are nonlinear interconnected cooperative systems with saturating interactions. These systems encode the possible outcomes of a decision process into different steady states of the dynamics. In particular, they are characterized by two main attractors in the positive and negative orthant, representing two choices of agreement among the agents, associated to the Perron-Frobenius eigenvector of the system. In this paper we give conditions for the appearance of other equilibria of mixed sign. The conditions are inspired by Perron-Frobenius theory and are related to the algebraic connectivity of the network. We also show how all these equilibria must be contained in a solid disk of radius given by the norm of the equilibrium point which is located in the positive orthant.

The problem of path planning for mobilesensors with the task of target monitoring is considered. A receding horizon optimal control approach based on the information filter is presented, where the limited field of view of the sensor can be modeled by introducing binary variables. The resulting nonlinear mixed integer problem to be solved in each sample, with no apparent tractable solution, is shown to be equivalent to a problem that robustly can be solved to global optimality using off-the-shelf optimization tools.

The kernel-based regularization method has two core issues: kernel design and hyperparameter estimation. In this paper, we focus on the second issue and study the properties of several hyperparameter estimators including the empirical Bayes (EB) estimator, two Steins unbiased risk estimators (SURE) (one related to impulse response reconstruction and the other related to output prediction) and their corresponding Oracle counterparts, with an emphasis on the asymptotic properties of these hyperparameter estimators. To this goal, we first derive and then rewrite the first order optimality conditions of these hyperparameter estimators, leading to several insights on these hyperparameter estimators. Then we show that as the number of data goes to infinity, the two SUREs converge to the best hyperparameter minimizing the corresponding mean square error, respectively, while the more widely used EB estimator converges to another best hyperparameter minimizing the expectation of the EB estimation criterion. This indicates that the two SUREs are asymptotically optimal in the corresponding MSE senses but the EB estimator is not. Surprisingly, the convergence rate of two SUREs is slower than that of the EB estimator, and moreover, unlike the two SUREs, the EB estimator is independent of the convergence rate of Phi(T)Phi/N to its limit, where Phi is the regression matrix and N is the number of data. A Monte Carlo simulation is provided to demonstrate the theoretical results. (C) 2018 Elsevier Ltd. All rights reserved.

There are two primary sources of sensor measurements for driver behavior profiling within insurance telematics and fleet management. The first is the on-board diagnostics system, typically found within most modern cars. The second is the global navigation satellite system, whose associated receivers commonly are embedded into smartphones or off-the-shelf telematics devices. In this paper, we present maximum likelihood and maximum a posteriori estimators for the problem of fusing speed measurements from these two sources to jointly estimate a vehicles speed and the scale factor of the wheel speed sensors. In addition, we analyze the performance of the estimators by use of the Cramer-Rao bound, and discuss the estimation of model parameters describing measurement errors and vehicle dynamics. Last, simulations and real-world data are used to show that the proposed estimators yield a substantial performance gain compared to when employing only one of the two measurement sources.

Introduction: The use of left ventricular assist device (LVAD) has grown rapidly. Adverse events do continue to occur. In recent years, analysis of LVAD sound recordings emerged as a means to monitor pump function and detect pump thrombosis. The aim of this study was to characterize the sounds from HeartMate II and to evaluate the use of handheld iOS devices for sound recordings. Method: Signal analysis of LVAD sound recordings, with dedicated recording equipment and iOS devices, was performed. Two LVADs running in mock loop circuits were compared to an implanted LVAD. Spectral analysis and parametric signal models were explored to quantify the sound and potentially detect changes in it. Results: The sound recordings of two LVADs in individual mock loop circuits and a third one implanted in a patient appeared to be similar. Qualitatively, sound characteristics were preserved following changes in pump speed. Recordings using dedicated equipment showed that HeartMate II sound comprises low-frequency components corresponding to pump impeller rotation, as well as high-frequency components due to a pulse width modulation of the electric power to the pump. These different signal components interact and result in a complicated frequency spectrum. The iPhone and iPod recordings could not reproduce the sounds as well as the dedicated equipment. In particular, lower frequencies were affected by outside disturbances. Discussion: This article outlines a systematic approach to LVAD sound analysis using signal processing methods to quantify and potentially detect changes, and describes some of the challenges, for example, with the use of inexpensive recording devices.

This note considers the identification of large-scale one-dimensional networks consisting of identical LTI dynamical systems. A subspace identification method is developed that only uses local input-output information and does not rely on knowledge about the local state interaction. The proposed identification method estimates the Markov parameters of a locally lifted system, following the state-space realization of a single subsystem. The Markov-parameter estimation is formulated as a rank minimization problem by exploiting the low-rank property and the two-layer Toeplitz structural property in the data equation, whereas the state-space realization of a single subsystem is formulated as a structured low-rank matrix-factorization problem. The effectiveness of the proposed identification method is demonstrated by simulation examples.

Time-correlated single-photon counting lidar provides very high-resolution range measurements, making the technology interesting for 3D imaging of objects behind foliage or other obscuration. We study six peak detection approaches and compare their performance from several perspectives: detection of double surfaces within the instantaneous field of view, range accuracy, performance under sparse sampling, and the number of outliers. The results presented are based on reference measurements of a characterization target. Special consideration is given to the possibility of resolving two surfaces closely separated in range within the field of view of a single pixel. An approach based on fitting a linear combination of impulse response functions to the collected data showed the best overall performance. (c) 2018 Society of Photo-Optical Instrumentation Engineers (SPIE).

This paper presents a new method to compute upper and lower bounds of any voltage or current of an arbitrary linear electric circuit model with uncertain parameters. The bounds are in the frequency domain, and when compared to a previously proposed method, this novel approach provides a higher level of guarantee. The reason is that the bounds are not only computed for a set of fixed frequencies but also computed to a set of intervals of frequencies. The details of the proposed approach, especially the equivalent uncertain element models, are given. Additionally, tests are performed on problems with low and high number of uncertain parameters. Contrary to the classical method of Monte Carlo, the results are not based on a random choice of parameters and do not depend on the number of iterations. It is shown on an example that the classical method of Monte Carlo needs a high number of iterations to reach results in agreement with the proposed method. Then, it leads to higher computation times of several orders of magnitude.

Nonlinear cooperative systems associated to vector fields that are concave or subhomogeneous describe well interconnected dynamics that are of key interest for communication, biological, economical, and neural network applications. For this class of positive systems, we provide conditions that guarantee existence, uniqueness and stability of strictly positive equilibria. These conditions can be formulated directly in terms of the spectral radius of the Jacobian of the system. If control inputs are available, then it is shown how to use state feedback to stabilize an equilibrium point in the interior of the positive orthant.

This study proposes a new model which is fully specified for automated seizure onset detection and seizure onset prediction based on electroencephalography (EEG) measurements. We processed two archetypal EEG databases, Freiburg (intracranial EEG) and CHB-MIT (scalp EEG), to find if our model could outperform the state-of-the art models. Four key components define our model: (1) multiscale principal component analysis for EEG de-noising, (2) EEG signal decomposition using either empirical mode decomposition, discrete wavelet transform or wavelet packet decomposition, (3) statistical measures to extract relevant features, (4) machine learning algorithms. Our model achieved overall accuracy of 100% in ictal vs. inter-ictal EEG for both databases. In seizure onset prediction, it could discriminate between inter-ictal, pre-ictal, and ictal EEG with the accuracy of 99.77%, and between inter-ictal and pre-ictal EEG states with the accuracy of 99.70%. The proposed model is general and should prove applicable to other classification tasks including detection and prediction regarding bio-signals such as EMG and ECG. (C) 2017 Elsevier Ltd. All rights reserved.

In a continuous-time nonlinear driftless control system, an involutive flow is a composition of input profiles that does not excite any Lie bracket. Such flow composition is trivial, as it corresponds to a "forth and back" cyclic motion obtained rewinding the system along the same path. The aim of this paper is to show that, on the contrary, when a (nonexact) discretization of the nonlinear driftless control system is steered along the same trivial input path, it produces a net motion, which is related to the gap between the discretization used and the exact discretization given by a Taylor expansion. These violations of involutivity can be used to provide an estimate of the local truncation error of numerical integration schemes. In the special case in which the state of the driftless control system admits a splitting into shape and phase variables, our result corresponds to saying that the geometric phases of the discretization need not obey an area rule, i.e., even zero-area cycles in shape space can lead to nontrivial geometric phases. (C) 2017 Elsevier B.V. All rights reserved.

The local approach to linear parameter varying (LPV) system identification consists in interpolating individually estimated local linear time invariant (LTI) models corresponding to fixed values of the scheduling variable. It is shown in this paper that, without any global structural assumption of the considered LPV system, individually estimated local state-space LTI models do not contain sufficient information for determining similarity transformations making them coherent. It is possible to estimate these similarity transformations from input-output data under appropriate excitation conditions. (C) 2017 Published by Elsevier Ltd.

The ensemble Kalman filter (EnKF) is a Monte Carlo-based implementation of the Kalman filter (KF) for extremely high-dimensional, possibly nonlinear, and non-Gaussian state estimation problems. Its ability to handle state dimensions in the order of millions has made the EnKF a popular algorithm in different geoscientific disciplines. Despite a similarly vital need for scalable algorithms in signal processing, e.g., to make sense of the ever increasing amount of sensor data, the EnKF is hardly discussed in our field. This self-contained review is aimed at signal processing researchers and provides all the knowledge to get started with the EnKF. The algorithm is derived in a KF framework, without the often encountered geoscientific terminology. Algorithmic challenges and required extensions of the EnKF are provided, as well as relations to sigma point KF and particle filters. The relevant EnKF literature is summarized in an extensive survey and unique simulation examples, including popular benchmark problems, complement the theory with practical insights. The signal processing perspective highlights new directions of research and facilitates the exchange of potentially beneficial ideas, both for the EnKF and high-dimensional nonlinear and non-Gaussian filtering in general.

Passive magnetic sensors measure the magnetic field density in three axes and are often integrated on a single chip. These low-cost sensors are widely used in car navigation as well as in battery powered navigation equipment such as smartphones as part of an electronic compass. We focus on a train localization application with multiple, exclusively onboard sensors and a track map. This approach is considered as a base technology for future railway applications such as collision avoidance systems or autonomous train driving. In this paper, we address the following question: how beneficial are passive magnetic measurements for train localization? We present and analyze measurements of two different magnetometers recorded on a regional train at regular passenger service. We show promising correlations of the measurements with the track positions and the traveled switch way. The processed data reveals that the railway environment has repeatable, location-dependent magnetic signatures. This is considered as a novel approach to train localization, as the use of these magnetic signals at first view is not obvious. The proposed methods based on passive magnetic measurements show a high potential to be integrated in new and existing train localization approaches.

Devising the planar routes of minimal length that are required to pass through predefined neighborhoods of target points plays an important role in reducing the missions operating cost. Two versions of the problem are considered. The first one assumes that the ordering of the targets is fixed a priori. In such a case, the optimal route is devised by solving a convex optimization problem formulated either as a second-order cone program or as a sum-of-squares optimization problem. Additional route properties, such as continuity and minimal curvature, are considered as well. The second version allows the ordering of the targets to be optimized to further reduce the route length. We show that such a problem can be solved by introducing additional binary variables, which allows the route to be designed using off-the-shelf mixed-integer solvers. A case study that shows that the proposed strategy is computationally tractable is presented.

This paper proposes a decentralized control strategy for the voltage regulation of islanded inverter-interfaced microgrids. We show that an inverter-interfaced microgrid under plug-and-play (PnP) functionality of distributed generations (DGs) can be cast as a linear time-invariant system subject to polytopic-type uncertainty. Then, by virtue of this novel description and use of the results from theory of robust control, the microgrid control system guarantees stability and a desired performance even in the case of PnP operation of DGs. The robust controller is a solution of a convex optimization problem. The main properties of the proposed controller are that: 1) it is fully decentralized and local controllers of DGs that use only local measurements; 2) the controller guarantees the stability of the overall system; 3) the controller allows PnP functionality of DGs in microgrids; and 4) the controller is robust against microgrid topology change. Various case studies, based on time-domain simulations in MATLAB/SimPowerSystems Toolbox, are carried out to evaluate the performance of the proposed control strategy in terms of voltage tracking, microgrid topology change, PnP capability features, and load changes.

A common issue with many system identification problems is that the true input to the system is unknown. This paper extends a previously presented indirect modelling framework that deals with identification of systems where the input is partially or fully unknown. In this framework, unknown inputs are eliminated by using additional measurements that directly or indirectly contain information about the unknown inputs. The resulting indirect predictor model is only dependent on known and measured signals and can be used to estimate the desired dynamics or properties. Since the input of the indirect model contains both known inputs and measurements that could all be correlated with the same disturbances as the output, estimation of the indirect model has similar challenges as a closed-loop estimation problem. In fact, due to the generality of the indirect modelling framework, it unifies a number of already existing system identification problems that are contained as special cases. For completeness, the paper is concluded with one method that can be used to estimate the indirect model as well as an experimental verification to show the applicability of the framework.

Most navigation systems today rely on global navigation satellite systems (gnss), including in cars. With support from odometry and inertial sensors, this is a sufficiently accurate and robust solution, but there are future demands. Autonomous cars require higher accuracy and integrity. Using the car as a sensor probe for road conditions in cloud-based services also sets other kind of requirements. The concept of the Internet of Things requires stand-alone solutions without access to vehicle data. Our vision is a future with both invehicle localization algorithms and after-market products, where the position is computed with high accuracy in gnss-denied environments. We present a localization approach based on a prior that vehicles spend the most time on the road, with the odometer as the primary input. When wheel speeds are not available, we present an approach solely based on inertial sensors, which also can be used as a speedometer. The map information is included in a Bayesian setting using the particle filter (PF) rather than standard map matching. In extensive experiments, the performance without gnss is shown to have basically the same quality as utilizing a gnss sensor. Several topics are treated: virtual measurements, dead reckoning, inertial sensor information, indoor positioning, off-road driving, and multilevel positioning.

There is a big trend nowadays to use event-triggered proximity report for indoor positioning. This paper presents a generic received-signal-strength (RSS) threshold optimization framework for generating informative proximity reports. The proposed framework contains five main building blocks, namely the deployment information, RSS model, positioning metric selection, optimization process and management. Among others, we focus on Gaussian process regression (GPR)-based RSS models and positioning metric computation. The optimal RSS threshold is found through minimizing the best achievable localization root-mean-square-error formulated with the aid of fundamental lower bound analysis. Computational complexity is compared for different RSS models and different fundamental lower bounds. The resulting optimal RSS threshold enables enhanced performance of new fashioned low-cost and low-complex proximity report-based positioning algorithms. The proposed framework is validated with real measurements collected in an office area where bluetooth-low-energy (BLE) beacons are deployed.

In this paper, we propose a distributed algorithm for solving coupled problems with chordal sparsity or an inherent tree structure which relies on primal–dual interior-point methods. We achieve this by distributing the computations at each iteration, using message-passing. In comparison to existing distributed algorithms for solving such problems, this algorithm requires far fewer iterations to converge to a solution with high accuracy. Furthermore, it is possible to compute an upper-bound for the number of required iterations which, unlike existing methods, only depends on the coupling structure in the problem. We illustrate the performance of our proposed method using a set of numerical examples.

This paper reviews the vast literature on static output feedback design for linear time-invariant systems including classical results and recent developments. In particular, we focus on static output feedback synthesis with performance specifications, structured static output feedback, and robustness. The paper provides a comprehensive review on existing design approaches including iterative linear matrix inequalities heuristics, linear matrix inequalities with rank constraints, methods with decoupled Lyapunov matrices, and non-Lyapunov-based approaches. We describe the main difficulties of dealing with static output feedback design and summarize the main features, advantages, and limitations of existing design methods. (C) 2016 International Federation of Automatic Control. Published by Elsevier Ltd. All rights reserved.

Todays 4G LTE systems bring unprecedented mobile broadband performance to over a billion of users across the globe. Recently, work on a 5G mobile communication system has begun, and next to a new 5G air interface, LTE will be an essential component. The evolution of LTE will therefore strive to meet 5G requirements and to address 5G use cases. In this article, we provide an overview of foreseen key technology areas and components for LTE Release 14, including latency reductions, enhancements for machine-type communication, operation in unlicensed spectrum, massive multi-antenna systems, broadcasting, positioning, and support for intelligent transportation systems.

Background: The mode of action of a drug on its targets can often be classified as being positive (activator, potentiator, agonist, etc.) or negative (inhibitor, blocker, antagonist, etc.). The signed edges of a drug-target network can be used to investigate the combined mechanisms of action of multiple drugs on the ensemble of common targets. Results: In this paper it is shown that for the signed human drug-target network the majority of drug pairs tend to have synergistic effects on the common targets, i.e., drug pairs tend to have modes of action with the same sign on most of the shared targets, especially for the principal pharmacological targets of a drug. Methods are proposed to compute this synergism, as well as to estimate the influence of the drugs on the side effect of another drug. Conclusions: Enriching a drug-target network with information of functional nature like the sign of the interactions allows to explore in a systematic way a series of network properties of key importance in the context of computational drug combinatorics.

Inspired by ideas taken from the machine learning literature, new regularization techniques have been recently introduced in linear system identification. In particular, all the adopted estimators solve a regularized least squares problem, differing in the nature of the penalty term assigned to the impulse response. Popular choices include atomic and nuclear norms (applied to Hankel matrices) as well as norms induced by the so called stable spline kernels. In this paper, a comparative study of estimators based on these different types of regularizers is reported. Our findings reveal that stable spline kernels outperform approaches based on atomic and nuclear norms since they suitably embed information on impulse response stability and smoothness. This point is illustrated using the Bayesian interpretation of regularization. We also design a new class of regularizers defined by "integral" versions of stable spline/TC kernels. Under quite realistic experimental conditions, the new estimators outperform classical prediction error methods also when the latter are equipped with an oracle for model order selection. (C) 2016 Elsevier Ltd. All rights reserved.

It is a well-known fact that externally positive linear systems may fail to have a minimal positive realization. In order to investigate these cases, we introduce the notion of minimal eventually positive realization, for which the state update matrix becomes positive after a certain power. Eventually positive realizations capture the idea that in the impulse response of an externally positive system the state of a minimal realization may fail to be positive, but only transiently. As a consequence, we show that in discrete-time it is possible to use downsampling to obtain minimal positive realizations matching decimated sequences of Markov coefficients of the impulse response. In continuous-time, instead, if the sampling time is chosen sufficiently long, a minimal eventually positive realization leads always to a sampled realization which is minimal and positive.

In this letter, numerical algorithms for computing the marginal version of the Bayesian Cramér–Rao bound (M-BCRB) for jump Markov nonlinear systems and jump Markov linear Gaussian systems are proposed. Benchmark examples for both systems illustrate that the M-BCRB is tighter than three other recently proposed BCRBs

In this work we present a practical calibration algorithm that calibrates a magnetometer using inertial sensors. The calibration corrects for magnetometer sensor errors, for the presence of magnetic disturbances and for misalignment between the magnetometer and the inertial sensor axes. It is based on a maximum likelihood formulation and is formulated as an offline method. It is shown to give good results using data from two different commercially available sensor units. Using the calibrated magnetometer measurements in combination with the inertial sensors to determine orientation, is shown to lead to significantly improved heading estimates.

Mathematical models of real life systems and processes are essential in today’s industrial work. To be able to construct such models is therefore a fundamental skill in modern engineering...

In this chapter we show that chordal structure can be used to devise efficient optimization methods for many common model predictive control problems. The chordal structure is used both for computing search directions efficiently as well as for distributing all the other computations in an interior-point method for solving the problem. The chordal structure can stem both from the sequential nature of the problem as well as from distributed formulations of the problem related to scenario trees or other formulations. The framework enables efficient parallel computations.

Comprehension of speech in noise is a challenge for hearing-impaired (HI) individuals. Electroencephalography (EEG) provides a tool to investigate the effect of different levels of signal-to-noise ratio (SNR) of the speech. Most studies with EEG have focused on spectral power in well-defined frequency bands such as alpha band. In this study, we investigate how local functional connectivity, i.e. functional connectivity within a localized region of the brain, is affected by two levels of SNR. Twenty-two HI participants performed a continuous speech in noise task at two different SNRs (+3 dB and +8 dB). The local connectivity within eight regions of interest was computed by using a multivariate phase synchrony measure on EEG data. The results showed that phase synchrony increased in the parietal and frontal area as a response to increasing SNR. We contend that local connectivity measures can be used to discriminate between speech-evoked EEG responses at different SNRs.

The unscented Kalman filter (UKF) is a method to solve nonlinear dynamic filtering problems, which internally uses the unscented transform (UT). The behavior of the UT is controlled by design parameters, seldom changed from the values suggested in early UT/UKF publications. Despite the knowledge that the UKF can perform poorly when the parameters are improperly chosen, there exist no wide spread intuitive guidelines for how to tune them. With an application relevant example, this paper shows that standard parameter values can be far from optimal. By analyzing how each parameter affects the resulting UT estimate, guidelines for how the parameter values should be chosen are developed. The guidelines are verified both in simulations and on real data collected in an underground mine. A strategy to automatically tune the parameters in a state estimation setting is presented, resulting in parameter values inline with developed guidelines.

A belief-space planning problem for GNSS-denied areas is studied where the location and number of landmarks available are unknown when performing the planning. To be able to plan an informative path in this situation, an algorithm using virtual landmarks to position the platform during the planning phase is studied.  The virtual landmarks are selected to capture the expected information available in different regions of the map, based on the beforehand known landmark density. The main contribution of this work is a better approximation of the obtained information from the virtual landmarks and a theoretical study of the properties of the approximation. Furthermore, the proposed approximation, in itself and its use in a path planner, is investigated with successful results.

An optical flow inspired magnetic-field based odometry estimation process is presented. The estimation process is based upon taking "image" like measurements of the magnetic-field using a magnetometer array. From the measurements a model of the local field is learned. Using the learned model the pose change that gives the smallest prediction error of the measurement at the next time instant is calculated. Two models for describing the magnetic-field are presented, and the performance of the odometry estimation process when using the two models is evaluated. The evaluation shows that at a high signal-to-noise ratio the pose change can be estimated with an error of only a few percentage of the true pose change. Further, the evaluation shows that the uncertainty of the estimate can be consistently estimated. Thus, the proposed odometry estimation process can be used to reduce the navigation error growth rate of, for example, inertial navigation systems by providing reliable odometry information when passing by magnetized objects.

Reinforcement learning techniques have successfully been applied to solve challenging problems. Among the more famous examples are playing games such as Go and real-time computer games such as StarCraft II. In addition, reinforcement learning has successfully been deployed in cyber-physical systems such as robots playing a curling-based game. These are all important and significant achievements indicating that the techniques can be of value for the aerospace industry. However, to use these techniques in the aerospace industry, very high requirements on verification and validation must be met. In this position paper, we outline four key problems for verification and validation of reinforcement learning techniques. Solving these are an important step towards enabling reinforcement learning techniques to be used in safety critical domains such as the aerospace industry.

Deep state space models (SSMs) are an actively researched model class for temporal models developed in the deep learning community which have a close connection to classic SSMs. The use of deep SSMs as a black-box identification model can describe a wide range of dynamics due to the flexibility of deep neural networks. Additionally, the probabilistic nature of the model class allows the uncertainty of the system to be modelled. In this work a deep SSM class and its parameter learning algorithm are explained in an effort to extend the toolbox of nonlinear identification methods with a deep learning based method. Six recent deep SSMs are evaluated in a first unified implementation on nonlinear system identification benchmarks. Copyright (C) 2021 The Authors.

The selections for the model orders and the number of controller parameters have not been discussed for many data-driven iterative learning control (ILC) methods. If they are not chosen carefully, the estimated model and designed controller will lead to either large variance or large bias. In this paper we try to use the kernel-based regularization method (KRM) to handle the model estimation problem and the controller design problem for unknown repetitive linear time-varying systems. In particular, we have used the diagonal correlated kernel and the marginal likelihood maximization method for the two problems. Numerical simulation results show that smaller mean square errors for each time instant are obtained by using the proposed ILC method in comparison with an existing data-driven ILC approach. Copyright (C) 2021 The Authors.

In this paper we present a distributed second-order augmented Lagrangian method for distributed model predictive control. We distribute the computations for search direction, step size, and termination criteria over what is known as the clique tree of the problem and calculate each of them using message passing. The algorithm converges to its centralized counterpart and it requires fewer communications between sub-systems as compared to algorithms such as the alternating direction method of multipliers. Results from a simulation study confirm the efficiency of the framework. Copyright (C) 2021 The Authors.

The use of the CDIO framework in the development of the ECIU University is presented. The paper discusses the relatively moderate adaptations and modifications of the CDIO Syllabus and Standards that are necessary to make the documents applicable also in this context. Since challenge-based learning (CBL) is central learning format in the ECIU University, special attention is given to the connections between CBL method, the conceive-design-implement-operate sequence and project-based learning, which is central in the CDIO framework. The paper presents both general aspects and examples of the applications and activities within ECIU University and Linköping University (LiU). The main messages of the paper are that the development of the ECIU University will benefit from applying the CDIO framework since it offers references for what an education should give, in terms of knowledge and skills, and how an education program should be designed. In addition, the components of the CDIO framework require a moderate amount of adaptation to be directly applicable. Examples of the ongoing implementation activities at LiU. 

Introduction 

Higher education institutions (HEIs) face challenges assessing the relevance of educational programmes. Upon graduation, the student should have acquired knowledge and understanding, competence and skills as well as good judgement and approaches to operate in a changing labour market. Ideas on new programmes and courses mainly emanate from research findings identified at the HEIs. Needs and expectations from external stakeholders have the potential to further contribute if room for collaboration is created. Extensive rapid societal changes increase this need for collaboration. 

The Standards and Guidelines for Quality Assurance in the European Higher Education Area (ESG) state that higher education aims to fulfil multiple purposes, including preparing students for active citizenship, future careers, personal development and create a broad, advanced knowledge base and stimulate research and innovation (ENQA, 2015). However, different stakeholders may have other priorities. Therefore, quality assurance needs to take these different perspectives into account. 

Relevance is considered an aspect of quality in higher education, but many HEIs in Sweden lack structures, processes, and methods for assessing relevance and involving stakeholders in these processes.

This project aimed to increase the knowledge and provide methods to systematically assess relevance and use university-industry collaborations as tools for educational development. 

Methods/Approach 

This paper is a summary of the project MERUT focused on methods for assessment of relevance in higher education. The project was carried out during 2017-2020 with financial support from Vinnova (The Swedish Innovation Agency) and involved seven Swedish HEIs. The connection to future career paths is often stated as the primary factor to describe the relevance of educational programmes and was selected as the focus of MERUT. Data have been collected using workshops, meetings, literature reviews, interview studies, and surveys. Important parts of the work have been interviews with external stakeholders in different labour-market areas who, in various ways, are involved in higher education, most often in advisory boards. Also, interviews with quality coordinators at university programmes as well as at faculty and university management levels at each HEI have been carried out.

The seven participating Swedish universities in MERUT have been Karolinska Institutet, Kristianstad University, KTH Royal Institute of Technology, Mälardalen University, Linköping University, Stockholm University, and Umeå University.

Results The project resulted in knowledge, tools, and methods to work systematically with relevance in higher education. The results can be summarized as follows:

·         Interviews with HEIs showed that they collaborate with external stakeholders in many ways, primarily around teaching and learning and to a lesser extent around programme management and quality assurance. 

·         Further, the involvement of stakeholders varied both between and within the universities (faculties, subject areas, levels of education). Therefore, it is difficult to evaluate, in a systematic way, how collaboration is included in the quality systems of the HEIs. 

·         A toolbox with methods for how to involve external stakeholders in the process of assessing and developing the relevance of a study programme. These methods can be used in continuous quality work, in major curriculum revisions as well as in the establishment of new programmes.

·         A checklist for external stakeholders' involvement in educational development to facilitate and clarify roles, structures, and tasks in connection with external stakeholders in both the development and operational phases of a study programme. 

The results show that there are many similarities between the HEIs in the study in terms of relevance assessment and dimensioning decisions. However, the potential of a systematic collaboration with stakeholders and society for relevance assessment and dimensioning of education is not yet fully being realised. 

Conclusion 

HEIs interact with external stakeholders in many ways within education. However, it is rare to find examples where external stakeholders are involved in the quality assurance process, at least not in a systematic way. The MERUT project has developed recommendations for collaboration perspectives and stakeholder participation in the governing of educational programmes. A systematic dialogue and interaction with stakeholders contribute to a mutual understanding of different stakeholder groups’ needs and expectations, and their view of quality of higher education. Furthermore, to consider relevance as a quality aspect creates a basis for a more methodical assessment process where external stakeholders can contribute in a clear role. MERUT has developed a toolbox and a checklist to facilitate such systematic interaction and collaboration with stakeholder groups. A conclusion of  this project is that a reciprocal, transparent, and systematic approach leads to a sustainable educational collaboration with improved quality and relevance of higher education.

It would constitute a large gain for society if the HEIs are able to systematically and by efficient processes take external stakeholders’ and societal needs and expectations into account when building comprehensive and systematic relevance assessment processes and in dimensioning of education.  

References

European Association for Quality Assurance in Higher Education (ENQA). (2015). Standards and Guidelines for Quality Assurance in the European Higher Education Area (ESG). Brussels, EURASHE.

In this paper, we study the influence of ill-conditioned regression matrix on two hyper-parameter estimation methods for the kernel-based regularization method: the empirical Bayes (EB) and the Steins unbiased risk estimator (SURE). First, we consider the convergence rate of the cost functions of EB and SURE, and we find that they have the same convergence rate but the influence of the ill-conditioned regression matrix on the scale factor are different: for upper bounds, the scale factor for SURE contains one more factor cond(Phi(T)Phi) than that of EB, where Phi is the regression matrix and cond(.) denotes the condition number of a matrix. This finding indicates that when Phi is ill-conditioned, i.e., cond(Phi(T)Phi) is large, the cost function of SURE converges slower than that of EB. Then we consider the convergence rate of the optimal hyper-parameters of EB and SURE, and we find that they are both asymptotically normally distributed and have the same convergence rate, but the influence of the ill-conditioned regression matrix on the scale factor are different. In particular, for the ridge regression case, we show that the optimal hyper-parameter of SURE converges slower than that of EB with a factor of 1/n(2), as cond(Phi(T)Phi) goes to infinity, where n is the FIR model order.

This work introduces antagonistic interactions into the so-called biased assimilation model of opinion dynamics, a nonlinear model expressing the bias of the agents towards their own opinions. In this model, opinions exchanged in a signed network are multiplied by a state dependent term having the bias as exponent. For small values of the bias parameters, while for structurally balanced networks polarization always occurs, we show that for structurally unbalanced networks also a state of indecision (corresponding to the centroid of the opinion hypercube) is a local attractor. When instead the biases are all large, the opinions usually polarize. In particular, for large enough biases if all agents take the same initial stance (i.e., are all on the same side of the indecision state), then the opinions polarize all to the same extreme value for both structurally balanced network and structurally unbalanced network.

We propose a sequence of pedagogical steps for introducing the Youla-Kucera parametrization, starting from the internal model principle, and introducing the control structures of disturbance observer and internal model control along the way. We provide some background on the concepts and a brief survey of their treatment in textbooks on control. Copyright (C) 2020 The Authors.

The task of maneuvering a multi-steered articulated vehicle in confined environments is difficult even for experienced drivers. In this work, we present an optimization-based trajectory planner targeting low-speed maneuvers in unstructured environments for multi-steered N-trailer vehicles, which are comprised of a car-like tractor and an arbitrary number of interconnected trailers with fixed or steerable wheels. The proposed trajectory planning framework is divided into two steps, where a lattice-based trajectory planner is used in a first step to compute a resolution optimal solution to a discretized version of the trajectory planning problem. The output from the lattice planner is then used in a second step to initialize an optimal control problem solver, which enables the framework to compute locally optimal trajectories that start at the vehicles initial state and reaches the goal state exactly. The performance of the proposed optimization-based trajectory planner is evaluated in a set of practically relevant scenarios for a multi-steered 3-trailer vehicle with a car-like tractor where the last trailer is steerable. Copyright (C) 2020 The Authors.

This paper presents an optimization-based receding horizon trajectory planning algorithm for dynamical systems operating in unstructured and cluttered environments. The proposed approach is a two-step procedure that uses a motion planning algorithm in a first step to efficiently find a feasible, but possibly suboptimal, nominal solution to the trajectory planning problem where in particular the combinatorial aspects of the problem are solved. The resulting nominal trajectory is then improved in a second optimization-based receding horizon planning step which performs local trajectory refinement over a sliding time window. In the second step, the nominal trajectory is used in a novel way to both represent a terminal manifold and obtain an upper bound on the cost-to-go online. This enables the possibility to provide theoretical guarantees in terms of recursive feasibility, objective function value, and convergence to the desired terminal state. The established theoretical guarantees and the performance of the proposed algorithm are verified in a set of challenging trajectory planning scenarios for a truck and trailer system. Copyright (C) 2020 The Authors.

Driving heavy-duty vehicles, such as buses and tractor-trailer vehicles, is a difficult task in comparison to passenger cars. Most research on motion planning for autonomous vehicles has focused on passenger vehicles, and many unique challenges associated with heavy-duty vehicles remain open. However, recent works have started to tackle the particular difficulties related to on-road motion planning for buses and tractor-trailer vehicles using numerical optimization approaches. In this work, we propose a framework to design an optimization objective to be used in motion planners. Based on geometric derivations, the method finds the optimal trade-off between the conflicting objectives of centering different axles of the vehicle in the lane. For the buses, we consider the front and rear axles trade-off, whereas for articulated vehicles, we consider the tractor and trailer rear axles trade-off. Our results show that the proposed design strategy produces planned paths that considerably improve the behavior of heavy-duty vehicles by keeping the whole vehicle body in the center of the lane.

Deep learning is a topic of considerable interest today. Since it deals with estimating - or learning - models, there are connections to the area of System Identification developed in the Automatic Control community. Such connections are explored and exploited in this contribution. It is stressed that common deep nets such as feedforward and cascadeforward nets are nonlinear ARX (NARX) models, and can thus be easily incorporated in System Identification code and practice. The case of LSTM nets is an example of NonLinear State-Space (NLSS) models. Copyright (C) 2020 The Authors.

For robust Gaussian process regression problems where the measurements are contaminated by outliers, a likelihood/measurement noise model with heavy-tailed distributions should be used to improve the prediction performance. In this paper, we propose to use a G-confluent distribution as the measurement noise model and a coordinate ascent variational inference method to infer the overall statistical model. In contrast with the commonly used Students t distribution, the G-confluent distribution can also be written as a Gaussian scale mixture, but its inverse scale follows a Beta distribution rather than a Gamma distribution, and its main advantage is that it is more flexible for modeling outliers while being equally suitable for variational inference. Numerical simulations based on benchmark data show that the G-confluent distribution performs better than or as well as the Students t distribution. Copyright (C) 2020 The Authors.

The Covid-19 pandemic has spawned numerous dynamic modeling attempts aimed at estimation, prediction, and ultimately control. The predictive power of these attempts has varied, and there remains a lack of consensus regarding the mechanisms of virus spread and the effectiveness of various non-pharmaceutical interventions that have been enforced regionally as well as nationally. Setting out in data available in the spring of 2020, and with a nowfamous model by Imperial College researchers as example, we employ an information-theoretical approach to shed light on why the predictive power of early modeling approaches have remained disappointingly poor. Copyright (C) 2020 The Authors.

This paper presents an optimization-based receding horizon trajectory planning algorithm for dynamical systems operating in unstructured and cluttered environments. The proposed approach is a two-step procedure that uses a motion planning algorithm in a first step to efficiently find a feasible, but possibly suboptimal, nominal solution to the trajectory planning problem where in particular the combinatorial aspects of the problem are solved. The resulting nominal trajectory is then improved in a second optimization-based receding horizon planning step which performs local trajectory refinement over a sliding time window. In the second step, the nominal trajectory is used in a novel way to both represent a terminal manifold and obtain an upper bound on the cost-to-go online. This enables the possibility to provide theoretical guarantees in terms of recursive feasibility, objective function value, and convergence to the desired terminal state. The established theoretical guarantees and the performance of the proposed algorithm are verified in a set of challenging trajectory planning scenarios for a truck and trailer system.   

This paper addresses the problem of change detection for a quadcopter in the presence ofwind disturbances. Different aspects of the quadcopter dynamics and various flight conditions have beeninvestigated. First, the wind is modeled using the Dryden wind model as a sum of a low-frequent and aturbulent part. Since the closed-loop control can compensate for system changes and disturbances andthe effect of the wind disturbance is significant, the residuals obtained from a standard simulation modelcan be misleading. Instead, a sensor-to-sensor submodel of the quadcopter is selected to detect a changein the payload using the Instrumental Variables (IV) cost function. It is shown that the mass variationcan be detected using the IV cost function in different flight scenarios.

An approach for belief space planning is presented, where knowledge about the landmark density is used as prior, instead of explicit landmark positions.

Having detailed maps of landmark positions  in a previously unvisited environment is considered unlikely in practice. Instead, it is argued that landmark densities should be used, as they could be estimated from other sources, such as ordinary maps or aerial imagery.

It is shown that it is possible to use virtual landmarks to approximate the landmark density to solve the presented problem. This approximation is also shown to give small errors during evaluation.

The approach is tested in a simulated environment, in conjunction with an extended information filter (EIF), where the computed path is shown to be superior compared to other alternative paths used as benchmarks.

An inference method for Gaussian process augmented state-space models are presented. This class of grey-box models enables domain knowledge to be incorporated in the inference process to guarantee a minimum of performance, still they are flexible enough to permit learning of partially unknown model dynamics and inputs. To facilitate online (recursive) inference of the model a sparse approximation of the Gaussian process based upon inducing points is presented. To illustrate the application of the model and the inference method, an example where it is used to track the position and learn the behavior of a set of cars passing through an intersection, is presented. Compared to the case when only the state-space model is used, the use of the augmented state-space model gives both a reduced estimation error and bias.

A wearable microphone array platform is used tolocalize stationary sound sources and amplify the sound inthe desired directions using several beamforming methods. Theplatform is equipped with inertial sensors and a magnetometerallowing predictions of source locations during orientationchanges and compensation for the displacement in the arrayconfiguration. The platform is modular, open and 3D printedto allow for easy reconfiguration of the array and for reuse inother applications, e.g., mobile robotics. The software componentsare based on open source. A new method for source localizationand signal reconstruction using Taylor expansion of the signals isproposed. This and various standard and non-standard Directionof Arrival (DOA) methods are evaluated in simulation andexperiments with the platform to track and reconstruct multipleand single sources. Results show that sound sources can belocalized and tracked robustly and accurately while rotating theplatform and that the proposed method outperforms standardmethods at reconstructing the signals.

Search-And-Rescue (SAR) is one of manyfields with applications benefiting from the increasingavailability of Unmanned Aerial Systems (UASs).  Most UAS applications rely on the UAS’s capabilityto carry a camera and stream video data for manualor automated processing. However, this relies onunobstructed views of the target, which limits the applicability of these systems.  In this paper, we instead describe the development and initial application testing of a system with a UAS-carried harmonic radar. This sensor is designed to detect the presence of Recco radar reflectors, commonly found integrated into alpine clothes and gear. The reflectors can be detected through vegetation and snow and is independent of many external factors such as lighting conditions. The paper describesthe system design and provides initial real-world results. The initial tests show fruitful results and opens up several avenues of continued research and development.

GNSS receivers are vulnerable to spoofing attacks, where false satellite signals are transmitted to trick the receiver to provide false position and/or time estimates. Novel algorithms are proposed for spoofing mitigation by exchanging double differences of pseudorange, or carrier phase, measurements between multiple GNSS receivers. In scenarios where the spoofing system utilizes a single transmit antenna, the pseudorange, and carrier phase, measurements that are associated with the spoofing signal can be detected and removed. Simulated meaconing attacks generated with a Spirent hardware simulator and measurements obtained with a modified version of GNSS-SDR are used to evaluate the proposed algorithms. Spoofing mitigation using pseudorange measurements is possible, for receivers that are separated at least five meters apart. With a receiver separation of 20 meters, the pseudorange double difference algorithm is able to correctly authenticate at least six of seven pseudoranges within 30 seconds. The carrier phase approach enables mitigation of spoofing signals at shorter receiver distances. However, this approach requires a more accurate time synchronization between the receivers.

In almost all parts of the world, breast cancer is one of the major causes of death among women. But at the same time, it is one of the most curable cancers if it is diagnosed at early stage. This paper tries to find a model that diagnose and classify breast cancer with high accuracy and help to both patients and doctors in the future. Here we develop a model using Normalized Multi Layer Perceptron Neural Network to classify breast cancer with high accuracy. The results achieved is very good (accuracy is 99.27%). It is very promising result compared to previous researches where Artificial Neural Networks were used. As benchmark test, Breast Cancer Wisconsin (Original) was used.

In this paper, we analyze a Linear Quadratic (LQ) control problem in terms of the average cost and the structure of the value function. We develop a completely model-free reinforcement learning algorithm to solve the LQ problem. Our algorithm is an off-policy routine where each policy is greedy with respect to all previous value functions. We prove that the algorithm produces stable policies given that the estimation errors remain small. Empirically, our algorithm outperforms the classical Q and off-policy learning routines.

Signed Laplacian matrices generally fail to be diagonally dominant and may fail to be stable. For both undirected and directed graphs, in this paper we present conditions guaranteeing the stability of signed Laplacians based on the property of eventual positivity, a Perron-Frobenius type of property for signed matrices. Our conditions are necessary and sufficient for undirected graphs, but only sufficient for digraphs, the gap between necessity and sufficiency being filled by matrices who have this Perron-Frobenius property on the right but not on the left side (i.e., on the transpose). An exception is given by weight balanced signed digraphs, where eventual positivity corresponds to positive semidefinitness of the symmetric part of the Laplacian. Analogous conditions are obtained for signed stochastic matrices.

We propose a generalization of the popular nonlinear ARX model structure by treating its parameters as varying over time. The parameters are considered generated by linear filters operating on the models regressors. The filters are expressed as a sum of atoms that are either sum of damped exponentials and sinusoids, or sinusoids with time varying frequencies. This form allows us to enforce stability of the parameter evolution as well as leverage the atomic norm minimization framework for inducing sparsity. It also facilitates easy incorporation of smoothness related priors that that making it possible to treat these models as nonlinear extensions of the familiar LPV models. (C) 2019, IFAC (International Federation of Automatic Control) Hosting by Elsevier Ltd. All rights reserved.

In urban environments, cellular network-based positioning of user equipment (UE) is a challenging task, especially in frequently occurring non-line-of-sight (NLOS) conditions. This paper investigates the use of two machine learning methods – neural networks and random forests – to estimate the position of UE in NLOS using best received reference signal beam power measurements. We evaluated the suggested positioning methods using data collected from a fifth-generation cellular network (5G) testbed provided by Ericsson. A statistical test to detect NLOS conditions with a probability of detection that is close to 90% is suggested. We show that knowledge of the antenna are crucial for accurate position estimation. In addition, our results show that even with a limited set of training data and one 5G transmission point, it is possible to position UE within 10 meters with 80% accuracy.

The CDIO Syllabus survey has successfully been applied to the Bachelor’s and Master’s programs in Experimental and Medical Biosciences, within the Faculty of Medicine and Health Sciences at Linköping University, Sweden. The programs are and have been, subject to considerable redesign with strong influence from the CDIO framework. One of the main drivers for the redesign is a shift concerning the main job market after graduation, from an academic career to industry and healthcare. One of the steps in the development process has been to carry out a Syllabus survey based on an adapted version of the CDIO Syllabus. The survey was sent out to students and to various categories of professionals, and in total 87 responses were received. The adapted version of the Syllabus and the design, execution, and outcome of the survey is presented.

The CDIO framework is an integrated and important part of the new quality assurance system within the Faculty of Science and Engineering at Linköping University. Both the CDIO Syllabus and the CDIO Standards are used extensively in the system. First, the paper presents the development and use of the second generation of course matrices (previously denoted ITU-matrices) and program matrices, which build upon an adapted and extended version of the CDIO Syllabus. The extension is made to also include bachelor’s and master’s program in subjects outside the engineering field. Second, the paper presents how the CDIO Standards are used in the quality reports, which are vital parts of the quality assurance systems. As a result, the CDIO framework is used for the design, management, and quality assurance of all education programs ( approximately 60 programs) within the Faculty of Science and Engineering at Linköping University.

In this paper, tightness relations (or inequalities) between Bayesian lower bounds (BLBs) on the mean-squared-error are derived which result from the marginalization of a joint probability density function (pdf) depending on both parameters of interest and extraneous or nuisance parameters. In particular,it is shown that for a large class of BLBs, the BLB derived from the marginal pdf is at least as tight as the corresponding BLB derived from the joint pdf. A Bayesian linear regression example is used to illustrate the tightness relations

This paper considers the problem of gathering information about features of interest in adversarial environments using mobile robots equipped with sensors. The problem is formulated as an informative path planning problem where the objective is to maximize the gathered information while minimizing the tracking performance of the adversarial observer. The optimization problem, that at first glance seems intractable to solve to global optimality, is shown to be equivalent to a mixed-integer semidefinite program that can be solved to global optimality using off-the-shelf optimization tools.

This paper addresses the problem of retrieving consistentestimates in a distributed network where the communication between the nodes is constrained such that only the diagonal elements of the covariance matrix are allowed to be exchanged. Several methods are developed for preserving and/or recovering consistency under the constraints imposed by the communication protocol. The proposed methods are used in conjunction with the covariance intersection method and the estimation performance is evaluated based on information usage and consistency. The results show that among the proposed methods, consistency can be preserved equally well at the transmitting node as at the receiving node.

This paper investigates the problem of controlling a complex network with reduced control energy. Two centrality measures are defined, one related to the energy that a control, placed on a node, can exert on the entire network, and the other related to the energy that the network exerts on a node. We show that by combining these two centrality measures conflicting control energy requirements, like minimizing the average energy needed to steer the state in any direction and the energy needed for the worst direction, can be simultaneously taken into account. From an algebraic point of view, the node ranking that we obtain from the combination of our centrality measures is related to the non-normality of the adjacency matrix of the graph.

This paper proposes a method to generate informative trajectories for a mobile sensor that tracks agile targets.The goal is to generate a sensor trajectory that maximizes the tracking performance, captured by a measure of the covariance matrix of the target state estimate. The considered problem is acombination of estimation and control, and is often referred to as informative path planning (IPP). When using nonlinear sensors, the tracking performance depends on the actual measurements, which are naturally unavailable in the planning stage.The planning problem hence becomes a stochastic optimization problem, where the expected tracking performance is used inthe objective function. The main contribution of this work is anapproximation of the problem based on deterministic sampling of the predicted target distribution. This is in contrast to prior work, where only the most likely target trajectory is considered.It is shown that the proposed method greatly improves the ability to track agile targets, compared to a baseline approach.   

By arranging a large number of inertial sensors in an array and fusing their measurements, it is possible to create inertial sensor assemblies with a high performance-to-price ratio. Recently, a maximum likelihood estimator for fusing inertial array measurements collected at a given sampling instance was developed. In this paper, the maximum likelihood estimator is extended by introducing a motion model and deriving a maximum a posteriori estimator that jointly estimates the array dynamics at multiple sampling instances. Simulation examples are used to demonstrate that the proposed sensor fusion method have the potential to yield significant improvements in estimation accuracy. Further, by including the motion model, we resolve the sign ambiguity of gyro-free implementations, and thereby open up for implementations based on accelerometer-only arrays.

We present a technique for kernel-based identification of Wiener systems. We model the impulse response of the linear block with a Gaussian process. The static nonlinearity is modeled with a combination of basis functions. The coefficients of the static nonlinearity are estimated, together with the hyperparameters of the covariance function of the Gaussian process model, using an iterative algorithm based on the expectation-maximization method combined with elliptical slice sampling to sample from the posterior distribution of the impulse response given the data. The same sampling method is then used to find the posterior-mean estimate of the impulse response. We test the proposed algorithm on a benchmark of randomly-generated Wiener systems.

In Massive MIMO, the pilot contamination effect reduces the spectral efficiency (SE) gains and superimposed pilot (SP) transmission has been proposed to mitigate this effect. SP is based on transmitting pilot and data symbols simultaneously to allow for longer pilots and no pilot overhead. This work studies the optimal power control strategies in the uplink of a Massive MIMO system with SP and detection based on maximum ratio combining The optimization objectives arc maximum product of SINRs and max-min fairness, and these problems are reformulated as geometric programs which allow for efficient implementations. The numerical results indicate that the SE gains from the optimal power control with respect to the heuristic statistical channel inversion power control, are more significant when the interference from pilot symbols is subtracted before data detection.

The main idea in this paper is to implement a distributed primal-dual interior-point algorithm for loosely coupled Quadratic Programming problems. We implement this in Julia and show how can we exploit parallelism in order to increase the computational speed. We investigate the performance of the algorithm on a Model Predictive Control problem.

In this work we consider a nonlinear interconnected system describing a decision-making process in a community of agents characterized by the coexistence of collaborative and antagonistic interactions. The resulting signed graph is in general not structurally balanced. It is shown in the paper that the decision-making process is affected by the frustration of the signed graph, in the sense that a nontrivial decision can be reached only if the social commitment of the agents is high enough to win the disorder introduced by the frustration in the network. The higher the frustration of the graph, the higher the commitment strength required from the agents.

We explore the problem of identification of LPV models when the scheduling variables are not known in advance and the model parameters exhibit a dynamic dependence on them. We consider an affine ARX model structure whose parameters vary with time. We solve for the models parameters and scheduling variables in two steps. In the first step, we use the measured input-output data to realize a parameter trajectory by solving a regularized Hankel matrix rank minimization problem. The regularization penalty is guided by the prior knowledge regarding the nature of systems time variation. In the second step, the scheduling variables are estimated as parameters of a sparse ARX structure relating the models parameters to the measured input-output variables. The effectiveness of the proposed approach is illustrated with two practical examples. (C) 2018, IFAC (International Federation of Automatic Control) Hosting by Elsevier Ltd. All rights reserved.

The demand for fuel-efficient transport solutions are steadily increasing with the goal of reducing environmental impact and increasing efficiency. Heavy-Duty Vehicle (HDV) platooning is a promising concept where multiple HDVs drive together in a convoy with small intervehicular spacing. By doing this, the aerodynamic drag is reduced which in turn lowers fuel consumption. We propose a novel Model Predictive Control (MPC) framework for longitudinal control of the follower vehicle in a platoon consisting of two HDVs when no vehicle-to-vehicle communication is available. In the framework, the preceding vehicle's velocity profile is predicted using artificial neural networks which uses a topographic map of the road as input and is trained offline using synthetic data. The gear shifting and mass of consumed fuel for the controlled follower vehicle is modeled and used within the MPC controller. The efficiency of the proposed framework is verified in simulation examples and is benchmarked with a currently available control solution.  

In this work, we study the multiple kernel based regularized system identification with the hyper-parameter estimated by using the Steins unbiased risk estimators (SURE). To approach the problem, a QR factorization is first employed to compute SUREs objective function and its gradient in an efficient and accurate way. Then we propose an algorithm to solve the SURE problem, which contains two parts: the outer optimization part and the inner optimization part. For the outer optimization part, the coordinate descent algorithm is used and for the inner optimization part, the projection gradient algorithm is used. Finally, the efficacy of the proposed algorithm is demonstrated by numerical simulations. (C) 2018, IFAC (International Federation of Automatic Control) Hosting by Elsevier Ltd. All rights reserved.

In this paper, we study the asymptotic properties of the generalized cross validation (GCV) hyperparameter estimator and establish its connection with the Steins unbiased risk estimators (SURE) as well as the mean squared error (MSE). It is shown that as the number of data goes to infinity, the GCV has the same asymptotic property as the SURE does and both of them converge to the best hyperparameter in the MSE sense. We illustrate the efficacy of the result by Monte Carlo simulations. (C) 2018, IFAC (International Federation of Automatic Control) Hosting by Elsevier Ltd. All rights reserved.

Because of the increased demand on fault detection, monitoring and predictive maintenance, online or recursive identification is playing a more important role in systems engineering. In the recent releases of System Identification Toolbox (TM) for MATLAB (R), this has been reflected in a more substantial support for online techniques. This contribution gives an account of these improvements. It covers the addition of nonlinear filtering algorithms, such as the extended Kalman filter, the unscented Kalman filter and particle filters. The traditional recursive estimation techniques for polynomial models have also been enhanced with a more versatile syntax. Several new Simulink (R) blocks have been developed to support Simulink (R) models with online estimation. (C) 2018, IFAC (International Federation of Automatic Control) Hosting by Elsevier Ltd. All rights reserved.

Using Maximum Likelihood (or Prediction Error) methods to identify linear state space model is a prime technique. The likelihood function is a nonconvex function and care must be exercised in the numerical maximization. Here the focus will be on affine parameterizations which allow some special techniques and algorithms. Three approaches to formulate and perform the maximization are described in this contribution: (1) The standard and well known Gauss Newton iterative search, (2) a scheme based on the EM (expectation-maximization) technique, which becomes especially simple in the affine parameterization case, and (3) a new approach based on lifting the problem to a higher dimension in the parameter space and introducing rank constraints. (C) 2018, IFAC (International Federation of Automatic Control) Hosting by Elsevier Ltd. All rights reserved.

This paper concerns the estimation of a dynamic model from two measured signals when it is not clear which signal should be used as input to the model. In this case, both a forward and an inverse model can be estimated. Here, a basic instrumental variable approach is used and it is shown that the forward and inverse model estimators give identical parameter estimates provided that corresponding model structures have been used. Furthermore, it is shown that this scenario occurs when properties of a quadcopter are estimated from accelerometer and gyro signals and, hence, that it does not matter which signal is used as input.

In this paper we show how to synthesize simple explicit MPC controllers based on approximate dynamic programming. Here, a given MPC optimization problem over a finite horizon is solved iteratively as a series of problems of size one. The optimal cost function of each subproblem is approximated by a quadratic function that serves as a cost-to-go function for the subsequent iteration. The approximation is designed in such a way that closed-loop stability and recursive feasibility is maintained. Specifically, we show how to employ sum-of-squares relaxations to enforce that the approximate cost-to-go function is bounded from below and from above for all points of its domain. By resorting to quadratic approximations, the complexity of the resulting explicit MPC controller is considerably reduced both in terms of memory as well as the on-line computations. The procedure is applied to control an inverted pendulum and experimental data are presented to demonstrate viability of such an approach.

In this paper, Bayesian lower bounds (BLBs) are obtained via a general form of the Pythagorean theorem where the inner product derives from the joint or the a-posteriori probability density function (pdf). When joint pdf is considered, the BLBs obtained encompass the Weiss-Weinstein family (WWF). When a-posteriori pdf is considered, by resorting to an embedding between two ad hoc subspaces, it is shown that any ”standard” BLBs of the WWF admits a ”tighter” form which upper bounds the ”standard” form. Interestingly enough, this latter result may explain why the ”standard” BLBs of the WWF are not always as tight as expected, as exemplified in the case of the Bayesian Cram´er-Rao Bound. As a consequence an updated definition of efficiency is proposed, as well as the introduction of an updated class of efficient estimators.

This paper presents a systematic approach for computing local solutions to motion planning problems in non-convex environments using numerical optimal control techniques. It extends the range of use of state-of-the-art numerical optimal control tools to problem classes where these tools have previously not been applicable. Today these problems are typically solved using motion planners based on randomized or graph search. The general principle is to define a homotopy that transforms, or preferably relaxes, the original problem to an easily solved problem. In this work, it is shown that by combining a Sequential Quadratic Programming (SQP) method with a homotopy approach that gradually transforms the problem from a relaxed one to the original one, practically relevant locally optimal solutions to the motion planning problem can be computed. The approach is demonstrated in motion planning problems in challenging 2D and 3D environments, where the presented method significantly outperforms both a state-of-the-art numerical optimal control method and a state-of-the-art open-source optimizing sampling-based planner commonly used as benchmark. 

An experimental comparison of two feed-forward based frictioncompensation methods is presented. The first method is based on theLuGre friction model, using identified friction model parameters, andthe second method is based on B-spline network, where the networkweights are learned from experiments. The methods are evaluated andcompared via experiments using a six axis industrial robot carryingout circular movements of different radii. The experiments show thatthe learning-based friction compensation gives an error reduction ofthe same magnitude as for the LuGre-based friction compensation.

The number of sensors used in tracking scenariosis constantly increasing, this puts high demands on the trackingmethods to handle these data streams. Central processing (ideallyoptimal) puts high demands on the central node, is sensitive toinaccurate sensor parameters, and suffers from the single pointof failure problem. Decentralizing the tracking can improve this,but may give considerable performance loss. The newly presentedinverse covariance intersection method, proven to be consistent,even under unknown track cross correlations, is benchmarkedagainst alternatives. Different track-to-track methods, includingsmoothed association over a window, are compared. A scenariowith objects tracked in multiple cameras, not necessarily opti-mized for tracking, are used to give realism to the evaluations.

In this paper, marginal versions of the Bayesian Bhattacharyya lower bound (BBLB), which is a tighter alternative to the classical Bayesian Cramer-Rao bound, for discrete-time filtering are proposed. Expressions for the second and third-order marginal BBLBs are obtained and it is shown how these can be approximately calculated using particle filtering. A simulation example shows that the proposed bounds predict the achievable performance of the filtering algorithms better.