Preferencje help
Widoczny [Schowaj] Abstrakt
Liczba wyników

Znaleziono wyników: 10

Liczba wyników na stronie
first rewind previous Strona / 1 next fast forward last
Wyniki wyszukiwania
Wyszukiwano:
w słowach kluczowych:  numerical stability
help Sortuj według:

help Ogranicz wyniki do:
first rewind previous Strona / 1 next fast forward last
EN
Further tests of EVT4 data acquisition system for electrical capacitance tomography are presented. The modular system, which can have up to 32 channels with an individual analogue to digital converter, was designed to ensure small uncertainty of capacitance measurement at high speed of imaging. The system’s performance in the context of 3D imaging was experimentally verified. In particular, we show that the measurement of changes in capacitance due to a small change of an electric permittivity distribution for the most distant electrodes in a suitably designed 3D sensor is possible using our system. Cross-plane measurements together with the measurements for the pairs of most distant electrodes are essential for accurate reconstruction of 3D distributions. Due to sensitivity of capacitance measurements obtained in the hardware, the measurements for all electrode pairs can be used in the inverse problem – the system of equations can be extended. Although the numerical condition number of a matrix of such a system is high, image reconstruction is possible from the data obtained in our system. The results of 3D image reconstruction for simple test objects are shown.
PL
Przedstawiono dalsze testy systemu akwizycji danych EVT4 do elektrycznej tomografii pojemnościowej. Modułowy system, który może mieć do 32 kanałów z indywidualnym przetwornikiem analogowo-cyfrowym, został zaprojektowany w celu zapewnienia małej niepewności pomiaru pojemności przy dużej prędkości obrazowania. Wydajność systemu w kontekście obrazowania 3D została zweryfikowana eksperymentalnie. W szczególności pokazujemy, że możliwy jest pomiar zmian pojemności wywoływanych niewielką zmianą rozkładu przenikalności elektrycznej dla najbardziej odległych elektrod odpowiednio zaprojektowanego czujnika 3D przy pomocy naszego systemu. Pomiary międzypłaszczyznowe wraz z pomiarami par elektrod najbardziej odległych są niezbędne do dokładnej rekonstrukcji rozkładów 3D. Ze względu na wrażliwość pomiarów pojemności uzyskanych w opracowanym urządzeniu, pomiary dla wszystkich par elektrod mogą być wykorzystane w problemie odwrotnym – układ równań może zostać rozszerzony. Chociaż współczynnik uwarunkowania numerycznego macierzy takiego układu jest wysoki, możliwa jest rekonstrukcja obrazu z wykorzystaniem danych uzyskiwanych w naszym systemie. Pokazane są wyniki rekonstrukcji obrazu 3D dla prostych obiektów testowych.
EN
This paper presents a programmable system-on-chip implementation to be used for acceleration of computations within hidden Markov models. The high level synthesis (HLS) and “divide-and-conquer” approaches are presented for parallelization of Baum-Welch and Viterbi algorithms. To avoid arithmetic underflows, all computations are performed within the logarithmic space. Additionally, in order to carry out computations efficiently – i.e. directly in an FPGA system or a processor cache – we postulate to reduce the floating-point representations of HMMs. We state and prove a lemma about the length of numerically unsafe sequences for such reduced precision models. Finally, special attention is devoted to the design of a multiple logarithm and exponent approximation unit (MLEAU). Using associative mapping, this unit allows for simultaneous conversions of multiple values and thereby compensates for computational efforts of logarithmic-space operations. Design evaluation reveals absolute stall delay occurring by multiple hardware conversions to logarithms and to exponents, and furthermore the experiments evaluation reveals HMMs computation boundaries related to their probabilities and floating-point representation. The performance differences at each stage of computation are summarized in performance comparison between hardware acceleration using MLEAU and typical software implementation on an ARM or Intel processor.
PL
Jedną z uniwersalnych, a zarazem prostych metod numerycznego rozwiązywania zadań nieustalonego przewodzenia ciepła w obudowie budynku jest metoda objętości skończonych. Jednak przy sformułowaniu jawnym konieczne jest dobranie kroku czasowego gwarantującego stabilność rozwiązania. W tym artykule wartość takiego kroku czasowego została określona z warunku zerowej różnicy pomiędzy teoretycznym a numerycznym rozwiązaniem dla danej geometrii ściany. Rozpatrzono wpływ współczynnika przewodzenia ciepła ściany na długość kroku. Warunki brzegowe przyjęte w artykule odpowiadały obciążeniom ogniowym ściany schronu.
EN
One of the most simplest and versatile methods for numeric solving of the problems of transient heat conduction in the building envelope is the finite volume method. In explicit scheme, however, there is a problem of assuming correct time step for solution stability. In this article, stable step value was determined from the condition of zero difference between the theoretical and numerical solution for the particular geometry of the wall. The influence of the thermal conductivity on the step length has been evaluated. Condition for the correctness of solutions obtained from the analysis in this article is sharper than the condition of equations stability according to the literature. For the boundary conditions, the adopted solution can be applied to the study of resistance of buildings shelters to long fire.
EN
In the paper the influence of the eigenvalues of the matrix on the stability of the numerical simulations of binary alloy solidification process was analyzed. One step Θ scheme time integration for the modeling was used. The occurrence of the equation clotting full, diagonal and lumped capacitance matrix was considered. In the explicit schema (conditionally stable) the major influence on the proper selection of the size of the critical time step was proved. The physically incorrect results are given when the maximum allowable time step is exceeded.
EN
A simple semi-recursive routine for nonlinearity recovery in Hammerstein systems is proposed. The identification scheme is based on the Haar wavelet kernel and possesses a simple and compact form. The convergence of the algorithm is established and the asymptotic rate of convergence (independent of the input density smoothness) is shown for piecewise-Lipschitz nonlinearities. The numerical stability of the algorithm is verified. Simulation experiments for a small and moderate number of input-output data are presented and discussed to illustrate the applicability of the routine.
EN
Supernodal technique, used in the field of sparse linear equation solvers, introduces dense kernel for speeding-up the computations. SuperLU solver is the most sophisticated and efficient example of this technique. The research reported here shows that supernodal solver is really fast, but this spectacular speed is obtained through abandoning of numerical stability. It has been shown that SuperLU solver fails in whole class of practical real-world problems. An attempt to characterize and diagnose this problems in terms of matrix parameters has been madę. Pointer solver, being stable numerically and fairly fast alternative for supernodal solver is presented quite comprehensively.
PL
Metoda superwęzłów stosowana przy rozwiązywaniu rzadkich układów równań algebraicznych polega na specjalnym traktowaniu przydiagonalnych fragmentów macierzy rzadkich. Fragmenty te traktowane są jako submacierze pełne, przez co uzyskuje się przyspieszenie obliczeń kosztem, jak podają niektórzy badacze, utraty stabilności numerycznej. Praca jest próbą analizy tego problemu w kategoriach parametrów macierzy rozwiązywanego układu równań. Przedstawia stabilną numerycznie i dość szybką alternatywę - solver ,,wskaźnikowy" (pointer solver).
7
EN
Laplace transforms in two variables can very be useful for solving certain partial differential equations, namely those describing transient behaviour of linear dynamical systems. In practice, it is often either too difficult or even impossible to obtain corresponding originals by analytic methods. In such cases methods that enable getting original numerically have to be applied. The 2D-NIL T method based on FFT, recently published and verified in Matlab language, seems be well usable. Its main advantage lies in high speed of calculation, however, it is necessary to connect it always with proper technigue of acceleration of the convergence to achieve required accuracy. It was shown either the epsilon or the quotient-difference algorithms are very suitable for this purpose. In the paper an error analysis, comparison and estimation of the optimal 2D-NILT parameters are newly presented.
PL
Dwuwymiarowa transformata Laplace'a (2D-NILT) jest niezwykle użyteczna przy rozwiązywaniu pewnych równań różniczkowych, w szczególności opisujących stany nieustalone systemów liniowych. W wielu praktycznych wypadkach otrzymanie odpowiadających rozwiązań metodami analitycznymi jest zbyt trudne lub nawet niemożliwe. W takich wypadkach niezbędna jest metoda numeryczna. Metoda 2D-NILT wykorzystująca FFT, ostatnio opublikowana i zaimplementowana w Matlabie wydaje się tu być bardzo użyteczna. Jej główną zaletą jest szybkość, uzyskiwana pod warunkiem zastosowania odpowiedniej metody przyspieszającej zbieżność. Jak pokazano, doskonale nadają się do tego: algorytm epsilon lub algorytm różnicy współczynnika. W pracy przedstawiono analizę dokładności i estymację optymalnych parametrów metody 2D-NILT.
EN
Propagation of a Brillouin precursor in a Lorentz dispersive medium is considered. The precursor is excited by a sine-modulated initial signal, with its envelope described by a hyperbolic tangent function. The purpose of the paper is to show how the rate of growth of the signal affects the form of the Brillouin precursor. Uniform asymptotic approach, pertinent to coalescing saddle points, is applied in the analysis. The results are illustrated with numerical examples.
EN
A one-dimensional electromagnetic problem of Sommerfeld precursor evolution, resulting from a finite rise-time signal excitation in a dispersive Lorentz medium is considered. The effect of the initial signal rate of growth as well as of the medium dumping on the precursor shape and its magnitude is discussed. The analysis applied is based on an approach employing uniform asymptotic expansions. In addition, new approximate formulas are given for the location of the distant saddle points which affect local frequency and dumping of the precursor. The results obtained are illustrated numerically and compared with the results known from the literature.
10
Content available remote Measuring and Maintaining Consistency: a Hybrid Ftf Algorithm
EN
Due to the versatility as well as its ease of implementation, the Fast Transversal Filters algorithm is attractive for many adaptive filtering applications. However, it is not widely used because of its undesirable tendency to diverge when operating in finite precision arithmetic. To compensate, modifications to the algorithm have been introduced that are either occasional (performed when a predefined condition(s) is violated) or structured as part of the normal update iteration. However, in neither case is any confidence explicitly given that the computed parameters are in fact close to the desired ones. Here, we introduce a time invariant parameter that provides the user with more flexibility in establishing confidence in the consistency of the updated filter parameters. Additionally, we provide evidence through the introduction of a hybrid FTF algorithm that when sufficient time is given prior to catastrophic divergence, the update parameters of the FTF algorithm can be adjusted so that consistency can be acquired and maintained.
first rewind previous Strona / 1 next fast forward last
JavaScript jest wyłączony w Twojej przeglądarce internetowej. Włącz go, a następnie odśwież stronę, aby móc w pełni z niej korzystać.