The transformed $l_1$ penalty (TL1) functions are a one parameter family of bilinear transformations composed with the absolute value function. When acting on vectors, the TL1 penalty interpolates $l_0$ and $l_1$ similar to $l_p$ norm, where $p$ is in $(0,1)$. In our companion paper, we showed that TL1 is a robust sparsity promoting penalty in compressed sensing (CS) problems for a broad range of incoherent and coherent sensing matrices. Here we develop an explicit fixed point representation for the TL1 regularized minimization problem. The TL1 thresholding functions are in closed form for all parameter values. In contrast, the $l_p$ thresholding functions ($p$ is in $[0,1]$) are in closed form only for $p=0;1;1=2;2=3$, known as hard, soft, half, and 2/3 thresholding respectively. The TL1 threshold values differ in subcritical (supercritical) parameter regime where the TL1 threshold functions are continuous (discontinuous) similar to soft-thresholding (half-thresholding) functions. We propose TL1 iterative thresholding algorithms and compare them with hard and half thresholding algorithms in CS test problems. For both incoherent and coherent sensing matrices, a proposed TL1 iterative thresholding algorithm with adaptive subcritical and supercritical thresholds (TL1IT-s1 for short), consistently performs the best in sparse signal recovery with and without measurement noise.

In this paper, a structure-preserving algorithm is developed for the computation of a semi-stabilizing solution of a Generalized Algebraic Riccati Equation (GARE). The semi-stabilizing solution of GAREs has been used to characterize the solvability of the (J, J ′ )-spectral factorization problem in control theory for general rational matrices which may have poles and zeros on the extended imaginary axis. The main difficulty in solving such a GARE lies in the fact that its associated Hamiltonian/skew-Hamiltonian pencil has eigenvalues on the extended imaginary axis. Consequently, it is not clear which eigenspace of the associated Hamiltonian/skew-Hamiltonian pencil can characterize the desired semi-stabilizing solution. That is, it is not clear which eigenvectors and principal vectors corresponding to the eigenvalues on the extended imaginary axis should be contained in the eigenspace that we wish to compute. Hence, the well-known generalized eigenspace approach for the classical algebraic Riccati equations cannot be employed directly. The proposed algorithm consists of a structure-preserving doubling algorithm (SDA) and a postprocessing procedure to determine the desired eigenvectors and principal vectors corresponding to the purely imaginary and infinite eigenvalues. Under mild assumptions, linear convergence of rate 1/2 for the SDA is proved. Numerical experiments illustrate that the proposed algorithm performs efficiently and reliably.

Thetransmissioneigenvalueproblem,besidesitscriticalroleininversescattering problems, deserves special interest of its own due to the fact that the corresponding differen- tial operator is neither elliptic nor self-adjoint. In this paper, we provide a spectral analysis and propose a novel iterative algorithm for the computation of a few positive real eigen- values and the corresponding eigenfunctions of the transmission eigenvalue problem. Based on approximation using continuous finite elements, we first derive an associated symmetric quadratic eigenvalue problem (QEP) for the transmission eigenvalue problem to eliminate the nonphysical zero eigenvalues while preserve all nonzero ones. In addition, the derived QEP enables us to consider more refined discretization to overcome the limitation on the number of degrees of freedom. We then transform the QEP to a parameterized symmet- ric definite generalized eigenvalue problem (GEP) and develop a secant-type iteration for solving the resulting GEPs. Moreover, we carry out spectral analysis for various existence intervals of desired positive real eigenvalues, since a few lowest positive real transmission eigenvalues are of practical interest in the estimation and the reconstruction of the index of refraction. Numerical experiments show that the proposed method can find those desired smallest positive real transmission eigenvalues accurately, efficiently, and robustly.

We consider the solution of the large-scale nonsymmetric algebraic Riccati equation XCX − XD − AX + B = 0, with M ≡ [D,−C;−B,A] ∈ R(n1+n2)×(n1+n2) being a nonsingular M-matrix. In addition, A and D are sparselike, with the products A−1u, A−⊤u, D−1v, and D−⊤v computable in O(n) complexity (with n = max{n1, n2}), for some vectors u and v, and B, C are low ranked. The structure-preserving doubling algorithms (SDA) by Guo, Lin, and Xu [Numer. Math., 103 (2006), pp. 392–412] is adapted, with the appropriate applications of the Sherman–Morrison– Woodbury formula and the sparse-plus-low-rank representations of various iterates. The resulting large-scale doubling algorithm has an O(n) computational complexity and memory requirement per iteration and converges essentially quadratically. A detailed error analysis, on the effects of truncation of iterates with an explicit forward error bound for the approximate solution from the SDA, and some numerical results will be presented.

Starting from the roots of the equation 3^x + 4^x = 5^x discussed in the Tenth Grade, we utilized the method, the algebraic extension of the rational number eld, to produce the ways to judge whether the roots of the basic exponential equation with a form as a^x+b^x = c^x are rational or not. For equations with more than two terms on the left side, as is the equation a_1^x+a_2^x+
···

+a_n^x = d^x, the determination of whether the root was irrational was comparatively dicult.Therefore, we provided a prevalent method for the examination of the root of a three-term equation as well as a conclusion that if the equation doesn't have integer roots, the roots won't be a rational number with a denominator of two. Finally, based on the method, the algebraic extension of the rational number eld, we concluded that under special occasions, the root of the equation can't be some rational numbers with certain denominators.

Pell equation is an important research object in elementary number theory of indefinite equation. its form is simple, but it is rich in nature. Many number theory problems can be transformed into the problem of Pell equation’s solvability. However, the previous methods in determining the Pell equation’s solvability are sophisticated for calculation, which leads to the lack of efficiency. This paper gives new and more widely used methods to determine the solvability of Pell equation, including several necessary conditions, sufficient conditions and necessary and sufficient conditions.

I heard the concept of the ”equal‐sum point” in a math summer camps and found it very interesting. After searching on the Internet,I found that the ”equal‐sum point” was associated with Soddy point. Soddy point is a magic point in plane geometry, which was found by British physicist, chemist Frederick.Soddy. Some research about the Soddy point has been done in foreign countries, but the properties are not comprehensive. In China, Teacher HuasongHuang raised the concept of the “equal‐sum point” and the “equal‐difference point” and drew some properties, I found that these two points were very similar to the Soddy point, but they did not link the two points with Soddy point. In this article ,I will connect the Soddy point with the “equal‐sum point” and the “equal‐difference point” and study their properties more deeply and find some new discoveries.

We study the super-resolution (SR) problem of recovering point sources consisting of a collection of isolated and suitably separated spikes from only the low frequency measurements. If the peak separation is above a factor in $(1, 2)$ of the Rayleigh length (physical resolution limit), L1 minimization is guaranteed to recover such sparse signals. However, below such critical length scale, especially the Rayleigh length, the $L_1$ certificate no longer exists. We show several local properties (local minimum, directional stationarity, and sparsity) of the limit points of minimizing two $L_1$ based nonconvex penalties, the difference of $L_1$ and $L_2$ norms $(L_{1-2})$ and capped $L_1 (CL_1)$, subject to the measurement constraints. In one and two dimensional numerical SR examples,the local optimal solutions from dierence of convex function algorithms outperform the global L1 solutions near or below Rayleigh length scales either in the accuracy of ground truth recovery or in nding a sparse solution satisfying the constraints more accurately.

Free kicks play an important role in a modern soccer match. A skillful soccer player usually uses the opportunity of a free kick to score a goal. There are mainly two ways to give a free kick, falling ball and banana ball. The last thesis analyzes the path of falling soccer in two-dimensional plane. This thesis mainly focuses on using geometrical ways to generate a graph in three-dimensional coordinate system to analyze the path of banana ball and find out the most suitable distance to give a free kick.

Based on Leslie Matrix, the paper tries to predict the population fluctuation of Nanjing by considering the regional distinction and latest policy impacts, and hereby comes up with four conclusions: a) The population of Nanjing will see a relatively rapid growth between 2013 and 2020; b) Given the situation, the elderly population of the city will experience a much faster increase in the next 40 years, causing an elderly bulge that not only poses heavy burden to social welfare but also gives rise to a disproportionate population structure; c) The proposed retirement age extension in theory, will reduce the pressure on pension system by retaining active labors, and thus alleviating the problems of aging population to some extent; d) The newly-issued Two-Child Policy will not change the trends in or transform the population structure, therefore its effects on population aging is rather limited.

This paper proposes an Integrated Obesity Index (IOI) based on height, weight, waist, and hip, physical indexes that can be accessed through regular physical examinations, on top of common obesity standards, including Body Mass Index (BMI), Waist‐Hip Ratio (WHR), and Body Fat Percentage (BFP). Integrated Obesity Index is aimed to be feasible and convenient under application, as well as comprehensive, compensating for the limitations of existing classifications of obesity. Thus, this paper proposes and examines the application formula of the IOI by establishing geometric models. Adopting the data from epidemiological survey, this papers brings physical indicators, gender, BFP and other variables into the model for data analysis, applies linear regression on the resulting images, and obtains the IOI cut‐off point. By conducting comprehensive data analyses to various obesity classifications, the paper demonstrates that the IOI has significant advantages over BMI, WHP, or BFP and will compensate for the deficiencies of the existing obesity classifications to a certain extent in application.

This article applies functional knowledge and computer software to describe the relation between people distribution and evacuation efficency and safety in high-rise buildings.The people distribution contains two parts:the allocation of people in different storeys and the distribution of people in a certain storey. Through this mathematical model,we attempt to find the optimum plan of people distribution and use it to guide the distribution of different functional areas in high-rise buildings.

It is well known that group is a set with an algebra operation, which is often denoted by G . A simple group is the group which has only the 2 trivial normal subgroups. The importance of the simple groups may analogue to the prime numbers in the number theory. The classification of the finite simple groups was the central problem of 20th century’s’ algebra. In this research report, we use Sylow theorem, Burnside theorem, group action and some other elementally group theory methods to obtain some criteria of simple group . By using these criteria we determine the possible simple groups whose order less than 700. Generalizes the results in the reference [1]，[2]， [3]，[4].We conclude that Theorem A group G of order less or equal than 700 could not be simple except for |G|∈{60,168,360,504,660 and all primes}.

Computational techniques derived from digital signal processing are playing a significant role in the security and digital copyrights of audio, video, and visual arts. In light of the quantum computing, the corresponding algorithms are becoming a new research direction in today’s high-technology world. The nature of quantum computer guarantees the security of quantum data, so a safe and effective quantum watermarking algorithm is in demand. Quantum watermarking is the technique that embeds the invisible quantum signal into quantum multimedia data for copyright protection. This research presents a quantum watermarking algorithm in discrete M-band wavelet domain. Different from most traditional algorithms, we propose a new algorithm which applies quantum watermarking in M-band wavelet domain. Assured by the Heisenberg uncertainty principle and quantum no-cloning theorem, the security of quantum watermark can reach a very high-level standard. In other words, this watermarking algorithm can defeat nearly all attackers, no matter using classical computer or quantum computer.

The aim of this study is to figure out the optimization of survival probability of the generals in the card game "Legends of the Three Kingdoms". First, based on combinatorics and probability theory, the PASCAL program was performed to analyze the probabilities of success when several popular generals perform their special skills in the game. The results are as follows: Second, seven typical generals, including Zhou Tai, Zhang Jiao, Zhen Ji, Lu Xun, Ma Chao, Xiahou Dun, and Guo Jia, were chosen, and their abilities including attacking and defense were arbitrarily quantified according to a model constructed in the present study. The scores are as follows: Card-drawing Changing Attacking Defense Blood-selling Zhen Ji 5.62 Zhou Tai 1.04 Ma Chao 1.68 Xiahou Dun 1.19 Guo Jia 3.36 Lu Xun 3.55 Zhang Jiao 2.55 The abilities of different generals were sorted, and by comparison with the statistical record online, the reliability of the model was proved. Thus, the irrationality of the game is identified in this study, and we recommend the game designer to modify the abilities of some generals and to obtain an optimized value according to the output of the model, by which the game would be more popular.

Cartographical projection is mathematical cartography. Map projection is the mathematical model of geoid. Gauss-Kruger projection is unable to realize seamless splicing between adjoining sheet maps. While Mercator projection, especially web Mercator projection, is able to realize seamless splicing but has larger measurement error. How to combine with the advantage of the two projections, apply the method of map projection transformation and realize accurate measurement under the framework of a global map, is the purpose of thisarticle. In order to demonstrate that the problem of measurement could be solved by map projection transformation, we studied the mathematical principle of Gauss - Kruger Projection and the Mercator Projection and designed the algorithm model based on characteristics of both projections. And established simulated data in ArcMap, calculated the error of the three indicators-the area, the length and the angle- under areas in different latitude(0°-4°、30°-34°、60°-64°)and longitude (120°-126°)range which validated the algorithm model. The conclusion suggests that it is feasible that using web Mercator Projection is capable to achieve a global map expression. Meanwhile, calculations of area, distance and angle using map projection transformation principle are feasible as well in Gauss - Kruger Projection. This algorithm combines the advantage of the two different kinds of projection, is able to satisfy uses' high demand and covers the shortage of measurement function of online map services in China. It has certain practical value.

We study homogenization of G-equation with a ow straining term (or the strain G-equation) in two dimensional periodic cellular ow. The strain G-equation is a highly non-coercive and non-convex level set Hamilton-Jacobi equation. The main objective is to investigate how the ow induced straining (the nonconvex term) in uences front propagation as the ow intensity A increases. Three distinct regimes are identified. When A is below the critical level, homogenization holds and the turbulent ame speed sT (effective Hamiltonian) is well-defined for any periodic ow with small divergence and is enhanced by the cellular ow as $s_T \ge O(A/logA)$. In the second regime where A is slightly above the critical value, homogenization breaks down, and $s_T$ is not well defined along any direction. Solutions become a mixture of fast moving part and a stagnant part. When $A$ is sufficiently large, the whole ame front ceases to propagate forward due to the flow induced straining. In particular, along directions $p = (1; 0)$ and $(0;1)$, $s_T$ is well-defined again with a value of zero (trapping). A partial homogenization result is also proved. If we consider a similar but relatively simpler Hamiltonian, the trapping occurs along all directions. The analysis is based on the two-player dierential game representation of solutions, selection of game strategies and trapping regions, and construction of connecting trajectories.

Dijkstra and Bellman-ford are the two basic methods to solve Single-source shortest paths problem，and SPFA is the optimized version of Bellman-ford. However, since the algorithm time complexity is 0 (kn), SPFA’s time efficiency has some instabilities. This paper will compare the time efficiency between Dijkstra optimized with heap and three kinds of SPFA in different density graphs. Meanwhile, random data is created to test optimized Dijkstra and 3 types of SPFA. Finally, according to the calculation result of two methods’ time efficiency, some suggestion is given for their application under different situation.

The problem we study in this paper originated from a simple high school math competition (AMC 12) question on a geometric probability model: “A frog makes 3 jumps, each exactly 1 meter long. The directions of the jumps are chosen independently at random. What is the probability that the frog’s final position is no more than 1 meter from its starting position?” [1] From here, we changed the number of jumps the frog made, the length of each jump and the dimensions of the space. Using recursion, we obtained the recursive formula of the probability of a frog landing x meters within its original spot after jumping m 1-meter jumps in an N -dimension space by integrating the probability distribution of the former step. We also suggested a concept “the intensity of the probability field” to describe the relative probability of each spot and gave an expression of the intensity of each spot of a 2-dimensional space when the frog makes m jumps.

The log-optimal strategy in gambling and asset allocation problem attempts to maximize the expectation of the logarithmic rate of return but not the gross wealth itself. This strategy has been shown to have long termsuperiority overother strategies in theory. Another remarkable virtue is it allows updating the real-time market information quickly which could increase of the logarithmic rate of return. We will prove these properties ourselves in a simple way. We further consider how to do asset allocation by applying the log-optimal strategy in practice. The distribution function of the returns is usually unknown and need to be estimated from the history in a real world. The estimation error will affect the asset allocation directly, but such effect mayresult in a "butterfly effect" which could bring an investment disaster. Therefore, it is an important issue to answer whether there is logoptimal strategy resisted "butterfly effect". We develop a new log-optimal strategy based on the allocation utility which is a quadratic approximation to the expectation of the logarithmic rate of return. We show that the upper bound of the allocation utility deviation can be controlled by the l_1-norm of portfolio and l_∞-norm of the bias of variance matrix estimator. This turns out that our proposed strategy is robust. We apply our strategy for NYSE data. The results showed that our strategy has high performance in return.

In this report, we have mainly developed an abstract framework of k bit (t,n) secret sharing framework for cloud storage. Such framework is clean and it can be implemented. A successful implementation of the framework would provide users with protection when the system is under the attack on its confidentiality,integrity and reliability. Furthermore, such system has its own encryption by using permutations and tailor made error detection, location and data rescue. We make use of Lagrange polynomials and take the advantages of the algebraic property “t distinct points on the plane can uniquely determine a polynomial function of degree t-1” to design a k bit (t,n) -secret sharing distributed storage. We employ the set with unique factorization property (UFP) so that we simply need to calculate the y intercept of a Lagrange polynomial and then use a look up table to recover a secret.Moreover, the set which has minimum UFP would help us to design storage with smallest containers. In addition to the algebraic methods, we can utilize the geometric facts that three non collinear points determine a unique circle and four non coplanar points determine a unique sphere to construct k bit (3,n) and k bit (4,n) secret sharing storage respectively. To generalize to arbitrary case, it is straight forward if we have defined the Haar measure on the higher dimensional unit sphere. The last method is an application of Chinese Remainder Theorem (CRT) and we have designed k bit (t,n) secret sharing distributed storage and one of the designs can produce containers with the half size of the original secret. However, such k is no longer unrestricted and it has to be chosen from an interval. We have developed a C program for implementing both algebraic, geometric k bit (3,n) secret sharing distributed storages as well as the CRT method. The performances of both algebraic and geometric designs are satisfactory in term of processing time and compressed container size. The container size is even half of the size of the original secret in the CRT case and it is also very speedy.

The diffusive transport of passive tracers or particles can be enhanced by incompressible, turbulent flow fields. Analyzing the effective behavior is a challenging problem with theoretical and practical importance in many areas of science and engineering, ranging from the transport of mass, heat, and pollutants in geophysical flows to sea ice dynamics and turbulent combustion. The long time, large scale behavior of such systems is equivalent to an enhanced diffusion process with an effective diffusivity tensor D*. Two different formulations of integral representations for D* were developed for the case of time-independent fluid velocity fields, involving spectral measures of bounded self-adjoint operators acting on vector fields and scalar fields, respectively. Here, we extend both of these approaches to the case of space-time periodic velocity fields, allowing for chaotic dynamics, providing rigorous integral representations for D* involving spectral measures of unbounded self-adjoint operators.We prove the different formulations are equivalent. Their correspondence follows from a one-to-one isometry between the underlying Hilbert spaces. We also develop a Fourier method for computing D*, which captures the phenomenon of residual diffusion related to Lagrangian chaos of a model flow. This is reflected in the spectral measure by a concentration of mass near the spectral origin.

In this paper, we prove that the ODE system, $ \dot{x}=sinz+cosy$ $ \dot{y}=sinx+cosz$ $ \dot{z}=siny+cosx$ whose right-hand side is the Arnold-Beltrami-Childress (ABC) flow with parameters A = B = C = 1, has periodic orbits on $(2 \pi T)^3$ with rotation vectors parallel to (1, 0, 0), (0, 1, 0), and (0, 0, 1). An application of this result is that the well-known G-equation model for turbulent combustion with this ABC flow on $R^3$ has a linear (i.e., maximal possible) flame speed enhancement rate as the amplitude of the flow grows.

No abstract uploaded!

We show the existence of complete, asymptotically flat Cauchy initial data for the vacuum Einstein field equations, free of trapped surfaces, whose future development must admit a trapped surface. Moreover, the datum is exactly a constant time slice in Minkowski space-time inside and exactly a constant time slice in Kerr space-time outside. The proof makes use of the full strength of Christodoulou’s work on the dynamical formation of black holes and Corvino-Schoen’s work on the construction of initial data sets.