Domain Coloring of Complex Functions

Total Page:16

File Type:pdf, Size:1020Kb

Domain Coloring of Complex Functions Domain Coloring of Complex Functions Konstantin Poelke and Konrad Polthier 1 Introduction 2 What is a Function? Let us briefly recap the definition of a Visualizing functions is an omnipresent function to fix terminology. A function f task in many sciences and almost every day consists of three parts: first, a set D of in- we are confronted with diagrams in news- put values, which is called the domain of the papers and magazines showing functions of function, second, a set Y called the range of all possible flavours. Usually such func- f and third, for every input value x ∈ D, tions are visualized by plotting their func- a unique value y ∈ Y , called the function tion graph inside an appropriate coordinate value of f at x, denoted f(x). The set Γ(f) system, with the probably most prominent of all pairs (a, f(a)), a ∈ D, is a subset of choice being the cartesian coordinate sys- the product set D × Y and called the func- tem. This allows us to get an overall im- tion graph of f. pression of the function’s behaviour as well One particular type of functions that are as to detect certain distinctive features such widely used in engineering and physics are as minimal or maximal points or points complex functions, i.e. functions f : D ⊆ where the direction of curvature changes. C → Y ⊆ C whose domain and range In particular, we can “see” the dependence are subsets of the complex numbers, and between input and output. However, this we will focus on complex functions in the technique is limited to three dimensions, following. We will identify C with the Eu- 2 simply because we do not know how to clidean plane R by the bijective assignment draw higher-dimensional cartesian coordi- x + yi 7→ (x,y), justifying the term complex nate systems. plane. Note that the function graph Γ(f) of a complex function lives inside the product space C × C ∼= R4. This article gives a short overview about the method of domain coloring for com- 3 Domain Coloring plex functions, which have four-dimensional In contrast to function graph plotting in- function graphs and therefore cannot be vi- side Euclidean two or three space, the tech- sualized traditionally. We discuss several nique of domain coloring does not need ad- color schemes focussing on various aspects ditional spatial dimensions for the range of complex functions, and provide Java-like but rather uses color dimensions to encode pseudocode examples explaining the cru- the function values. More precisely, given a cial ideas of the coloring algorithms to al- complex function f : D → Y , one defines low for easy reproduction. For a thorough a color scheme as a function col : Y → treatment of domain colorings from a more HSB, which assigns to every value y ∈ Y mathematical point of view see [1]. a color col(y), specified in HSB color space 1 for instance. In practice D is usually a rect- reducing the problem to the flat case in R2. 2 angular region of R . Then an initially cho- 4 Color Schemes sen resolution divides D into a discretized 4.1 Color Wheels domain Dh of pixels corresponding to small Every non-zero point z =(u,v) in the Eu- rectangles of fixed width and height. Ev- clidean plane has a unique representation in ery pixel i is identified with a point zi in D terms of polar coordinates as z = r·exp(iϕ), where the function is evaluated, for example if one requires the angle ϕ to lie in the in- the midpoint or a corner. Then one evalu- tervall [−π,π). ϕ is called the argument ates the pull-back function f ∗col := col◦f of z, denoted arg z, and r is the modu- at every point zi ∈ D and assigns the re- ∗ lus or absolute value |z|. The advantage sulting color value f col(zi) to the pixel i, of this representation is that the argument leading to a coloring of the whole domain can be directly identified with a hue value, Dh. This procedure is shown in Figure 1, which is usually given as an absolute an- Listing 1 gives an easy implementation. gle value between 0 and 360 or as a rela- z1 f col tive value between 0 and 1. We will make steady use of the function HSBColor(float hue, z3 float sat, float bri) which constructs an inte- z2 ger color representation from the specified Dh Y hue, saturation and brightness values, each Figure 1: Every pixel i in the discretized domain of them lying in the interval [0, 1]. Such Dh corresponds to a value zi ∈ D and is colored functions are usually contained in the stan- ∗col with the color f (zi). dard libraries of modern programming lan- Listing 1: Plotting procedure. The domain D is a guages, e.g. java.awt.Color.HSBtoRGB in Java or rectangular region specified by its lower left and up- colorsys.hsv_to_rgb in Python. per right corner, whose (x,y)-coordinates are given This leads to a first color scheme as shown in the array double[] dom of size 4. The resolution in Figure 2, using the argument arg z as its parameters determine the number of pixels of the only parameter. Listing 2 gives its imple- resulting image. mentation. Image plot(double[] dom, Fun f, Fun col, int xRes, int yRes){ // create new image object of size xRes*yRes Image im = new Image(xRes,yRes); // Increments in width and height direction double xInc = (dom[2]-dom[0])/xRes; double yInc = (dom[3]-dom[1])/yRes; for(y=0; y<yRes; y++){ for(x=0; x<xRes; x++){ // midpoint where f is evaluated z = (dom[0]+(x+0.5)*xInc, dom[1]+(y+0.5)*yInc); // assign computed color to pixel im.setColor(x,y,col(f(z))); }} return im; Figure 2: Hue-based color wheel and resulting do- } main coloring for z 7→ z2. The function duplicates We just mention for completeness that the argument of every point in the plane. As a re- this idea can easily be transferred to sur- sult the upper half plane is mapped to the whole faces D embedded in R3 by using local co- plane, as every color already appears on the upper half. ordinate charts and texture mappings and 2 Listing 2: A hue based color wheel using the argu- ment arg z as its only input parameter. int col(Point z){ //arg only defined for non-zero value if (z==0) return black; // h between 0 and 1 h = (arg(z)+π)/2π; return HSBColor(h, 1,1); } In case the argument is undefined, i.e. if z = 0, the function returns black. We will Figure 3: The discrete color wheel divides the plane neglect this form of easy exception handling into its quadrants. The image on the right shows a plot of the Joukowski function (z+1/z)/2. Inside a as well as the conversion of arg z into the circle of radius 1 the plane appears mirrored along relative parameter h ∈ [0, 1] in the follow- the horizontal real axis. ing. Instead of using the hue value as a contin- scheme to functions of a type particularly uous parameter, a discretized version of the important in algebra and complex geome- color wheel divides the plane into n equally- try. sized cells, where in each cell the color is constant. A natural choice is n = 4, divid- ing the Euclidean plane into its four quad- rants. This makes it easy to detect which points are mapped into which quadrant, see Figure 3 and Listing 3. Listing 3: A discretized color wheel implementation dividing the plane into four quadrants. int col(Point z){ // array holding colors for cells Figure 4: A simple gradient color scheme. The Color[] colors = {blue,red,yellow,green}; left picture shows the function f(z) = z5 − 1, // compute cell index whose zeros are exactly the fifth roots of unity. int i = floor(h*4) % 4; return colors[i]; The five endpoints are the zeros whereas the other } endpoints all meet at infinity, where f has a pole of order five. The horse-shoe-shaped image Alternatively one can use the argument on the right displays the meromorphic function 2 2 as a parameter for interpolating between fm(z)=(z − 1)(z + 1) /((z + i)(z − i) ). The two two or more colors leading to color gradi- endpoints are of order one and correspond to the ents wrapped around the origin. A sim- terms (z − 1) and (z + i). The two points where ple gradient color scheme between white black and white interchange are of order two, de- and black, discontinuous along the nega- termined by the remaining terms. tive real axis, is then given by a function Although the preceding color schemes blend(black,white, h) which accepts two colors only used the argument of a point and uses the relative hue value h as a weight z = r exp(iϕ), they already gave fairly good for a linear interpolation between the col- representations for complex functions. By ors’ RGB values such that blend(black,white,0) taking the modulus into account we can en- returns black and blend(black,white,1) returns rich these color schemes with several addi- white. Figure 4 shows an application of this tional features. As an example, Listing 4 3 marks zeros and poles as black and white this phenomenom, it is helpful to imple- spots to make it easy to detect these distin- ment grid-like color schemes such as rect- guished points of a complex function, im- angular grids or polar grids.
Recommended publications
  • GPU-Based Visualization of Domain-Coloured Algebraic Riemann Surfaces
    GPU-based visualization of domain-coloured algebraic Riemann surfaces Stefan Kranich∗ 11 November 2015 Abstract We examine an algorithm for the visualization of domain-coloured Riemann surfaces of plane algebraic curves. The approach faithfully reproduces the topology and the holomorphic structure of the Riemann surface. We discuss how the algorithm can be implemented efficiently in OpenGL with geometry shaders, and (less efficiently) even in WebGL with multiple render targets and floating point textures. While the generation of the surface takes noticeable time in both implementations, the visualization of a cached Riemann surface mesh is possible with interactive performance. This allows us to visually explore otherwise almost unimaginable mathematical objects. As examples, we look at the complex square root and the folium of Descartes. For the folium of Descartes, the visualization reveals features of the algebraic curve that are not obvious from its equation. 1 Introduction 1.1 Mathematical background The following basic example illustrates what we would like to visualize. Example 1.1. Let y be the square root of x, p y = x: If x is a non-negative real number, we typically define y as the non-negative real arXiv:1507.04571v3 [cs.GR] 10 Nov 2015 number whose square equals x, i.e. we always choose the non-negative solution of the equation y2 − x = 0 (1) p as y = x. For negative real numbers x, no real number y solves Equation 1. However, if we define the imaginary unit i as a number with the property that i2 = −1 then the square root of x becomes the purely imaginary number y = ipjxj: ∗Zentrum Mathematik (M10), Technische Universit¨atM¨unchen, 85747 Garching, Germany; E-mail address: [email protected] 1 Im Im Re Re Figure 1.2: When a complex number (black points) runs along a circle centred at the origin of the complex plane, its square roots (white points) move at half the angular velocity (left image).
    [Show full text]
  • Visualization of Complex Function Graphs in Augmented Reality
    M A G I S T E R A R B E I T Visualization of Complex Function Graphs in Augmented Reality ausgeführt am Institut für Softwaretechnik und Interaktive Systeme der Technischen Universität Wien unter der Anleitung von Univ.Ass. Mag. Dr. Hannes Kaufmann durch Robert Liebo Brahmsplatz 7/11 1040 Wien _________ ____________________________ Datum Unterschrift Abstract Understanding the properties of a function over complex numbers can be much more difficult than with a function over real numbers. This work provides one approach in the area of visualization and augmented reality to gain insight into these properties. The applied visualization techniques use the full palette of a 3D scene graph's basic elements, the complex function can be seen and understood through the location, the shape, the color and even the animation of a resulting visual object. The proper usage of these visual mappings provides an intuitive graphical representation of the function graph and reveals the important features of a specific function. Augmented reality (AR) combines the real world with virtual objects generated by a computer. Using multi user AR for mathematical visualization enables sophisticated educational solutions for studies dealing with complex functions. A software framework that has been implemented will be explained in detail, it is tailored to provide an optimal solution for complex function graph visualization, but shows as well an approach to visualize general data sets with more than 3 dimensions. The framework can be used in a variety of environments, a desktop setup and an immersive setup will be shown as examples. Finally some common tasks involving complex functions will be shown in connection with this framework as example usage possibilities.
    [Show full text]
  • Domain Coloring of Complex Functions
    Domain Coloring of Complex Functions Contents 1 Introduction Not only in mathematic we are confronted with graphs of all possible data sets. In economics and science we want to present the results by plotting their function graph inside a optimal coordinate system. In complex analysis we work with holomorphic functions which are complex differen- tiable. Like in real analysis we are interesting to plot such functions as graph inside a optimal coordinate systems. The graph of holomorphic functions is a subset of the four-dimensional coordinate system. But we are limited to three dimension, because we do not know how to draw objects in higher dimensional spaces. Therefore we are forced to construct a methode to visualize such functions. In this sec- tion we give a short overview of the domain coloring which allows us to visualize the graph of complex funtions by using colors. First we give a overview of complex analysis and holomorphic function. After we have presented the main idea behind the domain coloring we discuss a example to understand the main concept of domain coloring. 1 2 Introduction to Complex Analysis 2 Introduction to Complex Analysis In this section we give a short review of the main idea of complex analysis which we need to understand the methode of domain coloring. The reader should familiar with the construction of complex numbers and the representation of such numbers with polar coordinates. In complex anlysis we work with complex function which consist of three parts: first, a set D ⊂ C of input values, which is called the domain of the function, second, the range of f in C and third, for every input value x 2 D, a unique function value f(x) in the range of f.
    [Show full text]
  • Una Introducción Al Método De Dominio Colorado Con Geogebra Para
    101 Una introducción al método de dominio colorado con GeoGebra para la visualización y estudio de funciones complejas Uma introdução ao método do domínio colorado com GeoGebra para visualizar e estudar funções complexas An introduction the method domain coloring with GeoGebra for visualizing and studying complex JUAN CARLOS PONCE CAMPUZANO1 0000-0003-4402-1332 researchgate.net/profile/Juan_Ponce_Campuzano geogebra.org/u/jcponce http://dx.doi.org/10.23925/2237-9657.2020.v9i1p101-119 RESUMEN Existen diversos métodos para visualizar funciones complejas, tales como graficar por separado sus componentes reales e imaginarios, mapear o transformar una región, el método de superficies analíticas y el método de dominio coloreado. Este último es uno de los métodos más recientes y aprovecha ciertas características del color y su procesamiento digital. La idea básica es usar colores, luminosidad y sombras como dimensiones adicionales, y para visualizar números complejos se usa una función real que asocia a cada número complejo un color determinado. El plano complejo puede entonces visualizarse como una paleta de colores construida a partir del esquema HSV (del inglés Hue, Saturation, Value – Matiz, Saturación, Valor). Como resultado, el método de dominio coloreado permite visualizar ceros y polos de funciones, ramas de funciones multivaluadas, el comportamiento de singularidades aisladas, entre otras propiedades. Debido a las características de GeoGebra en cuanto a los colores dinámicos, es posible implementar en el software el método de dominio coloreado para visualizar y estudiar funciones complejas, lo cual se explica en detalle en el presente artículo. Palabras claves: funciones complejas, método de dominio coloreado, colores dinámicos. RESUMO Existem vários métodos para visualizar funções complexas, como plotar seus componentes reais e imaginários separadamente, mapear ou transformar uma região, o método de superfície analítica e o método de domínio colorido.
    [Show full text]
  • Arxiv:2002.05234V1 [Cs.GR] 12 Feb 2020 Have Emerged
    VISUALIZING MODULAR FORMS DAVID LOWRY-DUDA Abstract. We describe practical methods to visualize modular forms. We survey several current visualizations. We then introduce an approach that can take advantage of colormaps in python's matplotlib library and describe an implementation. 1. Introduction 1.1. Motivation. Graphs of real-valued functions are ubiquitous and com- monly serve as a source of mathematical insight. But graphs of complex functions are simultaneously less common and more challenging to make. The reason is that the graph of a complex function is naturally a surface in four dimensions, and there are not many intuitive embeddings available to capture this surface within a 2d plot. In this article, we examine different methods for visualizing plots of mod- ular forms on congruent subgroups of SL(2; Z). These forms are highly symmetric functions and we should expect their plots to capture many dis- tinctive, highly symmetric features. In addition, we wish to take advantage of the broader capabilities that exist in the python/SageMath data visualization ecosystem. There are a vast number of color choices and colormaps implemented in terms of python's matplotlib library [Hun07]. While many of these are purely aesthetic, some offer color choices friendly to color blind viewers. Further, some are designed with knowledge of color theory and human cognition to be perceptually uniform with respect to human vision. We describe this further in x3. 1.2. Broad Overview of Complex Function Plotting. Over the last 20 years, different approaches towards representing graphs of complex functions arXiv:2002.05234v1 [cs.GR] 12 Feb 2020 have emerged.
    [Show full text]
  • Real-Time Visualization of Geometric Singularities Master’S Thesis by Stefan Kranich
    TECHNISCHE UNIVERSITÄT MÜNCHEN Department of Mathematics Real-time Visualization of Geometric Singularities Master’s thesis by Stefan Kranich Supervisor: Prof. Dr. Dr. J¨urgen Richter-Gebert Advisor: Prof. Dr. Dr. J¨urgen Richter-Gebert Submission date: September 17, 2012 Abstract Dynamic geometry is a branch of geometry concerned with movement of elements in geometric constructions. If we move an element of a construction and update the positions of the other elements accordingly, we traverse a continuum of instances of the construction, its configuration space. In many geometric constructions, interesting effects occur in the neighbourhood of geometric singularities, i.e. around points in configuration space at which ambiguous geometric operations of the construction branch. For example, at a geometric singularity two points may exchange their role of being the upper/lower point of intersection of two circles. On the other hand, we frequently see ambiguous elements jump unexpectedly from one branch to another in the majority of dynamic geometry software. If we want to explain such effects, we must consider complex coordinates for the elements of geometric constructions. Consequently, what often hinders us to fully comprehend the dynamics of geometric constructions is that our view is limited to the real perspective. In particular, the two-dimensional real plane exposed in dynamic geometry software is really just a section of the multidimensional complex configuration space which determines the dynamic behaviour of a construction. The goal of this thesis is to contribute to a better understanding of dynamic geom- etry by making the complex configuration space of certain constructions accessible to visual exploration in real-time.
    [Show full text]
  • An Introduction to the Riemann Hypothesis
    An introduction to the Riemann hypothesis Author: Alexander Bielik [email protected] Supervisor: P¨arKurlberg SA104X { Degree Project in Engineering Physics Royal Institute of Technology (KTH) Department of Mathematics September 13, 2014 Abstract This paper exhibits the intertwinement between the prime numbers and the zeros of the Riemann zeta function, drawing upon existing literature by Davenport, Ahlfors, et al. We begin with the meromorphic continuation of the Riemann zeta function ζ and the gamma function Γ. We then derive a functional equation that relates these functions and formulate the Riemann hypothesis. We move on to the topic of finite-ordered functions and their Hadamard products. We show that the xi function ξ is of finite order, whence we obtain many useful properties. We then use these properties to find a zero-free region for ζ in the critical strip. We also determine the vertical distribution of the non-trivial zeros. We finally use Perron's formula to derive von Mangoldt's explicit formula, which is an approximation of the Cheby- shev function . Using this approximation, we prove the prime number theorem and conclude with an implication of the Riemann hypothesis. Contents Introduction 2 1 The statement of the Riemann hypothesis3 1.1 The Riemann zeta function ζ .........................................3 1.2 The gamma function Γ.............................................4 1.3 The functional equation............................................7 1.4 The critical strip................................................8 2 Zeros in the critical strip 10 2.1 Functions of finite order............................................ 10 2.2 The Hadamard product for functions of order 1............................... 11 2.3 Proving that ξ has order at most 1.....................................
    [Show full text]
  • Visualizing Complex Functions Using Gpus
    Visualizing Complex Functions Using GPUs Khaldoon Ghanem RWTH Aachen University German Research School for Simulation Sciences Wilhelm-Johnen-Straße 52428 Jülich E-mail: [email protected] Abstract: This document explains some common methods of visualizing complex functions and how to imple- ment them on the GPU. Using the fragment shader, we visualize complex functions in the complex plane with the domain coloring method. Then using the vertex shader, we visualize complex func- tions defined on a unit sphere like spherical harmonics. Finally, we redesign the marching tetrahedra algorithm to work on the GPGPU frameworks and use it for visualizing complex scaler fields in 3D space. 1 Introduction GPUs are becoming more and more attractive for solving computationally demanding problems. This is because they are cheaper than CPUs in two senses. First, they provide more performance for less cost .i.e. cheaper GFLOPs. Second, they are more energy efficient .i.e. cheaper running times. This reduced cost comes at the expense of less general purpose architecture and hence a different program- ming model. There are two main ways of programming GPUs. The first one is used in the graphics community using shading languages like HLSL, GLSL and Cg. In these languages, the programmer deals with vertices and fragments and processes them with the so called vertex and fragment shaders, respectively. Actu- ally, this has been the only way of programming GPUs for a while. Fortunately, in the recent years, frameworks for general programming have been developed like CUDA and OpenCL. They are much more suited for expressing the problem more abstractly in terms of threads.
    [Show full text]
  • Vizualiz´Acia Komplexn´Ych Funkci´I Pomocou Riemannov
    UNIVERZITA KOMENSKEHO´ V BRATISLAVE FAKULTA MATEMATIKY, FYZIKY A INFORMATIKY RNDr. Miroslava Val´ıkova´ Autorefer´atdizertaˇcnejpr´ace VIZUALIZACIA´ KOMPLEXNYCH´ FUNKCI´I POMOCOU RIEMANNOVYCH´ PLOCH^ na z´ıskanie akademick´ehotitulu philosophiae doctor v odbore doktorandsk´ehoˇst´udia 9.1.7 Geometria a topol´ogia Bratislava 2014 Dizertaˇcn´apr´acabola vypracovan´av dennej forme doktorandsk´ehoˇst´udiana Katedre algebry, geometrie a didaktiky matematiky Fakulty matematiky, fyziky a informatiky Univerzity Komensk´ehov Bratislave. Predkladateˇl: RNDr. Miroslava Val´ıkov´a Katedra algebry, geometrie a didaktiky matematiky FMFI UK, Mlynsk´adolina, 842 48 Bratislava Skoliteˇ ˇl: doc. RNDr. Pavel Chalmoviansk´y,PhD. Katedra algebry, geometrie a didaktiky matematiky FMFI UK, Mlynsk´adolina, 842 48 Bratislava Oponenti: ................................................................... ................................................................... ................................................................... ................................................................... ................................................................... ................................................................... Obhajoba dizertaˇcnejpr´acesa kon´a..................... o ............. h pred komisiou pre obhajobu dizertaˇcnejpr´acev odbore doktorandsk´ehoˇst´udiavy- menovanou predsedom odborovej komisie prof. RNDr. J. Korbaˇs,CSc. Geometria and topol´ogia{ 9.1.7 Geometria a topol´ogia na Fakulte matematiky, fyziky a informatiky Univerzity
    [Show full text]
  • Extending Riemann Mapping Capabilities for the Sage Mathematics Package
    Extending Riemann Mapping Capabilities for the Sage Mathematics Package Ethan Van Andel Calvin College 2011 1 !4#04'#5' This senior honor project focused on expanding and improving the riemann module that I developed as part of research project at Calvin College in the summer of 2009. riemann is part of the Sage software system and provides computation and visualization tools for Riemann mapping—an important component of complex analysis. The riemann module is the only publicly available implementation of Riemann mapping. Its primary use is educational and explorational, although it could be applied to scientific research as well. This project involved a great deal of mathematical and computer science research to develop and implement the techniques needed. I have attempted to provide a thorough explanation of the relevant mathematical background. Nevertheless, readers without a background in complex analysis may be lost at times. Generally, if the reader finds terminology that they want defined more thorougly, Wikipedia is an effective resource for filling that gap. In addition, since this project is highly interwoven with the 2009 research, some of the work described in this report was actually done two years ago. I have tried to keep the focus on the work for this project, discussing previous work only when it is relevant to an understanding of the whole. To avoid unnecessary complications, strict distinctions will not be made as to the exact place of every bit of work. As a rough guide, the work that was done for this project includes some general and usability improvements, substantial optimization, the exterior Riemann map, the Ahlfors spiderweb, and analytic error testing.
    [Show full text]
  • Domain Coloring on the Riemann Sphere María De Los Ángeles Sandoval-Romero Antonio Hernández-Garduño
    The Mathematica® Journal Domain Coloring on the Riemann Sphere María de los Ángeles Sandoval-Romero Antonio Hernández-Garduño Domain coloring is a technique for constructing a tractable visual object of the graph of a complex function. The package complexVisualize.m improves on existing domain coloring techniques by rendering a global picture on the Riemann sphere (the compactification of the complex plane). Additionally, the package allows dynamic visualization of families of Möbius transformations. In this article we discuss the implementation of the package and illustrate its usage with some examples. ■ 1. Introduction Domain coloring is a technique that uses a color spectrum to compensate for the missing dimension in the graph of a complex function. It was first introduced in [1], with a detailed description of the complex-plane case in [2]. A general definition is found in [3]. More precisely, consider a function f : U → V between two sets U, V (for example, two complex manifolds). Choose a “color" function, κ : V → HSB, where HSB denotes the Hue-Saturation-Brightness standard color space. Next, for any z ∈ U, compute f(z) and evaluate the resulting color κ ◦ f(z), assigning this color to the preimage z. The final colored domain has all the information (through the color) needed to say where the point z ∈ U gets mapped. Of course, the effectiveness of domain coloring depends on the choice of an adequate color scheme κ. Geometrically, the HSB color space is identified with an inverted solid cone. In cylin- drical coordinates, HSB is parametrized by (s b cos θ, s b sin θ, b), where θ ∈ [0, 2 π], s ∈ [0, 1], and b ∈ [0, 1].
    [Show full text]
  • Complex Representations of Functions
    6 Complex Representations of Functions “He is not a true man of science who does not bring some sympathy to his studies, and expect to learn something by behavior as well as by application. It is childish to rest in the discovery of mere coincidences, or of partial and extraneous laws. The study of geometry is a petty and idle exercise of the mind, if it is applied to no larger system than the starry one. Mathematics should be mixed not only with physics but with ethics; that is mixed mathematics. The fact which interests us most is the life of the naturalist. The purest science is still biographical.” Henry David Thoreau (1817-1862) 6.1 Complex Representations of Waves We have seen that we can seek the frequency content of a function f (t) defined on an interval [0, T] by looking for the Fourier coefficients in the Fourier series expansion ¥ a0 2pnt 2pnt f (t) = + ∑ an cos + bn sin . 2 n=1 T T The coefficients take forms like 1 Euler’s formula can be obtained using x 2 Z T 2pnt the Maclaurin expansion of e : an = f (t) cos dt. T T ¥ xn 1 xn 0 ex = = 1 + x + x2 + ··· + + ··· . ∑ n! 2! n! However, trigonometric functions can be written in a complex ex- n=0 ponential form. This is based on Euler’s formula (or, Euler’s identity):1 We formally set x = iq Then, ¥ n iq (iq) eiq = cos q + i sin q. e = ∑ n=0 n! (iq)2 (iq)3 (iq)4 The complex conjugate is found by replacing i with −i to obtain = 1 + iq + + + + ··· 2! 3! 4! −iq (q)2 (q)3 (q)4 e = cos q − i sin q.
    [Show full text]