本文利用相空间重构技术和混沌理论讨论了开都河日径流的混沌性质。通过日径流时间序列的功率谱分析,从定性角度讨论了日径流时间序列的混沌特征。进一步根据互信息量法得到相空间重构的延时,再根据Cao方法得到相空间重构的嵌入维数。利用Matlab软件计算得到相空间重构的延时和最佳嵌入维数分别为τ=6,m=14。这样将一维的开都河日径流时间序列重构成14维的相空间。通过最小数据量法计算出开都河日径流时间序列最大Lyapunov指数。利用最大Lyapunov指数对开都河日径流时间序列进行定量混沌分析。最后通过二阶Volterra自适应一步模型进行模拟。结果表明:开都河日径流时间序列的功率谱是连续的,功率谱呈现随频率增高而以指数方式递减趋势,区别于具有离散尖峰谱特征的周期时间序列和具有连续的、频率和振幅不相关谱特征的随机时间序列。这从定性角度表明开都河日径流时间序列具有混沌特征。通过计算得到开都河日径流时间序列的最大Lyapunov指数0〈λmax=0.0097〈1,从定量角度表明开都河日径流时间序列具有较弱的混沌特征。利用二阶Volterra自适应一步模型模拟得到相关系数和相对均方根误差分别0.9376和0.2390。这说明利用Volterra自适应模型模拟效果较好。
The runoff process is very complex. It is influenced by various environmental factors, such as precipitation, temperature, topographic conditions and the types of land utilization. In this paper, Chaos Theory, a mathematical sub-discipline that studies complex systems and the phase space reconstruction technology were used to investigate the chaotic characteristics of daily runoff process in Kaidu River Basin. By power spectrum analysis of the daily runoff time series, the chaos characteristics were discussed from qualitative perspective. The time delay of phase space reconstruction was calculated by the mutual information method and the optimal embedding dimension was determined by the Cao method. With Matlab software, the time delay and optimal embedding dimension were calculated as τ =6 and m=14 respectively. Thus, the one-dimensional daily runoff time series in the basin was successfully reconstructed into a multi-dimensional phase space. Furthermore, the chaos regarding the daily runoff time series in Kaidu River Basin was quantitatively analyzed by the maximum Lyapunov index. In the end, the simulation was implemented by the second order Volterra adaptive one- step model. The results indicate that the power spectrum of the daily runoff time series in Kaidu River Basin is continuous and it decreases exponentially with the increase of frequency. Therefore, from a qualitative perspective, it shows that the daily time series in Kaidu River Basin have chaotic characteristics. The maximum Lyapunov index of the daily time series is 0λ max =0.00971, also indicating chaotic characteristics from the quantitative view. The correlation coefficient and the relative root mean square errors (RRMSE) obtained by the simulation are 0.9376 and 0.2390 respectively, which suggests a satisfactory simulation result of the Volterra adaptive model. In addition, graphical comparison also indicates that the Volterra adaptive model can simulate the variation pattern of the daily runoff time series well. In detail, the proportion