Posted by: louisyangliu | March 5, 2009

Some Thoughts on Infinity Norm

The infinity norm may be one of the usual models to study about in different situations in mathematics, not the continuous ones, but the discrete or singular ones.

We know that ||(x,y)||_{\infty}:=max(|x|,|y|), for any (x,y)\in \mathbb{R}^{2}, can be expressed as an integral on S^{1} by spreading a measure on it, which is

||(x,y)||_{\infty}=c\int_{S^{1}}|x\cos\theta+y\sin\theta|(\delta_{\frac{\pi}{4}}(\theta)+\delta_{\frac{3\pi}{4}}(\theta))d\theta (1)

for some constant c because max(|x|,|y|)=\frac{1}{2}(|x+y|+|x-y|), where

\delta_{\frac{\pi}{4}}(\theta)=\begin{cases}<br /> +\infty, & \theta=\frac{\pi}{4}\\<br /> 0, & \theta\ne\frac{\pi}{4}.\end{cases} (2)

satisfying \int_{-\infty}^{\infty}\delta_{\frac{\pi}{4}}(\theta)\, dx=1 and

\delta_{\frac{3\pi}{4}}(\theta)=\begin{cases}<br /> +\infty, & \theta=\frac{3\pi}{4}\\<br /> 0, & \theta\ne\frac{3\pi}{4}.\end{cases} (3)

satisfying \int_{-\infty}^{\infty}\delta_{\frac{3\pi}{4}}(\theta)\, dx=1 are modified Dirac delta functions.

One might also be able to put an measure on S^{3} for the complex norm ||(z,w)||_{\infty}=max(|z|,|w|), (z,w)\in \mathbb{C}^{2}, which can be written in real norm as ||(x,y,u,v)||:=max(\sqrt{x^{2}+y^{2}},\sqrt{u^{2}+v^{2}}). A problem is what function f would satisfy

||(x,y)||_{\infty}=\int_{S^{3}}|x\cos\xi_{1}\cos\eta+y\cos\xi_{2}\sin\eta|f(\xi_{1},\xi_{2},\eta)d\xi_{1}d\xi_{2}d\eta (4)

if one parametrizes S^{3} by Hopf coordinates

z=e^{i\xi_{1}}\sin\eta,\, w=e^{i\xi_{2}}\cos\eta. (5)

One can show that the function f is actually some function independent of \xi_{1} and \xi_{2} by the invariance of the complex norm under U(1)\times U(1) action, so it can be denoted as f(\eta). But it is unknown yet what exactly the function f(\eta) is, because the invariance yields that f(\eta)=\delta_{\frac{\pi}{4}}(\eta)+\delta_{\frac{3\pi}{4}}(\eta), but this doesn’t produce an appropriate function in the general expression for the complex infinity norm. If one looks back the case of real infinity norm, the way of expressing the infinity norm is by making use of the absolute values on sum and difference smartly to sift out the maximum from two number. But in the complex case the double integral on the torus

\left\{ (\frac{\sqrt{2}}{2}e^{i\xi_{1}},\frac{\sqrt{2}}{2}e^{i\xi_{2}}):\xi_{1},\xi_{2}\in[0,2\pi]\right\} (6)

eliminates the feasibility of sifting out the maximal modulus.

Posted by: louisyangliu | October 5, 2008

On Stokes’ Theorem

Stokes’ theorem plays an important role in computing integrals of differential forms on manifolds, and its consequence, the Cauchy integral formula, is quite fundamental in complex analysis in one variable as well.

Suppose M is a compact 2-dimensional submanifold with smooth boundary in the plane \mathbb{R}^{2}. From symplectic geometry, the natural symplectic form on T^{*}\mathbb{R}^{2} can be written as \omega=d\alpha, where \alpha=\xi dx+\eta dy. If one computes the integral,

\int_{D^{*}(\partial M)}\alpha\wedge\omega=\int_{D^{*}(\partial M)}\xi dx\wedge dy\wedge d\eta+\int_{D^{*}(\partial M)}\eta dy\wedge dx\wedge d\xi(1)

where D^{*}(\partial M) denotes the codisk bundle of \partial M, since \partial M is an one dimensional curve in \mathbb{R}^{2}, it makes the integral containing dx\wedge dy measuring the perturbations of base points of cotangent vectors in D^{*}(\partial M) be zero. Hence each term of the right hand side of (1) is zero. So one has that \int_{D^{*}(\partial M)}\alpha\wedge\omega=0, which is also true in general when M is a compact hypersurface contained in some affine n-1-plane in \mathbb{R}^{n}.

By the way, Gowers wrote a post, one way of looking at Cauchy theorem, which is very interesting and tells us that the Cauchy’s theorem is the natural generalization in 2-dimension of the statement f'=0\Longrightarrow f is a constant for f:\mathbb{R\rightarrow R}, that, we know, is a direct consequence of the second fundamnetal theorem of calculus which is the special case of the Stokes’ theorem in one dimension. So in this sense, Stokes’ theorem is the “generator” of all the theorems and statements of this kind.

Posted by: louisyangliu | July 9, 2008

Talk on an Aggregation Model of Marine Paticles

I’ll talk about an aggregation model of marine particles in VIGRE graduate student seminar this Friday. This is one part of joint work with Adrian Burd. I’ll talk about the moments for log-normal distribution, its application on size distribution of marine particles, the multi-modes model, and some results from numerical experiments. Here is our Slides for the talk. By the way, the picture by Adrian in the slides, which I named as “sunrise or sunset” , is actually sunset, he took it from the Gulf of Mexico.

Posted by: louisyangliu | July 1, 2008

Infinitesimal Calculus, Rocket, and Phoenix Mars Lander

As a fan of the Phoenix Mars Lander, I recently paid close attention to the news about it since it successfully landed around the north pole of the Mars. So far the good news from Mars to human beings are the finding of ice water, friendly soil environment for some plants, and other things.

The course of development in space technology is a great history. From Newton’s invention of infinitesimal calculus in mathematics in 1665 to his three laws of motion; from the first manned aeroplane by Wright Brothers in 1903, the theoretic foundation of which is the lift equation which basically says the lift force is directly proportional to the square of velocity of object and the area of its wings, L=kv^2A, and had appeared in Newton’s “The Mathematical Principles of Natural Philosophy”, to rockets launching and sending satellites to space; from the first manned landing on the moon in 1969 to the Phoenix Mars Lander nowadays, and so on, what great strides human have made.

In MATLAB, one can obtain solutions to differential equations in symbolic or numerical manner. For simple first order ODEs, one can use the command “dsolve” to solve solvable equations analytically. In particular, it can be used to obtain anti-derivatives. Numerically, one can assign numerical initial conditions, using “ode15″ for some stiff ODEs, where the stiffness is determined by the Jacobian of the ODE system, or other similar commands, depending on the type of equation. The steps usually includes defining functions/equations in M-files, setting numerical domains and step sizes for basic variables and initial conditions, applying suitable solver codes, and graphing the numerical results, and so forth.

For instance, one can build a model for simulating the coagulation of particles in aqua-environment. By the central limit theorem in classic probability theory, the distribution of particles in different sizes can be considered as  a log-normal distribution, or a composition of log-normal distributions

N_i(D) = \frac{m_i}{D \sigma_i \sqrt{2 \pi}}e^{-\frac{(\ln (D) - \mu_i)^2}{2\sigma_i^2}},     (1)

where D is the diameter of particle in ocean environment, \mu_i and \sigma_i are parameters of different modes of log-normal distributions. By the theory from particle physics, the interaction of different particles has a coagulation kernel, such as Brownian kernel

\beta_{Br}(D,\tilde{D})=12\pi(D+\tilde{D})      (2)

for two particles with diameters D and \tilde{D}. Now let us consider a single-mode model based on Koziol-Leighton’s result on aerosol dynamics, as the following ODE for this modeling,

\frac{dm(t)}{dt}=\frac{1}{2\mu e^{\frac{\sigma^2}{2}}}\int_{0}^{\infty}\int_{0}^{\infty}\beta_{Br}(D,\tilde{D})N(D,t)N(\tilde{D},t)(\tilde{(D}^{3}+D^{3})^\frac{1}{3}-D)dDd\tilde{D},

where N (D,t)=\frac{m(t)}{D \sigma \sqrt{2 \pi}}e^{-\frac{(\ln (D) - \mu )^2}{2\sigma^2}} and similarly for N (\tilde{D},t). Given initial data to m, we can solve out m(t) by solver “ode23s” or “ode15s”, and obtain a simple simulation for particles coagulation.

Posted by: louisyangliu | May 18, 2008

Earthquake Relief and Prayer

A strong earthquake hit eastern Sichuan of China six days ago. Thousands of lives were lost in the quake, many of them were elementary or high school students under seventeen. Let’s pray with Dion’s song “The Prayer” for the people suffered from the disaster. Please lend a helping hand, donations can be made at UNICEF or Red Cross.

I was attending a conference on geometric group theory, geometric analysis, and mapping class group in the Japan-U.S. Mathematics Institute (JAMI) of JHU at Baltimore last weekend. JAMI with a history of about twenty years holds conference every year for promoting communications and cooperations of American mathematicians and Japanese mathematical society.

Mapping class group \mathcal{M}, the group of orientation-preserving self-homeomorphisms of a topological space, becomes another interesting algebraic object in recent decades besides homology groups as algebraic invariants for topological spaces. Elements in mapping class group induce actions on homology groups of a space, the kernel called Torelli group \mathcal{T} is important for understanding mapping class group because of the following short exact sequence

1 \to \mathcal{T}(\Sigma_{g}) \to \mathcal{M}(\Sigma_{g}) \to Sp(2g, \mathbf{\mathbb{Z})} \to 1    (1)

where Sp(2g,\mathbb{Z}) is the symplectic group of degree 2g over \mathbb{Z}, and \Sigma_{g} is a genus g surface. There are many interesting investigations on \mathcal{T} and the so-called Johnson homomorphism that gives the free part of the abelianization of \mathcal{T}.

Posted by: louisyangliu | April 7, 2008

A Petition to Help Support Maths at USQ

Please go to Terry Tao’s Blog, and give your endorsement to his petition to help support maths at University of Southern Queensland (USQ) by Signing here.

Older Posts »

Categories