CN104463382B - 气候变化预测以及趋势分析方法 - Google Patents
气候变化预测以及趋势分析方法 Download PDFInfo
- Publication number
- CN104463382B CN104463382B CN201510003459.1A CN201510003459A CN104463382B CN 104463382 B CN104463382 B CN 104463382B CN 201510003459 A CN201510003459 A CN 201510003459A CN 104463382 B CN104463382 B CN 104463382B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msubsup
- msub
- msup
- mfrac
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Strategic Management (AREA)
- Development Economics (AREA)
- Economics (AREA)
- Entrepreneurship & Innovation (AREA)
- Educational Administration (AREA)
- Operations Research (AREA)
- Marketing (AREA)
- Game Theory and Decision Science (AREA)
- Quality & Reliability (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Information Retrieval, Db Structures And Fs Structures Therefor (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
Abstract
本发明实施例提供一种气候变化预测以及趋势分析方法。本发明方法,包括:处理器接收气候变量数据以及与所述气候变量数据对应的时间值;所述处理器采用模糊聚类算法计算所述气候变量数据的渐进变化与跳跃变化;所述处理器根据所述渐进变化与所述跳跃变化对气候变化进行预测以及趋势分析。本发明实施例准确体现气候变化中跳跃变化,有利于进一步准确预测资源的合理利用及区域可持续发展,有助于分析对应该跳跃变化的成因。
Description
技术领域
本发明实施例涉及计算机计算领域,尤其涉及一种气候变化预测以及趋势分析方法。
背景技术
近年来,洪涝、干旱和冰雨等极端气候引起了全世界各地人们的关注。造成该现象的原因有很多。例如,冷暖气流的急剧变化等自然因素,同时,人类对森林的滥砍滥伐、汽车的尾气排放等人类活动的人为因素也占很大的比例。
现有技术对此采用混沌-人工神经网络方法对气候变化进行预测以及趋势分析,该方法是结合混沌理论和人工神经网络,互相抵消两者之间的不足之处。分析流域气温、降雨和径流量的变化过程。进而对于气候的变化预测以及趋势分析。
该种方法对于气候的变化体现了气候变化的整体趋势,但是无法准确地体现气候跳跃变化。
发明内容
本发明实施例提供一种气候变化预测以及趋势分析方法,以克服现有中无法准确地对于气候变化中的跳跃变化作出预测以及趋势分析的问题。
本实施例提供了一种气候变化预测以及趋势分析方法,包括:
处理器接收气候变量数据以及与所述气候变量数据对应的时间值;
所述处理器采用模糊聚类算法计算所述气候变量数据的渐进变化与跳跃变化;
所述处理器根据所述渐进变化与所述跳跃变化对气候变化进行预测以及趋势分析。
进一步地,所述处理器采用模糊聚类算法计算所述气候变量数据的渐进变化与跳跃变化,包括:
步骤一、处理器初始化所述模糊聚类算法参量,设置所述模糊聚类算法终止阈值、初始分割阶数和最小分割阶数;
步骤二、所述处理器求得隶属度、混合概率,并判断所述混合概率是否大于0,若是,则直接更新聚类参数以及所述隶属度,若否,则降低分割阶数后更新聚类参数以及所述隶属度;
步骤三、所述处理器判断相邻两次迭代过程的模糊划分矩阵的相对误差是否大于终止阈值,若否,则返回步骤二,若是,则根据最小信息长度准则选择最优分割阶数;
步骤四、所述处理器将所述最优分割阶数与最优聚类参数代入模糊聚类算法得出所述气候变量数据的渐进变化与跳跃变化。
进一步地,所述处理器初始化所述模糊聚类算法参量,包括:
处理器采用分段常函数得到在时间坐标集合上kmax个初始模糊聚类中心为
其中,N为气候变量数据的个数,m为模糊加权指数且m=2,ti为气候变量数据对应的时间,p=int(N/kmax),当k=kmax,令kp=N;
处理器将时间坐标集合上的初始偏差定义为
其中,N为气候变量数据的个数,m为模糊加权指数且m=2,ti为气候变量数据对应的时间,μk为在时间坐标集合上的初始模糊聚类中心,sk(ti)为如上所述分段常函数;
所述处理器从气候变量数据集合中随机选取kmax个数据作为初始中心,得到
vk∈{xi,1≤i≤N},1≤k≤kmax (3)
其中,N为气候变量数据的个数,xi为气候变量数据;
所述处理器令所述气候变量数据集合中的初始协方差与单位矩阵成正比,得到
其中,为kmax阶单位矩阵, 为气候变量数据的平均值;
所述处理器令初始混合概率αk与所述kmax成反比,得到
αk=1/kmax,1≤k≤kmax (5)
其中,kmax为初始聚类数。
进一步地,所述处理器将所述最优分割阶数代入模糊聚类算法得出所述气候变量数据的渐进变化与跳跃变化,包括:
处理器采用高斯函数得到在时间坐标集合上kmax个模糊分割βk(ti),公式为
其中ti为气候变量数据对应的时间,μk为在时间坐标集合上的初始模糊聚类中心,μk为时间坐标集合上的初始偏差;
进一步地,所述处理器求得隶属度、混合概率,包括:
采用公式
计算隶属度,其中,所述为第k个聚类的混合概率,p(t|ηk)为在时间坐标集合上的第k个高斯函数,p(x|ηk)为在气候变量数据集合上第k个高斯函数,m为模糊加权指数且m=2;
采用公式
计算混合概率,其中,所述ui,k为第i个数据属于第k个聚类的隶属度,c为每个聚类成员所含参数的维数;
所述则更新聚类参数以及所述隶属度,包括:
处理器采用期望值最大算法将所述聚类参数以及隶属度更新为
其中,所述ui,k为第i个数据属于第k个聚类的隶属度,xi为气候变量数据,ti为气候变量数据xi对应的时间,μk与分别为在时间坐标集合上的第k个高斯函数p(t|ηk)的中心和偏差,vk与Σk分别为在气候变量数据集合上第k个高斯函数p(x|ηk)的中心和方差,m为模糊加权指数且m=2。
进一步地,所述根据最小信息长度准则选择最优分割阶数,包括:
处理器根据最小信息长度准则,采用公式
求得最小信息长度,其中,所述c为每个聚类成员所含参数的个数,αk为第k个聚类的混合概率,N为气候变量数据的个数,knz为当前聚类个数,p(t|ηk)为在时间坐标集合上的第k个高斯函数,p(x|ηk)为在气候变量数据集合上第k个高斯函数。
进一步地,所述步骤三之后,还包括:
处理器判断若相邻两次迭代过程的模糊划分矩阵的相对误差低于终止阈值时,删除具有最低概率的聚类成员并返回所述步骤二,查找更小的信息长度值。
本发明实施例处理器接收气候变量数据以及与所述气候变量数据对应的时间值,采用模糊聚类算法计算所述气候变量数据的渐进变化与跳跃变化,并根据所述渐进变化与所述跳跃变化对齐气候变化进行预测以及趋势分析,准确地体现了气候变化中的跳跃变化,有利于进一步准确预测资源的合理利用及区域可持续发展,有助于分析对应该跳跃变化的成因。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作一简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明气候变化预测以及趋势分析方法流程图;
图2为本发明模糊聚类算法的流程图;
图3为本发明仿真数据示意图;
图4为本发明仿真数据的初始模糊分割示意图;
图5为本发明仿真数据的最小信息长度准则的演化示意图;
图6为本发明仿真数据的最优模糊分割结果示意图;
图7为本发明仿真数据的最优隶属度结果示意图;
图8为本发明仿真数据的最优明确分割结果示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
图1为本发明气候变化预测以及趋势分析方法流程图,如图1所示,本实施例方法,包括:
步骤101、处理器接收气候变量数据以及与所述气候变量数据对应的时间值;
步骤102、所述处理器采用模糊聚类算法计算所述气候变量数据的渐进变化与跳跃变化;
步骤103、所述处理器根据所述渐进变化与所述跳跃变化对气候变化进行预测以及趋势分析。
具体来说,处理器接收气候变量数据,同时还接收与该气候变量数据相对应的时间值。例如,气候变化数据为花园口水温站提供的径流数据,并且该径流数据还对应标明年份的时间值。处理器根据该气候变量数据以及时间值,采用模糊聚类算法计算气候变量数据的渐进变化与跳跃变化。从而预测气候变化与分析气候变化的趋势。
假设数据可以被具有参数η={η1,…,ηK}的多元(包括时间坐标)高斯混合分布有效建模,其中xi=[x1,i,x2,i,…,xd,i]T(1≤i≤N),ti是数据xi对应的时间点,K是聚类的数目,ηk是第k(1≤k≤K)个聚类参数构成的向量。
因为时间变量ti独立于变量xi,所以数据zi到第k个聚类原型的距离测度D2(zi,ηk)由下式表示:
其中ηk表示第k个聚类原型的参数,表示第k个聚类的混合概率。p(ti|ηk)表示在时间坐标集合上时间ti属于第k个聚类的概率,且
其中μk与分别表示在时间坐标集合上的第k个高斯函数p(t|ηk)的中心和方差,p(xi|ηk)表示在数据集合上数据xi属于第k个聚类的概率,且
其中vk与Σk分别表示在数据集合上第k个高斯函数p(x|ηk)的中心和方差。第k个聚类的原型参数K个聚类的所有参数的集合定义为η≡{η1,…,ηK}。此时,概率密度函数定义为:
在时间坐标集合上时间点ti属于第k个聚类的隶属度βk(ti)∈[0,1]表示时间点ti属于第k个时间序列片段区间的程度,即
本发明实施例处理器接收气候变量数据以及与所述气候变量数据对应的时间值,采用模糊聚类算法计算所述气候变量数据的渐进变化与跳跃变化,并根据所述渐进变化与所述跳跃变化对齐气候变化进行预测以及趋势分析,准确地体现了气候变化中的跳跃变化,有利于进一步准确预测资源的合理利用及区域可持续发展,有助于分析对应该跳跃变化的成因。
图2为本发明模糊聚类算法的流程图,如图2所示,本实施例方法,包括:
步骤201、处理器初始化所述模糊聚类算法参量,设置所述模糊聚类算法终止阈值、初始分割阶数和最小分割阶数;
具体来说,处理器采用分段常函数得到在时间坐标集合上kmax个初始模糊聚类中心为
其中,N为气候变量数据的个数,m为模糊加权指数且m=2,ti为气候变量数据对应的时间,p=int(N/kmax),当k=kmax,令kp=N;
处理器将时间坐标集合上的初始偏差定义为
其中,N为气候变量数据的个数,m为模糊加权指数且m=2,ti为气候变量数据对应的时间,μk为在时间坐标集合上的初始模糊聚类中心,sk(ti)为如上所述分段常函数;
所述处理器从气候变量数据集合中随机选取kmax个数据作为初始中心,得到
vk∈{xi,1≤i≤N},1≤k≤kmax (8)
其中,N为气候变量数据的个数,xi为气候变量数据;
所述处理器令所述气候变量数据集合中的初始协方差与单位矩阵成正比,得到
其中,为kmax阶单位矩阵, 为气候变量数据的平均值;
所述处理器令初始混合概率αk与所述kmax成反比,得到
αk=1/kmax,1≤k≤kmax (10)
其中,kmax为初始聚类数。
步骤202、所述处理器求得隶属度、混合概率,并判断所述混合概率是否大于0,若是,则直接更新聚类参数以及所述隶属度,若否,则降低分割阶数后更新聚类参数以及所述隶属度;
具体来说,采用公式
计算隶属度,其中,所述为第k个聚类的混合概率,p(t|ηk)为在时间坐标集合上的第k个高斯函数,p(x|ηk)为在气候变量数据集合上第k个高斯函数,m为模糊加权指数且m=2;
采用公式
计算混合概率,其中,所述ui,k为第i个数据属于第k个聚类的隶属度,c为每个聚类成员所含参数的维数;
所述则更新聚类参数以及所述隶属度,包括:
处理器采用期望值最大算法将所述聚类参数以及隶属度更新为
其中,所述ui,k为第i个数据属于第k个聚类的隶属度,xi为气候变量数据,ti为气候变量数据xi对应的时间,μk与分别为在时间坐标集合上的第k个高斯函数p(t|ηk)的中心和偏差,vk与Σk分别为在气候变量数据集合上第k个高斯函数p(x|ηk)的中心和方差,m为模糊加权指数且m=2。
步骤203、所述处理器判断相邻两次迭代过程的模糊划分矩阵的相对误差是否大于终止阈值,若否,则返回步骤二,若是,则根据最小信息长度准则选择最优分割阶数;
步骤204、所述处理器将所述最优分割阶数与最优聚类参数代入模糊聚类算法得出所述气候变量数据的渐进变化与跳跃变化。
具体来说,基于模糊聚类的最小信息长度准则:
且
其中c是每个聚类成员所含参数的个数(即ηk的维数),N是数据集合Z所含数据的总数。当αk等于零,MML的惩罚项将趋于-∞。通过只编码概率不为零的聚类成员参数来解决这个问题。令表示具有非零概率的聚类成员的数目,则
所以,基于模糊聚类的MML准则为
模糊聚类删减策略具体为:
EM算法最小化代价函数有下面的M步骤:
对于k:αk>0,更新参数可得
其中ui,k在E步骤的更新公式如下
在删除运算中采用针对混合模型的分量形式的EM算法(CEM2)。这个算法并不同时更新所有的聚类成员参数,而是序贯更新聚类成员的参数:更新α1和η1,计算
更新α2和η2,重新计算等等。当删除一个聚类成员时(αk(l+1)=0),CEM2算法重新分配概率密度来增加其余聚类成员的生存机会。这样初始聚类的数目可以任意大,而不产生任何问题。CEM2算法的计算负担仅稍重于EM算法。
图3为本发明仿真数据示意图,图3所示为黄河流域花园口水文站1919-2004年年径流数据。
图4为本发明仿真数据的初始模糊分割示意图。
具体来说,处理器采用分段常函数得到在时间坐标集合上kmax个初始模糊聚类中心为
其中,N为气候变量数据的个数,m为模糊加权指数且m=2,ti为气候变量数据对应的时间,p=int(N/kmax),当k=kmax,令kp=N;
处理器将时间坐标集合上的初始偏差定义为
其中,N为气候变量数据的个数,m为模糊加权指数且m=2,ti为气候变量数据对应的时间,μk为在时间坐标集合上的初始模糊聚类中心,sk(ti)为如上所述分段常函数;
所述处理器从气候变量数据集合中随机选取kmax个数据作为初始中心,得到
vk∈{xi,1≤i≤N},1≤k≤kmax (30)
其中,N为气候变量数据的个数,xi为气候变量数据;
所述处理器令所述气候变量数据集合中的初始协方差与单位矩阵成正比,得到
其中,为kmax阶单位矩阵, 为气候变量数据的平均值;
所述处理器令初始混合概率αk与所述kmax成反比,得到
αk=1/kmax,1≤k≤kmax (32)
其中,kmax为初始聚类数。
进一步地,所述处理器将所述初始模糊聚类参数代入模糊聚类算法得出所述气候变量数据的初始模糊分割,包括:
处理器采用高斯函数得到在时间坐标集合上kmax个模糊分割βk(ti),公式为
其中ti为气候变量数据对应的时间,μk为在时间坐标集合上的初始模糊聚类中心,为时间坐标集合上的初始偏差。
图5为本发明仿真数据的最小信息长度准则的演化示意图。
具体来说,基于模糊聚类的最小信息长度准则:
且
其中c是每个聚类成员所含参数的个数(即ηk的维数),N是数据集合Z所含数据的总数。当αk等于零,MML的惩罚项将趋于-∞。通过只编码概率不为零的聚类成员参数来解决这个问题。令表示具有非零概率的聚类成员的数目,则
所以,基于模糊聚类的MML准则为
图6为本发明仿真数据的最优模糊分割结果示意图,如图6所示,本实施例方法将所述最优分割阶数与最优聚类参数代入模糊聚类算法得出最优模糊分割结果。
具体来说,处理器采用高斯函数得到在时间坐标集合上kopt个模糊分割βk(ti),公式为
其中ti为气候变量数据对应的时间,μk为在时间坐标集合上的初始模糊聚类中心,为时间坐标集合上的初始偏差。
图7为本发明仿真数据的最优隶属度结果示意图,如图7所示,本实施例方法将所述最优分割阶数与最优聚类参数代入模糊聚类算法得出最优隶属度结果。
具体来说,采用公式
计算隶属度,其中,所述为第k个聚类的混合概率,p(t|ηk)为在时间坐标集合上的第k个高斯函数,p(x|ηk)为在气候变量数据集合上第k个高斯函数,m为模糊加权指数且m=2;
图8为本发明仿真数据的最优明确分割结果示意图,如图8所示,本实施例方法将所述最优分割阶数与最优聚类参数代入模糊聚类算法得出最优明确分割结果。
具体来说,采用公式
计算最优明确分割结果。
本实施例中,时间序列模糊分割算法的具体为:
(1)输入数据集合设置终止阈值ε,初始分割阶数kmax,最小分割阶数kmin。
(2)初始聚类参数其中其中μk与分别表示在时间坐标集合上的第k个高斯函数p(t|ηk)的中心和方差,vk与Σk分别表示在数据集合上第k个高斯函数p(x|ηk)的中心和方差,
(3)计算隶属度
(4)更新混合概率
(5)如果更新聚类参数和隶属度
(6)如果降低当前分割阶数knz,令knz=knz-1。如果knz<kmin,算法结束;令k=k+1,如果k<kmax,回到步骤(3)。
(7)计算最小信息长度准则
(8)如果||U(l+1)-U(l)||≥ε,令l=l+1,回到第(3)步。
(9)如果MessLen(η(l+1),Z)≤MessLenmin,令
MessLenmin=MessLen(η(l+1),Z) (52)
ηopt=η(l+1) (53)
(10)令k*=argmink{αk>0},knz=knz-1。如果knz≥kmin,回到步骤(3)。
本发明实施例处理器接收气候变量数据以及与所述气候变量数据对应的时间值,采用模糊聚类算法计算所述气候变量数据的渐进变化与跳跃变化,并根据所述渐进变化与所述跳跃变化对齐气候变化进行预测以及趋势分析,准确地体现了气候变化中的跳跃变化,有利于进一步准确预测资源的合理利用及区域可持续发展,有助于分析对应该跳跃变化的成因。
最后应说明的是:以上各实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述各实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分或者全部技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的范围。
Claims (2)
1.一种气候变化预测以及趋势分析方法,其特征在于,包括:
S1处理器接收气候变量数据以及与所述气候变量数据对应的时间值;
S2所述处理器采用模糊聚类算法计算所述气候变量数据的渐进变化与跳跃变化,其步骤包括:
步骤一、处理器初始化所述模糊聚类算法参量,设置所述模糊聚类算法终止阈值、初始分割阶数和最小分割阶数,所述处理器初始化所述模糊聚类算法参量步骤包括:
处理器采用分段常函数得到在时间坐标集合上kmax个初始模糊聚类中心为
<mrow>
<msub>
<mi>&mu;</mi>
<mi>k</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</msubsup>
<msup>
<mrow>
<mo>&lsqb;</mo>
<msub>
<mi>s</mi>
<mi>k</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>t</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mi>m</mi>
</msup>
<msub>
<mi>t</mi>
<mi>i</mi>
</msub>
</mrow>
<mrow>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</msubsup>
<msup>
<mrow>
<mo>&lsqb;</mo>
<msub>
<mi>s</mi>
<mi>k</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>t</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mi>m</mi>
</msup>
</mrow>
</mfrac>
<mo>,</mo>
<mn>1</mn>
<mo>&le;</mo>
<mi>k</mi>
<mo>&le;</mo>
<msub>
<mi>k</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,N为气候变量数据的个数,m为模糊加权指数且m=2,ti为气候变量数据对应的时间,p=int(N/kmax),当k=kmax,令kp=N;
处理器将时间坐标集合上的初始偏差定义为
<mrow>
<msubsup>
<mi>&sigma;</mi>
<mi>k</mi>
<mn>2</mn>
</msubsup>
<mo>=</mo>
<mfrac>
<mrow>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</msubsup>
<msup>
<mrow>
<mo>&lsqb;</mo>
<msub>
<mi>s</mi>
<mi>k</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>t</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mi>m</mi>
</msup>
<mrow>
<mo>(</mo>
<msub>
<mi>t</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msub>
<mi>&mu;</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>t</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msub>
<mi>&mu;</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
<mi>T</mi>
</msup>
</mrow>
<mrow>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</msubsup>
<msup>
<mrow>
<mo>&lsqb;</mo>
<msub>
<mi>s</mi>
<mi>k</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>t</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mi>m</mi>
</msup>
</mrow>
</mfrac>
<mo>,</mo>
<mn>1</mn>
<mo>&le;</mo>
<mi>k</mi>
<mo>&le;</mo>
<msub>
<mi>k</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,N为气候变量数据的个数,m为模糊加权指数且m=2,ti为气候变量数据对应的时间,μk为在时间坐标集合上的初始模糊聚类中心,sk(ti)为如上所述分段常函数;
所述处理器从气候变量数据集合中随机选取kmax个数据作为初始中心,得到
vk∈{xi,1≤i≤N},1≤k≤kmax (3)
其中,N为气候变量数据的个数,xi为气候变量数据;
所述处理器令所述气候变量数据集合中的初始协方差与单位矩阵成正比,得到
<mrow>
<msub>
<mi>&Sigma;</mi>
<mi>k</mi>
</msub>
<mo>=</mo>
<msup>
<mi>&sigma;</mi>
<mn>2</mn>
</msup>
<msub>
<mi>I</mi>
<mrow>
<msub>
<mi>k</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>&times;</mo>
<msub>
<mi>k</mi>
<mi>max</mi>
</msub>
</mrow>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,为kmax阶单位矩阵, 为气候变量数据的平均值;
所述处理器令初始混合概率αk与所述kmax成反比,得到
αk=1/kmax,1≤k≤kmax (5)
其中,kmax为初始聚类数;
步骤二、所述处理器求得隶属度、混合概率,并判断所述混合概率是否大于0,若是,则直接更新聚类参数以及所述隶属度,若否,则降低分割阶数后更新聚类参数以及所述隶属度,所述处理器求得隶属度、混合概率,包括:
采用公式
<mrow>
<msubsup>
<mi>u</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>k</mi>
</mrow>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>=</mo>
<mfrac>
<msup>
<mrow>
<mo>&lsqb;</mo>
<msubsup>
<mi>&alpha;</mi>
<mi>k</mi>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mi>p</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>t</mi>
<mi>i</mi>
</msub>
<mo>|</mo>
<msubsup>
<mi>&eta;</mi>
<mi>k</mi>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
<mi>p</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>|</mo>
<msubsup>
<mi>&eta;</mi>
<mi>k</mi>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<mi>m</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mrow>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<msub>
<mi>k</mi>
<mi>max</mi>
</msub>
</msubsup>
<msup>
<mrow>
<mo>&lsqb;</mo>
<msubsup>
<mi>&alpha;</mi>
<mi>j</mi>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mi>p</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>t</mi>
<mi>i</mi>
</msub>
<mo>|</mo>
<msubsup>
<mi>&eta;</mi>
<mi>j</mi>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
<mi>p</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>|</mo>
<msubsup>
<mi>&eta;</mi>
<mi>j</mi>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<mi>m</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
计算隶属度,其中,所述kmax为初始聚类数,αk=p(ηk)≥0为第k个聚类的混合概率,且满足条件为在时间坐标集合上的第k个高斯函数,p(x|ηk)为在气候变量数据集合上第k个高斯函数,m为模糊加权指数且m=2;
采用公式
<mrow>
<msubsup>
<mi>&alpha;</mi>
<mi>k</mi>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msubsup>
<mo>=</mo>
<mfrac>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
<mo>{</mo>
<mn>0</mn>
<mo>,</mo>
<mrow>
<mo>(</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msubsup>
<mi>u</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>k</mi>
</mrow>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mfrac>
<mi>c</mi>
<mn>2</mn>
</mfrac>
<mo>}</mo>
</mrow>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>k</mi>
</munderover>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
<mo>{</mo>
<mn>0</mn>
<mo>,</mo>
<mrow>
<mo>(</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msubsup>
<mi>u</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>k</mi>
</mrow>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mfrac>
<mi>c</mi>
<mn>2</mn>
</mfrac>
<mo>}</mo>
</mrow>
</mfrac>
<mo>,</mo>
<mo>{</mo>
<msubsup>
<mi>&alpha;</mi>
<mn>1</mn>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msubsup>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<msubsup>
<mi>&alpha;</mi>
<msub>
<mi>k</mi>
<mi>max</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msubsup>
<mo>}</mo>
<mo>=</mo>
<mfrac>
<mrow>
<mo>{</mo>
<msubsup>
<mi>&alpha;</mi>
<mn>1</mn>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msubsup>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<msubsup>
<mi>&alpha;</mi>
<msub>
<mi>k</mi>
<mi>max</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msubsup>
<mo>}</mo>
</mrow>
<mrow>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<msub>
<mi>k</mi>
<mi>max</mi>
</msub>
</msubsup>
<msubsup>
<mi>&alpha;</mi>
<mi>k</mi>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msubsup>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
计算混合概率,其中,所述kmax为初始聚类数,ui,k为第i个数据属于第k个聚类的隶属度,c为每个聚类成员所含参数的维数,N为气候变量数据的个数;
所述处理器采用期望值最大算法将所述聚类参数μk、vk、Σk以及隶属度ui,k更新为
<mrow>
<msubsup>
<mi>&mu;</mi>
<mi>k</mi>
<mrow>
<mo>&prime;</mo>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
</msubsup>
<mo>=</mo>
<mfrac>
<mrow>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</msubsup>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>u</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>k</mi>
</mrow>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
<mi>m</mi>
</msup>
<msub>
<mi>t</mi>
<mi>i</mi>
</msub>
</mrow>
<mrow>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</msubsup>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>u</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>k</mi>
</mrow>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
<mi>m</mi>
</msup>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>&sigma;</mi>
<mi>k</mi>
<mrow>
<mo>&prime;</mo>
<mn>2</mn>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msup>
<mo>=</mo>
<mfrac>
<mrow>
<msubsup>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</msubsup>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>u</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>k</mi>
</mrow>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
<mi>m</mi>
</msup>
<mrow>
<mo>(</mo>
<msub>
<mi>t</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msubsup>
<mi>&mu;</mi>
<mi>k</mi>
<mrow>
<mo>&prime;</mo>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>)</mo>
</mrow>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>t</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msubsup>
<mi>&mu;</mi>
<mi>k</mi>
<mrow>
<mo>&prime;</mo>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>)</mo>
</mrow>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
<mi>T</mi>
</msup>
</mrow>
<mrow>
<msubsup>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</msubsup>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>u</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>k</mi>
</mrow>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
<mi>m</mi>
</msup>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msubsup>
<mi>v</mi>
<mi>k</mi>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msubsup>
<mo>=</mo>
<mfrac>
<mrow>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</msubsup>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>u</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>k</mi>
</mrow>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
<mi>m</mi>
</msup>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
</mrow>
<mrow>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</msubsup>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>u</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>k</mi>
</mrow>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
<mi>m</mi>
</msup>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>10</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msubsup>
<mo>&Sigma;</mo>
<mi>k</mi>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msubsup>
<mo>=</mo>
<mfrac>
<mrow>
<msubsup>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</msubsup>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>u</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>k</mi>
</mrow>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
<mi>m</mi>
</msup>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msubsup>
<mi>v</mi>
<mi>k</mi>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msubsup>
<mi>v</mi>
<mi>k</mi>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
<mi>T</mi>
</msup>
</mrow>
<mrow>
<msubsup>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</msubsup>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>u</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>k</mi>
</mrow>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
<mi>m</mi>
</msup>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>11</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msubsup>
<mi>u</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>k</mi>
</mrow>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msubsup>
<mo>=</mo>
<mfrac>
<msup>
<mrow>
<mo>&lsqb;</mo>
<msubsup>
<mi>&alpha;</mi>
<mi>k</mi>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msubsup>
<mi>p</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>t</mi>
<mi>i</mi>
</msub>
<mo>|</mo>
<msubsup>
<mi>&eta;</mi>
<mi>k</mi>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
<mi>p</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>|</mo>
<msubsup>
<mi>&eta;</mi>
<mi>k</mi>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<mi>m</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mrow>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<msub>
<mi>k</mi>
<mi>max</mi>
</msub>
</msubsup>
<msup>
<mrow>
<mo>&lsqb;</mo>
<msubsup>
<mi>&alpha;</mi>
<mi>j</mi>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msubsup>
<mi>p</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>t</mi>
<mi>i</mi>
</msub>
<mo>|</mo>
<msubsup>
<mi>&eta;</mi>
<mi>j</mi>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
<mi>p</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>|</mo>
<msubsup>
<mi>&eta;</mi>
<mi>j</mi>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<mi>m</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>12</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,所述ui,k为第i个数据属于第k个聚类的隶属度,xi为气候变量数据,ti为气候变量数据xi对应的时间,μ′k与分别为在时间坐标集合上的第k个高斯函数p(t|ηk)的中心和偏差,vk与Σk分别为在气候变量数据集合上第k个高斯函数p(x|ηk)的中心和方差,m为模糊加权指数且m=2;
步骤三、所述处理器判断相邻两次迭代过程的模糊划分矩阵的相对误差是否大于终止阈值,若否,则返回步骤二,若是,则根据最小信息长度准则选择最优分割阶数,所述根据最小信息长度准则选择最优分割阶数,包括:处理器根据最小信息长度准则,采用公式
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mi>M</mi>
<mi>e</mi>
<mi>s</mi>
<mi>s</mi>
<mi>L</mi>
<mi>e</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<msup>
<mi>&eta;</mi>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msup>
<mo>,</mo>
<mi>Z</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mo>-</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<mi>log</mi>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>:</mo>
<msub>
<mi>&alpha;</mi>
<mi>k</mi>
</msub>
<mo>></mo>
<mn>0</mn>
</mrow>
</munder>
<msub>
<mi>&alpha;</mi>
<mi>k</mi>
</msub>
<mi>p</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>t</mi>
<mi>i</mi>
</msub>
<mo>|</mo>
<msubsup>
<mi>&eta;</mi>
<mi>k</mi>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
<mi>p</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>|</mo>
<msubsup>
<mi>&eta;</mi>
<mi>k</mi>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<mfrac>
<mi>c</mi>
<mn>4</mn>
</mfrac>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>:</mo>
<msub>
<mi>&alpha;</mi>
<mi>k</mi>
</msub>
<mo>></mo>
<mn>0</mn>
</mrow>
</munder>
<mi>log</mi>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>N&alpha;</mi>
<mi>k</mi>
</msub>
</mrow>
<mn>12</mn>
</mfrac>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mfrac>
<msub>
<mi>k</mi>
<mrow>
<mi>n</mi>
<mi>z</mi>
</mrow>
</msub>
<mn>2</mn>
</mfrac>
<mi>log</mi>
<mfrac>
<mi>N</mi>
<mn>12</mn>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<msub>
<mi>k</mi>
<mrow>
<mi>n</mi>
<mi>z</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>c</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
<mn>2</mn>
</mfrac>
<mo>.</mo>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>13</mn>
<mo>)</mo>
</mrow>
</mrow>
求得最小信息长度,其中,所述c为每个聚类成员所含参数的个数,αk为第k个聚类的混合概率,N为气候变量数据的个数,knz为当前聚类个数,p(t|ηk)为在时间坐标集合上的第k个高斯函数,p(x|ηk)为在气候变量数据集合上第k个高斯函数;
步骤四、所述处理器将所述最优分割阶数与最优聚类参数代入模糊聚类算法得出所述气候变量数据的渐进变化与跳跃变化,包括:处理器采用高斯函数得到在时间坐标集合上kmax个模糊分割βk(ti),公式为
<mrow>
<msub>
<mi>&beta;</mi>
<mi>k</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>t</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>A</mi>
<mi>k</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>t</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<msub>
<mi>k</mi>
<mi>max</mi>
</msub>
</msubsup>
<msub>
<mi>A</mi>
<mi>j</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>t</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>,</mo>
<mn>1</mn>
<mo>&le;</mo>
<mi>k</mi>
<mo>&le;</mo>
<msub>
<mi>k</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>14</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,k为当前模糊聚类数目,ti为气候变量数据对应的时间,μk为在时间坐标集合上的初始模糊聚类中心,为时间坐标集合上的初始偏差;
S3所述处理器根据所述渐进变化与所述跳跃变化对气候变化进行预测以及趋势分析。
2.根据权利要求1所述的方法,其特征在于,所述步骤三之后,还包括:
处理器判断若相邻两次迭代过程的模糊划分矩阵的相对误差低于终止阈值时,删除具有最低概率的聚类成员并返回所述步骤二,查找更小的信息长度值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510003459.1A CN104463382B (zh) | 2015-01-05 | 2015-01-05 | 气候变化预测以及趋势分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510003459.1A CN104463382B (zh) | 2015-01-05 | 2015-01-05 | 气候变化预测以及趋势分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104463382A CN104463382A (zh) | 2015-03-25 |
CN104463382B true CN104463382B (zh) | 2017-11-07 |
Family
ID=52909386
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510003459.1A Expired - Fee Related CN104463382B (zh) | 2015-01-05 | 2015-01-05 | 气候变化预测以及趋势分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104463382B (zh) |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103258142A (zh) * | 2013-06-04 | 2013-08-21 | 哈尔滨工业大学 | 基于天气模式的风电场短期风速预测方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA2446853A1 (en) * | 2001-05-15 | 2002-11-21 | Psychogenics Inc. | Systems and methods for monitoring behavior informatics |
-
2015
- 2015-01-05 CN CN201510003459.1A patent/CN104463382B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103258142A (zh) * | 2013-06-04 | 2013-08-21 | 哈尔滨工业大学 | 基于天气模式的风电场短期风速预测方法 |
Non-Patent Citations (3)
Title |
---|
春播期天气模糊聚类分析及判别预报;郑玉兰;《山东气象》;19880828;第14-15页 * |
景德镇市天气状况的模糊聚类分析;石卫华;《改革与开放》;20100925;第117、119页 * |
用模糊聚类分析方法作春季连阴雨预报;邹浩 等;《气象学报》;19830302;第41卷(第2期);第231-235页 * |
Also Published As
Publication number | Publication date |
---|---|
CN104463382A (zh) | 2015-03-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105426970B (zh) | 一种基于离散动态贝叶斯网络的气象威胁评估方法 | |
Wu et al. | Evolving RBF neural networks for rainfall prediction using hybrid particle swarm optimization and genetic algorithm | |
CN109146162B (zh) | 一种基于集成循环神经网络的概率风速预测方法 | |
CN107704875A (zh) | 基于改进ihcmac神经网络的建筑负荷预测方法和装置 | |
WO2021027153A1 (zh) | 交通流数据分析模型的构建方法和装置 | |
CN107705556A (zh) | 一种基于支持向量机和bp神经网络结合的交通流预测方法 | |
CN110071502A (zh) | 一种短期电力负荷预测的计算方法 | |
WO2021196743A1 (zh) | 热带气旋强度预报信息的生成方法及系统 | |
CN105225000B (zh) | 一种基于模糊序优化的风功率概率模型非参数核密度估计方法 | |
CN108280998A (zh) | 基于历史数据动态选择的短时交通流预测方法 | |
CN113537626B (zh) | 一种聚合信息差异的神经网络结合时间序列预测的方法 | |
CN112396152A (zh) | 一种基于cs-lstm的洪水预报方法 | |
CN108984830A (zh) | 一种基于模糊网络分析的建筑能效评价方法及装置 | |
CN106682764A (zh) | 一种基于并行预测策略的公共楼宇日前空调负荷预测方法 | |
CN113780684A (zh) | 一种基于lstm神经网络的智慧楼宇用户用能行为预测方法 | |
CN113361785A (zh) | 配电网短期负荷预测方法、装置、终端及存储介质 | |
CN116468138A (zh) | 空调负荷预测方法、系统、电子设备及计算机存储介质 | |
CN115915364A (zh) | 通信基站的节能方法、装置、计算机可读介质及设备 | |
Guo et al. | A combined model based on sparrow search optimized BP neural network and Markov chain for precipitation prediction in Zhengzhou City, China | |
CN108537370A (zh) | 基于混合智能降维算法的特大流域水电站群优化调度方法 | |
CN109426901A (zh) | 一种中长期用电预测方法及装置 | |
CN115238948A (zh) | 小水电发电量预测方法及装置 | |
CN114580762A (zh) | 一种基于XGBoost的水文预报误差校正方法 | |
CN104463382B (zh) | 气候变化预测以及趋势分析方法 | |
Zhang et al. | A modified random forest based on kappa measure and binary artificial bee colony algorithm |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20171107 Termination date: 20190105 |
|
CF01 | Termination of patent right due to non-payment of annual fee |