研究論文

2017年

1. Isao Nakazawa and Ken Umeno, "Fundamental study on almost periodic frequency arrangements for super-multi-access radio communication systems", to appear in IEICE Communications Express (ComEx) (2017).
2, Hirofumi Tsuda and Ken Umeno, "Weyl Spreading Sequence Optimizing CDMA", to
appear in IEICE Transactions on Communications (2017).
3. Hirofumi Tsuda and Ken Umeno, "Non-Linear Programming: Maximize SINR for Designing Spreading Sequence", to
appear in IEEE Transactions on Communications (2017).
4. Atsushi Iwasaki and Ken Umeno, "Further Improving Security of Vector Stream Cipher", Nonlinear Theory and Its Applications, IEICE, vol. 8, No.3, (2017) pp. 215-223.
5. Hiroki Okada and Ken Umeno, "Randomness Evaluation with the Discrete Fourier Transform Test Based on Exact Analysis of the Reference Distribution", IEEE Transactions on Information Forensics and Security, vol. 12, No.5, (2017) pp. 1218-1226.
6. Ken-ichi Okubo and Ken Umeno, "New Chaos Indicators for Systems with Extremely Small Lyapunov Exponents", in Chaos, Complexity and Transport,Edited by Xavier Leoncini, Christophe Eloy and Gwenn Boedec, (World Scientific,
2017), pp. 185-203.
7. Takuya Iwata and Ken Umeno,"Preseismic ionospheric anomalies detected before the 2016 Kumamoto earthquake",Journal of Geophysical Research, vol. 122 (2017) pp.3602-3616, doi:10.1002/2017JA023921
8. Atsushi Iwasaki and Ken Umeno,"One-stroke polynomials over a ring of modulo 2^w", JSIAM Letters, Vol. 9 (2017) pp.5-9, doi:10.14495/jsiaml.9.5

2016年

1. Hirofumi Tsuda and Ken Umeno, "Orthogonal basis spreading sequence for optimal CDMA", JSIAM Letters, vol. 8 (2016) pp. 77-80, doi:10.14495/jsiaml.8.77
2. Iwata Takuya and Ken Umeno, "Correlation Analysis for Preseismic Total Electron Content Anomalies around the 2011 Tohoku-Oki Earthquake",Journal of Geophysical Research, ,vol. 121(2016) pp.8969-8984, doi: 10.1002/2016JA023036
3. Ken Umeno and Ken-ichi Okubo,Exact Lyapunov exponents of the generalized Boole transformations, Prog. Theor. Exp. Phys. 021A01, pp.1-10, doi: 10.1093/ptep/ptv195
4. Atsushi Iwasaki and Ken Umeno, "Improving security of Vector Stream Cipher", Nonlinear Theory and Its Applications, IEICE, vol. 7 (2016) pp. 30-37, doi: 10.1587/nolta.7.30
5. Ken Umeno, "Ergodic transformatons on R preserving Cauchy laws",Nonlinear Theory and Its Applications, IEICE, vol. 7 (Invited Paper) (2016) pp.14-20, doi:10.1587/nolta.7.14
6. Shin-itiro Goto, "Contact geometric descriptions of vector fields on dually flat spaces and their applications in electric circuit models and nonequilibrium statistical mechanics",Journal of Mathematical Physics, vol. 57, (2016)102702, doi: http://dx.doi.org/10.1063/1.4964751

2015年

1. H. Okada, K. Umeno, D. Handoko, M. Ihsan, R. Maharani, P. Nursetia, H. Rosyad, Warsito, "Brain Electrical Capacitance Volume Tomography Signals Analysis with Moving Maximum Lyapunov Exponents", Advanced Science, Engineering and Medicine, vol. 7, No. 10 (2015), pp. 897-899.
2. 大久保　健一,　梅野　健,弱カオス系のカオス性判定について, 日本応用数理学会論文誌, vol.25, No.3(2015), pp. 165-190.
3. 岡田 大樹、梅野　健、"新たな非線形時系列解析の手法　—移動最大リアプノフ指数線によるカオス解析ー", レーザー研究, vol. 43, No. 6(2015), pp.
4. Chen-An Yang, Kung Yao, Ken Umeno, and Ezio Biglieri, "Superefficient Monte Carlo Simulations" in Simulation Technologies in Networking and Communications Selecting the Best Tool for Test, Edited by Al-Sakib Khan Pathan, Mohammad Mostafa Monowar Shafiullah Khan (CRC Press, 2015) pp.69-91.

2014年

1. Ryo Takahashi and Ken Umeno, "Performance Evaluation of CDMA Using Chaotic Spreading Sequence with Constant Power in Indoor Power Line Fading Channels", IEICE Transactions on Fundamentals of Electronics, Communications and Computer Sciences,E67-A,No. 7 (2014), pp.1619-1622, doi: 10.1587/transfun.E97.A.1619

2013年

1. K. Umeno and M. H. Kao, "Chaos Theory as the answer to limited spectrum?", ITU News, No.10, 2013.
2. Chen-An Yang, Kung Yao, Ken Umeno and Ezio Biglieri, "Using Deterministic Chaos for Superefficient Monte Carlo Simulations", IEEE Circuits and Systems Magazine, 13 (2013), pp.26 - 35.
3. Shun Ogawa, Spectral and formal stability criteria of spatially inhomogeneous solutions to the Vlasov equation for the Hamiltonian mean-field model, Phys. Rev. E 87, 062107 (2013), arXiv:1301.1130
4. 梅野健　情報の統計力学　数理科学 No. 600 (2013), pp. 35-41
5. K. Umeno and A.-H. Sato, "Chaotic Method for Generating q-Gaussian Random Variables", IEEE Transactions on Information Theory, 59 (2013), pp.3199 - 3209.

修士論文

＊修士課程学生は修了のために、修士論文を作成します。
ここ数年の修士論文の題目および摘要は、
以下のリストの氏名をクリックしてご覧下さい。

特別研究報告書

＊４回生は卒業のために、特別研究報告書を作成します。
ここ数年間の特別研究報告書の題目および摘要は、
以下のリストの氏名をクリックしてご覧下さい。

有井 伴樹

Study on CDMA Systems with Primitive Root Codes (原始根符号を用いたCDMAシステムの研究)

In this thesis, we study on CDMA systems with primitive root codes. First, we explain primitive roots modulo a prime number $p$, and introduce the concept of safe prime numbers. Second, using the properties of primitive roots, we construct synchronous CDMA systems with primitive root codes which are complex number sequences. In this paper, we construct the system with the following situations: (1) The same primitive root $q_1$ is assigned to every user. (2) Different primitive roots are assigned to each user. (3) Some users are assigned by the same primitive root $q_1$, and other users are assigned by another primitive root $q_2$. Third, we analyse theoretically Bit Error Rate (BER) of the systems and effects of a safe prime number. Fourth, we simulate numerically BER of the systems. Finally, we compare the theoretical and numerical analysis of BER and evaluate the performance of the systems.

ページトップへ

入江 哲史

Study on the high-precision Monte-Carlo computation using random numbers with nonuniform density (非一様乱数を使用する高精度モンテカルロ計算についての研究)

In this paper, a generalized method and several practical methods of the high-precision Monte-Carlo computation using random numbers with nonuniform density which decrease are proposed. Also the results of comparison between the proposed methods and previous methods, which use random or quasi-random number sequences with uniform density, are shown for applications to several multidimensional integral computation. As a result, the superiority of this methods is demonstrated.

ページトップへ

久世 友博

Robustness of interdependent networks with degree-correlated inter-connections (次数相関を持って複数のネットワークが相互につながっている場合のロバスト性)

Modern systems are constructed with multiple networks that are connected to each other. For example, electrical systems are constructed with the power grids and their communication support systems. In such a interdependent network, failure of nodes in one constituent network leads nodes in the other network to fail. This happens recursively and leads to a cascade of failures. It is known that the interdependent networks with random inter-connections have weaker robustness, tolerance to failures, than the individual networks. However, if the interdependent networks have degree correlations between the networks constructing them, the robustness of the interdependent networks may be changed. Since actual interdependent networks have some correlations, we investigate the effects of the correlations on the networks.
A group that nodes connected with several links is called a cluster and if the number of nodes in the cluster is large, the network is robust. We perform simulations for various ratios of the initial failure of nodes and evaluate the cluster sizes after the cascade of failures. We show that when a interdependent network has a positive degree correlation between two networks that construct it, it has the stronger robustness than that for the networks with no degree correlation. Moreover, as the result of the numerical simulation this system shows a percolation phase transition and the threshold is approximately a linear function of the correlation coefficient. Then, we show not only the numerical simulation results but theoretical ones for the robustness of the interdependent networks. The theory can be applied to the interdependent networks with any degree distributions and any inter-correlations. The theoretical results correspond to the numerical simulation results mainly at any case.

ページトップへ

ページトップへ

ページトップへ

家治川 博

Linear Regression Analysis of Foreign Exchanges with a Method of Segmenting Time Series Based on the Likelihood-Ratio Test

There has been a lot of researches to analyze time series data. On the other hand, a variation of information which we can obtain is getting wider and wider these days because of the development of information technology. Therefore, it is important to deal with various information for a purpose to analyze financial time series data more accurately.
In this paper, we conduct a linear regression analysis and set Google search queries as explanatory variables of the model in order to analyze financial time series data with taking various information into consideration. Additionally, we also set the interbank exchange frequencies and the volatility calculated by GARCH(1,1) model as explanatory variables.
Since financial time series are not generally modeled as a stationary process, but modeled as a non-stationary process, we assume that the non-stationary time series consists of several stationary segments with different properties. In order to discriminate the joint of these distributions, we conduct a likelihood- ratio test. According to this test, the point which maximizes the likelihood- ratio between the null model (homogeneous disturbance distribution) and the alternative model (a mixture of two different normal distributions) is regarded as the most possible combination point. Since likelihood usually includes sample errors, we evaluate the significance level of it by using the bootstrap method. We employ the bootstrap distribution as a discriminant measure to divide the time series into two segments at an adequate point, recursively.
We apply this method to the data of foreign exchange market and make an analysis of them in terms of segmenting points and regression coefficients. We compare our proposed method with the ordinary linear regression method visually and numerically. We conclude that Google search term which segments the time series at a point different from other search term has an additional relation with the explained variable. Our proposed method can detect the data point at which the tendency of the time series changes and can analyze the time series better than the ordinary method.

ページトップへ

守田 悠三

Analysis of Foreign Exchange Rates Based on Parametric Risk Assessment Procedures withq-Gaussian and Pearson type IV Distributions

Recently, it has been much easier for individuals to buy and sell foreign currencies in the market and it is becoming more important to understand the foreign currency risks to hold foreign currencies safely. Exchange rates sometimes fluctuate unpredictably and it could cause loss of the deposit. Especially, it is well-known that the volatilities of the price fluctuate de- pending on the time period, which seems to cause fat-tailedness observed in actual data. Therefore, it is important to regard its fat-tailedness when we estimate the risk from historical data.
In this paper, we introduce parametric risk assessment procedures in the foreign exchange market. In order to consider the ruin probability when we have some deposit, we assume two types of distributions, the q-Gaussian and the Pearson type IV distributions that log-returns in foreign exchange rates obey.
We perform parameter estimation of q-Gaussian distribution for 30 currency pairs with a maximum likelihood method. The parameter q is estimated in the range of from 1.3 to 1.7, and we confirm that the empirical distributions of the market data have fat-tails. To check whether the estimated parameters are statistically significant, we calculate p-values of two types of statistical test, Kolmogorov-Smirnov (KS) test and Anderson-Darling (AD) test. We reveal that all p-values in KS test are larger than 0.1, although p-values of 9 currency pairs are less than 0.1. This means that the log-returns obey the q-Gaussian as a whole, but do not when we focus on tails of the distributions.
We also perform parameter estimation of the Pearson type IV distribution, which has skewness, for 30 currency pairs. In this case, the average p-values in both KS and AD test are better than those of the q-Gaussian and the p-values in AD test for only 5 pairs are less than 0.1. This means that the Pearson type IV is better-fitted to the market data. Therefore, we reveal that the model with skewness is preferred for risk estimation. We calculate 1% Value-at-Risk (VaR) in both cases. The difference of the VaR between the q-Gaussian and the Pearson type IV is about 10%. This indicates that VaR with the q-Gaussian could cause underestimation of the risk.

ページトップへ

吉村 玄太

Capacity-Approaching LDPC Codes Constructed from Extended Protographs

Low-density parity-check (LDPC) codes are in a class of the most powerful error correcting codes available today. In this paper, to improve the decoding performance we introduce a new class of LDPC codes constructed from a template called an extended protograph, which belongs to the superset of the popular protographs. We exploit the extrinsic information transfer (EXIT) chart and asymptotic ensemble weight enumerator techniques to predict the capacity-approaching performance of extended protograph code ensembles. EXIT charts compute the iterative decoding threshold of a code ensemble, which contributes to the waterfall region performance for the decoding error rate. Asymptotic ensemble weight enumerators estimate whether or not the minimum distance of code ensemble increases linearly with code length, which affects the error floor performance for the decoding error rate. Taking account of these conflicting performance indices, we apply the simulated annealing (SA) algorithm to find extended protograph based codes which have low iterative decoding threshold and the linear minimum distance growth property. Finally, the performance of optimized extended protograph based codes over the binary-input additive white Gaussian noise (BIAWGN) channel are compared with that of existing protograph based codes and the quasi-cyclic (QC) code adopted by IEEE 802.16e (WiMAX) standard.

ページトップへ

池田 和雄

Opinion propagation using partisan voter　models on several networks

The partisan voter model is the one of　models which treat the social opinion dynamics. We use them in order to treat opinion dynamics on several networks. This model is a modified version of the voter model which describes the evolution to consensus in the networks of nodes, voters, that possess respective opinions within a discrete set. In the partisan voter model, each node has also an innate and fixed liking for one opinion. This liking determines the probability that each node changes its opinion. We calculate the convergence time toward consensus of opinions in several networks by using the partisan voter model. The convergence time in the scale-free networks, the BA model, is smaller than that in the complete graph. With the modified BA model we investigate the dependence on the exponent of the degree distribution, that is, on the number and size of hubs on the network. If the exponent is small, the number of hubs is small and the size of hubs is large. From our simulation, for the small exponent, the convergence time is shorter than that for the large exponent. The correlation of degrees is controlled by rewires of the links in the network initially created by the BA model. For assortative networks, we confirm with our simulation that convergence time is, in general, longer than that for uncorrelated networks. On the other hand, for disassortative networks, it is shorter than that for uncorrelated one in general.

ページトップへ

野田 実

Japanese hotel statistics in terms of regional room capacities

In this paper, in order to understand the dependence of stay capacity on regionality, we propose a method to determine districts depending on the number of rooms and to classify its districts. We empirically analyze the geographical positions and the number of rooms about 2,881 Japanese hotels which have 582,898 rooms in total. Firstly, we conduct a clustering analysis of the regional stay capacity by centroid method. Secondly, we introduce the maximum entropy principle in order to divide regional areas into some levels. It may be concluded that the rank size distribution for the number of rooms in the cluster is fitted with a power-law function and that the scaling exponent is dependent on the number of clusters.

ページトップへ

森岡 篤

Optimization of routing strategies for data transfer in peer-to-peer networks

Recently, peer-to-peer file-sharing systems have become familiar and the information traffic in the networks is increasing. Therefore it causes various traffic problems in peer-to-peer networks. In this paper, we model some features of the peer-to-peer networks, and investigate the traffic problems. Peer-to-peer networks have two notable characters. One is that each peer frequently searches for a file and download it from a peer who has the requested file. To decide whether a peer has the requested file or not in modeling of the search and download process, we introduce a file-parameter Pj , which expresses the normalized amount of files stored in peer j. It is assumed that if Pj is large, peer j has many files and can meet other peers' requests with high probability. The other character is that peers leave and join into the network repeatedly. Many researchers address traffic problems of data transfer in computer communication networks. To our knowledge, however, no reports focus on those in peer-to-peer networks whose topology changes with time. For routing paths of data transfer, generally, the shortest paths are used in usual computer networks. In this paper, we introduce a new optimal routing strategy which uses weights of peers to avoid traffic congestion. We find that the new routing strategy is superior to the shortest path strategy in terms of data traveling time when many peers join in the data transfer.

ページトップへ

永田 啓悟

Analysis of image encryption schemes using chaotic maps

Along with the development of the information and telecommunications networks, various kinds of encryption schemes have been developed by many researchers. Some of them pay attention to chaotic encryption. They utilize chaotic properties such as initial sensitivity, synchronization, and so on, for enhancing security and effect of encryption. In this paper, we consider image encryption schemes using baker's map, which is one of chaotic maps. Image encryption scheme is divided into two phases, "permutation" and "diffusion". Permutation is encryption on positions of pixels, and diffusion is that on gray values. Both of them are necessary for secure encryption. In some previous studies, baker's map is used for permutation. With repeated application of this map, the pixels are permutated, and the original image comes to look featureless and random like noise. We can decrypt it easily if and only if the widths of the rectangles, called keys, are known because baker's map is reversible. First, we analyze the properties of baker's map and improve it to surmount the weakness. Secondly, we apply baker's map to diffusion. Finally, we evaluate the security of the proposed encryption scheme, and prove the usefulness of it.

ページトップへ

中本 武志

Estimating the tail index of distributions: Case study on the foreign exchange market

We study an unconditional distribution derived from Alfarano-Lux model which has the two parameters: the herding propensity and the autonomous switching tendency. One of the parameters, the herding propensity characterizes the tail shape of its distribution, but does not accord with the well-known Hill's estimator in case of foreign exchange market data. In this paper, we transform the probability density function of Alfarano-Lux model to the expanded form and obtain the analytical form of its cumulative distribution function. Additionally, we explain the reason why the difference between the estimates of Alfarano-Lux model and Hill's estimator is generated, and conduct Kolmogorov-Smirnov test to measure the goodness-of-fit. As an application of Alfarano-Lux model to the foreign exchange market data, we measure the fluctuation risk for some currency pairs.

ページトップへ

西岡 謙太

Log returns of stock values and q-Gaussian distributions: Application to the risk-assessment

We assess a risk of financial time series with the distribution estimated from them. Because it is not easy to infer its tail shape due to a lack of data in a practical manner.(1) We introduce Value at Risk(VaR) to a risk measure and compare it with variance under the q-Gaussian assumption.(2) We examine performance of the maximum likelihood estimator with the q-Gaussian log-likelihood function.(3) By using the distribution estimates, we verify errors of the VaR to estimated one. Finally we conduct an empirical analysis on log-returns of a stock traded in the Tokyo Stock Exchange.

ページトップへ

大谷 卓也

Sequential associative memories on complex networks

The models of associative memories in neural networks, such as Hopfield model, are quite tractable mathematically, and applied to many areas. One of them is a model of sequential associative memories. Most of the models, however, including sequential associative memories, are usually applied to all-to-all networks. Since the studies of neural networks originate in brain science, it is natural to apply them to more realistic networks. Many researchers have been interested in complex networks, which is said to describe properties of real networks in society, computers, and neurons, that is, the small-world and scale-free properties. For example, random graph, the Watts-Strogatz model, the Barabási-Albert model, are the most famous complex networks which have some interesting properties. These days, the studies of complex networks are developed drastically. In this paper, the model of sequential associative memories is applied to these complex networks. The main subject is how the network topology affects the temperature dependence of the retrieval performance. Computer simulation reveals that the temperature dependence varies with the network topology, and it also comes out that the local performance on each node is increasing monotonically with its degree. Modifying the existing theory to be suitable for complex networks, we have been obtained a new approximated equation for overlap using a mean-field approximation and the the central limit theorem. This approximated equation can describe the temperature dependence of performance with any networks whose degree distributions are known. Numerical solution of the equation on the Barabási-Albert model is compatible with simulation results. On the other networks, the numerical solutions do not agree well with the results. Loops of networks are considered as one of the possible causes for the inagreement because the loops could disturb the independence of node states which is essential to apply the central limit theorem. As the Barabási-Albert model is known to have a very small number of loops than the other complex networks, the above compatibility between simulation and theoretical results is consistent with the loop effects. In order to determine whether the loops are the cause of the error, in the last of this paper, we present new approximated simultaneous equations which evaluate the effects of loops. As far as the numerical solution on random graph, it is likely that the loops have little effect on the performance of sequential associative memories on complex networks.

ページトップへ

砂川 敦志

Active random walkers: a simple model fro catalyst dynamics

Some problems related to active random walkers are studied with use of a model for a chemical system which consists of catalysts(C), products(P) and reactants(R). For numerical experiments we use a Monte Carlo method to follow motion of C- and P-particles. For analytical considerations a Langevin model for particle dynamics is introduced, which is converted to a coupled diffusion equation, whose linear stability is studied in relation to numerical experiments. We assume either attractive or repulsive interactions between C-particles and P-particles and the catalysts C are able to change environment locally by producing P-particles from R-particles. P-particles and C-particles can diffuse in the system and P-particles are characterized by a decay constant, which prevents the number of P-particles from increasing indefinitely. Our main results include

◦1) In case of attractive interaction between C- and P-particles, the density fields of C-particles and P-particles tend to a stationary bound state, which consists of a bump of C-particles and a bump of P-particles occupying the same small region (in a long time limit).
◦2) In case of repulsive interaction between C- and P-particles, we observe (irregular) density waves for both C-particles and P-particles, in which a region of high P-particles density corresponds to a region of low C-particles density, and vice versa.
◦3) When the amount of R-particles, available for the chemical reaction, is finite we observe that a ring, inside of which R-particles are consumed completely, expands outward.
◦4) Linear stability analyses give us information on the characteristic length of the stationary state and this turned out to be consistent with our findings (1) and (2) above.

ページトップへ

君塚 誠

Properties of coupled double well systems with delay and noise

The interplay between noise and delay in physical as well as biological systems gives rise to a lot of interesting phenomena and is gathering interest of many researchers. An especially hot area in this connection is the semiconductor laser emitter (VCSEL in short), which is modeled by a Brownian particle in a double-well potential under delayed force. If there were no delayed force, the system would show a simple barrier-crossing or hopping, which means for VCSEL the transition from the vertically polarized state of light to the horizontally polarized one and vice versa. In the presence of delay, which represents the delayed feedback for VCSEL, this hopping is modified due to strong correlation between (x(t)) and (x(t-τ)) , with (x(t)) and (τ) denoting position of the Brownian particle at time t and the delay time, respectively. In this work we consider a model which consists of (N)-Brownian particles put in a double-well potential and investigate some theoretical problems such as the positional stationary state distribution function (Pss) and the time correlation function(TCF). This model may be considered as a model for VCSEL in which (N) laser emitters are connected in cascade. At the moment we have no experimental results available for comparison with our theory. Thus we performed ourselves computer experiments to obtain some physical quantities of interest. For theoretical analysis of Brownian dynamics in a double well potential, we employ a two-state approximation, which replaces the continuous position x(t) with s(t), which takes only two values +1 ( if (x(t)>0) ) and -1 ( if (x(t)<0) ). This greatly simplifies our problem and at the same time gives much insight into our problem. Some results obtained in this work include, 1) the positional stationary state distribution function (Pss) for the Brownian particle model is calculated (computer) experimentally and its (τ) and feedback strength ε dependences are clarified. These dependences are theoretically studied for the case (N=2) based on the two-state approximation, which results in (Pss) well correlated with our experiments. It is remarked that this (τ) dependence is absent for the case (N=1, 2) the time correlation function is also calculated numerically by solving the Langevin equation and theoretically based on the two-state approximation.

ページトップへ

近藤 健夫

Signal response in scale-free network of bistable units

Information processing in biological systems, such as brains or cell membranes, has several features in common: (i) it is an analog information processing; (ii) scale-free networks play an important role; and (iii) each information processing unit possesses strong nonlinearity. We studied the efficiency of information processing in a scale-free network, consisting of many interacting nonlinear units with double-well potential, both theoretically and (computer) experimentally. We introduced the gain (G) to quantify the efficiency, which is defined as the ratio of the output strength aL of the maximum response unit L to the periodic input signal strength (A), i.e. (G≡aL/A) . In the previous work, Acebron et al. calculated (G) as a function of coupling (or interaction) strength (λ) i.e. &math (G(λ)) , by simulations and they found that (G(λ)) showed plateau behavior in some range of (λ) They tried to understand their numerical results with emphasis put on a hub, which is a typical feature of the scale-free network. That is, they analyzed dynamics based on a simplified model (starlike network). However, the analysis contains some problems: (i) it is limited to a very small (λ) region; and (ii) the overall role of the scale-free network for information processing is not touched upon. In this study, we develop a formalism, in which the gain (G) is studied based on a one-body problem, which turns out to be a good approximation to replace the dynamics of the complex network. We analyzed this model, and obtained four theoretical predictions: (i) the value (Gplateau) of the plateau height of (G(λ)) (ii) the initial transient behavior of (G(λ)) in a small (λ) region, (Gtrans(λ)) with finiteness of the system taken into account, (iii) the asymptotic (G) value, (Gsync) after full synchronization is achieved for large (λ) and (iv) the critical value sync) , beyond which all units are fully synchronized and the gain (G); becomes (Gsync) We compared these theoretical predictions with computer experiments and confirmed that our theory reproduces experimental results at least semi-quantitatively. From this we may say that we could understand (G(λ)) in the whole coupling strength range by revealing the role of the complex network in information processing.

ページトップへ

西村 麻衣子

Scaling analysis on quotation activities in the foreign exchange market: Empirical investigation and stochastic modeling

We investigate quotation activities in the foreign exchange market both empirically and theoretically. We found the scaling relationship between means of the number of quotations during window lengths and their standard deviations. We confirm that the scaling exponent temporally changes from 0.8 to 0.9 depending on observation days and that it tends to be unity as the time window length increases. We also estimate a cross-correlation coefficients and relative frequencies of the quotation activities. As a result, it is found that the scaling exponent and the cross-correlation coefficients show a significant correspondence relationship. Therefore, it is concluded that scaling analysis would be one of adequate ways to grasp the total market state. Besides, we propose a stochastic model to understand the market participants' activities and conduct a theoretical analysis for the proposed model with parameter fitting by empirical data. Comparing empirical results with theoretical ones, we examine the adequency of the model. Consequently, we found that the fluctuation of the probability with which market participants decide their attitudes and the weight of each currency pair have an important role to determine the value of a scaling exponent under the condition that the probability is homogeneous for every market participant.

ページトップへ

ページトップへ

岩間 心平

Effects of time delay on nonlinear stochastic systems

Considerable attention is paid to some stochastic systems, whose dynamics is determined by the present state x(t) and the state x(t-T) in the past with r(> 0) denoting the delay time. Usually this delay is ascribed to the finite speed of information transmission and it is rather natural that effects of delay are intensively studied mainly for biological systems, e.g. as models to describe postural sway, visual feedback, and brain activity, to mention a few. Effects of delay are also studied for chemical, physical, and engineering systems. So long as we know, there have been no systematic studies on nonlinear systems where effects of delay play essential roles to determine system properties. Also from the viewpoint of methodology we have no reliable method, which can be applied for large delay time and strong nonlinearity. From these points stochastic systems with delay are offering many problems, interesting both from mathematical and physical viewpoints. Recent progress in understanding the delayed system may be partly due to the advent of the Fokker-Planck equation for the (one-body) distribution function p(x, t). However this equation is not a closed one for the distribution function p(x, t) in the sense that it contains additional 'collision' term, expressed in terms of two-body conditional distribution function p(y, t-T | x, t) for x(t-T) = y when x(t) = x. Developing some approximation schemes to make it a closed one, we give detailed discussions on the range of validity of each approximation scheme for the delay Fokker-Planck equation. Our strategy is as follows: First we note that the delay Fokker-Planck equation has two important parameters, the delay time T and the strength E of a delay term. We propose three approximation schemes, namely, which are supposed to be applicable for (i) small T, (ii) small E and (iii) large T regions. By applying these schemes to a double-well potential system and comparing theoretical predictions with (numerical) experiments, we estimate how our proposed schemes work to this problem.

ページトップへ

川本 大樹

Efficient packet routing strategy in complex networks

We investigate new packet routing strategies which mitigate traffic congestion on complex networks. Instead of using shortest paths, we propose efficient paths which avoid hubs on scale-free networks with a weight of each node. Firstly, we compare the routing strategy using degree-based weights with that using betweenness-based weights on two types of scale-free networks. The strategy using degree-based weights is more efficient than that using betweenness-based weights on scale-free networks generated by a preferential attachment. On the other hand, the ascendancy of the strategy using degree-based weights over that using betweenness-based weights is reversed on scale-free networks composed by taking into account the distance between nodes. Next, we consider the heuristic algorithm which improves step by step routing properties on congestion by using the information of betweenness of each node in every step. We propose new heuristic algorithm which balances traffic on networks by achieving minimization of the maximum betweenness in the much smaller number of iteration steps.

ページトップへ

酒井 洋

Discrimination and Cluster methods for multi-dimensional time series of the Foreign Exchange Market based on Spectral Distances

In this thesis, the similarities between multi-dimensional time series extracted from high frequency financial data of the foreign exchange market are measured and the hierarchical clustering is performed based on them. We introduce two methods to calculate their similarities based on spectral distances. One is the sum of the Kullback-Leibler divergence between two normalized power spectra of time series of each elements for different observations, and the other is that between two largest eigenvalues of cross spectral matrix computed from multiple time series. In order to verify adequateness of these methods, we introduce the agent-based model of the foreign exchange market in which N market participants exchange M currency pairs, and perform numerical simulation and the hierarchical clustering with pseudo price movements obtained from it. As the results, it is found that the tendency that time series are unified sequentially as the difference of parameters became large, and that results obtained by means of two methods are almost similar. Finally applying this procedure to actual tick data we confirm that this can extract meaningful information from a large amount of data. The day when the rate movements of several currency pairs are greatly abnormal has a tendency that that on the day belongs to a separated cluster from those on other days. Therefore, it is concluded that this procedure is applicable to automatic extraction of meaningful information about markets from enomarous amounts of financial data.

ページトップへ

多羅尾 光記

Information Flow and Causality in Coupled Systems

In studies of a complex system, one of the major concerns is the detection and quantification of "causal interdependencies" among many dynamical subsystem, which constitute the complex system. For systems studied via physics or chemistry, this interdependency may be called "coupling" or "interaction". The purpose of this thesis is to study both physical and man-made systems on equal footing, using the concept of entropy transfer (or information flow) between subsystems. From physics we know that heat flows from a high to a low temperature region and sound wave propagates through a system with a constant velocity. These may be regarded as a kind of information transfer and there have been many works to understand these phenomena from the viewpoint of information flow. Based on the recent advance in methodology to quantify entropy transfer, as developed by Schreiber, we consider three systems, (1) a linearly coupled Langevin system, (2) FitzHugh?-Nagumo neural networks and (3) foreign exchange market. For each system we quantitatively analyze the system dynamics, especially the interrelation among subsystems, based on entropy transfer. The characteristics of these systems is that they have many units that are coupled with each other, unidirectionally or bidirectionally. In the linear system composed of two Brownian particles, we derive the transfer entropy rate in a theoretical way by taking the limit dt→0 where dt represents the sampling time. And we compare it with the entropy production rate. In the non-linear system that consists of many neurons, the information flow in three types of network, (i) linear array type, (ii) triangular type and (iii) a star-like one is analyzed. For a foreign exchange market, we treat the thirteen currency pairs. By calculating the entropy transfer among currency pairs, we observed some rules, which are rather vague but suggestive and seem to be consistent with our intuition.

ページトップへ

小崎 元也

Application of the Beck model to stock markets: Value-at-Risk and portfolio risk assessment

We have applied the Beck model, developed for mechanical systems that exhibit scaling properties, to stock markets. Our study revealed that the Beck model elucidates properties of stock market returns and is applicable to practical use such as the Value-at-Risk estimation and the portfolio analysis. We have performed empirical analysis with daily/intraday data of the S&P 500 index returns and found that the volatility fluctuation of real markets iw well consistent with the assumptions of the Beck model: the volatility fluctuates in much larger time scale than return itself and the inverse of variance, or "inverse temperature," beta obeys Gamma-distribution. As predicted by the Beck model. the method of Value-at-Risk (VaR), one of the most significant indicator in risk management, is studied for q-Gaussian distribution. Our proposed method enables the VaR estimation in consideration of tail risk, which is underestimated by the variance-covariance method. A framework of portfolio risk assessment under the existence of tail risk is considered. We have proposed a multi-asset model with a single volatility fluctuation shared by all assets, named the single beta model, and empirically examined the agreement between the model an and an imaginary portfolio with Dow Jones indices. It turns out that the single beta model gives good approximation to portfolios composed of the assets with non-Gaussian and correlated returns.

ページトップへ

新谷 幸平

Empirical analysis and numerical simulation of foreign exchange market dynamics

In this thesis, the dynamics of the foreign currency markets is analyzed with the tick frequency, which is the trace of market participatnts' activities. With the double-threshold agent model, the relation between the tick frequency time series of two agent groups is discussed. Empirical analysis with actual tick data indicates that the similarity of the tick frequency time series is getting similar as time passes and the similarity network among currency pairs changes by the areas where the markets are active. It means that in the present financial markets, market participants become more homogeneous, and the risks in the financial markets must be calibrated depending on areas, respectively. Alternative analysis indicates that the similarity between the best ask price and the similarity between the tick frequency time series influence each other constantly and weakly, and that the co-movement of them quickly and largely changes. Finally, a way of constructing the real-time market monitoring system is discussed.

ページトップへ

ページトップへ

徐 普永

Self-tuning of activation energy in a two-state system

Many researchers have observed, in various biological and physical systems, that noise can enhance responses(i.e. output signals) of nonlinear systems to a weak periodic driving force(i.e. input signals) in a positive way, thus transfering more information than a case without noise. The phenomenon is closely related to stochastic resonance(SR), which states that there exists optimal noise strength T* where best information transfer is achieved. However it seems that not much attention has been paid to a situation where noise is rather weak and information transfer is very limited. In order to solve this problem, we attempt to apply an idea of self-tuning (ST), first proposed to explain high sensitivity of auditory systems of animals. That is, we consider a simple adaptation process of an activation energy(to be regarded as a threshold), which turns out to work well for a two-state system (TSS) in a weak noise region. Improvement of information processing ability of the TSS in a strong noise region was made possible by analytically studying the adaptation equation. We apply our ST method to a double-well potential system (DWPS), showing that the ST method works well also for the DWPS. Some quantities, not directly obtainable theoretically such as the first-passage-time distribution function, are calculated based on Monte-Carlo simulations. Finally we consider ST and SR from the viewpoint of 'energy transfer' from input signals to reservoirs. It is shown from numerical experiments that this energy transfer shows similar behavior as the information transfer.

ページトップへ

長沼 佑樹

Packet routing strategy using neural networks on complex networks

We investigate routing strategies on complex networks. Firstly, using neural networks, we introduce a routing strategy where path lengths and queue lenghts are taken into account within a framework of statistical physics. The performance in this strategy becomes more efficient from improvement of the distance term. At the same time, we analyze how the properties of networks influence the performance of this strategy. Secondly, we propose a routing strategy where connection weights in neural networks are adjusted by local information. We also confirm how the distance term and the properties of networks influence the performance of this adjustive strategy.

ページトップへ

松田 雄馬

Synchronization of a randomly coupled map model for neural networks

Neurons are known to oscillate synchronously with each other to achieve various functions. Studies on randomly coupled Hodgkin-Huxley neuron models have show that neurons synchronize their activity even in noisy environments. Synchronization is reported to occur when the size of networks is large enough even for sparsely connected neurons. Because fluctuation is neglected in large size networks, synchronization is controlled only by average states of neurons. Although it is widely observed in nature, especially in nervous systems, most of the mechanism is still unknown. In this paper, the condition for synchronization on neural networks is investigated with numerical simulations of a one-dimensional map neuron model. Then, the synchronization in noisy fields is theoretically understood throughout the distribution of neuron states formulated by inductive statistics and Gaussian mixture models.

ページトップへ

ページトップへ

岡田 大樹

離散フーリエ変換検定法とその改良についての研究

また，DFT検定の問題点の考察を通して，DFT検定とは検定対象が多少異なるが，周期的特徴を検出する独自の検定法を考案した．この検定法について，少数のフーリエ係数のみが大きくなるような乱数の検定においては離散フーリエ変換検定法より周期的特徴の検出能力が優れることが確認できた．
それの他に，離散フーリエ変換をより一般化した変換について考察し，それをもとにして作る新たな離散変換を離散フーリエ変換と置き換えるという変更を離散フーリエ変換検定法に加えることで，チェビシェフ写像などを生成法の基にするカオスな擬似乱数系列等を検定可能とする，擬似乱数系列の持つカオス性（＝非乱数性）を棄却検定する検定法を構成した．

ページトップへ

長谷川 史晃

嗅覚系モデルにみられる神経集団システムの研究

また，力学系が振動するためのパラメータの条件を求め，これが正しいことを示すために，力学系が振動する時のリミットサイクルの内部の面積を解析的に求め，求めたパラメータの条件のもとで，面積が存在することを導き出した．

ページトップへ

亀山 慎吾

原子時系発生システムの安定化と高精度化に関する研究

そして、最後に精度の観点から現在米国のNISTで用いられている標準時の時系アルゴリズムにおいて、指数フィルタと呼ばれるものについて議論した。また、これにLSFを導入することで負の相関をもたせ、より高精度なシステムになる可能性があることを提案した。この提案は今後の研究の課題の1つとしてこれから調査していかなければならないことである。

ページトップへ

猶原 僚也

パワー一定カオス拡散符号を用いたCDMA通信システムの性能評価：ルベーグスペクトラムフィルタの適用

また、自己相関特性を制御し、他のユーザからの干渉雑音を抑制するルベーグスペクトラムフィルタ（Lebesgue-Spectrum-Filter）という線形フィルタを導入することによって、ビット誤り率やSNRが低下したり、SIR（信号対干渉比）が上昇することをシミュレーションによって示す。具体的には、従来のシステムより同じビット誤り率でユーザが15\%増加できてそれが非同期CDMAシステムでは解析的に最適である。

ページトップへ

ページトップへ

ページトップへ

ページトップへ

ページトップへ

ページトップへ

ページトップへ

ページトップへ

吉村 玄太

適応型ネットワークにおける潜伏期を持つ感染症のダイナミクス

ネットワークのダイナミクス（トポロジーつまり接続状態の変化）とネットワーク上のダイナミクス（ノー ド状態の変化）が相互作用するようなネットワークを適応型ネットワークと呼ぶ。本論文では潜伏期を持つ感染症が適応型ネットワークにおいてどのような定常状態に達するかを明らかにする。まず適応型ネットワークにおける潜伏期を持つ感染症のダイナミクスを表現するため、Susceptible-Infected-Susceptibleモデル（SISモデル）を発展させたSusceptible-Latent-Infected-Susceptibleモデル（SLISモデル）を導入する。次にネットワークのノード状態及びトポロジーを表す平均場量を用いて近似し、導入した SLISモデルを支配する連立常微分方程式を導出する。この方程式を解くことで、SLISモデルの定常状態における感染ノード密度を理論的に求める。また SLISモデルを計算機上で実装し、シミュレーションを行うことで定常状態における感染ノード密度を数値的に求める。感染確率pに対して、無病状態から流行状態への転移が起きる侵入閾値 pinv 、及び流行状態から無病状態への転移が起きる持続閾値 pper が定義できる。従来のSISモデルでは、リンクの張り替えにより感染症との接触を避ける防衛の効果によって、多くの場合に pper < pinv となり、2つの閾値の間では履歴に応じて無病状態と流行状態の2つの安定な状態をとる、いわゆる双安定性を有する。これに対して本研究で導入したSLISモデルでは、潜伏期が長くなるほど防衛の効果が薄れ、 pinv, pper が減少していくこと、双安定性が失われることを、得られた解析及びシミュレーションの結果から示す。た pinv, pper の、潜伏期の長さを特徴付けるxに対する依存性を示す。さらに潜伏期における自覚症状の有無により定常状態、特に pper や流行状態の感染ノード密度に変化が生じることを示す。

ページトップへ

ページトップへ

ページトップへ

ページトップへ

ページトップへ

ページトップへ

ページトップへ

ページトップへ

木津 幸子

時間周期外力の働く系における揺らぎ定理

エントロピー生成あるいはエントロピー生成率は、非平衡な過程を特徴付ける最も基本的かつ一般的な物理 量である。1990年代初頭に見出された「揺らぎ定理」の出現によりエントロピー生成は益々その重要性を増しつつある。観測時間τでのエントロピー生成を Δ とすると、「揺らぎ定理」はその確率分布 p(Δ) p(Δ)/p(-Δ)=exp(Δ) を満足することを主張する。小さな体系や短い観測時間では揺らぎが大きく、熱力学第2法則に反してエントロピー生成が負になることもありうるが、揺らぎ定理はその確率を定量的に議論することを可能にした。
ただし、確率分布 p(Δ) から直ちに熱力学第2法則 ⟨Δ⟩= ∫ dΔp(Δ)Δ > 0 が導かれるので、揺らぎ定理は熱力学第2法則を、揺らぎを含めた形に一般化（拡張）したものであるといえ る。揺らぎ定理研究の初期には、主として決定論的な力学系が扱われたが、次第に確率系も研究の対象になり、現在では非決定論的な、すなわちストキヤス ティックな力学系も多くの研究者の関心を集めている。また、「揺らぎ定理」の多くのものでは観測時間Tが十分大きいときに成立することを主張するが、現在、この制限を緩める方向でいくつかの研究がなされている。本研究においては、ポテンシャルV(x)の中を運動する（過減衰）ブラウン粒子に周期外の時間周期外力が作用する系をランジェバン方程式で記述し、 i）観測時間丁において、上に述べた揺らぎ定理が厳密に成立することを示し、 ii）これを以下のポテンシャルをもつ3つの系に適用した。
a）調和振動子系 V(x)=x2/2
b）2次一4次系 V(x)=(x+1)2(x-1)2
c）区分線形系 V(x)=x(2f) (0 ≤ x < f) , V(x)=(1-x)/2(1-f) (f ≤ x < 1)

ページトップへ

ページトップへ

武田 隆之

社会ネットワークにおけるランダム性を含んだ意見形成モデルについて

あるネットワークを構成する者の間で周囲の意見を取り入れ、周囲の意見の平均値と自己の意見の積に雑音を加えたものを元に自己の意見の更新がなされるモデルを考える。一定時間後にそのネットワーク全体の意見がどのような結論状態 に達するかはネットワークの構造と雑音の強度に依存する。いま初期状態をランダムに決定し、雑音の最大値を指定した上で、十分なステップ数だけ意見の更新を行った状態を結論状態とみなす。リンク数とリンク先が固定された規則的なネットワーク、リンク数は固定でありながら単位時間ごとにリンク先がランダムに決定されるネットワーク、リンク数の分布がべき乗になりリンク先が固定されたネットワークに対して結論状態を求める。平均場近似を用いて理論値を導出し、各ネットワークの結論状態と比較する。また、一定の割合で自己の意見を固持するようなモデルを導入し、固持する割合による結論状態の差異を調べ、ネットーク内で固持する割合が一定であるモデルとそうでないモデルの間での比較を行い、ネットワークの構造ごとの結論状態の差異について考察する。

ページトップへ

ページトップへ

平松 将

多変量自己回帰モデルを用いたスペクトル距離の分析

データに基づき体系の背後に横たわるダイナミックスを定量化することによって、 現象の理解を深める目的で、本研究では、時系列の緩和の程度を定量化するスペ クトラル・エントロピーに基づいたスペクトル距離の評価方法を提案する。そして、 多変量自己回帰モデルによって人工的に生成した時系列に対して、提案するスペクトル距離の定量化手法が妥当であることを示す。

ページトップへ

ページトップへ