Ce mémoire traite de l’harmonicitié d’instruments de musique à travers la géométrie spectrale. Nous y présentons, en premier lieu, les résultats connus concernant la corde de guitare, le tambour circulaire et puis le tabla ; le premier est harmonique, le deuxième ne l’est pas et puis le dernier s’en approche. Le cas de la timbale est ce qui constitue la majeure partie de notre travail. L’ingénieur-physicien Robert E. Davis en avait déjà étudié la quasi-harmonicité et nous faisons ici une relecture mathématique de sa démarche. En alliant les méthodes analytiques et numériques, nous montrons que la caisse de résonance de la timbale permet à la fois d’ajuster les fréquences de vibration de la forme ω_(i1) , avec 1 ≤ i ≤ 5, afin qu’elles s’approchent du rapport idéal 2 : 3 : 4 : 5 : 6, et elle permet aussi d’étouffer certains autres modes dissonants. Pour ce faire, nous élaborons un modèle simplifié de timbale cylindrique basé sur la physique et sur ce que propose Davis dans sa thèse. Ce modèle nous fournit un système d’équations divisé en trois parties : la vibration de la peau et la pression à l’intérieur et à l’extérieur de la timbale. Nous utilisons la méthode des fonctions de Green pour trouver les expressions des deux pressions. Nous nous servons de celles-ci ainsi que d’un développement en série de Fourier-Bessel modifiée pour résoudre les équations de la vibration de la peau. La résolution de ces équations se ramène finalement à celle d’un système matriciel infini dont nous faisons l’analyse numériquement. À l’aide de Mathématica et de ce système matriciel, nous trouvons les fréquences de vibration de la timbale, ce qui nous permet d’analyser l’harmonicité de l’instrument. Grâce à une mesure de dissonance, nous optimisons l’harmonicité de la timbale en fonction du rayon du cylindre, de sa hauteur et de la tension. / This thesis deals with the harmonicity of musical instruments through spectral geometry. First, we present the known results concerning the guitar string, the circular drum and the tabla ; the first is harmonic, the second is not, and the last is somewhere in between. The case of the timpani constitutes the major part of our work. The physicist-engineer Robert E. Davis had already studied its quasi-harmonicity and here we undergo a mathematical proofreading of his approach. By combining analytical and numerical methods, we show that the sound box of the timpani allows an adjustement of the vibration frequencies of the form ω_(i1) , with 1 ≤ i ≤ 5, so that they get close to the ideal 2 : 3 : 4 : 5 : 6 ratio, while it also stifles some other dissonant modes. To do so, we develop a simplified model of a cylindrical timpani based on physics and on what Davis suggests in his thesis. This model provides a system of equations divided into three parts : the vibration of the skin and the pressure inside and outside the timpani. We use the method of Green’s functions to find the expressions of the pressures. We use these together with a modified Fourier-Bessel series development to solve the equations of the vibration of the skin. In the end, the solving of these equations is reduced to an infinite matrix system that we analyze numerically. Using Mathematica and this matrix system, we find the vibrational frequencies of the timpani, which allows us to analyze the harmonicity of the instrument. Thanks to a measure of dissonance, we optimize the harmonicity of different timpani models with different cylinder radii, heights and tensions.

<p>Predmet istraživanja u doktorskoj disertaciji je metoda za konstrukciju<br />lokalizacionih skupova za spektar i pseudospektar nelinearnih problema<br />karakterističnih korena bazirana na Geršgorinovoj teoremi i njenim<br />generalizacijama koja koristi osobine poznatih podklasa H-matrica.<br />Navedena tvrđenja i primeri rasvetljavaju odnose između navedenih<br />lokalizacionih skupova, što je posebno značajno za primenu u praksi.<br />Sadržaj ovog rada time predstavlja polaznu tačku za dublja istraživanja na<br />temu konstrukcije lokalizacionih skupova za spektar i pseudospektar<br />nelinearnih problema karakterističnih korena Geršgorinovog tipa.</p> / <p>The subject of research in the doctoral dissertation is a method for constructing<br />spectra and pseudospectra localization sets for nonlinear eigenvalue problems<br />based on Ger&scaron;gorin theorem and its generalizations, that uses the properties of<br />well-known subclasses of H-matrices. Theorems and examples given in this<br />paper are showing relations between stated localization sets, which is very<br />important for practical applications. Therefore, the content of this paper represent<br />the starting point for deeper explorations on the subject of constructing spectra<br />and pseudospectra localization sets for Ger&scaron;gorin type nonlinear eigenvalue<br />problems.</p>

[pt] Neste trabalho, são apresentadas algumas aplicações da otimização topológica para problemas de autovalor onde o principal objetivo é maximizar um determinado autovalor, como por exemplo uma frequência natural de vibração ou uma carga crítica linearizada, usando elementos finitos poligonais em domínios bidimensionais arbitrários. A otimização topológica tem sido comumente utilizada para minimizar a flexibilidade de estruturas sujeitas a restrições de volume. A ideia desta técnica é distribuir uma certa quantidade de material em uma estrutura, sujeita a carregamentos e condições de contorno, visando maximizar a sua rigidez. Neste trabalho, o objetivo é obter uma distribuição ótima de material de maneira a maximizar uma determinada frequência natural (para mantê-la afastada da frequência de excitação externa, por exemplo) ou maximizar a menor carga crítica linearizada (para garantir um nível mais elevado de estabilidade da estrutura). Malhas poligonais construídas usando diagramas de Voronoi são empregadas na solução do problema de otimização topológica. As variáveis de projeto, i.e. as densidades do material, utilizadas no processo de otimização, são associadas a cada elemento poligonal da malha. Vários exemplos de otimização topológica, tanto para problemas de frequências naturais de vibração quanto para cargas críticas linearizadas, são apresentados para demonstrar a funcionalidade e a aplicabilidade da metodologia proposta. / [en] In this work, we present some applications of topology optimization for eigenvalue problems where the main goal is to maximize a specified eigenvalue, such as a natural frequency or a linearized buckling load using polygonal finite elements in arbitrary two-dimensional domains. Topology optimization has commonly been used to minimize the compliance of structures subjected to volume constraints. The idea is to distribute a certain amount of material in a given design domain subjected to a set of loads and boundary conditions such that to maximize its stiffness. In this work, the objective is to obtain the optimal material distribution in order to maximize the fundamental natural frequency (e.g. to keep it away from an external excitation frequency) or to maximize the lowest critical buckling load (e.g. to ensure a higher level of stability of the structures). We employ unstructured polygonal meshes constructed using Voronoi tessellations for the solution of the structural topology optimization problems. The design variables, i.e. material densities, used in the optimization scheme, are associated with each polygonal element in the mesh. We present several topology optimization examples for both eigenfrequency and buckling problems in order to demonstrate the functionality and applicability of the proposed methodology.

This thesis presents a means of enhancing the iterative calculation techniques used in electronic structure calculations, particularly Kohn-Sham DFT. Based on the subspace iteration method of the FEAST eigenvalue solving algorithm, this nonlinear FEAST algorithm (NLFEAST) improves the convergence rate of traditional iterative methods and dramatically improves their robustness. A description of the algorithm is given, along with the results of numerical experiments that demonstrate its effectiveness and offer insight into the factors that determine how well it performs.

The famous Moore's law states: Since the invention of the integrated circuit, the number of transistors that can be placed on an integrated circuit has increased exponentially, doubling approximately every two years. As a result of the downscaling of the size of the transistor, quantum effects have become increasingly important while affecting significantly the device performances. Nowadays, at the nanometer scale, inter-atomic interactions and quantum mechanical properties need to be studied extensively. Device and material simulations are important to achieve these goals because they are flexible and less expensive than experiments. They are also important for designing and characterizing new generation of electronic device such as silicon nanowire or carbon nanotube (CNT) transistors. Several modeling methods have been developed and applied to electronic structure calculations, such as: Hartree-Fock, density functional theory (DFT), empirical tight-binding, etc. For transport simulations, most of the device community focuses on studying the stationary problem for obtaining characteristics such as I-V curves. The non-equilibrium transport problem is then often addressed by solving a multitude of time-independent Schrodinger-type equation for all possible energies. On the other hand, for many other electronic applications including high-frequency electronics response (e.g. when a time-dependent potential is applied to the system), the description of the system behavior necessitate insights on the time dependent electron dynamics. To address this problem, it is then necessary to solve a time-dependent Schrodinger-type equation. In this thesis, we will focus on solving time-dependent problems with application to CNTs. We will be identifying all the numerical difficulties and propose new effective modeling and numerical schemes to address the current limitations in time-dependent quantum simulations. we will point out that two numerical errors may occur: an integration error and the anti-commutation issue error; the direct computation above being mathematically equivalent to performing the integration of the time dependent Hamiltonian using a rectangle numerical quadrature formula along the total simulation times. After careful study and many numerical experiments, we found that the Gaussian quadrature scheme provides a good trade off between computational consumption and numerically accuracy, meanwhile unitary, stability and time reversal properties are well preserved. The new Gaussian quadrature integration scheme uses (i) much fewer points in time to approximate the integral of the Hamiltonian, (ii) ordered exponential to factorize the time evolution operator, (iii) FEM discretize techniques (iv) and at last, the FEAST eigenvalue solver to diagonalize and solve each exponential.

Max-plus代数は，実数全体に無限小元を付加した集合に，加法として最大値をとる演算，乗法として通常の加法を考えた代数系である．本論文では，max-plus線形方程式に対するCramerの公式の類似物を用いて，線形方程式の解空間の基底が構成できることを示した．さらに固有値問題に関連して，max-plus行列の固有ベクトルの概念を2通りの観点から拡張した． / The max-plus algebra is the semiring with addition "max" and multiplication "+". In the present thesis, the author gives a combinatorial characterization of solutions of linear systems in terms of the max-plus Cramer's rule. Further, the author extends the concept of eigenvectors of max-plus matrices from two different perspectives. / 博士(理学) / Doctor of Philosophy in Science / 同志社大学 / Doshisha University

Optical testing in adverse environments, ophthalmology and applications where characterization by curvature is leveraged all have a common goal: accurately estimate wavefront shape. This dissertation investigates wavefront sensing techniques as applied to optical testing based on gradient and curvature measurements. Wavefront sensing involves the ability to accurately estimate shape over any aperture geometry, which requires establishing a sampling grid and estimation scheme, quantifying estimation errors caused by measurement noise propagation, and designing an instrument with sufficient accuracy and sensitivity for the application. Starting with gradient-based wavefront sensing, a zonal least-squares wavefront estimation algorithm for any irregular pupil shape and size is presented, for which the normal matrix equation sets share a pre-defined matrix. A Gerchberg–Saxton iterative method is employed to reduce the deviation errors in the estimated wavefront caused by the pre-defined matrix across discontinuous boundary. The results show that the RMS deviation error of the estimated wavefront from the original wavefront can be less than λ/130~ λ/150 (for λ equals 632.8nm) after about twelve iterations and less than λ/100 after as few as four iterations. The presented approach to handling irregular pupil shapes applies equally well to wavefront estimation from curvature data. A defining characteristic for a wavefront estimation algorithm is its error propagation behavior. The error propagation coefficient can be formulated as a function of the eigenvalues of the wavefront estimation-related matrices, and such functions are established for each of the basic estimation geometries (i.e. Fried, Hudgin and Southwell) with a serial numbering scheme, where a square sampling grid array is sequentially indexed row by row. The results show that with the wavefront piston-value fixed, the odd-number grid sizes yield lower error propagation than the even-number grid sizes for all geometries. The Fried geometry either allows sub-sized wavefront estimations within the testing domain or yields a two-rank deficient estimation matrix over the full aperture; but the latter usually suffers from high error propagation and the waffle mode problem. Hudgin geometry offers an error propagator between those of the Southwell and the Fried geometries. For both wavefront gradient-based and wavefront difference-based estimations, the Southwell geometry is shown to offer the lowest error propagation with the minimum-norm least-squares solution. Noll’s theoretical result, which was extensively used as a reference in the previous literature for error propagation estimate, corresponds to the Southwell geometry with an odd-number grid size. For curvature-based wavefront sensing, a concept for a differential Shack-Hartmann (DSH) curvature sensor is proposed. This curvature sensor is derived from the basic Shack-Hartmann sensor with the collimated beam split into three output channels, along each of which a lenslet array is located. Three Hartmann grid arrays are generated by three lenslet arrays. Two of the lenslets shear in two perpendicular directions relative to the third one. By quantitatively comparing the Shack-Hartmann grid coordinates of the three channels, the differentials of the wavefront slope at each Shack-Hartmann grid point can be obtained, so the Laplacian curvatures and twist terms will be available. The acquisition of the twist terms using a Hartmann-based sensor allows us to uniquely determine the principal curvatures and directions more accurately than prior methods. Measurement of local curvatures as opposed to slopes is unique because curvature is intrinsic to the wavefront under test, and it is an absolute as opposed to a relative measurement. A zonal least-squares-based wavefront estimation algorithm was developed to estimate the wavefront shape from the Laplacian curvature data, and validated. An implementation of the DSH curvature sensor is proposed and an experimental system for this implementation was initiated. The DSH curvature sensor shares the important features of both the Shack-Hartmann slope sensor and Roddier’s curvature sensor. It is a two-dimensional parallel curvature sensor. Because it is a curvature sensor, it provides absolute measurements which are thus insensitive to vibrations, tip/tilts, and whole body movements. Because it is a two-dimensional sensor, it does not suffer from other sources of errors, such as scanning noise. Combined with sufficient sampling and a zonal wavefront estimation algorithm, both low and mid frequencies of the wavefront may be recovered. Notice that the DSH curvature sensor operates at the pupil of the system under test, therefore the difficulty associated with operation close to the caustic zone is avoided. Finally, the DSH-curvature-sensor-based wavefront estimation does not suffer from the 2π-ambiguity problem, so potentially both small and large aberrations may be measured.

No description available.

One of the problems with the modern radio communication is the lack of availableradio frequencies. Recent studies have shown that, while the available licensed radiospectrum becomes more occupied, the assigned spectrum is significantly underutilized.To alleviate the situation, cognitive radio (CR) technology has been proposedto provide an opportunistic access to the licensed spectrum areas. Secondary CRsystems need to cyclically detect the presence of a primary user by continuouslysensing the spectrum area of interest. Radiowave propagation effects like fading andshadowing often complicate sensing of spectrum holes. When spectrum sensing isperformed in a cooperative manner, then the resulting sensing performance can beimproved and stabilized. In this thesis, two fully distributed and adaptive cooperative Primary User (PU)detection solutions for CR networks are studied. In the first part of this thesis we study a distributed energy detection schemewithout using any fusion center. Due to reduced communication such a topologyis more energy efficient. We propose the usage of distributed, diffusion least meansquare (LMS) type of power estimation algorithms with different network topologies.We analyze the resulting energy detection performance by using a commonframework and verify the theoretical findings through simulations. In the second part of this thesis we propose a fully distributed detection scheme,based on the largest eigenvalue of adaptively estimated correlation matrices, assumingthat the primary user signal is temporally correlated. Different forms of diffusionLMS algorithms are used for estimating and averaging the correlation matrices overthe CR network. The resulting detection performance is analyzed using a commonframework. In order to obtain analytic results on the detection performance, theadaptive correlation matrix estimates are approximated by a Wishart distribution.The theoretical findings are verified through simulations. / <p>QC 20170908</p>

In this thesis, a method to calculate eigenpairs is implemented for the Multipsi library. While the standard implemtentations use the Davidson method with Rayleigh-Ritz extraction to calculate the eigenpairs with the lowest eigenvalues, the new method uses the harmonic Davidson method with the harmonic Rayleigh-Ritz extraction to calculate eigenpairs with eigenvalues near a chosen target. This is done for Configuration Interaction calculations and for Multiconfigurational methods. From calculations, it seems the new addition to the Multipsi library is worth investigating further as convergence for difficult systems with a lot of near-degeneracy was improved. / I denna avhandling implementeras en metod för att beräkna egenpar för Multipsi-biblioteket. Medan standardimplementeringarna använder Davidson-metoden med Rayleigh-Ritz-extraktion för att beräkna egenparen med de lägsta egenvärdena, använder den nya metoden den harmoniska Davidson-metoden med den harmoniska Rayleigh-Ritz-extraktionen för att beräkna egenparen med egenvärden nära ett valt mål. Detta görs för konfigurationsinteraktionsberäkningar och för multikonfigurationsmetoder. Utifrån beräkningarna verkar det nya tillskottet till Multipsi-biblioteket vara värt att undersöka vidare eftersom konvergensen för svåra system med mycket nära degenerering förbättrades.