On the multidimensional distribution of the subset sum generator of pseudorandom numbers

We show that for a random choice of the parameters, the subset sum pseudorandom number generator produces a sequence of uniformly and independently distributed pseudorandom numbers. The result can be useful for both cryptographic and quasi-Monte Carlo applications and relies on bounds of exponential sums.

     

Arithmetic functions with linear recurrence sequences

We obtain asymptotic formulas for all the moments of certain arithmetic functions with linear recurrence sequences. We also apply our results to obtain asymptotic formulas for some mean values related to average orders of elements in finite fields.     

Computing Aviation Sparing Policies: Solving a Large Nonlinear Integer Program

Deployed US Navy aircraft carriers must stock a large number of spare parts to support the various types of aircraft embarked on the ship. The sparing policy determines the spares that will be stocked on the ship to keep the embarked aircraft ready to fly. Given a fleet of ten or more aircraft carriers and a cost of approximately 50 million dollars per carrier plus the cost of spares maintained in warehouses in the United States, the sparing problem constitutes a significant portion of the Navy’s resources. The objective of this work is to find a minimum-cost sparing policy that meets the readiness requirements of the embarked aircraft. This is a very large, nonlinear, integer optimization problem. The cost function is piecewise linear and convex while the constraint mapping is highly nonlinear. The distinguishing characteristics of this problem from an optimization viewpoint are that a large number of decision variables are required to be integer and that the nonlinear constraint functions are essentially “black box” functions; that is, they are very difficult (and expensive) to evaluate and their derivatives are not available. Moreover, they are not convex. Integer programming problems with a large number of variables are difficult to solve in general and most successful approaches to solving nonlinear integer problems have involved linear approximation and relaxation techniques that, because of the complexity of the constraint functions, are inappropriate for attacking this problem. We instead employ a pattern search method to each iteration of an interior point-type algorithm to solve the relaxed version of the problem. From the solution found by the pattern search on each interior point iteration, we begin another pattern search on the integer lattice to find a good integer solution. The best integer solution found across all interations is returned as the optimal solution. The pattern searches are distributed across a local area network of non-dedicated, heterogeneous computers in an office environment, thus, drastically reducing the time required to find the solution.     

Atmospheric aerosol limb scanning based on the lunar eclipses photometry

The work is devoted to the analysis of the surface photometric observations of two total lunar eclipses in 2004. The lunar surface relative brightness distribution inside the umbra was used to retrieve the vertical distribution of aerosol extinction of the solar radiation expanding by a tangent path and its dependence on the location at the limb of the Earth. The upper altitude of troposphere aerosol layer was estimated for different latitude zones. The correlation between additional aerosol extinction in the upper troposphere and cyclones was investigated.     

Synthesis of a new tricyclic 3-(tetrazol-5-yl)pyridine system from 2-(azidomethyl)nicotinonitriles

2-Methyl-3-cyanopyridines were converted into the corresponding 2-azidomethyl derivatives, which then underwent an intramolecular cycloaddition reaction. A novel heterocyclic system containing a 3-(tetrazol-5-yl)pyridine unit was obtained in this way.     

Support vector machine via nonlinear rescaling method

In this paper we construct the linear support vector machine (SVM) based on the nonlinear rescaling (NR) methodology (see [Polyak in Math Program 54:177–222, 1992; Polyak in Math Program Ser A 92:197–235, 2002; Polyak and Teboulle in Math Program 76:265–284, 1997] and references therein). The formulation of the linear SVM based on the NR method leads to an algorithm which reduces the number of support vectors without compromising the classification performance compared to the linear soft-margin SVM formulation. The NR algorithm computes both the primal and the dual approximation at each step. The dual variables associated with the given data-set provide important information about each data point and play the key role in selecting the set of support vectors. Experimental results on ten benchmark classification problems show that the NR formulation is feasible. The quality of discrimination, in most instances, is comparable to the linear soft-margin SVM while the number of support vectors in several instances were substantially reduced.     

Fractional Order Disturbance Observer for Robust Vibration Suppression

For the first time, the fractional order disturbance observer (FO-DOB) is proposed for vibration suppression applications such as hard disk drive servo control. It has been discovered in a recently published US patent application (US20010036026) that there is a tradeoff between phase margin loss and strength of the low frequency vibration suppression. Given the required cutoff frequency of the low pass filter, also known as the Q-filter, it turns out that the relative degree of the Q-filter is the major tuning knob for this tradeoff. The solution in US20010036026 was based on an integer order Q-filter with a variable relative degree. This actually motivated the use of a fractional order Q-filter. The fractional order disturbance observer is based on the fractional order Q-filter. The implementation issue is also discussed. The nice point of this paper is that the traditional DOB is extended to the fractional order DOB with the advantage that the FO-DOB design is now no longer conservative nor aggressive, i.e., given the cutoff frequency and the desired phase margin, we can uniquely determine the fractional order of the low pass filter.