Adaptive Markov chain Monte Carlo and Bayesian filtering for state space models

This thesis is concerned with the state and parameter estimation in state space models. The estimation of states and parameters is an important task when mathematical modeling is applied to many different application areas such as the global positioning systems, target tracking, navigation, brain imaging, spread of infectious diseases, biological processes, telecommunications, audio signal processing, stochastic optimal control, machine learning, and physical systems. In Bayesian settings, the estimation of states or parameters amounts to computation of the posterior probability density function. Except for a very restricted number of models, it is impossible to compute this density function in a closed form. Hence, we need approximation methods. A state estimation problem involves estimating the states (latent variables) that are not directly observed in the output of the system. In this thesis, we use the Kalman filter, extended Kalman filter, Gauss–Hermite filters, and particle filters to estimate the states based on available measurements. Among these filters, particle filters are numerical methods for approximating the filtering distributions of non-linear non-Gaussian state space models via Monte Carlo. The performance of a particle filter heavily depends on the chosen importance distribution. For instance, inappropriate choice of the importance distribution can lead to the failure of convergence of the particle filter algorithm. In this thesis, we analyze the theoretical Lᵖ particle filter convergence with general importance distributions, where p ≥2 is an integer. A parameter estimation problem is considered with inferring the model parameters from measurements. For high-dimensional complex models, estimation of parameters can be done by Markov chain Monte Carlo (MCMC) methods. In its operation, the MCMC method requires the unnormalized posterior distribution of the parameters and a proposal distribution. In this thesis, we show how the posterior density function of the parameters of a state space model can be computed by filtering based methods, where the states are integrated out. This type of computation is then applied to estimate parameters of stochastic differential equations. Furthermore, we compute the partial derivatives of the log-posterior density function and use the hybrid Monte Carlo and scaled conjugate gradient methods to infer the parameters of stochastic differential equations. The computational efficiency of MCMC methods is highly depend on the chosen proposal distribution. A commonly used proposal distribution is Gaussian. In this kind of proposal, the covariance matrix must be well tuned. To tune it, adaptive MCMC methods can be used. In this thesis, we propose a new way of updating the covariance matrix using the variational Bayesian adaptive Kalman filter algorithm.

Spin Polarization and Normal Reflection at Nanoscale Point Contacts

This thesis presents point-contact measurements between superconductors (Nb, Ta, Sn,Al, Zn) and ferromagnets (Co, Fe, Ni) as well as non-magnetic metals (Ag, Au, Cu, Pt).The point contacts were fabricated using the shear method. The differential resistanceof the contacts was measured either in liquid He at 4.2 K or in vacuum in a dilutionrefrigerator at varying temperature down to 0.1 K. The contact properties were investigatedas function of size and temperature. The measured Andreev-reflection spectrawere analysed in the framework of the BTK model – a three parameter model that describescurrent transport across a superconductor - normal conductor interface. Theoriginal BTK model was modified to include the effects of spin polarization or finitelifetime of the Cooper pairs. Our polarization values for the ferromagnets at 4.2 K agree with the literature data, but the analysis was ambiguous because the experimental spectra both with ferromagnets and non-magnets could be described equally well either with spin polarization or finite lifetime effects in the BTK model. With the polarization model the Z parametervaries from almost 0 to 0.8 while the lifetime model produces Z values close to 0.5. Measurements at lower temperatures partly lift this ambiguity because the magnitude of thermal broadening is small enough to separate lifetime broadening from the polarization. The reduced magnitude of the superconducting anomalies for Zn-Fe contacts required an additional modification of the BTK model which was implemented as a scaling factor. Adding this parameter led to reduced polarization values. However, reliable data is difficult to obtain because different parameter sets produce almost identical spectra.

Valmistavan tehtaan tuotannon tehostaminen reaaliaikaisen tiedon avulla

Teollinen toimintaympäristö on muuttunut globalisaation seurauksena. Sisäisten ja ulkoisten asiakkaiden vaatimukset ovat entistä haastavampia toimitusajan, laadun ja kustannusten suhteen. Yritysten on automatisoitava toimintoja, jotta voidaan saavuttaa kilpailuetua ja parantaa yrityksen tehokkuutta. Diplomityön tarkoitus on kirjallisuus- ja vertailututkimuksen yhteenvetona koostaa asiakasyrityksen tarpeisiin soveltuva toimintomalli reaaliaikaisen tiedon hyödyntä-miseen valmistavan tehtaan tuotannon tehostamisessa. Kirjallisuustutkimusten avulla pyritään analysoimaan valmistavan teollisuuden erityispiirteet toimintaympäristön osalta, kartoittamaan soveltuvin teknologia haastavaan toimintaympäristöön sekä löytää keinot kuinka reaaliaikaista tietoa hyödynnetään tuotannon tehokkuuden toteutumisessa. Vertailututkimuksessa analysoidaan kahta reaaliaikaista tietoa toiminnoissaan hyödyntävää yritystä, joiden tuloksia verrataan kirjallisuustutkimuksen aineistoon. Työn soveltavassa osuudessa luotiin toimintomalli asiakasyrityksen tarpeisiin, huomioiden tehtaan toiminnot kokonaisuutena. Toimintomallin avulla yritys pystyy toteuttamaan reaaliaikaisen tiedon keräämisen sekä hyödyntämisen yrityksen toi-minnoista. Toimintomallin toimivuus todettiin soveltamalla pilot-projekti asiakasyrityksen tuotantoon ja analysoitiin saavutettuja tuloksia haastattelujen sekä kvalitatiivisten tulosten perusteella. Työn tutkimukset osoittivat haasteet valmistavan teollisuuden toimintaympäristön ja automaatioteknologian soveltamisessa. Yritysten tulee suunnitteella toimintamallin kokonaisuus ja pohtia järjestelmän toimivuutta ongelmatilanteissa. Reaaliaikaisen tiedonkeruun soveltaminen onnistuneesti mahdollistaa hyötyjä monella tapaa. Merkittävimmät ovat tehokkuuden kehittyminen toimintojen nopeutumisena sekä resurssisäästöt toimintojen automatisoinnin johdosta. Reaaliaikaista tietoa voidaan hyödyntää myös suorituskyvyn johtamisessa mittariston avulla.

Novel virtual environment and real-time simulation based methods for improving life-cycle efficiency of non-road mobile machinery

Virtual environments and real-time simulators (VERS) are becoming more and more important tools in research and development (R&D) process of non-road mobile machinery (NRMM). The virtual prototyping techniques enable faster and more cost-efficient development of machines compared to use of real life prototypes. High energy efficiency has become an important topic in the world of NRMM because of environmental and economic demands. The objective of this thesis is to develop VERS based methods for research and development of NRMM. A process using VERS for assessing effects of human operators on the life-cycle efficiency of NRMM was developed. Human in the loop simulations are ran using an underground mining loader to study the developed process. The simulations were ran in the virtual environment of the Laboratory of Intelligent Machines of Lappeenranta University of Technology. A physically adequate real-time simulation model of NRMM was shown to be reliable and cost effective in testing of hardware components by the means of hardware-in-the-loop (HIL) simulations. A control interface connecting integrated electro-hydraulic energy converter (IEHEC) with virtual simulation model of log crane was developed. IEHEC consists of a hydraulic pump-motor and an integrated electrical permanent magnet synchronous motorgenerator. The results show that state of the art real-time NRMM simulators are capable to solve factors related to energy consumption and productivity of the NRMM. A significant variation between the test drivers is found. The results show that VERS can be used for assessing human effects on the life-cycle efficiency of NRMM. HIL simulation responses compared to that achieved with conventional simulation method demonstrate the advances and drawbacks of various possible interfaces between the simulator and hardware part of the system under study. Novel ideas for arranging the interface are successfully tested and compared with the more traditional one. The proposed process for assessing the effects of operators on the life-cycle efficiency will be applied for wider group of operators in the future. Driving styles of the operators can be analysed statistically from sufficient large result data. The statistical analysis can find the most life-cycle efficient driving style for the specific environment and machinery. The proposed control interface for HIL simulation need to be further studied. The robustness and the adaptation of the interface in different situations must be verified. The future work will also include studying the suitability of the IEHEC for different working machines using the proposed HIL simulation method.

Tuotekustannuslaskennan nykytila ja kustannusseurannan kehittäminen mallimoottoriajattelun avulla

Tutkimuksen tarkoituksena on tutkia globaalin konsernin yhden liiketoimintayksikön tuotekustannuslaskennan nykytilaa. Lisäksi tutkimuksessa selvitetään, miten tuotekohtaista kustannusseurantaa voidaan kehittää mallimoottoriajatuksen avulla. Tutkimus on toteutettu laadullisena case-tutkimuksena yhden organisaation tietojen pohjalta. Teoriaosuuden lähdeaineistot koostuvat pääosin kustannuslaskennan ja -johtamisen perusteoksista ja tieteellisistä artikkeleista. Empiriaosuuden tiedot pohjautuvat haastatteluihin, tietojärjestelmiin ja tutustumiseen organisaatioon. Tutkimuksessa selvisi, että liiketoimintayksikkö ei tällä hetkellä seuraa tuotekohtaisia kustannuksia yksittäisten tuotteiden tasolla. Kustannusseuranta tapahtuu sen sijaan suurempien kokonaisuuksien keskimääräisten kustannuksien tasolla. Tuotekustannuslaskenta on toteutettu perinteiseksi menetelmäksi luokiteltavalla laskentatavalla, jossa välilliset kustannukset kohdistetaan yleiskustannuslisäprosenttien avulla. Tutkimuksen perusteella yleiskustannuksien kohdistamisperusteissa on havaittavissa viitteitä kustannuksien vääristymisestä. Tuotetason kustannuksien seurantaan kehitettiin mallimoottoriajatukseen pohjautuva kustannusmalli, jonka avulla seurataan tarkasti valikoitujen tuotteiden kustannuksien kehittymistä sekä kustannusrakennetta. Mallin avulla voidaan lisätä tuotetason kustannustietoisuutta liiketoimintayksikössä sekä tehdä havaintoja tuotekohtaisten kustannuksien kehityssuunnasta. Mallin kustannustietona käytetään olemassa olevan kustannuslaskentajärjestelmän tietoja. Tästä johtuen mallin kustannustiedoissa on havaittavissa myös viitteitä kustannuksien vääristymisestä.

Analytic evaluation of three-phase short circuit demagnetization and hysteresis loss risk in rotor-surface-magnet Permanent-Magnet Synchronous Machines

Permanent magnet synchronous machines (PMSM) have become widely used in applications because of high efficiency compared to synchronous machines with exciting winding or to induction motors. This feature of PMSM is achieved through the using the permanent magnets (PM) as the main excitation source. The magnetic properties of the PM have significant influence on all the PMSM characteristics. Recent observations of the PM material properties when used in rotating machines revealed that in all PMSMs the magnets do not necessarily operate in the second quadrant of the demagnetization curve which makes the magnets prone to hysteresis losses. Moreover, still no good analytical approach has not been derived for the magnetic flux density distribution along the PM during the different short circuits faults. The main task of this thesis is to derive simple analytical tool which can predict magnetic flux density distribution along the rotor-surface mounted PM in two cases: during normal operating mode and in the worst moment of time from the PM’s point of view of the three phase symmetrical short circuit. The surface mounted PMSMs were selected because of their prevalence and relatively simple construction. The proposed model is based on the combination of two theories: the theory of the magnetic circuit and space vector theory. The comparison of the results in case of the normal operating mode obtained from finite element software with the results calculated with the proposed model shows good accuracy of model in the parts of the PM which are most of all prone to hysteresis losses. The comparison of the results for three phase symmetrical short circuit revealed significant inaccuracy of the proposed model compared with results from finite element software. The analysis of the inaccuracy reasons was provided. The impact on the model of the Carter factor theory and assumption that air have permeability of the PM were analyzed. The propositions for the further model development are presented.

Virtuaaliprototyypin käyttö roottorin dynamiikan analysoinnissa

Työn tavoitteena oli selvittää kaupallisen dynamiikansimulointiohjelmiston so-veltuvuus roottoridynamiikan analysointiin. Työssä keskityttiin erityisesti rootto-rin dynamiikkaan vaikuttavien epäideaalisuuksien mallintamiseen. Simulointitu-losten tarkkuutta selvitettiin mittauksilla. Lisäksi vertailtiin yleiskäyttöisen dyna-miikan simulointiohjelmiston ja roottoridynamiikan erikoisohjelmiston teoriaa. Tutkittava roottori oli paperikoneen putkitela. Telan joustavuus kuvattiin ele-menttimenetelmällä ratkaistujen moodien avulla. Elementtimallissa huomioitiin telan vaipan seinämänpaksuusvaihtelu, joka vaikuttaa telan massa- ja jäykkyysja-kaumaan. Dynamiikkaohjelmistossa mallinnettiin telan tuennasta tulevat herätteet. Dynamiikkaohjelmistona käytettiin ADAMS:ia ja FEM-ohjelmana ANSYS:stä. Tuloksista havaittiin käytetyn menetelmän soveltuvan roottoridynamiikan ana-lysointiin ja roottorin epäideaalisuuksien mallintamiseen. Simulointimallilla saa-tiin esille murtolukukriittiset pyörimisnopeudet ja telan kriittinen pyörimisnopeus vastasi hyvin mittaustuloksia.

Development and Effect Analysis of the Asteri Consultative Auditing Process - Safety and Security Management in Educational Institutions

The Finnish legislation requires for a safe and secure learning environment. However, the comprehensive, risk based safety and security management (SSM) and the management commitment in the implementation and development of the SSM are not mentioned in the legislation. Multiple institutions, operators and researchers have studied and developed safety and security in educational institutions over the past decade. Typically the approach has been fragmented and without bringing up the importance of the comprehensive SSM. The development needs of the safety and security operations in universities have been studied. However, in universities of applied sciences (UASs) and in elementary schools (ESs), the performance level, strengths and weaknesses of the comprehensive SSM have not been studied. The objective of this study was to develop the comprehensive, risk based SSM of educational institutions by developing the new Asteri consultative auditing process and study its effects on auditees. Furthermore, the performance level in the comprehensive SSM in UASs and ESs were studied using Asteri and the TUTOR model developed by the Keski-Uusimaa Department for Rescue Services. In addition, strengths, development needs and differences were identified. In total, 76 educational institutions were audited between the years 2011 and 2014. The study is based on logical empiricism, and an observational applied research design was used. Auditing, observation and an electronic survey were used for data collection. Statistical analysis was used to analyze the collected information. In addition, thematic analysis was used to analyze the development areas of the organizations mentioned by the respondents in the survey. As one of the main contributions, this research presents the new Asteri consultative auditing process. Organizations with low performance levels on the audited subject benefit the most from the Asteri consultative auditing process. Asteri may be usable in many different types of audits, not only in SSM audits. As a new result, this study provides new knowledge on attitudes related to auditing. According to the research findings, auditing may generate negative attitudes and the auditor should take them into account when planning and preparing for audits. Negative attitudes can be compensated by producing added value, objectivity and positivity for the audit and, thus, improve the positive effects of auditing on knowledge and skills. Moreover, as the results of this study shows, auditing safety and security issues do not increase feelings of insecurity, but rather increase feelings of safety and security when using the new Asteri consultative auditing process with the TUTOR model. The results showed that the SSM in the audited UASs was statistically significantly more advanced than that in the audited ESs. However, there is still room for improvement in the ESs and the UASs as the approach to the SSM was fragmented. It can be assumed that the majority of Finnish UASs and ESs do not likely meet the basic level of the comprehensive, risk based the SSM.