本发明公开了一种混合系统的混沌状态检测方法及系统。该方法包括获取混合系统混合过程的视频；提取视频中的多帧图像对应的混合瞬态值；将多个混合瞬态值和多个混合瞬态值对应的时刻作为基本元素，生成时间序列；对时间序列进行去噪处理；采用改进的0‑1测试法，检测时间序列是否混沌。采用本发明的混沌状态检测方法及系统，能够在高噪音情况下精确地判断出混合物是否处于混沌状态，进而判断混合物是否混合均匀，快速精确地找到停止搅拌的时机。
1.一种混合系统的混沌状态检测方法，其特征在于，所述方法，包括：获取混合系统混合过程的视频；
提取所述视频中的多帧图像对应的混合瞬态值；所述混合瞬态值是根据所述混合系统的物性参数生成的；
将多个所述混合瞬态值和多个所述混合瞬态值对应的时刻作为基本元素，生成时间序列；
采用改进的0‑1测试法，检测所述时间序列是否混沌；所述改进的0‑1测试法是在0‑1测试法计算均方位移的公式中加入噪音干扰项；
所述采用改进的0‑1测试法，检测所述时间序列是否混沌，具体包括：计算所述时间序列的第一转换变量和第二转换变量；
根据所述第一转换变量、所述第二转换变量和所述噪音干扰项，计算所述时间序列的均方位移；
计算所述均方位移的增长率，得到所述时间序列对应的增长率序列；
判断所述增长率序列的中位值是否大于预设截止值，得到第二判断结果；
若所述第二判断结果为是，则所述时间序列处于混沌状态；
若所述第二判断结果为否，则所述时间序列处于非混沌状态；
所述计算所述时间序列的第一转换变量和第二转换变量，具体包括：采用如下公式计算所述第一转换变量和所述第二转换变量：式中，φ(j)表示降采样处理后的时间序列，j表示采样时间，j＝1,2,...,N，N表示降采样处理后的时间序列的长度，n表示数量， pc(n)表示第n个第一转换变量，qc(n)表示第n个第二转换变量，c表示抑制参数，c∈(0,2π)；
所述根据所述第一转换变量、所述第二转换变量和所述噪音干扰项，计算所述时间序列的均方位移，具体包括：
采用如下公式计算所述均方位移：式中，Mc(n)表示均方位移，ηn表示随机变量，ηn∈[‑0.5,0.5]，σ表示噪音值，σ∈[0,1]，pc(j)表示第j个第一转换变量，pc(j+n)表示第j+n个的第一转换变量，qc(j)表示第j个第二转换变量，qc(j+n)表示第j+n个第二转换变量。
2.根据权利要求1所述的混合系统的混沌状态检测方法，其特征在于，在所述对所述时间序列进行去噪处理，之后还包括：判断所述时间序列是否过采样，得到第一判断结果；
若所述第一判断结果为是，则对所述时间序列进行降采样处理；
若所述第一判断结果为否，则执行步骤“采用改进的0‑1测试法，检测所述时间序列是否混沌”。
3.根据权利要求1所述的混合系统的混沌状态检测方法，其特征在于，所述对所述时间序列进行去噪处理，具体包括：
采用平均滤波法或小波去噪法或施赖伯去噪法对所述时间序列进行去噪处理。
4.根据权利要求1所述的混合系统的混沌状态检测方法，其特征在于，所述计算所述均方位移的增长率，具体包括：
式中，Kc表示均方位移的增长率，ξ表示向量集合， ξ(j)表示向量集合ξ中的第j个数值，Δ(j)表示均方位移集合Δ中的第j个数值，Δ表示均方位移集合，cov(ξ,Δ)表示向量集合ξ和均方位移集合Δ的协方差，s表示向量集合ξ的长度，var(ξ)表示向量集合ξ的方差，var(Δ)表示均方位移集合Δ的方差，表示向量集合ξ的平均数， 表示均方位移集合Δ的平均数。
5.一种混合系统的混沌状态检测系统，其特征在于，所述系统，包括：视频获取模块，用于获取混合系统混合过程的视频；
混合瞬态值提取模块，用于提取所述视频中的多帧图像对应的混合瞬态值；所述混合瞬态值是根据所述混合系统的物性参数生成的；
时间序列生成模块，用于将多个所述混合瞬态值和多个所述混合瞬态值对应的时刻作为基本元素，生成时间序列；
混沌序列检测模块，用于采用改进的0‑1测试法，检测所述时间序列是否混沌；所述改进的0‑1测试法是在0‑1测试法计算均方位移的公式中加入噪音干扰项。
6.根据权利要求5所述的混合系统的混沌状态检测系统，其特征在于，所述系统，还包括：
过采样判断模块，用于判断所述时间序列是否过采样，得到第一判断结果；
若所述第一判断结果为是，则执行降采样子模块；若所述第一判断结果为否，则执行所述混沌序列检测模块；
所述混沌序列检测模块，具体包括：转换变量计算子模块，用于计算所述时间序列的第一转换变量和第二转换变量；
均方位移计算子模块，用于根据所述第一转换变量、所述第二转换变量和所述噪音干扰项，计算所述时间序列的均方位移；
增长率序列生成子模块，用于计算所述均方位移的增长率，得到所述时间序列对应的增长率序列；
混沌状态判断子模块，用于判断所述增长率序列的中位值是否大于预设截止值，得到第二判断结果；若所述第二判断结果为是，则执行混沌状态确定单元；若所述第二判断结果为否，则执行非混沌状态确定单元；
混沌状态确定单元，用于确定所述时间序列处于混沌状态；
非混沌状态确定单元，用于确定所述时间序列处于非混沌状态；
所述转换变量计算子模块，还用于采用如下公式计算所述第一转换变量和所述第二转换变量：
式中，φ(j)表示降采样处理后的时间序列，j表示采样时间，j＝1,2,...,N，N表示降采样处理后的时间序列的长度，n表示数量， pc(n)表示第n个第一转换变量，qc(n)表示第n个第二转换变量，c表示抑制参数，c∈(0,2π)；
所述均方位移计算子模块，用于采用如下公式计算所述均方位移：式中，Mc(n)表示均方位移，ηn表示随机变量，ηn∈[‑0.5,0.5]，σ表示噪音值，σ∈[0,1]，pc(j)表示第j个第一转换变量，pc(j+n)表示第j+n个的第一转换变量，qc(j)表示第j个第二转换变量，qc(j+n)表示第j+n个第二转换变量。
一种混合系统的混沌状态检测方法及系统
技术领域
[0001] 本发明涉及冶金化学混合技术领域，特别是涉及一种混合系统的混沌状态检测方法及系统。
背景技术
[0002] 混合是化工、冶金、材料、制药、造纸、建筑、生物和环境等众多领域中必不可少的操作单元之一，在物理反应过程中搅拌混合主要是为了获得反应物混合的均匀度，而在化
学反应参加的搅拌过程中则是为了优化反应进程。因此，多相搅拌混合效果的评价对搅拌
反应器的优化设计及反应过程的优化起着重要作用。描述和评价搅拌体系宏观混合效果的
方法众多，主要方法包括对搅拌功率、循环流量、排出流量及混匀时间等参数的测量。在冶
金方面，多相混合物(即混合系统)处于混沌状态，则代表混合物混合均匀，需停止对混合系
统的搅拌，其中，电导率法、光学法、热测法等方法是常用的测量混匀时间的技术。但是，上
述方法只能量化两相流混合瞬态的混合效果，若要进一步分析多相流混合瞬态随时间演化
趋势以发掘随机过程中的规律，则需要表征时间序列的混沌特性。
[0003] 一般地，混沌识别方法有庞加莱映像法、相图法、李雅普诺夫(Lyapunov)Lyapunov指数法和柯尔莫哥洛夫(Kolmogorov)熵等方法，但这些混沌识别方法受噪声等因素的干
扰，存在识别精度不高的问题，这导致实验或生产中，技术人员不能准确识别停止搅拌的时
机，造成了混合不均匀或者资源的浪费。因此，削弱噪音干扰对混合物状态识别的影响，准
确快速地找到停止搅拌的时机，具有理论和实用的双重价值。
发明内容
[0004] 本发明的目的是提供一种混合系统的混沌状态检测方法及系统，能够在高噪音情况下，精确地判断出混合物开始处于混沌状态的时间。
[0005] 为实现上述目的，本发明提供了如下方案：
[0006] 一种混合系统的混沌状态检测方法，包括：
[0007] 获取混合系统混合过程的视频；
[0008] 提取所述视频中的多帧图像对应的混合瞬态值；所述混合瞬态值是根据所述混合系统的物性参数生成的；
[0009] 将多个所述混合瞬态值和多个所述混合瞬态值对应的时刻作为基本元素，生成时间序列；
[0010] 对所述时间序列进行去噪处理；
[0011] 采用改进的0‑1测试法，检测所述时间序列是否混沌；所述改进的0‑1测试法是在0‑1测试法计算均方位移的公式中加入噪音干扰项。
[0012] 可选的，在所述对所述时间序列进行去噪处理，之后还包括：
[0013] 判断所述时间序列是否过采样，得到第一判断结果；
[0014] 若所述第一判断结果为是，则对所述时间序列进行降采样处理；
[0015] 若所述第一判断结果为否，则执行步骤“采用改进的0‑1测试法，检测所述时间序列是否混沌”。
[0016] 可选的，所述对所述时间序列进行去噪处理，具体包括：
[0017] 采用平均滤波法或小波去噪法或施赖伯去噪法对所述时间序列进行去噪处理。
[0018] 可选的，所述采用改进的0‑1测试法，检测所述时间序列是否混沌，具体包括：
[0019] 计算所述时间序列的第一转换变量和第二转换变量；
[0020] 根据所述第一转换变量、所述第二转换变量和所述噪音干扰项，计算所述时间序列的均方位移；
[0021] 计算所述均方位移的增长率，得到所述时间序列对应的增长率序列；
[0022] 判断所述增长率序列的中位值是否大于预设截止值，得到第二判断结果；
[0023] 若所述第二判断结果为是，则所述时间序列处于混沌状态；
[0024] 若所述第二判断结果为否，则所述时间序列处于非混沌状态。
[0025] 可选的，所述计算所述时间序列的第一转换变量和第二转换变量，具体包括：
[0026] 采用如下公式计算所述第一转换变量和所述第二转换变量：
[0027]
[0028]
[0029] 式中，φ(j)表示降采样处理后的时间序列，j表示采样时间，j＝1,2,...,N，N表示降采样处理后的时间序列的长度，n表示数量， pc(n)表示第n个第一转换变
量，qc(n)表示第n个第二转换变量，c表示抑制参数，c∈(0,2π)。
[0030] 可选的，所述根据所述第一转换变量、所述第二转换变量和所述噪音干扰项，计算所述时间序列的均方位移，具体包括：
[0031] 采用如下公式计算所述均方位移：
[0032]
[0033] 式中，Mc(n)表示均方位移，ηn表示随机变量，ηn∈[‑0.5,0.5]，σ表示噪音值，σ∈[0,1]，pc(j)表示第j个第一转换变量，pc(j+n)表示第j+n个的第一转换变量，qc(j)表示第j
个第二转换变量，qc(j+n)表示第j+n个第二转换变量。
[0034] 可选的，所述计算所述均方位移的增长率，具体包括：
[0035] 采用如下公式计算所述均方位移的增长率：
[0036]
[0037] 其中，
[0038]
[0039]
[0040]
[0041]
[0042] 式中，Kc表示均方位移的增长率，ξ表示向量集合， ξ(j)表示向量集合ξ中的第j个数值，Δ(j)表示均方位移集合Δ中的第j个数值，Δ表示均方位移集合，
cov(ξ,Δ)表示向量集合ξ和均方位移集合Δ的协方
差，s表示向量集合ξ的长度，var(ξ)表示向量集合ξ的方差，var(Δ)表示均方位移集合Δ的
方差，表示向量集合ξ的平均数，表示均方位移集合Δ的平均数。
[0043] 一种混合系统的混沌状态检测系统，包括：
[0044] 视频获取模块，用于获取混合系统混合过程的视频；
[0045] 混合瞬态值提取模块，用于提取所述视频中的多帧图像对应的混合瞬态值；所述混合瞬态值是根据所述混合系统的物性参数生成的；
[0046] 时间序列生成模块，用于将多个所述混合瞬态值和多个所述混合瞬态值对应的时刻作为基本元素，生成时间序列；
[0047] 去噪模块，用于对所述时间序列进行去噪处理；
[0048] 混沌序列检测模块，用于采用改进的0‑1测试法，检测所述时间序列是否混沌；所述改进的0‑1测试法是在0‑1测试法计算均方位移的公式中加入噪音干扰项。
[0049] 可选的，所述系统，还包括：
[0050] 过采样判断模块，用于判断所述时间序列是否过采样，得到第一判断结果；
[0051] 若所述第一判断结果为是，则执行降采样子模块；若所述第一判断结果为否，则执行所述混沌序列检测模块；
[0052] 降采样子模块，用于对所述时间序列进行降采样处理。
[0053] 可选的，所述混沌序列检测模块，具体包括：
[0054] 转换变量计算子模块，用于计算所述时间序列的第一转换变量和第二转换变量；
[0055] 均方位移计算子模块，用于根据所述第一转换变量、所述第二转换变量和所述噪音干扰项，计算所述时间序列的均方位移；
[0056] 增长率序列生成子模块，用于计算所述均方位移的增长率，得到所述时间序列对应的增长率序列；
[0057] 混沌状态判断子模块，用于判断所述增长率序列的中位值是否大于预设截止值，得到第二判断结果；若所述第二判断结果为是，则执行混沌状态确定单元；若所述第二判断
结果为否，则执行非混沌状态确定单元；
[0058] 混沌状态确定单元，用于确定所述时间序列处于混沌状态；
[0059] 非混沌状态确定单元，用于确定所述时间序列处于非混沌状态。
[0060] 与现有技术相比，本发明的有益效果是：
[0061] 本发明提出了一种混合系统的混沌状态检测方法及系统，方法包括获取混合系统混合过程的视频；提取视频中的多帧图像对应的混合瞬态值；将多个混合瞬态值和多个混
合瞬态值对应的时刻作为基本元素，生成时间序列；对时间序列进行去噪处理；采用改进的
0‑1测试法，检测时间序列是否混沌；改进的0‑1测试法是在0‑1测试法计算均方位移的公式
中加入噪音干扰项。本发明通过在0‑1测试检测中加入噪音值参数，能够在高噪音情况下精
确地判断出混合系统(即混合物)是否处于混沌状态，进而判断混合物是否混合均匀，快速
精确地找到停止搅拌的时机。
附图说明
[0062] 为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施
例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图
获得其他的附图。
[0063] 图1为本发明实施例中所提供的混合系统的混沌状态检测方法的流程图；
[0064] 图2为本发明实施例中所提供的混合系统的混沌状态检测系统的结构示意图；
[0065] 图3为本发明实施例中所提供的顶吹气液两相流搅拌实验监测图。
具体实施方式
[0066] 下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于
本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他
实施例，都属于本发明保护的范围。
[0067] 本发明的目的是提供一种混合系统的混沌状态检测方法及系统，能够在高噪音情况下，精确地判断出混合物开始处于混沌状态的时间。
[0068] 为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本发明作进一步详细的说明。
[0069] 实施例
[0070] 图1为本发明实施例中所提供的混合系统的混沌状态检测方法的流程图；如图1所示，本发明提供了一种混合系统的混沌状态检测方法，包括：
[0071] 步骤101：获取混合系统混合过程的视频。
[0072] 步骤102：提取视频中的多帧图像对应的混合瞬态值。混合瞬态值是根据混合系统的物性参数生成的。
[0073] 步骤103：将多个混合瞬态值和多个混合瞬态值对应的时刻作为基本元素，生成时间序列。
[0074] 步骤104：对时间序列进行去噪处理。
[0075] 步骤104具体包括：采用平均滤波法或小波去噪法或施赖伯去噪法对时间序列进行去噪处理，对时间序列进行去噪处理时优先使用施赖伯去噪法。
[0076] 步骤104之后还包括：
[0077] 判断时间序列是否过采样，得到第一判断结果；若第一判断结果为是，则对时间序列进行降采样处理；若第一判断结果为否，则执行步骤“采用改进的0‑1测试法，检测时间序
列是否混沌”。
[0078] 具体地，对时间序列进行降采样处理的方法有两种，第一种降采样处理的方法是隔固定频率抽取一个信号的向下采样法，第二种降采样处理的方法截取固定长度的局部信
号，取最大值和最小值的向下采样法；本发明对时间序列进行降采样处理时优先使用第二
种降采样处理的方法。
[0079] 步骤105：采用改进的0‑1测试法，检测时间序列是否混沌。改进的0‑1测试法是在0‑1测试法计算均方位移的公式中加入噪音干扰项。检测出当前时间序列为混沌时间序列，
证明混合物处于混沌状态，此时为停止搅拌的时机。
[0080] 步骤105具体包括：
[0081] 步骤1051：计算时间序列的第一转换变量和第二转换变量。
[0082] 步骤1051具体包括：
[0083] 采用如下公式计算第一转换变量和第二转换变量：
[0084]
[0085]
[0086] 式中，φ(j)表示降采样处理后的时间序列，j表示采样时间，j＝1,2,...,N，N表示降采样处理后的时间序列的长度，n表示数量， pc(n)表示第n个第一转换变
量，qc(n)表示第n个第二转换变量，c表示抑制参数，c∈(0,2π)。
[0087] 步骤1052：根据第一转换变量、第二转换变量和噪音干扰项，计算时间序列的均方位移。
[0088] 步骤1052具体包括：
[0089] 采用如下公式计算均方位移：
[0090]
[0091] 式中，Mc(n)表示第n个均方位移，ηn表示随机变量，ηn∈[‑0.5,0.5]，σ表示噪音值，σ∈[0,1]，pc(j)表示第j个第一转换变量，pc(j+n)表示第j+n个的第一转换变量，qc(j)表示
第j个第二转换变量，qc(j+n)表示第j+n个第二转换变量。
[0092] 步骤1053：计算均方位移的增长率，得到时间序列对应的增长率序列。
[0093] 步骤1053采用如下公式计算均方位移的增长率：
[0094]
[0095] 其中，
[0096]
[0097]
[0098]
[0099]
[0100] 式中，Kc表示均方位移的增长率，ξ表示向量集合， ξ(j)表示向量集合ξ中的第j个数值，Δ(j)表示均方位移集合Δ中的第j个数值，Δ表示均方位移集合，
cov(ξ,Δ)表示向量集合ξ和均方位移集合Δ的协方
差，s表示向量集合ξ的长度，var(ξ)表示向量集合ξ的方差，var(Δ)表示均方位移集合Δ的
方差，表示向量集合ξ的平均数，表示均方位移集合Δ的平均数。
[0101] 步骤1054：判断增长率序列的中位值是否大于预设截止值，得到第二判断结果；若第二判断结果为是，则时间序列处于混沌状态；若第二判断结果为否，则时间序列处于非混
沌状态。
[0102] 图2为本发明实施例中所提供的混合系统的混沌状态检测系统的结构示意图；如图2所示，本发明提供了一种混合系统的混沌状态检测系统，包括：
[0103] 视频获取模块201，用于获取混合系统混合过程的视频。
[0104] 混合瞬态值提取模块202，用于提取视频中的多帧图像对应的混合瞬态值；混合瞬态值是根据混合系统的物性参数生成的。
[0105] 时间序列生成模块203，用于将多个混合瞬态值和多个混合瞬态值对应的时刻作为基本元素，生成时间序列。
[0106] 去噪模块204，用于对时间序列进行去噪处理。去噪模块204采用平均滤波法或小波去噪法或施赖伯去噪法对时间序列进行去噪处理。
[0107] 混沌序列检测模块205，用于采用改进的0‑1测试法，检测时间序列是否混沌；改进的0‑1测试法是在0‑1测试法计算均方位移的公式中加入噪音干扰项。
[0108] 混沌序列检测模块205，具体包括：
[0109] 转换变量计算子模块，用于计算时间序列的第一转换变量和第二转换变量。
[0110] 转换变量计算子模块，具体包括：
[0111] 采用如下公式计算第一转换变量和第二转换变量：
[0112]
[0113]
[0114] 式中，φ(j)表示降采样处理后的时间序列，j表示采样时间，j＝1,2,...,N，N表示降采样处理后的时间序列的长度，n表示数量， pc(n)表示第n个第一转换变
量，qc(n)表示第n个第二转换变量，c表示抑制参数，c∈(0,2π)。
[0115] 均方位移计算子模块，用于根据第一转换变量、第二转换变量和噪音干扰项，计算时间序列的均方位移。
[0116] 均方位移计算子模块，具体包括：
[0117] 采用如下公式计算均方位移：
[0118]
[0119] 式中，Mc(n)表示第n个均方位移，ηn表示随机变量，ηn∈[‑0.5,0.5]，σ表示噪音值，σ∈[0,1]，pc(j)表示第j个第一转换变量，pc(j+n)表示第j+n个的第一转换变量，qc(j)表示
第j个第二转换变量，qc(j+n)表示第j+n个第二转换变量。
[0120] 增长率序列生成子模块，用于计算均方位移的增长率，得到时间序列对应的增长率序列。
[0121] 增长率序列生成子模块，具体包括：
[0122] 采用如下公式计算均方位移的增长率：
[0123]
[0124] 其中，
[0125]
[0126]
[0127]
[0128]
[0129] 式中，Kc表示均方位移的增长率，ξ表示向量集合， ξ(j)表示向量集合ξ中的第j个数值，Δ(j)表示均方位移集合Δ中的第j个数值，Δ表示均方位移集合，
cov(ξ,Δ)表示向量集合ξ和均方位移集合Δ的协方
差，s表示向量集合ξ的长度，var(ξ)表示向量集合ξ的方差，var(Δ)表示均方位移集合Δ的
方差，表示向量集合ξ的平均数，表示均方位移集合Δ的平均数。
[0130] 混沌状态判断子模块，用于判断增长率序列的中位值是否大于预设截止值，得到第二判断结果；若第二判断结果为是，则执行混沌状态确定单元；若第二判断结果为否，则
执行非混沌状态确定单元；
[0131] 混沌状态确定单元，用于确定时间序列处于混沌状态；
[0132] 非混沌状态确定单元，用于确定时间序列处于非混沌状态。
[0133] 本发明提供的混合系统的混沌状态检测系统，还包括过采样判断模块，用于判断时间序列是否过采样，得到第一判断结果；若第一判断结果为是，则执行降采样子模块；若
第一判断结果为否，则执行混沌序列检测模块；降采样子模块，用于对时间序列进行降采样
处理。
[0134] 具体地，图3为本发明实施例中所提供的顶吹气液两相流搅拌实验监测图如图3所示，在顶吹气液两相流搅拌实验过程中，按不同的埋管深度和不同的气体流量，设置了五组
工况，进行混合物混沌状态判断，从高速摄像机拍摄的混合物视频中提取多帧图像对应的
混合瞬态值；将多个混合瞬态值和多个混合瞬态值对应的时刻作为基本元素，生成时间序
列，对时间序列进行分析，去噪算法选用施赖伯去噪法，降采样方法选用取信号局部极小值
和极大值的向下采样方法，在输入5组工况的贝蒂数时间序列(即将多个混合瞬态值和多个
混合瞬态值对应的时刻作为基本元素，生成的时间序列)后混沌决策树代码运行结果如表1
所示，表1中L1‑L5表示不同的工况，百分数表示不同的噪音值，不同工况和噪音值对应的数
值为增长率。具体的，L1的埋管深度为5cm，气体流量为2000L/h；L2的埋管深度为6cm，气体
流量为1000L/h；L3的埋管深度为6cm，气体流量为2000L/h；L4的埋管深度为7cm，气体流量
为1000L/h；L5的埋管深度为7cm，气体流量为2000L/h。
[0135] 表1不同噪音值与增长率值的对应关系表
[0136]
[0137] 使用增长率值判断系统是否混沌，当增长率值大于或者等于截断值(截止值取0.85)时，当前时间序列为混沌时间序列，此时气液两相流混合物处于混沌状态，此时停止
搅拌；当增长率值小于截断值时，当前时间序列为非混沌时间序列，此时气液两相流混合物
处于非混沌状态，此时继续搅拌直至检测气液两相流混合物处于混沌状态。
[0138] 本发明提供的混合系统的混沌状态检测方法及系统，除了时间序列得提取之外，该方法不需要来自用户的任何输入，另外，本发明在0‑1测试中加入噪音项，对具有高噪音
的混合系统混沌状态的判别具有高度的准确性。
[0139] 本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。对于实施例公开的系统
而言，由于其与实施例公开的方法相对应，所以描述的比较简单，相关之处参见方法部分说
明即可。
[0140] 本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据
本发明的思想，在具体实施方式及应用范围上均会有改变之处。综上，本说明书内容不应理
解为对本发明的限制。
