Program
March 8 (Monday)
Below timetable is Taiwan time. Add 1 hour for Japanese side.
8:50-9:00 Opening
9:00-9:55 Session 1 Chair: Yue-Hong Li (National Cheng Kung University)
9:00-9:15 Kai-Yang Liu (National Yang Ming Chiao Tung University, M2) Lyapunov functions for some epidemic models with delay
9:15-9:30
Kaito Fujihara (Shimane University, M2)
Mathematical study on hidden emotions and body patterns observed in Cuttlefish
9:30-9:45
Tsuyoshi Ishizone (Meiji University, M2)
A method to infer latent dynamics model by combining deep neural networks with ensemble Kalman filter
9:45-9:55 Chun-Sheng Chen (National Central University, B)
The data forecast of COVID-19 by method of machine learning
10:20-11:20 Session 2 Chair: Masakazu Kuze (Hiroshima University) 10:20-10:35 Yusuke Fujita (Hiroshima University, M2)
Why do two-dimensional dunes always merge
10:35-10:50 Yu-Jie Ho (National Yang Ming Chiao Tung University, M1) The Second Laplacian Eigenvalue
10:50-11:05 Ryu Fujiwara (Meiji University, M2)
Discontinuous steady states of the nonlocal Allen-Cahn equation 11:05-11:20 Ya-Chi Chu (National Cheng Kung University, M2)
The convexity for the joint range of two quadratic functions 11:20-11:25 Group Photo
13:00-13:55 Session 3 Chair: Shun Ito (Meiji University) 13:00-13:15 Naoki Kondo (Meiji University, M2)
Mathematical modeling and analysis of filamentous active matter 13:15-13:30 Huan-Chi Chang (National Chung Cheng University, M3)
Using machine learning methods to predict RNA secondary structure
13:30-13:45
Koutaro Tanaka (Shimane University, M1)
Synchronous solutions and their stability in coupled salt-water oscillators
13:45-13:55
Yen-Jia Chen (National Center for Theoretical Sciences, B) Mathematical models of COVID-19 with diffusion effects and their data forecast
14:15-15:10 Session 4 Chair: Chih-Chiang Huang (National Taiwan University) 14:15-14:30 Ryosuke Sakai (Meiji University, M2)
On Pehlivan’s four-scroll chaotic system 14:30-14:45
Jin-Zhi Phoong (National Taiwan University, M2)
On synchronization analysis of complex coupled Kuramoto oscillators
14:45-15:00
Ryota Kubohara (Shimane University, M1)
Periodic solutions for three species Lotka-Volterra competitive equations
15:00-15:10 Jian-Yu Chen (National Central University, B) Combine visual cryptography with RSA algorithm
15:30-16:25 Session 5 Chair: Shun-Chieh Wang (National Taiwan University) 15:30-15:45 Yuto Sakai (Ryukoku University, M1)
Comparison principle for time fractional diffusion equations 15:45-15:55 Yueh-Tzu Hung (National Yang Ming Chiao Tung University, B)
Turing instability for a class of 3-component systems 15:55-16:10
Yuto Kikuchi (Shimane University, M1)
A consideration of a difference on two types of mutual interactions between individuals in flocks with numerical simulation
16:10-16:25 Hao-Yuan Chang (National Cheng Kung University, M1) Near rings and BIBD
March 9 (Tuesday)
Below timetable is Taiwan time. Add 1 hour for Japanese side.
8:50-9:00 Opening
9:00-10:00 Session 6 Chair: Pu Zhao Kow (National Taiwan University)
9:00-9:15
Kazuya Okamoto (Musashino University, M2)
A nonlinear difference equation with bistability as a new traffic flow model
9:15-9:30 Chou- Kao (National Cheng Kung University, M2) Semi-Lagrangian schemes for level set equation 9:30-9:45
Kazuhito Toyota (Future University Hakodate, M1)
Proposal and simulation of a reservoir computing by Belousov- Zhabotinsky reaction to generate a sine wave
9:45-10:00 Zhi-Hao Shi (National Taiwan University, M2) Traveling wave solution for a stage structure model
10:20-11:30 Session 7 Chair: Kazuki Shigyou (Hiroshima University)
10:20-10:35 Tzu-Hsien Young (National Yang Ming Chiao Tung University, M2) Detection of change points for Weibull distributed time series data 10:35-10:50
Ya-Chi Wang (National Yang Ming Chiao Tung University, M2) Image inpainting algorithm based on partial differential equation approach
10:50-11:10
Yue-Hong Li (National Cheng Kung University, D1)
The asymptotic expansion of the trace of heat kernel on S2 under S1-action
11:10-11:30 Shun Ito (Meiji University, D1)
Control of spatio-temporal chaos of Ginzburg-Landau equation
13:00-14:05 Session 8 Chair: Yu-Shuo Chen (Tamkang University)
13:00-13:20
Pu-Zhao Kow (National Taiwan University, D3)
The Lewy-Stampacchia inequality for the fractional Laplacian and its application to anomalous unidirectional diffusion equations
13:20-13:40
Masakazu Kuze (Hiroshima University, D2)
Chemical oscillations and waves on microbead in Belousov- Zhabotinsky reaction
13:40-14:00
Kazuki Shigyou (Hiroshima University, PD)
Development of light field microscope for whole brain imaging of C.
elegans 14:00-14:05 Group Photo
16:00-16:30 Awards session
Lyapunov functions for some epidemic models with delay
Kai-Yang Liu
National Yang Ming Chiao Tung University Email: peter85912@gmail.com
As far as the current situation, disease transmission is an important issue. Predic- tive mathematical models for epidemics are fundamental to understand if the disease will break out or die out. In this talk, we will focus on the global stability of the solution of various models, such as SIS, SEIR, SEIS, SIQR, and SEIQR models, with distribution delay based on the construction of various Lyapunov functions.
Mathematical Study on Hidden Emotions and Body Patterns Observed in Cuttlefish
Kaito Fujihara1†, Kazuyuki Nakamura2, Yurika Nakai3, Haruhiko Yasumuro3, Kohei Okamoto3, Yuzuru Ikeda3, Mayuko Iwamoto4
1 Shimane University, 2 Meiji University, 3 University of the Ryukyus, 4 Doshisha University
n19m014@matsu.shimane-u.ac.jp
Cuttlefish has a high intelligence and visual ability, and can change “body pattern”
rapidly, which includes body colors, patterns and textures. It has been suggested that this specific ability has a role not only as camouflage, but also as a communication tool [1]. Although their body patterns had been classified by visual confirmation [2], it is considered still insufficient because cuttlefish’s body patterns are more variety and complicated than human eye can distinguish. The difficulty of classification makes it difficult for us to understand a regularity and a meaning of cuttlefish’s body patterns in communication.
The purpose of this research is to figure out the communication system by body patterns in cuttlefish. In this research, we used the experimental data in a situation where two cuttlefishes seem to be communicating. For a classification in detail, we extracted a series of image data of the body patterns from the videos, and classified them 30 categories using machine learning (ML). Then, to estimate changes in emotion behind the body pattern expressed, we applied a Hidden Markov Model (HMM) to the time-series data on the categories of body patterns. In this presentation, we will discuss an emotion estimation of the cuttlefish based on psychological models, using the estimated results of HMM and the behavioral data.
References
[1] Hanlon, R. T. & Messenger, J. B. “Cephalopod behaviour (2nd ed.)”, Cambridge University Press, Cambridge (2018)
[2] Nakajima, R. & Ikeda, Y. “A catalog of the chromatic, postural, and locomotor behaviors of the pharaoh cuttlefish (Sepia pharaonis) from Okinawa Island, Japan”, Marine Biodiversity, 47, 735-753 (2017)
A Method to Infer Latent Dynamics Model by Combining Deep Neural Networks with Ensemble Kalman Filter
Tsuyoshi Ishizone (1), Tomoyuki Higuchi (2), Kazuyuki Nakamura (1,3) 1: Meiji University; 2: Chuo University; 3: JST, PRESTO
tsuyoshi.ishizone@gmail.com
Inference of latent dynamics model plays an important role in time-series analysis such as prediction, imputation, and control for time-series data. To infer the model, FIVO and its variants combine deep neural networks (DNNs) with sequential Monte Carlo (SMC) and provides the state-of-the-art predictive performance [1, 2]. However, these methods have mainly two weaknesses: less particle diversity and biased gradient estimates. The ensemble Kalman filter (EnKF) which is widely used in the meteorology has good characteristics to overcome these drawbacks. In our presentation, we will review the characteristics of the FIVO and discuss the applicability of the EnKF to infer the nonlinear dynamics.
References
[1] C. J. Maddison, D. Lawson, G. Tucker, N. Heess, M. Norouzi, A. Mnih, A. Doucet, and Y. W. Teh. “Filtering Variational Objectives”, in NeurIPS, 2017.
[2] A. K. Moretti, Z. Wang, L. Wu, I. Drori, and I. Pe’er. “Variational Objectives for Markovian Dynamics with Backwards Simulation”, in European Conference on Artificial Intelligence, 2020.
The Data Forecast of COVID-19 by method of machine learning
Chun-Sheng Chen National Central University E-mail: jasonk123121898@gmail.com
By using the data from US, South Korea, Brazil, India, Russia, Italy and also adjusted the related data parameters in numerical simulations to generate forecast and the effect of prediction of the epidemic situation for each country. Meanwhile to utilize US data to compare SQIARD with SIARD, and show the effects of predictions.
In this lecture, i will use mathematical models and data analysis methods to analyze the proportion of asymptomatic infections in various countries and predict the future trend of the epidemic. Then one of way is to rewrite the differential equations as discrete form with the consideration in [1]. Therefore, we establish two mathematical models, SQIARD and SIARD model, to simulate the COVID-19 epidemic. Obviously the SIARD model is only a simplified form of the SQIARD model. The difference between this two models are one include the parameters “number of people to be screened for the epidemic” and another one exclude it. It is very hard to get the data for this parameter but fortunately US government has complete data on the number of people in quarantine. This data can be used to generate the model parameters. In order to get the related model parameters, to remove the parameter Q(t) first and then use SIARD model to conduct the prediction data of US, South Korea, Brazil, India, Russia and Italy. The prediction of the SIARD model will be created first and then use US’s data which contains Q(t) to show the effect of prediction of the SQIARD model.
Since the most serious epidemic of COVID-19 was in 2020, so we take training data:
5/11 ∼ 8/19 and validation data: 8/20 ∼ 9/8. After training our model, i use this model to predict future: 9/9 ∼ 9/28. Finally, i use this model on date: 1/11 to predict 1/11 ∼ 1/31, until on date 1/31, i can use the really data to compare the predict data obtained by our model.
References
[1] Yi-Cheng Chen, Ping-En Luy, Graduate Student Member, IEEE, Cheng-Shang Chang, Fellow, IEEE, and Tzu-Hsuan Liux. A Time-dependent SIR model for COVID-19 with Undetectable Infected Persons. IEEE (2020)
Why do two-dimensional dunes always merge?
Yusuke Fujita1, Hiraku Nishimori2, Makoto Iima1
1Program of Mathematical and Life Sciences, Graduate School of Integrated Sciences for Life, Hiroshima University
2Meiji Institute for Advanced Study of Mathematical Sciences, Meiji University e-mail: m195052@hiroshima-u.ac.jp
Sand dunes are mounds generated through the transport of sand caused by the wind. We focus on barchan dunes and transverse dunes, both of which are generated under the unidirec- tional flow condition. In the vertical cross-section along the flow direction, the barchan dunes and transverse dunes are simplified by a triangular object in two-dimensional space.
Dunes move and interact (one another). Barchan dunes, dunes with a three-dimensional crescent-like structure, show a variety of interactions; merging, splitting, and ejection, during the collision [1]. On the other hand, the transverse dunes, dunes with quasi two-dimensional structure, always merge during the collision [2]. Clearly, the difference of collision behavior comes from the detailed flow pattern depending on the detailed shape of the dunes, however, the mechanism has not been clarified.
In this study, we analyzed flow structure around two-dimensional dune models, triangular cross-sections, by direct numerical simulation. When two cross-sections are placed at particular intervals, a dead-water region was observed between them, where the flow velocity was much slower than that in the surrounding area (Fig. 1). Further we analyze sand transport based on the flow velocity distribution, which leads to the conclusion that two-dimensional dunes always merge [3].
Figure 1: Snapshots of flow speed around dune models (Re = 1000). Height of model is h = 1.
Left: peak-to-peak distance is D = 6. Vortices are generated from both tops of the dunes.
Right: peak-to-peak distance is D = 3.9. Dead-water regions (blue) are observed between the dunes and downstream.
References
[1] A. Katsuki, M. Kikuchi, H. Nishimori, N. Endo, and K. Taniguchi, Cellular Model for Sand dunes with Saltation, Avalanche and Strong Erosion: Collisional Simulation of Barchans, Earth Surf. Process. Landforms 36, 372, (2011)
[2] A. D. Ferreira and M. R. M. Fino, A Wind Tunnel Study of Wind Erosion and Profile Reshaping of Transverse Sand Piles in Tandem, Geomorphology 139-140, 230, (2012) [3] Y. Fujita, H. Nishimori, and M. Iima, Dead-Water Region around Two-Dimensional Sand-
dune Models, J. Phys. Soc. Jpn. 89, 063901, (2020).
The Second Laplacian Eigenvalue
Yu-Jie Ho
National Yang Ming Chiao Tung University E-mail: yujie.am05@nctu.edu.tw
The second smallest normalized Laplacian eigenvalue is closely related to the con- nectivity of a graph. The Cheeger inequality states that graphs with a small second eigenvalue can be easily divided into two, and graphs with a large second eigenvalue are highly connected. In the presentation, I will show the application of graphs with a small second eigenvalue to image partition, and the application of graphs with a large second eigenvalue to computer networking.
Discontinuous steady states of the nonlocal Allen-Cahn equation
Ryu Fujiwara
Meiji University, Graduate School of Advanced Mathematical Sciences cs191009@meiji.ac.jp
Since Turing stated in 1953 that reaction diffusion equations exhibit complicated pat- tern due to the diffusion induced instability, a variety of reaction diffusion equations have been proposed and analyzed. There is much research about the structure of the solution of reaction diffusion equations, in particular, the Allen-Cahn equation. For example, Allen-Cahn equation has a travelling wave solution [2]. Moreover, the con- vexity of the domain is essential for the uniform solutions to be only stable steady states of the Allen-Cahn equation [5, 3].
Stable nonuniform steady states can appear in reaction diffusion equations with non- local effects, called nonlocal reaction diffusion equations. Bates and Chmaj analyzed the Allen-Cahn equation with nonlocal diffusion term defined by the convolution, and showed there exist stable and discontinuous steady states [1]. Kaliuzhyni-Verbovedsky and Medvedev in [4] formulated a nonlocal reaction-diffusion equation as
∂u
∂t(x, t) = D[Lu](x, t) + f (u(x, t)), [Lu](x) = Z 1
0
W (x, y)u(y)dy −
Z 1 0
W (x, y)dy
u(x), where W is a measurable function on [0, 1]2, f is a general function in R and D > 0.
We call the equation the nonlocal Allen-Cahn equation if f is bistable.
There are two goals for this research. One is to show the existence of the stable and discontinuous steady states of the nonlocal Allen-Cahn equation. We also show numerical results, which imply the existence of such solutions even if the sufficient conditions do not hold. The second goal is to show the existence of the solution globally in time for a general function f by proving that the nonlocal operator L is a sectorial operator.
References
[1] P. W. Bates and A. Chmaj, An integrodifferential model for phase transitions:
Stationary solutions in higher space dimensions, J. Statistic. Phys., 95 (1999), 1119–1139.
[2] P. C. Fife and J. B. McLeod, The approach of solutions of non-linear diffusion equations to travelling front solutions, Arch. Rat. Mech. Anal., 65 (1977), 335–
361.
[3] J. K. Hale and J. M. Vegas, A nonlinear parabolic equation with varying domain, Arch. Rational Mech. Anal., 84 (1984), 99–123.
[4] D. Kaliuzhyni-Verbovetskyi and G. S. Medvedev, The semilinear heat equation on sparse random graphs, SIAM J. Math. Anal., 49.2 (2017), 1333–1355.
[5] H. Matano, Aysmptotic behaviour and stability of solutions of semilinear elliptic equations, Publ. RIMS. Kyoto Univ., 15 (1979), 401–454.
The Convexity for The Joint Range of Two Quadratic Functions
Ya-Chi Chu
National Cheng Kung University Email: ycchu.academic@gmail.com
Given n × n symmetric matrices A and B, Dines in 1941 proved that the joint range set {(xTAx, xTBx)| x ∈ Rn} is always convex. Our work is concerned with non- homogeneous extension of the Dines theorem for the range set R(f, g) = {(f (x), g(x)) | x ∈ Rn}, f (x) = xTAx + 2aTx + a0 and g(x) = xTBx + 2bTx + b0. In this talk, we will illustrate that R(f, g) is convex if, and only if, any pair of level sets, {x ∈ Rn|f (x) = α}
and {x ∈ Rn|g(x) = β}, do not separate each other. With the novel geometric concept about separation, we also provide a polynomial-time procedure to practically check whether a given R(f, g) is convex or not.
This is a joint work with Professors Huu-Quang Nguyen and Ruey-Lin Sheu.
Mathematical modeling and analysis of filamentous active matter
Naoki Kondo (1), Elliott Ginder (2)
Graduate School of Advanced Mathematical Science, Meiji University 1 : cs191007@meiji.ac.jp
2 : eginder@meiji.ac.jp
Loosely defined, active matter refers to agents that displaying a function of self motion; such as the motion of flocks of birds, as well as camphor scrapings atop water.
Motivated by experiments from S. Nakata et al. [1], we consider the following model equation:
τ ut = ∆u − u + δΓ t > 0, γ(t, s) ∈ Γ, (x, y) ∈ R2
˙γ(t, s) = −∇u(γ(t, s)) + F t > 0, s ∈ [0, L(γ)), γ ∈ R2 u(0, x, y) = 0 (x, y) ∈ R2
lim
s→L(γ)
γ(t, s) = γ(t, 0) t > 0
γ(0, s) = γ0(s) s ∈ [0, L(γ))
(1)
The above model includes a reaction-diffusion equation and a Newtonian equation for expressing the motion of filamentous active matter which causes the Marangoni convection. Given initial and boundary conditions, u denotes the concentration of a chemical substance within the domain R2, γ(t, s) = (γ1(t, s), γ2(t, s)) designates the shape of the active matter, δΓ is a delta function, L(γ) is a length of Γ, s ∈ (0, L(γ)) is an arc-length parameter, F (Γ) is an outer force acting on the active matter, and τ is a nonnegative constants.
In this talk, we will discuss the stability analysis of the model equation, using the method of adiabatic elimination. In particular, we derive analytic expressions for ring- shaped stationary solutions in terms modified Bessel functions and show the results of numerical simulations.
References
[1] S. Nakata, K. Kayahara, M. Kuze, E. Ginder, M. Nagayama, and H. Nishimori, Soft Matter, Synchronization of self-propelled soft pendulums 2018, 14, pp. 3791–3798
Using machine learning methods to predict RNA secondary structure
Huan-Chi Chang
National Chung-Cheng University Email: mia701107@gmail.com
Machine learning has been widely used in many fields recently, including biolog- ical information. RNA plays an important role in the processes of genetic coding, translation, and gene expression. To understand the functionality of RNA, scientists need to know its secondary structure formed by the pairing of nucleotides. The struc- tures can be obtained by biological experimental technology directly, however, our research intends to use machine learning methods to predict RNA secondary structure through Neural Network, Random Forest, Extreme Gradient Boosting (XGBoosting), LightGBM (Light Gradient Boosting Machine, LGBM), etc. Comparing the predicted results with the experimental laboratory structure, calculating the average prediction accuracy, and exploring the reasons for the poor cases, we try to find feasible correction methods to improve the accuracy of predicting RNA secondary structure by machine learning methods.
Synchronous solutions and their stability in coupled salt-water oscillators
Koutaro Tanaka† Yasuhisa Saito Shimane University
†n20m008@matsu.shimane-u.ac.jp
A salt-water oscillator is a vibration phenomenon in which salt-water outflow and water inflow occur alternately at regular intervals by placing a cup with a small hole in the bottom containing salt-water in a container filled with water. In this talk, we will present stability results on synchronous solutions in coupled salt-water oscillators where two salt-water oscillators are placed close to each other, expressed by following equations
d2x1 dt2 − ε
(
α − β dx1 dt
2) dx1
dt + ω2x1 = δ
a11x1+ a12x2+ b11dx1
dt + b12dx2 dt
d2x2 dt2 − ε
(
α − β dx2 dt
2) dx2
dt + ω2x2 = δ
a21x1+ a22x2+ b21
dx1 dt + b22
dx2 dt
where x1, x2 represent the height difference between water surface and salt-water sur- face, and α, β, ω, ε, δ are positive constants with 0 < δ ε 1. Our results include mathematically new findings not found in [1].
References
[1] K. Yoshikawa, H. Kawakami, Nonlinear dynamics for chemical reaction sys- tems(Japanese), JSIAM Vol.4 No.3 Sept.(1994), pp. 238–258
Mathematical models of COVID-19 with diffusion effects and their data forecast
Yen-Jia, Chen
National Taiwan Normal University Email: rex19981002@gmail.com
In this paper, we introduced a SIARD model with diffusion. We discussed that there is a weak solution in the model, and then discussed the condition of basic reproduction number and disease free equilibrium of the model. Finally, we show the numerical simulation on SIARD diffusion model. After that, we train the related data parameters, in our numerical simulations, we respectively conduct the forecast of the data of US, South Korea, Brazil, India, Russia and Italy, and the effect of prediction of the epidemic situation in each country. Furthermore, we apply US data to compare SQIARD with SIARD, and display the effects of predictions.
On Pehlivan’s four-scroll chaotic system
Sakai Ryosuke
Graduate School of Advanced Mathematical Sciences, Meiji University cs191003@meiji.ac.jp
Since 1970, chaotic dynamical system have been attracting attention from many re- searchers. Lorenz equation[1] is one of the typical example of such systems. It has a double-scroll chaotic attractor called the Lorenz attractor. In recent years, the subject of interest has shifted to the study of multi-scroll attractors, which have more complex topological structures. In fact, there have been reports about the existence of 2, 3 and 4 scroll attractors. Therefore, it is important to elucidate the structure of the chaotic attractor of such chaotic systems. The purpose of this paper is to discuss the dynamics of Pehlivan’s model[2], which has a four-scroll attractor(see Figure 1) with a smooth function. It has been suggested by the previous research[3] that Pehlivan model has an invariable plane when b = 1 and its vector field affects the structure of the four- scroll attractor. This study focused on Pehlivan model when b=1. In this research, we revealed the dynamics when b = 1 and found the geometric structure of the attractor from the analysisby Poincare map when b is close to 1.
Figure 1: four-scroll attractor
References
[1] E. N. Lorenz, “Deterministic Nonperiodic Flow”, Atmospheric Sciences, 20, 130- 141, (1963)
[2] I. Pehlivan, I. M. Moroz, S. Vaidyanathan, “Analysis, synchronization and circuit design of a novel butterfly attractor”, Journal of Sound and Vibration, 333(20), 5077-5096, (2014)
[3] T. Miyaji, H. Okamoto, A. D. D. Craik, “A Four-Leaf Chaotic Attractor of a Three Dimensional Dynamical System’,International Journal of Bifurcation and Chaos, 25(2), (2015)
On Synchronization Analysis of Complex Coupled Kuramoto Oscillators
Jin-Zhi Phoong National Taiwan University
phoong1205@gmail.com
This paper studies the synchronization problem emerging in a network of coupled oscillators which can be found in various natural phenomena. Here, the complex sce- nario has been considered. The sufficient condition for the frequency synchronization has been established for the uniform case.
Periodic Solutions for Three Species Lotka-Volterra Competitive Equations
Ryota Kubohara† Yasuhisa Saito
Graduate School of Natural Science and Technology, Shimane University
†n20m005@matsu.shimane-u.ac.jp
There have been many studies on three species Lotka-Volterra competitive systems dX
dt = (ε1− X − β10 Y − r10 Z) X dY
dt = (ε2− r20 X − Y − β20 Z) Y (2) dZ
dt = (ε3− β30 X − r03Y − Z) Z,
where εi, βi0 and r0i (i = 1, 2, 3) are positive constants (for example, [1]). For the cyclic competition case βi0 = 0 (i = 1, 2, 3), the existence of periodic solutions was proven by [2]. In this talk, we try to generalize the result to non-cyclic competition cases (2) by some new method making use of several quadratic curves.
References
[1] M. L. Zeeman, Hopf bifurcations in competitive three-dimentional Lotka-Volterra systems, Dynamics and Stability of Systems 8(3), pp. 189–216
[2] H. Nakajima, Stability and periodicity for ecological system model (Japanese), Busseikenkyu 29(6), pp. 345–387
Combine Visual Cryptography with RSA Algorithm
Jian-Yu-Chen
National Central University Email : jianyu890619@gmail.com
In today’s age of advanced Internet, when we send an important image to others, it is difficult to guarantee that it will be transmitted with security. Thus, the encryption of the image is very important. Traditional visual cryptography used to extend to the k out of n scheme where a secret image is encrypted into n shares but only k shares are needed for decryption where
k ≤ n
However, this method often accompanied by image contrast reduction, and when putting it into computer, it might be cracked since it’s density is not uniform. Thus, we construct a method to encrypt the image by combining with RSA Algorithm [1] . Thank to the fact that finding the factors of a large composite number is difficult, RSA Algorithm is still not easy to crack. By encrypting the pixels of the image by RSA Algorithm to increase the security of encryption, this method attain the result what we are expected.
References
[1] The RSA Cryptosystem: History, Algorithm, Primes, Michael Calderbank, 2007
Comparison principle for time fractional diffusion equations
Yuto Sakai Ryukoku University t20m005@mail.ryukoku.ac.jp Introduction
We consider the following a time fractional diffusion equation (c
0Dα0u = uxx+ f (u), x ∈ R, t > 0, u(x, 0) = u0(x) ≥ 0, x ∈ R,
where 0 < α < 1 and c0D0αu is a Caputo derivative, which is defined as in
c
0Dα0u(x, t) := 1 Γ(1 − α)
Z t 0
(t − s)−α∂u(x, s)
∂s ds.
Here f (u) is a globally lipschitz continuous function and u0is a continuous function. we will give a proof of the comparison principle for the time fractional diffusion equation.
In order to prove it, we use the Mittag-Leffler function, which is one of generalization of exponential functions.
References
[1] S.D.Eidelman, A.N.Kochubei, Cauchy problem for fractional diffusion equations, J. Differential Equations 199 (2004) 211–255.
[2] J.Kemppainen, J.Siljander, R.Zacher, Representation of solutions and large-time behavior for fully nonlocal diffussion equations, J. Differential Equations 263
(2017),149–201.
Turing Instability for a Class of 3-component Systems
Yueh-Tzu Hung
National Yang Ming Chiao Tung University Email: kscity806.am06@nctu.edu.tw
We investigate the Turing instability for three-component systems with cross- diffusion. The main goal is to understand the effect of cross-diffusion on the stationary instability (S-instability) and wave stability (W-instability), respectively. Some suffi- cient and/or necessary conditions for Turing instability to occur are provided. Numer- ical simulations and applications are given to illustrate our theoretical results.
A consideration of a difference on two types of mutual interaction between individuals in flocks with numerical
simulation
Yuto Kikuchi1∗ , Mayuko Iwamoto2 1 Shimane University 2 Doshisha University
*n20m004@matsu.shimane-u.ac.jp
A behavior of self-propelled particles with rules of mutual interactions, which is a model for a flock of organisms [1], is greatly affected by kinds of the rule of inter- actions. Previous study in a jackdaw flock has reported that there are two types of interactions [2]: “metric interaction” that takes the average of vectors in all particles which exist within a certain range, and “topological interaction” that takes the average of vectors in a specific number of particles. The former is considered to be used in a situation for avoiding predators, while the latter is used for group flight. These simple rules of interaction help us to understand a flock-like behavior, but careful studies are required when we consider a mechanism of real organisms’ flock behavior. For ex- ample, it is not reasonable that in “metric interaction” individuals who exist slightly outside of the specified distance would not interact, and in “topological interaction”
individuals would interact if it is exactly within the specified number. In addition, it seems to be difficult for organisms to switch adequately between two rules in response to a situation, so we consider to be unlikely that organisms have a specific number for interaction, and have multiple rules for interaction in their environment.
The purpose of this study is to suggest that an idea of a plausible interaction rule which realize both flock behaviors observed in metric and topological interactions. In this study, we will propose a new model of the interaction using flock’s density: the smaller the density, the larger the distance over which the certain interact range. In this presentation, as a first step for building a model, numerical simulations with each model by two kinds of interactions are performed to verify each flock behavior and discuss the features of each interaction.
References
[1] Tamas Vicsek, et al., Novel Type of Phase Transition in a System of Self-Driven Particle, Physical Review Letters , 75(6), 1226-1229 (1995)
[2] Hangjian Ling, et al., Behavioural plasticity and the transition to order in jackdaw flocks, Nature Communications, 10, 5174 (2019)
Near rings and BIBD
Hao-Yuan Chang
National Cheng Kung University E-mail: L16091051@gs.ncku.edu.tw
A right near ring is a triple (N, +, ∗) where N is a set and + and ∗ are binary operations on N such that
1. (N, +) is a group (not necessarily ablelian);
2. multiplication is associative;
3. (a + b)c = ac + bc for all a, b, c ∈ N ;
I will introduce some interesting fact about near rings, and a special case called planar near rings. In 1968, Ferrero found a way to produce many examples of planner near rings, and we can produce balanced incomplete block designs (BIBDs) from them.
BIBD is a subject in combinatory which applied in experimental design, so the result helps giving more examples. Moreover, we can design a code base on those block designs made by planar near rings.
A nonlinear difference equation with bistability as a new traffic flow model
Kazuya Okamoto
Department of Mathematical Engineering, Graduate School of Engineering, Musashino University, 3-3-3 Ariake, Koto-ku, Tokyo 135-8181, Japan
poulkun67@gmail.com Akiyasu Tomoeda
Department of Mathematical Engineering, Faculty of Engineering, Musashino University, 3-3-3 Ariake, Koto-ku, Tokyo 135-8181, Japan
Macroscopic traffic flow models have the usefulness of reducing the cost of large-scale numerical simulation on traffic networks. This is because the macroscopic traffic flow model does not focus on each car in mathematically expressing the movement of the car, but on how dense the cars are in a lattice. By expressing this, it is not necessary to calculate the number of cars. In this session, we introduce a new macroscopic traffic flow model composed by the difference equation and linear analysis for it. This model is a reference model [1] with the past time added, which represents the inertia of the car. We show that the parameters of our model have the potential to eliminate traffic jams from the perspective of stability. It is understood that changes in average density change the stability of homogeneous flow, resulting in traffic jams. In other words, the homogeneous flow is stable when the average density is low, but when it exceeds the critical density, it becomes unstable and cannot be maintained, and transitions to inhomogeneous flow, which is observed as traffic jams. Actually, the conditions under which the homogeneous flow becomes unstable are derived from the linear stability analysis, and we consider the conditions under which traffic jams occur based on the unstable region. In addition, the phase diagram of the model was obtained using numerical simulation, and it was shown that the model has the property of bistability, which is very important in the traffic flow models [2]. Bistability is the property that two solutions are stable with respect to the same mean density. We explain that the model shows a density region where both homogeneous and inhomogeneous flows are stable. Then it shows numerically that the inhomogeneous flow is traveling in the opposite direction of the traveling direction of the car, like the traffic jams seen in the actual traffic flow.
References
[1] K. Higashi, J. Satsuma and T. Tokihiro, Rule 184 fuzzy cellular automaton as a mathematical model for traffic flow, 11 Apr 2020, arXiv:2004.05294
[2] M. Treiber, A. Kesting, Traffic flow dynamics: Data, Models and Simulations., (2013) Springer.
Semi-Lagrangian Schemes for Level Set Equation
Chou Kao
National Cheng Kung University E-mail: L16071158@gs.ncku.edu.tw
Semi-Lagrangian (SL) scheme is widely applied in computation of first order equation where the solution is discretized by Lagrangian formulation in time over a fixed grid in space. In this talk, I will discuss SL scheme for one dimensional scalar convection equation and Eikonal equation, where the solution curve moves in given velocity and in normal direction respectively. WENO scheme as a higher order interpolation will be mentioned. Then I will extend to SL approximation of level set equation in two dimensional space where a quasi-fifth order interpolation technique is proposed.
Proposal and simulation of a reservoir computing by Belousov-Zhabotinsky reaction to generate a sine wave
Kazuhito Toyota
Complex Information Science Field, Future University Hakodate g2120030@fun.ac.jp
The BZ(Belousov-Zhabotinsky) reaction is a chemical reaction in which the redox state changes slowly and cyclically. The BZ reaction is known to produce complex structures such as oscillatory rhythms, excitable patterns and chaos. It also has the property of changing its period under the influence of light stimulation. The BZ reac- tion is often used in experiments because it can represent complex chemical reaction systems with few reagents [1].
One of the recurrent neural networks that is well known as machine learning is a method called the ”reservoir computing” [2]. It is often characterized by the use of physical phenomena in the reservoir laid in the middle layer of a machine learning. In the past, other researches of the reservoir computing have been implemented using a bucket of water or a soft material that resembles an octopus foot as a reservoir [3].
However, there have not been many implementations yet due to problems such as the time and effort required for preparation. Conversely, the BZ reaction can be easily experimented and various chemical reaction systems can be represented by adjusting the reagents. Therefore, in this study, we try to propose the reservoir computing driven by BZ reaction. We regard the BZ reaction as a reservoir.
We investigate whether it is theoretically and numerically possible to perform a reservoir computing by BZ reaction. We compute the partial differential equation model of the BZ reaction numerically in a two-dimensional domain and placed some units to distill the time series data on that domain. We assume that the target output is set as a sine wave, and the back projection from the output is given by the light stimulation in the model. Varying the parameter of strength for the light in the model as a sine wave, we verify whether to output a sine wave similar to the light stimulation by the linear sum of the time series of the units. As a result, it is found that the sine wave could be output roughly.
References
[1] Hidetoshi Miike, Yoshihiro Mori, Tomohiko Yamaguchi (1997). Chemistry of non- equilibrium systems III Dynamics of reaction and diffusion systems (Hiheikoukei no Kagaku III Hannou · Kakusankei no Dainamikusu), Koudansha.
[2] H. Jaeger (2002). A tutorial on training recurrent neural networks, covering BPPT, RTRL, EKF and the ”echo state network” approach, GMD Report 159, German National Research Center for Information Technology, pp1-46.
[3] Kouhei Nakajima (2020). Physical reservoir computing -an introductory perspec- tive, Japanese Journal of Applied Physics 59, 060501, pp. 1-13.
Traveling wave solution for a stage structure model
Zhi-Hao,Shi
National Taiwan University Email: n12603579table@gmail.com
In the first part of this short talk, I will introduce some preliminary researches on a non-time-delayed two age stage structure model in mathematical ecology (in brief, stage structure model). Its characteristic is, different from the single stage model that uses a scalar function (or constant value) to represent its growth rate (or to say: fitness), the staged structure model uses a 2*2 matrix function (or constant matrix) to describe the succession of sub-populations in different life history stages of the total population.
Many variants of this model were first modeled between the late 20th century to the beginning of the 21st century. Also, some studies of parabolic problems have been carried out since early 21st century, for example, classical papers are [1] [2] [3].
What’s more, the minimum speed of the traveling wave fronts of this class of model was not determined until 2020 by Huang et al. [4] . Besides, the existence of the traveling wave and the linear determinancy of the wave velocity have not been further investigated. In the rest of this short talk, I will talk about some progress on the existence of the traveling wave solution of a special case based on the recent studies by Prof. Chen and me.
References
[1] Schreiber, S., Rudolf, V. H. (2008). Crossing habitat boundaries: coupling dynamics of ecosystems through complex life cycles. Ecology letters, 11(6), 576-587.
[2] Bouguima, S. M., Fekih, S., Hennaoui, W. (2008). Spacial structure in a juvenile- adult model. Nonlinear Analysis: Real World Applications, 9(3), 1184-1201.
[3] Bouguima, S. M., Mehlia, F. Z. (2012). Asymptotic behavior of an age-structured population model with diffusion. Journal of Applied Analysis and Computation, 2(4), 351-362.
[4] Huang, Q., Zhang, Y. (2020). Spread rates of a juvenile-adult population in constant and temporally variable environments. Theoretical Ecology, 1-16.
Detection of Change Points for Weibull Distributed Time Series Data
Tzu-Hsien Young
National Yang Ming Chiao Tung University E-mail: yhps910454@gmail.com
Wind power is a valuable sustainable energy that many countries have been agreesively developing. In energy resource assessments, understanding the regional wind properties is a crucial problem. A common way of the assessments is to model the wind speed data with a two-parameter Weibull distribution, in which the scale and shape parameter are estimated from the measurement data. This process suc- ceeds when the data follow a unique specific Weibull distribution. However, if there are change points, at which the wind speed distribution changes, the changes in wind properties cause the statistical error margins in the energy resource assessments. To address the issue more accurately, one approach is to identify the change points and perform the assessments on the dissected sectional data. For this purpose, we propose to use the prune exact linear time algorithm [1], combining with detection of changes of Weibull distribution parameters. The PELT algorithm has an excellent linear com- plexity, O(n), which is particularly suitable for large amount of time series data. The change point detection is based on the maximum-likelihood method. To demonstrate the method, we first construct thousands synthesized wind speed time series, contain- ing two sets of Weibull distributions, i.e. with one deliberately added change point.
The detection result indicates that the present method resolves the change points more accurately than the classical change point detection scheme with normal distributions.
Finally, we apply the method to the real Fuhai offshore wind mast data, taken between 2017/11/18 and 2018/10/26. The present method identified two significant change points occurring on 2018/02/26 and 2018/09/23, and they agree precisely with the transitional periods between the winter and summer monsoon seasons in Taiwan.
This research is supported in parts under MOST, 109-2218-E-001 -003 and is a joint work with Professors Pi-Wen Tsai, Chih-Yu Kuo and Wen-Yi Chang.
References
[1] R. Killick and P. Fearnhead and I. A. Eckley, Optimal Detection of Changepoints With a Linear Computational Cost, J. Am. Stat. Assoc., 107 (500), 1590-1598
Image Inpainting Algorithm Based on Partial Differential Equation Approach
Ya-Chi Wang
National Yang Ming Chiao Tung University E-mail: happy099088@gmail.com
In this talk, we numerically study the algorithms for image inpainting based on energy minimization and the low dimensional manifold. Li and Yao [1] have presented two mathematical models based on an energy function for image inpainting . Numerical schemes are introduced to solve the corresponding Euler-Lagrange equations. As for Low Dimensional Manifold Model (LDMM) which was studied for image processing in [2], LDMM uses the dimension of the patch manifold and the Euler-Lagrange equation is obtained by applying weighted nonlocal Laplacian (WNLL).
In this talk, numerical experiments are presented for removing the unwanted region by using Li-Yao model and LDMM ,respectively. This is a joint work with Professor Ming-Cheng Shiue.
References
[1] Li, S.J., Yao, Z.A., Low Dimensional Manifold Model for Image Processing, Imaging Sci. J., 2013, 61(3), pp. 292–300
[2] S. Oscher, Z. Shi, and W. Zhu, Low Dimensional Manifold Model for Image Pro- cessing, SIAM J. Imaging Sciences, 2017, 10, pp. 1669–1690
The Asymptotic Expansion of the Trace of Heat Kernel on S
2under S
1-action
Yue-Hong Li
National Cheng Kung University Email: L18091023@gs.ncku.edu.tw
Given a n-dimensional compact manifold M and a compact group G acting on M of isometries. We denote M the Laplacian on M and Eλ the eigenspace with respect to λ, an eigenvalue of − M. Consider (ρ, V ) a finite dimensional irreducible representation of G. Jochen Br¨uning and Ernst Heintze defined a function and show that in [2]
Lρ(t) :=X
λ≥0
e−λtdimHomG(V, Eλ) ∼ t−n2
∞
X
j=0
ajtj2, t → 0+
if G is of ranking 1.
In this talk, we restrict our attention on a specific case, M = S2 and G = S1 and introduce the method mentioned in [3] to achieve the coefficients of the expansion via studying the asymptotic expansion of the following double integral
Z 0
Z f (δ,y) 0
e−x2y2t xαyβdxdy
As a result, we want to find how does the geometry information(e.g. metric) affect on the coefficients aj in Lρ(t).
References
[1] S. Minakshisundaram and ˚A. Pleijel, Some Properties of The Eigenfunctions of The Laplace-Operator on Riemannian Manifolds, Canadian J.Math (1949), 242- 256.
[2] Jochen Br¨uning and Ernst Heintze, The Asymptotic Expansion of Mi- nakshisundaram -Pleijel in Tne Equivariant Case, Duke Mathematical Journal (1984).
[3] Jochen Br¨uning, On The Asymptotic Expansion of Some Integrals, Arch. Math (1984), 253-259.
[4] R. Wong, Asymptotic Approxiamtions of integrals, SIAM(2001).
Control of spatio-temporal chaos of Ginzburg-Landau equation
Shun Ito Meiji University catastrovanity@gmail.com
Ginzburg-Landau equation has two types-behavior: one is spatio-temporal chaos, the other is a limit cycle oscillation that is a spatially uniform periodic solution. Spatio- temporal chaos is located to near a fixed point inside a limit cycle on phase space. We numerically discovered that spatio-temporal chaos shifted to an outside of a limit cycle by a perturbation like a short time impulse converges on a limit cycle oscillation. Then, we prove analytically that an initial condition at an outside of limit cycle converges on a limit cycle oscillation.
The Lewy-Stampacchia Inequality for the Fractional Laplacian and Its Application to Anomalous Unidirectional
Diffusion Equations
Pu-Zhao Kow
National Taiwan University Email: d07221005@ntu.edu.tw
We consider a Lewy-Stampacchia-type inequality for the fractional Laplacian on a bounded domain in Euclidean space. Using this inequality, we can show the well- posedness of fractional-type anomalous unidirectional diffusion equations. This study is an extension of the work by Akagi-Kimura (2019) for the standard Laplacian. However, there exist several difficulties due to the nonlocal feature of the fractional Laplacian. We overcome those difficulties employing the Caffarelli-Silvestre extension of the fractional Laplacian.
This talk is prepared based on [1], which is a joint work with Masato Kimura.
References
[1] Pu-Zhao Kow and Masato Kimura, The Lewy-Stampacchia Inequality for Frac- tional Laplacian and Its Application to Anomalous Unidirectional Diffusion Equations. https://arxiv.org/pdf/1909.00588.pdf(arXiv:1909.00588).
Chemical oscillations and waves on microbead in Belousov-Zhabotinsky reaction
Masakazu Kuze
Graduate school of Integrated Science for Life, Hiroshima University kuzem1234@hiroshima-u.ac.jp
The Belousov−Zhabotinsky (BZ) reaction, which self-organizes spatio-temporal patterns based on the autocatalytic reaction, is a widely studied experimental system. The metal catalyst of the BZ reaction was loaded into a cation-exchange resin bead. We call this bead the BZ bead. Global oscillations (GO), which are the uniform oscillations in the entire bead, and traveling waves (TW), which are the propagating chemical waves from one end of the bead, were observed when BZ bead was immersed into a BZ solution without a catalyst (Fig. 1a).[1] Previous study has reported that GO and TW could be selectively and reversibly generated by positive and negative values of the electrical potential, E, respectively. In this study, we found the hysteresis on the switching of the spatio-temporal patterns depending on the scanning direction of E. In addition, we discuss our results in the context of an increased concentrations of the activator, HBrO2, or the inhibitor, Br−, to clarify the mechanism of switching between GO and TW.
The switching from GO to TW occurred near –0.2 V in negative scan, i.e., from E = +1.0 V to −1.0 V (Fig. 1b). On the other hand, when E was scanned from −1.0 V to +1.0 V (positive scan), the switching from TW to GO occurred near +0.2 V. In the negative scan, GO were maintained because oscillations in the bead were strongly affected by sufficient concentration of the inhibitor which accumulated near the positive electrode.
On the other hand, TW was maintained at a small value of E due to the electrochemical production of the activator in the positive scan.
Figure 1. (a) Snapshots of GO and TW for the single bead. Dotted circles correspond to the edges of the beads. Time interval: 2 s. (b) The generation point of the oscillations, Rl, at d = 0.63 mm of the BZ beads at the scanning of the electrical potential (E). Gray and black circles correspond to the negative and positive scan, respectively.
References
[1] M. Kuze, H. Kitahata, O. Steinbock, S. Nakata, Distinguishing the Dynamic Fingerprints of Two- and Three Dimensional Chemical Waves in Microbeads., J. Phys. Chem. A 122(8), pp. 1967–1971
Development of Light field microscope for whole brain imaging of C. elegans Kazuki Shigyou1, Haruka Maeoka1, Reiji Yagi2,
Kazuki Shimabukuro2, Shin Usuki3, Takuma Sugi1
1 Graduate School of Integrated Sciences f or Lif e, Hiroshima U niversity, Higashihiroshima, Japan
2Department of Mechanical Engineering, Shizuoka University, Hamamatsu, Japan
3 Graduate School of Science and Technology, ShizuokaUniversity, Hamamatsu, Japan kshigyo@hiroshima-u.ac.jp
Living organisms with nervous systems perform very sophisticated information processing, despite the fact that their brain volume is much smaller than that of latest computing units. For example, flies can instantly identify the spatial information of an oncoming human hand and, through fast feedback processing, can switch to avoidance behavior faster than humans can catch them [1]. In addition, even Caenorhabditis elegans (nematode), which have only 302 neurons, 1 mm long, and 50 to 100 m thick can memorize temperature [2] and odor, and dynamically change their behavior based on these memories. If these brain functions can be implemented as a device, it will contribute to realizing high-speed collision avoidance systems in cars and the development of advanced odor sensors. In order to understand and apply these function of the brain, we need to observe and quantify the whole brain activity. However, the time scale of the neural activity is milli-second order and the spatial scale is micro-meter to centi-meter scale with a complex 3D structure. Activities of the nervous system are difficult to observe with conventional observation equipment such as a confocal scanning laser microscope, because it has not enough temporal resolution. To solve this problem, Light field microscopes were developed by Leboy group in 2006 [3]. Light field microscopes can capture a three- dimensional structure in a single camera shot. Although this makes it possible to capture these 3D structures and nervous activity on milli-second order limited by camera frame rate, the spatial resolution of conventional light field microscopy is not enough to observe neural cells. Therefore, we developed a high-spatial resolution light field microscope. In this presentation, we will explain the principle of light field microscopes and observation results of the microscope we developed.
Figure 1: Conventional light field microscope image(Left). Our result(Right). Red scale bar is 100um.
References
[1] Trevor J. Wardill, Katie Knowles, Laura Barlow, Gervasio Tapia, Karin Nordstr¨om, Robert M. Olberg, and Paloma T. Gonzalez-Bellido. The Killer Fly Hunger Games: Target Size and Speed Predict Decision to Pursuit. Brain, Behavior and Evolution, 86(1):28–37, 2015.
[2] Takuma Sugi, Yukuo Nishida, and Ikue Mori. Regulation of behavioral plasticity by systemic temperature signaling in Caenorhabditis elegans. Nature Neuroscience, 14(8):984–992, 2011.
[3] Marc Levoy, Ren Ng, Andrew Adams, Matthew Footer, and Mark Horowitz. Light field microscopy. ACM Transactions on Graphics, 25(3):924, 2006.