Sensitivity analysis in optimized parametric curve fitting

Purpose – Curve fitting from unordered noisy point samples is needed for surface reconstruction in many applications -- In the literature, several approaches have been proposed to solve this problem -- However, previous works lack formal characterization of the curve fitting problem and assessment on the effect of several parameters (i.e. scalars that remain constant in the optimization problem), such as control points number (m), curve degree (b), knot vector composition (U), norm degree (k), and point sample size (r) on the optimized curve reconstruction measured by a penalty function (f) -- The paper aims to discuss these issues -- Design/methodology/approach - A numerical sensitivity analysis of the effect of m, b, k and r on f and a characterization of the fitting procedure from the mathematical viewpoint are performed -- Also, the spectral (frequency) analysis of the derivative of the angle of the fitted curve with respect to u as a means to detect spurious curls and peaks is explored -- Findings - It is more effective to find optimum values for m than k or b in order to obtain good results because the topological faithfulness of the resulting curve strongly depends on m -- Furthermore, when an exaggerate number of control points is used the resulting curve presents spurious curls and peaks -- The authors were able to detect the presence of such spurious features with spectral analysis -- Also, the authors found that the method for curve fitting is robust to significant decimation of the point sample -- Research limitations/implications - The authors have addressed important voids of previous works in this field -- The authors determined, among the curve fitting parameters m, b and k, which of them influenced the most the results and how -- Also, the authors performed a characterization of the curve fitting problem from the optimization perspective -- And finally, the authors devised a method to detect spurious features in the fitting curve -- Practical implications – This paper provides a methodology to select the important tuning parameters in a formal manner -- Originality/value - Up to the best of the knowledge, no previous work has been conducted in the formal mathematical evaluation of the sensitivity of the goodness of the curve fit with respect to different possible tuning parameters (curve degree, number of control points, norm degree, etc.)

Characterization of micro-bubble size distribution and flow configuration in DAF contact zone by a non-intrusive image analysis system and tracer tests

This paper aims to investigate the influence of some dissolved air flotation (DAF) process variables (specifically: the hydraulic detention time in the contact zone and the supplied dissolved air concentration) and the pH values, as pretreatment chemical variables, on the micro-bubble size distribution (BSD) in a DAF contact zone. This work was carried out in a pilot plant where bubbles were measured by an appropriate non-intrusive image acquisition system. The results show that the obtained diameter ranges were in agreement with values reported in the literature (10-100mm), quite independently of the investigated conditions. The linear average diameter varied from 20 to 30mm, or equivalently, the Sauter (d(3,2)) diameter varied from 40 to 50mm. In all investigated conditions, D(50) was between 75% and 95%. The BSD might present different profile (with a bimodal curve trend), however, when analyzing the volumetric frequency distribution (in some cases with the appearance of peaks in diameters ranging from 90-100mm). Regarding volumetric frequency analysis, all the investigated parameters can modify the BSD in DAF contact zone after the release point, thus potentially causing changes in DAF kinetics. This finding prompts further research in order to verify the effect of these BSD changes on solid particle removal efficiency by DAF.

Evaluation of flow regimes in a semi-cylindrical spouted bed through statistical, mutual information, spectral and Hurst`s analysis

Hydrodynamic studies were conducted in a semi-cylindrical spouted bed column of diameter 150 mm, height 1000 mm, conical base included angle of 60 degrees and inlet orifice diameter 25 mm. Pressure transducers at several axial positions were used to obtain pressure fluctuation time series with 1.2 and 2.4 mm glass beads at U/U-ms from 0.3 to 1.6, and static bed depths from 150 to 600 mm. The conditions covered several flow regimes (fixed bed, incipient spouting, stable spouting, pulsating spouting, slugging, bubble spouting and fluidization). Images of the system dynamics were also acquired through the transparent walls with a digital camera. The data were analyzed via statistical, mutual information theory, spectral and Hurst`s Rescaled Range methods to assess the potential of these methods to characterize the spouting quality. The results indicate that these methods have potential for monitoring spouted bed operation.

A longitudinal analysis of cytotoxic T lymphocyte precursor frequencies to the hepatitis B virus in chronically infected patients

Individuals with acute hepatitis B virus (HBV) infection characteristically mount a strong, multispecific cytotoxic T lymphocyte (CTL) response that is effective in eradicating virus. In contrast, this response in chronic carriers is usually weak or undetectable. Since it is generally acknowledged that HBV pathogenesis is immune-mediated, the occurrence of episodes of active liver disease in many carriers suggests that these individuals can mount active CTL responses to HBV. To see whether the detection of circulating CTLs is related to these flare episodes, we have determined the CTL precursor (CTLp) frequencies to HLA-A2-restricted viral peptides in seven patients over a 12-24-month period of their disease. Limiting dilution analyses (LDA) were performed longitudinally to five epitopes comprising the viral capsid (HBc), envelope (HBs) and polymerase (pol) proteins. Assays were performed against a mixture of peptides, or against each individual peptide, to measure overall CTL activity and the multispecificity of the responses, respectively. Since two of the patients were treated with recombinant human interleukin-12 (rHuIL-12) at the time, with one individual achieving complete disease remission a year later after being treated with interferon-alpha, we were also able to examine the effects of these cytokines on HBV cytotoxicity. Our results indicate that weak but detectable CTL responses do occur in chronic carriers which are generally associated with disease flares, although CTLps were also seen occasionally during minimal disease activity. The range of specificities varied between individuals and within each individual during the course of the disease. Finally, we also provide evidence that CTL reactivity is stimulated following treatment with certain cytokines, but is dependent on the time of administration.

Time-frequency and time-scale characterisation of the beat-by-beat high-resolution electrocardiogram

Proceedings of the Sixth Portuguese Conference on Bioemedical Engineering faro, Portugal

Quantifying consistency of Event Related Potentials across subjects based on single-trial topographic analysis

Introduction: Non-invasive brain imaging techniques often contrast experimental conditions across a cohort of participants, obfuscating distinctions in individual performance and brain mechanisms that are better characterised by the inter-trial variability. To overcome such limitations, we developed topographic analysis methods for single-trial EEG data [1]. So far this was typically based on time-frequency analysis of single-electrode data or single independent components. The method's efficacy is demonstrated for event-related responses to environmental sounds, hitherto studied at an average event-related potential (ERP) level. Methods: Nine healthy subjects participated to the experiment. Auditory meaningful sounds of common objects were used for a target detection task [2]. On each block, subjects were asked to discriminate target sounds, which were living or man-made auditory objects. Continuous 64-channel EEG was acquired during the task. Two datasets were considered for each subject including single-trial of the two conditions, living and man-made. The analysis comprised two steps. In the first part, a mixture of Gaussians analysis [3] provided representative topographies for each subject. In the second step, conditional probabilities for each Gaussian provided statistical inference on the structure of these topographies across trials, time, and experimental conditions. Similar analysis was conducted at group-level. Results: Results show that the occurrence of each map is structured in time and consistent across trials both at the single-subject and at group level. Conducting separate analyses of ERPs at single-subject and group levels, we could quantify the consistency of identified topographies and their time course of activation within and across participants as well as experimental conditions. A general agreement was found with previous analysis at average ERP level. Conclusions: This novel approach to single-trial analysis promises to have impact on several domains. In clinical research, it gives the possibility to statistically evaluate single-subject data, an essential tool for analysing patients with specific deficits and impairments and their deviation from normative standards. In cognitive neuroscience, it provides a novel tool for understanding behaviour and brain activity interdependencies at both single-subject and at group levels. In basic neurophysiology, it provides a new representation of ERPs and promises to cast light on the mechanisms of its generation and inter-individual variability.

Positive time-frequency distributions based on joint marginal constraints

This correspondence studies the formulation of members ofthe Cohen-Posch class of positive time-frequency energy distributions.Minimization of cross-entropy measures with respect to different priorsand the case of no prior or maximum entropy were considered. It isconcluded that, in general, the information provided by the classicalmarginal constraints is very limited, and thus, the final distributionheavily depends on the prior distribution. To overcome this limitation,joint time and frequency marginals are derived based on a "directioninvariance" criterion on the time-frequency plane that are directly relatedto the fractional Fourier transform.

Wigner higher-order spectra: definition, properties, computation and application to transient signal analysis

The Wigner higher order moment spectra (WHOS)are defined as extensions of the Wigner-Ville distribution (WD)to higher order moment spectra domains. A general class oftime-frequency higher order moment spectra is also defined interms of arbitrary higher order moments of the signal as generalizations of the Cohen’s general class of time-frequency representations. The properties of the general class of time-frequency higher order moment spectra can be related to theproperties of WHOS which are, in fact, extensions of the properties of the WD. Discrete time and frequency Wigner higherorder moment spectra (DTF-WHOS) distributions are introduced for signal processing applications and are shown to beimplemented with two FFT-based algorithms. One applicationis presented where the Wigner bispectrum (WB), which is aWHOS in the third-order moment domain, is utilized for thedetection of transient signals embedded in noise. The WB iscompared with the WD in terms of simulation examples andanalysis of real sonar data. It is shown that better detectionschemes can be derived, in low signal-to-noise ratio, when theWB is applied.

Immersioultraäänitutkimuksen kehittäminen

Teräksenvalmistajilta edellytetään jatkuvasti panostusta laadun ja laadunvarmistuksen kehittämiseen. Teräksen laatu ja puhtaus korostuvat varsinkin silloin, kun terästä käytetään vaativiin käyttökohteisiin, kuten autoteollisuuden tarpeisiin. Ultraäänitarkastusmenetelmää käytetään laadun-varmistuksessa teräksen sisävikojen etsimiseen. Ultraäänitarkastuksessa lähetetään suuritaajuuksista ääntä kappaleeseen. Ääni etenee materiaalissa ja heijastuu erilaisista epäjatkuvuuskohdista. Luotaimeen palaavaa ääntä analysoimalla saadaan tietoa teräksestä ja sen sisävioista. Ultraäänitarkastuksen ongelmana on vian tyypin määrittäminen hankaluus sekä herkkyys tutkittavan materiaalin aineominaisuuksille. Työn tavoitteena oli immersioultraäänitarkastuksen kehittäminen sovellettuna teräksenvalmistajan tarpeisiin. Materiaalin aineominaisuuksista tutkittiin seostuksen vaikutusta. Teräslajit tarkastettiin valssitilaisena, karkaistuna ja normalisoituna. Lisäksi tutkittiin kappaleen pinnankarheuden ja -muodon vaikutusta tarkastukseen. Vikatyyppien tunnistamisen mahdollisuuksia ultra-äänitarkastuksessa selvitettiin käyttäen FFT- taajuusanalyysiä. Erilailla lämpökäsitellyillä teräslajeilla näytti eniten tuloksiin vaikuttavan terästen raekoko. Valssitilaisilla teräksillä raekoko on suuri, jolloin ultraääni vaimenee voimakkaasti teräksessä. Huomattavaa kuitenkin oli, että mikäli lämpö-käsittelyillä ei teräksen raekokoa saada juuri pienennettyä, ei myöskään vaimeneminen vähene. Tämän vuoksi lämpökäsittely ei välttämättä ole aina tarpeellinen valmisteltaessa ultraääninäytteitä. Fourier’n taajuusanalyysissä huomattiin olevan eroavaisuuksia verrattaessa huokosista ja sulkeumista palaavien kaikujen taajuusspektrejä. Näiden tulosten perusteella näyttäisi olevan mahdollista käyttää FFT- menetelmää ultraääni-tarkastuksessa vikojen luokitteluun.

Offset-print mottling measurement using machine vision

Mottling is one of the key defects in offset-printing. Mottling can be defined as unwanted unevenness of print. In this work, diameter of a mottle spot is defined between 0.5-10.0 mm. There are several types of mottling, but the reason behind the problem is still not fully understood. Several commercial machine vision products for the evaluation of print unevenness have been presented. Two of these methods used in these products have been implemented in this thesis. The one is the cluster method and the other is the band-pass method. The properties of human vision system have been taken into account in the implementation of these two methods. An index produced by the cluster method is a weighted sum of the number of found spots, and an index produced by band-pass method is a weighted sum of coefficients of variations of gray-levels for each spatial band. Both methods produce larger indices for visually poor samples, so they can discern good samples from the poor ones. The difference between the indices for good and poor samples is slightly larger produced by the cluster method. 11 However, without the samples evaluated by human experts, the goodness of these results is still questionable. This comparison will be left to the next phase of the project.

Teollisuusprosessien signaalien automaattinen analysointi

Diplomityössä kehitettiin menetelmiä teollisuusprosessien signaalien automaattiseen havainnointiin ja luotiin työkalu tulosten esittämiseen. Työn tarkoituksena on nopeuttaa ja helpottaa prosessin ongelmien ratkaisua luokittelemalla signaalit matemaattisten menetelmien avulla. Koska prosessin mittaussignaalit ovat pääasiassa stokastisia, eli niitä ei voida etukäteen ennustaa, käsitellään signaaleita tilastomatemaattisin keinoin. Työstä rajattiin mittaushistorian käyttö, joten värähtelyiden tunnistus toteutettiin taajuusanalyysin avulla. Korrelaation avulla löydetään samankaltaiset signaalit. Testeissä todettiin, että työssä kehitetyt havainnoinnit toimivat eri näytteenottotaajuuksilla ja työkalun suoritusnopeus todettiin hyväksi. Lopuksi esiteltiin todellinen teollisuusprosessin ongelma ja siihen mahdollisia ratkaisuja.

Yrittäjyyskasvatus yhteistyöverkostoissa

Tutkimus käsittelee Yrittäjyyskasvatuksen Mittariston -projektia, jossa tutkimuskohteena on peruskoulun ensimmäisen asteen luokan- ja aineenopettajien näkemys ja kokemus yrittäjyyskasvatuksen verkostoyhteistyöstä. Tutkimuksen tarkoituksena oli selvittää miten hyvin opettajat tuntevat verkostoyhteistyötä, mikä on heidän tietämyksensä yrittäjyyskasvatuksesta ja kuinka tämä näkyy heidän työssään ja opetuksessaan. Tutkimuksen otos on 450 opettajaa. Tulokset analysoitiin SPSS-tilastomenetelmäohjelmalla. Tilastollisina tutkimusmenetelminä käytettiin jakaumien frekvenssianalyysiä, Faktori­analyysin Pääkomponenttianalyysiä ja Kaksisuuntaista varianssianalyysia (Anova). Tutkimuksen johtopäätöksenä voidaan todeta, että opettajien tiedot yhteistyö-verkostojen tarjoamista palveluista ovat hyvin hajanaiset. Ongelma jatkuu helposti niin kauan kunnes opettajien koulutusohjelmaan tuodaan lisää yrittäjyyskasvatus- ja yrittäjyysopintoja. Tämä pitäisi huomioida myös tulevissa opetussuunnitelmissa. Tämän tutkimuksen tavoitteena oli tuoda esille Yrittäjyyskasvatuksen mittariston tulosten kautta yrittäjyyskasvatuksen nykytila, tuoda ratkaisuja ehdotusten kautta opetukseen ja herättää keskustelua yrittäjyyskasvatuksen parantamiseksi.

La survie implantaire suite à une réhabilitation par mise en charge immédiate d’une prothèse totale mandibulaire reliée à deux implants non jumelés : une étude pilote expérimentale

Problématique : Les connaissances théoriques et pratiques au sujet de la mise en charge immédiate d’implants non jumelés chez les édentés sont limitées. Objectifs : Cette étude avait pour but de : (1) déterminer le taux de survie implantaire de 2 implants non jumelés supportant une prothèse totale mandibulaire suite à une mise en charge immédiate, (2) évaluer les changements des niveaux osseux et de stabilité implantaire survenus sur ces 2 implants durant une période de 4 mois et les comparer à un implant témoin, et (3) décrire les complications cliniques associées à ce mode de mise en charge. Méthodologie : Chez 18 individus édentés (âge moyen de 62±7 ans), cette étude de phase I avec un design pré/post a évalué les résultats cliniques suivant la mise en charge immédiate (<48 heures) de 2 implants non jumelés par une prothèse totale mandibulaire. À l’aide de radiographies périapicales, de sondages osseux et d’analyses de la fréquence en résonnance, les niveaux osseux péri-implantaires (en mm) et les niveaux de stabilité implantairte (en ISQ) de ces 2 implants insérés dans la région parasymphysaire ont été évalués à la chirurgie (T0) et au suivi de 4 mois (T1). Un implant non submergé et sans mise en charge inséré dans la région de la symphyse mandibulaire a été utilisé comme témoin. Les données ont été analysées avec des modèles mixtes linéaires, la méthode de Tukey ajustée, l’analyse de variance de Friedman et des tests de rang signés de Wilcoxon. Résultats : De T0 à T1, 3 implants mis en charge immédiatement ont échoué chez 2 patients. Le taux de survie implantaire obtenu était donc de 91,7% (33/36) et, par patient, de 88,9% (16/18). Aucun implant témoin n’a échoué. Les changements osseux documentés radiologiquement et par sondage autour des implants mis en charge immédiatement étaient, respectivement, de -0,2 ± 0,3 mm et de -0,5 ± 0,6 mm. Les pertes d’os de support implantaire n’ont pas été démontrées statistiquement différentes entre les implants avec mise en charge immédiate et les témoins. Les niveaux moyens de stabilité implantaire ont augmenté de 5 ISQ indépendamment de la mise en charge. Les niveaux moyens d’ISQ n’ont pas été démontrés statistiquement différents entre les implants avec mise en charge immédiate et les témoins à T0 ou T1. Cinq des 18 patients n’ont expérimenté aucune complication clinique, alors que 9 en ont eu au moins deux. Hormis les échecs implantaires, aucune de ces complications n’a entraîné de changements au protocole. Conclusion : Les résultats à court terme suggèrent que : (1) le taux de survie implantaire suite au protocole immédiat est similaire à ceux rapportés lors d’un protocole de mise en charge conventionnel (2) les changements d’os de support implantaire et de stabilité ne sont pas différents comparativement à ceux d’un implant témoin, (3) un niveau élevé d’expérience clinique et chirurgicale est requis pour effectuer les procédures et pour gérer les complications associées. Ces résultats préliminaires devraient être confirmés dans une étude clinique de phase II.

Detection of Parkinson’s disease tremor-onset with a frequency analysing artificial neural network

Deep Brain Stimulator devices are becoming widely used for therapeutic benefits in movement disorders such as Parkinson's disease. Prolonging the battery life span of such devices could dramatically reduce the risks and accumulative costs associated with surgical replacement. This paper demonstrates how an artificial neural network can be trained using pre-processing frequency analysis of deep brain electrode recordings to detect the onset of tremor in Parkinsonian patients. Implementing this solution into an 'intelligent' neurostimulator device will remove the need for continuous stimulation currently used, and open up the possibility of demand-driven stimulation. Such a methodology could potentially decrease the power consumption of a deep brain pulse generator.

Functional timing analysis of VLSI circuits containing complex gates

The recent advances in CMOS technology have allowed for the fabrication of transistors with submicronic dimensions, making possible the integration of tens of millions devices in a single chip that can be used to build very complex electronic systems. Such increase in complexity of designs has originated a need for more efficient verification tools that could incorporate more appropriate physical and computational models. Timing verification targets at determining whether the timing constraints imposed to the design may be satisfied or not. It can be performed by using circuit simulation or by timing analysis. Although simulation tends to furnish the most accurate estimates, it presents the drawback of being stimuli dependent. Hence, in order to ensure that the critical situation is taken into account, one must exercise all possible input patterns. Obviously, this is not possible to accomplish due to the high complexity of current designs. To circumvent this problem, designers must rely on timing analysis. Timing analysis is an input-independent verification approach that models each combinational block of a circuit as a direct acyclic graph, which is used to estimate the critical delay. First timing analysis tools used only the circuit topology information to estimate circuit delay, thus being referred to as topological timing analyzers. However, such method may result in too pessimistic delay estimates, since the longest paths in the graph may not be able to propagate a transition, that is, may be false. Functional timing analysis, in turn, considers not only circuit topology, but also the temporal and functional relations between circuit elements. Functional timing analysis tools may differ by three aspects: the set of sensitization conditions necessary to declare a path as sensitizable (i.e., the so-called path sensitization criterion), the number of paths simultaneously handled and the method used to determine whether sensitization conditions are satisfiable or not. Currently, the two most efficient approaches test the sensitizability of entire sets of paths at a time: one is based on automatic test pattern generation (ATPG) techniques and the other translates the timing analysis problem into a satisfiability (SAT) problem. Although timing analysis has been exhaustively studied in the last fifteen years, some specific topics have not received the required attention yet. One such topic is the applicability of functional timing analysis to circuits containing complex gates. This is the basic concern of this thesis. In addition, and as a necessary step to settle the scenario, a detailed and systematic study on functional timing analysis is also presented.