A search technique for global optimization in a chaotic environment

We describe a new algorithm which uses the trajectories of a discrete dynamical system to sample the domain of an unconstrained objective function in search of global minima. The algorithm is unusually adept at avoiding nonoptimal local minima and successfully converging to a global minimum. Trajectories generated by the algorithm for objective functions with many local minima exhibit chaotic behavior, in the sense that they are extremely sensitive to changes in initial conditions and system parameters. In this context, chaos seems to have a beneficial effect: failure to converge to a global minimum from a given initial point can often be rectified by making arbitrarily small changes in the system parameters.     

The Kissing Triangles: The Aesthetics of Mathematical Discovery

Papert's (1978) appeal to reconsider the power and possibilities of the aesthetic in mathematics learning is often ignored in mathematics education research. This paper begins with the premise, put forth by Dewey (1934), that the aesthetic structures many dimensions of inquiry and experience. In the same way that using particular paintings, musical compositions, or even everyday experiences has been instrumental to attempts by philosophers to understand the aesthetic dimensions of meaning and experience in artistic domains, I propose that analysing a particular encounter with mathematics may help reveal the nature and role of the often nebulous responses of elegance, beauty, and `fit' to which mathematicians lay claim in their mathematical activity. To achieve this, I draw on and adapt the defining features of the aesthetic character of experience set forth by the aesthetician Beardsley (1982). This, in turn, sheds light on the role thataesthetics can play in mathematical inquiry and experience, and provides initial categories and conjectures that can be used to investigate the potential roles of aesthetics in mathematics learning contexts.This revised version was published online in September 2005 with corrections to the Cover Date.     

ON HYPERBOLIC TIME DISCOUNTING IN EXHAUSTIBLE RESOURCE MODELS: AN APPLICATION TO WORLD OIL RESOURCES

ABSTRACT. Recent research on discounting in long term economic models involves hyperbolic discounting, in which the marginal discount rate shrinks as time passes. To investigate hyperbolic discounting and exhaustible resource allocation, this work develops a discrete‐time world oil model and model solution procedure, then uses the model to examine the consequences of adopting conventional (constant annual) discounting when hyperbolic discounting is appropriate, of adopting one hyperbolic discount rate path when a different hyperbolic path is appropriate, and of adopting hyperbolic discounting when conventional discounting is appropriate. Five conventional and two hyperbolic discount rate paths are considered. One hyperbolic path is that used by Nordhaus and Boyer [2000]; the other is that recommended by Weitzman [2001]. The generality of the findings is also assessed.     

Index Heuristics for Multiclass M/G/1 Systems with Nonpreemptive Service and Convex Holding Costs

We consider the optimal service control of a multiclass M/G/1 queueing system in which customers are served nonpreemptively and the system cost rate is additive across classes and increasing convex in the numbers present in each class. Following Whittle's approach to a class of restless bandit problems, we develop a Langrangian relaxation of the service control problem which serves to motivate the development of a class of index heuristics. The index for a particular customer class is characterised as a fair charge for service of that class. The paper develops these indices and reports an extensive numerical investigation which exhibits strong performance of the index heuristics for both discounted and average costs.     

二维椭圆型方程反问题中优化算法的比较

近年来，决定椭圆型方程系数反问题在地磁、地球物理、冶金和生物等实际问题上有着广泛的应用．本文讨论了二维的决定椭圆型方程系数反问题的数值求解方法．由误差平方和最小原则，这个反问题可化为一个变分问题，并进一步离散化为一个最优化问题，其目标函数依赖于要决定的方程系数．本文着重考察非线性共轭梯度法在此最优化问题数值计算中的表现，并与拟牛顿法作为对比．为了提高算法的效率我们适当选择加快收敛速度的预处理矩阵．同时还考察了线搜索方法的不同对优化算法的影响．数值实验的结果表明，非线性共轭梯度法在这类大规模优化问题中相对于拟牛顿法更有效．     

Polyelectrolyte complexes. VI. Polycation structure,polyanion molar mass,and polyion concentration effects on complex nanoparticles based on poly(sodium 2‐acrylamido‐2‐methylpropanesulfonate)

The formation and characterization of some interpolyelectrolyte complex (IPEC) nanoparticles based on poly(sodium 2‐acrylamido‐2‐methylpropanesulfonate) (NaPAMPS), as a function of the polycation structure, polyanion molar mass, and polyion concentration, were followed in this work. Poly(diallyldimethylammonium chloride) and two polycations (PCs) containing (N,N‐dimethyl‐2‐hydroxypropyleneammonium chloride) units in the backbone (PCA₅ and PCA₅D₁) were used as starting polyions. The complex stoichiometry, (n^?/n⁺)_iso, was pointed out by optical density at 500 nm (OD₅₀₀), polyelectrolyte titration, and dynamic light scattering. IPEC nanoparticle sizes were influenced by the polycation structure and polyanion molar mass only before the complex stoichiometry, which was higher for the more hydrophilic polycations (PCA₅ and PCA₅D₁) and for a higher NaPAMPS molar mass, and were almost independent of these factors after that, at a flow rate of the added polyion of about 0.28 mL × (mL PC)^?1 × h^?1. The IPEC nanoparticle sizes remained almost constant for more than 2 weeks, both before and after the complex stoichiometry, at low concentrations of polyions. NIPECs as stable colloidal dispersions with positive charges in excess were prepared at a ratio between charges (n^?/n⁺) of 0.7, and their storage colloidal stability, as a function of the polycation structure and polyion concentration (from 0.8 to ca. 7.8 mmol/L), was demonstrated. © 2004 Wiley Periodicals, Inc. J Polym Sci Part A: Polym Chem 42: 2495–2505, 2004     

Viscoelastic behavior of polypropylene/nitrile rubber thermoplastic elastomer blends: Application of Kerner's models for reactively compatibilized and dynamically vulcanized systems

The viscoelastic properties of binary blends of nitrile rubber (NBR) and isotactic polypropylene (PP) of different compositions have been calculated with mean‐field theories developed by Kerner. The phase morphology and geometry have been assumed, and experimental data for the component polymers over a wide temperature range have been used. Hashin's elastic–viscoelastic analogy principle is used in applying Kerner's theory of elastic systems for viscoelastic materials, namely, polymer blends. The two theoretical models used are the discrete particle model (which assumes one component as dispersed inclusions in the matrix of the other) and the polyaggregate model (in which no matrix phase but a cocontinuous structure of the two is postulated). A solution method for the coupled equations of the polyaggregate model, considering Poisson's ratio as a complex parameter, is deduced. The viscoelastic properties are determined in terms of the small‐strain dynamic storage modulus and loss tangent with a Rheovibron DDV viscoelastometer for the blends and the component polymers. Theoretical calculations are compared with the experimental small‐strain dynamic mechanical properties of the blends and their morphological characterizations. Predictions are also compared with the experimental mechanical properties of compatibilized and dynamically cured 70/30 PP/NBR blends. The results computed with the discrete particle model with PP as the matrix compare well with the experimental results for 30/70, 70/30, and 50/50 PP/NBR blends. For 70/30 and 50/50 blends, these predictions are supported by scanning electron microscopy (SEM) investigations. However, for 30/70 blends, the predictions are not in agreement with SEM results, which reveal a cocontinuous blend of the two. Predictions of the discrete particle model are poor with NBR as the matrix for all three volume fractions. A closer agreement of the predicted results for a 70/30 PP/NBR blend and the properties of a 1% maleic anhydride modified PP or 3% phenolic‐modified PP compatibilized 70/30 PP/NBR blend in the lower temperature zone has been observed. This may be explained by improved interfacial adhesion and stable phase morphology. A mixed‐cure dynamically vulcanized system gave a better agreement with the predictions with PP as the matrix than the peroxide, sulfur, and unvulcanized systems. © 2004 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 42: 1417–1432, 2004     

Relationship between the dynamic softening transition and wet sliding friction of elastomer compounds

For properly chosen elastomer compounds, thermorheological characterization is combined with an examination of the variation of the wet sliding friction with temperature. A conceptual argument leads to the assumption that the wet sliding friction should maximize at the energy dissipation peak associated with the dynamic softening transition at a characteristic frequency determined by the sliding speed and the effective smallest surface asperity scale. The dynamic softening transition is characterized with the peak in tan δ/G′ⁿ, where tan δ is the loss tangent, G′ is the elastic modulus, and n is a constant between 0 and 1. The William–Landel–Ferry transform is uncritically applied for extrapolating the position of the peak in tan δ/G′ⁿ at high frequencies. Even based on the criterion of tan δ, the results obtained on a concrete surface indicate that the effective smallest asperity scale is of order of 100 μm. © 2004 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 42: 2467–2478, 2004     

Searching in trees

In (Discrete Math. 17 (1977)181) Rivest introduced the search complexity of binary trees and proved that among all binary trees with a fixed search complexity the smallest ones are the so-called Fibonacci trees. This result is extended for q-trees. The structure of the smallest q-trees is again Fibonacci-like but more complicated than in the binary case. In addition an upper bound for the asymptotic growth of these trees is given.