• 四川郎酒股份有限公司获第十二届人民企业社会责任奖年度环保奖 2019-05-13
  • 银保监会新规剑指大企业多头融资和过度融资 2019-05-12
  • 韩国再提4国联合申办世界杯 中国网友无视:我们自己来 2019-05-11
  • 中国人为什么一定要买房? 2019-05-11
  • 十九大精神进校园:风正扬帆当有为 勇做时代弄潮儿 2019-05-10
  • 粽叶飘香幸福邻里——廊坊市举办“我们的节日·端午”主题活动 2019-05-09
  • 太原设禁鸣路段 设备在测试中 2019-05-09
  • 拜耳医药保健有限公司获第十二届人民企业社会责任奖年度企业奖 2019-05-08
  • “港独”没出路!“梁天琦们”该醒醒了 2019-05-07
  • 陈卫平:中国文化内涵包含三方面 文化复兴表现在其中 2019-05-06
  • 人民日报客户端辟谣:“合成军装照”产品请放心使用 2019-05-05
  • 【十九大·理论新视野】为什么要“建设现代化经济体系”?   2019-05-04
  • 聚焦2017年乌鲁木齐市老城区改造提升工程 2019-05-04
  • 【专家谈】上合组织——构建区域命运共同体的有力实践者 2019-05-03
  • 【华商侃车NO.192】 亲!楼市火爆,别忘了买车位啊! 2019-05-03
    • / 23
    • 下载费用:30 金币  

    重庆时时彩后二杀5码: 一种均匀分布的波动数据预测方法.pdf

    关 键 词:
    一种 均匀分布 波动 数据 预测 方法
      专利查询网所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
    摘要
    申请专利号:

    CN201110057078.3

    申请日:

    2011.03.09

    公开号:

    CN102156802A

    公开日:

    2011.08.17

    当前法律状态:

    授权

    有效性:

    有权

    法律详情: 授权|||实质审查的生效IPC(主分类):G06F 19/00申请日:20110309|||专利申请权的转移IPC(主分类):G06F 19/00变更事项:申请人变更前权利人:北京航空航天大学变更后权利人:北京航空航天大学变更事项:地址变更前权利人:100191 北京市海淀区学院路37号变更后权利人:100191 北京市海淀区学院路37号变更事项:共同申请人变更后权利人:中航航空电子有限公司登记生效日:20111221|||著录事项变更IPC(主分类):G06F 19/00变更事项:发明人变更前:冀利刚 张叔农 周梦 张红 陈本刚 冯畅变更后:张叔农 冀利刚 周梦 张红 陈本刚 冯畅|||公开
    IPC分类号: G06F19/00 主分类号: G06F19/00
    申请人: 北京航空航天大学
    发明人: 冀利刚; 张叔农; 周梦; 张红; 陈本刚; 冯畅
    地址: 100191 北京市海淀区学院路37号
    优先权:
    专利代理机构: 北京慧泉知识产权代理有限公司 11232 代理人: 王顺荣;唐爱华
    PDF完整版下载: PDF下载
    法律状态
    申请(专利)号:

    CN201110057078.3

    授权公告号:

    ||||||||||||

    法律状态公告日:

    2016.04.27|||2013.09.11|||2012.02.01|||2012.02.01|||2011.08.17

    法律状态类型:

    授权|||实质审查的生效|||专利申请权、专利权的转移|||著录事项变更|||公开

    摘要

    本发明一种均匀分布的波动数据预测方法,它有六大步骤:步骤一、去均值;步骤二、找大数;步骤三、趋势处理;步骤四、灰色GM(1,1)预测;步骤五、还原数据;步骤六、误差评估。本发明是建立在灰色GM(1,1)模型基础上,通过对输入数据做均值、三角变换处理,使其具有光滑性增长特点,然后再利用GM(1,1)进行预测,最后还原即可得到预测数据。本发明构思科学,计算简单,工作量小,预测精度高。它在系统可靠性分析技术领域里具有较好的实用价值和广阔的应用前景。

    权利要求书

    1.一种均匀分布的波动数据预测方法,其特征在于:该方法具体步骤如下:设原始的n个均匀分布波动数据组成的序列为X(0)={x(0)(1),x(0)(2),…,x(0)(n)},用其预测m个数据,得到Y(0)={y(0)(1),…y(0)(n),y(0)(n+1),…,y(0)(n+m)},式中x(0)(k)、y(0)(k)分别表示序列X(0)、Y(0)中第k个数据,并且还用{y(0)(1),y(0)(2),…y(0)(n)}检测预测误差,{y(0)(n+1),…,y(0)(n+m)}即为预测m数据步骤一:去均值:对于给定原始数据序列X(0),先计算其均值u,再用序列X(0)中的每一个数减去均值u,即通过计算r(0)(k)=x(0)(k)-u??(k=1,…,n)得到序列R(0)={r(0)(1),r(0)(2),…,r(0)(n)};步骤二:找大数:在序列R(0)中找出绝对值最大的数,记为r(0)max;步骤三:趋势处理:通过三角变换,即计算将序列R(0)处理成具有增长趋势的序列R(1)={r(1)(1),r(1)2,…,r(1)(n)},式中arcsin为反正弦函数,π为圆周率,取π≈3.14;步骤四:灰色GM(1,1)预测:利用GM(1,1)模型对序列R(1)进行预测,得到序列R(2)={r(2)(1),r(2)(2),…,r(2)(m+n)};步骤五:还原数据:根据得到的序列R(2),最后再还原到预测数据Y(0);其中式中sin是正弦函数;步骤六:误差评估:定义误差率序列为ER,其中比较最大的误差率ermax,若|ermax|≤15%预测效果较高,否则预测效欠佳。2.根据权利要求1所述的一种均匀分布的波动数据预测方法,其特征在于:步骤四中所述的利用GM(1,1)模型对序列R(1)进行预测,得到序列R(2)={r(2)(1),r(2)(2),…,r(2)(m+n)}后的具体实施过程如下:①累加:对序列R(1)进行累加,得到序列R(11)={r(11)(1),r(11)(2),…,r(11)(n)},其中②计算紧邻均值:设Z(1)={z(1)(2),…,z(1)(n)}为R(11)的紧邻均值生成序列,其中然后再根据紧邻均值序列Z(1)与序列R(1)计算:B=-z(1)(2)1-z(1)(3)1......-z(1)(n)1,]]>γN=r(1)(2)r(1)(3)...r(1)(n);]]>③估计参数:用最小二乘法计算GM(1,1)模型r(1)(k)+a·z(1)(k)=b中参数a,b的值,即计算④计算拟合累加值:根据GM(1,1)模型及估计参数值计算预测的累加值R(12)={r(12)(1),r(12)(2),…,r(12)(n+m)},其中r(12)(1)=r(1)(1),式中e为自然对数底数,取e≈2.718;⑤累减还原:对计算出的拟合累加值进行累减还原,得到序列R(2),其中r(2)(1)=r(12)(1),r(2)(k)=r(12)(k)-r(12)(k-1)(k=2,3,…,m+n),完成GM(1,1)模型预测过程。

    说明书

    一种均匀分布的波动数据预测方法

    (一)技术领域:

    本发明涉及一种数据预测方法,尤其涉及一种均匀分布的波动数据预测方法。它用于业界均匀分布的波动数据预测,属于系统可靠性分析技术领域。

    (二)背景技术:

    在日常的生产生活中,我们碰到的很多数据具有均匀分布、呈现波动形式的情况,例如内燃机缸盖振动信号即为波动数据,该振荡数据可以有效反映内燃机运行状况。若能对此类数据提前预估,为尽早发现问题提供依据对我们的生产、生活都具有重要意义。

    目前在数据预测方面虽然方法有很多,但是适合均匀分布的波动数据预测却是一片空白。例如曲线拟合方法虽然原理简单、实施方便,但预测误差通常较大,而且很难找到一条与此类数据相吻合的曲线。时间序列模型其预测方法是建立在平稳随机过程假设的基础上的,对于实际运行设备的状态参数,一股得不到线性平稳的时间序列,尤其故障发生时更是如此,因而时间序列模型的应用范围受到很大的限制??ǘ?Kalman)滤波器有其局限性,只能用于线性系统。人工神经网络模型的隐蔽单元数目难以确定,且学习速度较慢,容易陷入局部极小点,迭代次数多,学习时间长,使得人工神经网络的实际应用受到限制。模糊模型虽适用于对复杂系统的预测,但由于模糊预测系统还处于研究阶段,有些问题还需着重解决,如静态知识库无法反映设备的失效过程,没有实时控制的特性,从而削弱了实用性。

    灰色模型是一种近年发展起来的预测方法,由我国邓聚龙教授于上世纪80年代提出,该模型具有所需样本少、预测精度高、运算量小等特点。但是灰色模型在运用中也有局限,若数据本身呈一定的增长或下降趋势,其预测精度一股较高,若数据本身具有一定的光滑性,预测精度会更好。而直接将灰色模型运用到均匀分布的波动数据预测,效果并不是很好,因此我们需要寻找一种计算简单,预测精度高的算法来处理均匀分布的波动数据预测。

    (三)发明内容:

    1、目的:针对以上均匀分布的波动数据预测中存在的问题,本发明的目的是提供了一种均匀分布的波动数据预测方法,它是基于灰色GM(1,1)模型的改进算法,构造了一个可以用于解决此类问题的预测器,从而,克服了现有技术的不足。

    2、技术方案:本发明是建立在灰色GM(1,1)模型基础上,通过对输入数据做均值、三角变换处理,使其具有光滑性增长特点,然后再利用GM(1,1)进行预测,最后还原即可得到预测数据。本方法中包含数学运算,因此直接编写计算机程序来实现,编程语言、版本均不限,C、C++、Java、Matlab等均可。

    本发明是一种均匀分布的波动数据预测方法,该方法具体步骤如下:

    设原始的n个均匀分布波动数据组成的序列为X(0)={x(0)(1),x(0)(2),…,x(0)(n)},用其预测m个数据,得到Y(0)={y(0)(1),…y(0)(n),y(0)(n+1),…,y(0)(n+m)},式中x(0)(k)、y(0)(k)分别表示序列X(0)、Y(0)中第k个数据。并且还可以用{y(0)(1),y(0)(2),…y(0)(n)}检测预测误差,{y(0)(n+1),…,y(0)(n+m)}即为预测m数据。

    步骤一:去均值:对于给定原始数据序列X(0),先计算其均值u,再用序列X(0)中的每一个数减去均值u,即通过计算r(0)(k)=x(0)(k)-u??(k=1,…,n)得到序列R(0)={r(0)(1),r(0)(2),…,r(0)(n)};

    步骤二:找大数:在序列R(0)中找出绝对值最大的数,记为r(0)max;

    步骤三:趋势处理:通过三角变换,即计算

    将序列R(0)处理成具有增长趋势的序列R(1)={r(1)(1),r(1)2,…,r(1)(n)},式中arcsin为反正弦函数,π为圆周率,取π≈3.14;

    步骤四:灰色GM(1,1)预测:利用GM(1,1)模型对序列R(1)进行预测,得到序列R(2)={r(2)(1),r(2)(2),…,r(2)(m+n)}。其具体实施过程如下:

    ①累加:对序列R(1)进行累加,得到序列R(11)={r(11)(1),r(11)(2),…,r(11)(n)},其中

    ②计算紧邻均值:设Z(1)={z(1)(2),…,z(1)(n)}为R(11)的紧邻均值生成序列,其中然后再根据紧邻均值序列Z(1)与序列R(1)计算:

    B=-z(1)(2)1-z(1)(3)1......-z(1)(n)1,]]>γN=r(1)(2)r(1)(3)...r(1)(n);]]>

    ③估计参数:用最小二乘法计算GM(1,1)模型r(1)(k)+a·z(1)(k)=b中参数a,b的值,即计算

    ④计算拟合累加值:根据GM(1,1)模型及估计参数值计算预测的累加值R(12)={r(12)(1),r(12)(2),…,r(12)(n+m)},其中r(12)(1)=r(1)(1),式中e为自然对数底数,取e≈2.718;

    ⑤累减还原:对计算出的拟合累加值进行累减还原,得到序列R(2),其中r(2)(1)=r(12)(1),r(2)(k)=r(12)(k)-r(12)(k-1)(k=2,3,…,m+n),完成GM(1,1)模型预测过程。

    步骤五:还原数据:根据得到的序列R(2),最后再还原到预测数据Y(0)。其中式中sin是正弦函数。

    步骤六:误差评估:定义误差率序列为ER,其中比较最大的误差率ermax,若|ermax|≤15%预测效果较高,否则预测效欠佳。

    其中,在步骤一中所述“去均值”的计算流程见图1;

    其中,在步骤二中所述“找大数”的计算流程见图2;

    其中,在步骤三中所述“趋势处理”的计算流程见图3;

    其中,在步骤四的第①小步中所述“累加”的计算流程见图4;

    其中,在步骤四的第②小步中所述“计算紧邻均值”的计算流程见图5;

    其中,在步骤四的第③小步中所述“估计参数”的计算流程见图6;

    其中,在步骤四的第④小步中所述“计算拟合累加值”的计算流程见图7;

    其中,在步骤四的第⑤小步中所述“累减还原”的计算流程见图8;

    其中,在步骤五中所述“还原数据”的计算流程见图9;

    其中,在步骤六中所述“误差评估”的计算流程见图10;

    其中,步骤一、二、三、五是本发明的创新点,步骤四、六是已存在的算法。

    3、优点及功效:

    (1)本发明是基于灰色GM(1,1)模型所做的预测,因此具有它本身算法简单,运算量小等优点。

    (2)本发明的创新之处在于对原始波动数据的均值变换、三角变换处理,使处理后的数据具有光滑性增长趋势,满足了灰色GM(1,1)模型对数据的要求,在很大程度上提高了预测精度,并且保证了预测出的数据仍然具有波动性。

    (四)附图说明:

    图1是“去均值”的计算流程图;

    图2是“找大数”的计算流程图;

    图3是“趋势处理”的计算流程图;

    图4是GM(1,1)预测中“累加”的计算流程图;

    图5是GM(1,1)预测中“计算紧邻均值”的计算流程图;

    图6是GM(1,1)预测中“估计参数”的计算流程图;

    图7是GM(1,1)预测中“计算拟合累加值”的计算流程图;

    图8是GM(1,1)预测中“累减还原”的计算流程图;

    图9是“还原数据”的计算流程图;

    图10是“误差评估”的计算流程图;

    图11是某些波动数据的预测效果图;

    图12是本发明流程框图。

    图中标号及符号说明如下:

    X(0)表示原始波动数据,即X(0)={x(0)(1),x(0)(2),…,x(0)(n)},x(0)(i)表示序列X(0)中第i个数;

    n表示序列X(0)中数据个数;

    sum表示临时变量,用来计算序列X(0)的和;

    u表示序列X(0)的算术平均值;

    R(0)表示序列X(0)去均值后的新序列,r(0)(i)表示该序列中第i个数;

    r(0)max表示序列R(0)中绝对值最大的数;

    *表示乘法;

    /表示除法;

    arcsin表示反正弦函数;

    π表示圆周率,取π≈3.141592;

    R(1)表示由序列R(0)趋势处理后的新序列,r(1)(i)表示该序列中第i个数;

    R(11)表示由序列R(1)累加后的新序列,r(11)(i)表示该序列中第i个数;

    B表示一个n-1行2列的矩阵,由序列R(11)的紧邻均值生成;

    γN表示一个n-1维的列向量,由序列R(1)生成

    D表示一个2行2列的矩阵,D=BTB;

    C表示一个2维列向量,C=BTγN;

    P表示一个2行2列的矩阵,P=D-1=(BTB)-1;

    表示模型r(1)(k)+a·z(1)(k)=b中a、b的拟合值;

    e表示自然对数底数,e≈2.718;

    R(12)表示由序列R(0)拟合的累加值序列,r(12)(i)表示该序列中第i个数;

    m表示要预测的数据个数;

    R(2)表示由序列R(0)预测的新序列,r(2)(i)表示该序列中第i个数;

    sin表示正弦函数;

    Y(0)表示由原始序列X(0)预测的新序列,y(0)(i)表示该序列中第i个数;

    ER表示预测Y(0)前n数的相对误差序列,er(i)表示该序列中第i个数,ermax表示最大的相对误差。

    o表示原数据点;

    *表示新方法预测的数据点。

    (五)具体实施方式:

    参见图12,本发明一种均匀分布的波动数据预测方法,该方法具体步骤如下:

    设原始的n个均匀分布波动数据组成的序列为X(0)={x(0)(1),x(0)(2),…,x(0)(n)},用其预测m个数据,得到Y(0)={y(0)(1),…y(0)(n),y(0)(n+1),…,y(0)(n+m)},式中x(0)(k)、y(0)(k)分别表示序列X(0)、Y(0)中第k个数据。并且还可以用{y(0)(1),y(0)(2),…y(0)(n)}检测预测误差,{y(0)(n+1),…,y(0)(n+m)}即为预测m数据

    步骤一:去均值:对于给定原始数据序列X(0),先计算其均值u,再用序列X(0)中的每一个数减去均值u,即通过计算r(0)(k)=x(0)(k)-u??(k=1,…,n)得到序列R(0)={r(0)(1),r(0)(2),…,r(0)(n)};

    步骤二:找大数:在序列R(0)中找出绝对值最大的数,记为r(0)max;

    步骤三:趋势处理:通过三角变换,即计算将序列R(0)处理成具有增长趋势的序列R(1)={r(1)(1),r(1)2,…,r(1)(n)},式中arcsin为反正弦函数,π为圆周率,取π≈3.14;

    步骤四:灰色GM(1,1)预测:利用GM(1,1)模型对序列R(1)进行预测,得到序列R(2)={r(2)(1),r(2)(2),…,r(2)(m+n)}。其具体实施过程如下:

    ①累加:对序列R(1)进行累加,得到序列R(11)={r(11)(1),r(11)(2),…,r(11)(n)},其中

    ②计算紧邻均值:设Z(1)={z(1)(2),…,z(1)(n)}为R(11)的紧邻均值生成序列,其中然后再根据紧邻均值序列Z(1)与序列R(1)

    计算:

    B=-z(1)(2)1-z(1)(3)1......-z(1)(n)1,]]>γN=r(1)(2)r(1)(3)...r(1)(n);]]>

    ③估计参数:用最小二乘法计算GM(1,1)模型r(1)(k)+a·z(1)(k)=b中参数a,b的值,即计算

    ④计算拟合累加值:根据GM(1,1)模型及估计参数值计算预测的累加值R(12)={r(12)(1),r(12)(2),…,r(12)(n+m)},其中r(12)(1)=r(1)(1),式中e为自然对数底数,取e≈2.718;

    ⑤累减还原:对计算出的拟合累加值进行累减还原,得到序列R(2),其中r(2)(1)=r(12)(1),r(2)(k)=r(12)(k)-r(12)(k-1)(k=2,3,…,m+n),完成GM(1,1)模型预测过程。

    步骤五:还原数据:根据得到的序列R(2),最后再还原到预测数据Y(0)。其中式中sin是正弦函数。

    步骤六:误差评估:定义误差率序列为ER,其中比较最大的误差率ermax,若|ermax|≤15%预测效果较高,否则预测效欠佳。

    该预测方法是用来解决均匀分布的波动数据的预测问题,接下来将结合一个实例,详细的阐述该方法的具体实施方式:

    数据准备:首先准备一组波动数据,该数据是围绕某一值均匀分布的??梢酝ü杉?,也可通过计算机随机生成。本实例采用的20数据如下:

    3.7831?4.1926?4.1273?3.9068?4.0874?4.0550?4.1456?4.2302?3.9872?3.8689

    3.8765?3.8748?4.1952?3.9246?3.9407?4.2236?3.9586?4.1253?4.1282?3.9913

    利用该组数据,预测其后10个数据。

    数据预测:

    (1)去均值:由序列知n=20,首先计算得该序列的均值u=4.0311,再计算r(0)(k)=x(0)(k)-u(k=1,…,20),得到R(0)序列:

    -0.2481?0.1615??0.0962?-0.1244?0.0563??0.0238?0.1144??0.1990?-0.0439?-0.1623

    -0.1547?-0.1564?0.1641?-0.1065?-0.0904?0.1924?-0.0725?0.0942?0.0970??-0.0398

    这一步骤可以通过两个简单循环操作完成,具体过程参见图1;

    (2)找大数:在序列R(0)中找到得绝对值最大的数为r(0)max=-0.2481,这一步骤可通过一个循环操作完成,具体过程参见图2;

    (3)增长性处理:利用序列R(0),计算得到具有增长趋势的序列R(1):

    7.8540???11.8576?18.4514?25.6578?31.1871?37.6029?43.5029?49.3345??56.7266

    63.5447??69.7881?76.0803?80.9588?88.4084?94.6208?99.6432?107.1108?112.7080

    118.9786?125.8249

    这一步骤可通过一个循环操作完成,具体过程参见图3;

    (4)灰色GM(1,1)预测:

    ①累加:利用序列R(1),计算得到累加的序列R(11):

    7.8540???19.7116??38.1630???63.8208???95.0079??132.6108?176.1138?225.4483

    282.1749?345.7196?415.5077??491.5880??572.5468?660.9552?755.5760

    855.2192?962.3300?1075.0380?1194.0166?1319.8415

    这一步骤可通过一个循环操作完成,具体过程参见图4;

    ②计算B矩阵与r向量:计算(n=20),得到的两个向量分别为:

    B=-13.78281-28.93731-50.99191-79.41441-113.80941-154.36231-200.78111-253.81161-313.94731-380.61371-453.54791-532.06741-616.75101-708.26561-805.39761-908.77461-1018.68401-1134.52731-1256.92911,]]>γN=11.857618.451425.657831.187137.602943.502949.334556.726663.544769.788176.080380.958888.408494.620899.6432107.1108112.7080118.9786125.8249;]]>

    这一步骤可通过一个循环操作完成,具体过程参见图5;

    ③估计参数:计算得到估计参数值:这一步骤的实现,若是用Matlab计算可以通过调用命令inv来求逆,再计算矩阵乘法即可。若是用C等编程语言,需要用两个循环过程完成,第一个循环先计算BTγN,第二个循环计算BTB,求逆可以直接采用伴随余子式计算,最后直接用二阶方阵的乘法定义计算参数即可,具体过程参见图6;

    ④计算拟合累加值:通过计算r(12)(1)=r(1)(1),(k=1,…,29),得到拟合累加序列R(12):

    7.8540????33.2415???65.3677???100.4032??138.6115??180.2799??225.7218??275.2789

    329.3238??388.2631??452.5398??522.6373??599.0829??682.4514??773.3696

    872.5214??980.6523??1098.5754?1227.1776?1367.4259?1520.375??1687.1751

    1869.0806?2067.4593?2283.8030?2519.7387?2777.0406?3057.6436?3363.6578

    3697.3845

    这一步骤可通过一个循环操作完成,具体过程参见图7;

    ⑤累减还原:利用E(12)计算e(2)(1)=e(12)(1),e(2)(k)=e(12)(k)-e(12)(k-1)(k=2,3,…,m+n),得到序列E(2):

    7.8540???25.3875??32.1262??35.0355??38.2083??41.6684??45.4419?49.5571

    54.0449??58.9392??64.2767??70.0976??76.4456??83.3684??90.9182?99.1517

    108.1309?117.9231?128.6022?140.2483?152.9491?166.8001?181.9054

    198.3787?216.3437?235.9357?257.3019?280.6030?306.0142?333.7267

    这一步骤可通过一个循环操作完成,具体过程参见图8;

    (5)还原数据:计算得到最后预测序列Y(0):

    3.7831?3.9686?3.8694?4.1453?3.9102?4.2138?3.7846?4.1926?4.1789?3.8618

    3.7850?3.8248?3.8163?3.7847?3.9848?4.2747?3.7910?4.2776?3.9811?3.8075

    3.8239?4.1034?4.1061?4.1409?3.9286?4.1083?4.1065?4.2401?4.2688?3.8680

    这一步骤可通过一个循环操作完成,具体过程参见图9;

    误差评估:

    计算得到误差率序列ER:

    0?5.57%?6.40%?-5.91%?4.39%??-3.94%?8.95%?0.93%??-4.75%?0.17%?2.27%

    1.24%???9.40%?3.47%??-1.09%?-1.27%?4.16%?-3.78%?3.65%??4.56%

    这一步骤可通过一个循环操作完成,具体过程参见图10,通过比较发现,绝对值最大的

    误差率为9.40%,预测效果较好,预测效果图可参见图11。

    根据上述算法,可以编制计算机程序,打包成一个预测器,用于均匀分布的波动数据的预测。以下是Matlab程序代码:

    关于本文
    本文标题:一种均匀分布的波动数据预测方法.pdf
    链接地址://www.4mum.com.cn/p-5867926.html
    关于我们 - 网站声明 - 网站地图 - 资源地图 - 友情链接 - 网站客服 - 联系我们

    [email protected] 2017-2018 www.4mum.com.cn网站版权所有
    经营许可证编号:粤ICP备17046363号-1 
     


    收起
    展开
  • 四川郎酒股份有限公司获第十二届人民企业社会责任奖年度环保奖 2019-05-13
  • 银保监会新规剑指大企业多头融资和过度融资 2019-05-12
  • 韩国再提4国联合申办世界杯 中国网友无视:我们自己来 2019-05-11
  • 中国人为什么一定要买房? 2019-05-11
  • 十九大精神进校园:风正扬帆当有为 勇做时代弄潮儿 2019-05-10
  • 粽叶飘香幸福邻里——廊坊市举办“我们的节日·端午”主题活动 2019-05-09
  • 太原设禁鸣路段 设备在测试中 2019-05-09
  • 拜耳医药保健有限公司获第十二届人民企业社会责任奖年度企业奖 2019-05-08
  • “港独”没出路!“梁天琦们”该醒醒了 2019-05-07
  • 陈卫平:中国文化内涵包含三方面 文化复兴表现在其中 2019-05-06
  • 人民日报客户端辟谣:“合成军装照”产品请放心使用 2019-05-05
  • 【十九大·理论新视野】为什么要“建设现代化经济体系”?   2019-05-04
  • 聚焦2017年乌鲁木齐市老城区改造提升工程 2019-05-04
  • 【专家谈】上合组织——构建区域命运共同体的有力实践者 2019-05-03
  • 【华商侃车NO.192】 亲!楼市火爆,别忘了买车位啊! 2019-05-03
  • 河南11选5分布走势图 北京十一选五开奖彩票控 广东11选5开奖信息 开元棋牌官网 体育彩票11选5开奖结果 中国人民银行股票行情大盘走势 香港正版彩票资料大全 gptjw股票推荐网 吉祥棋牌 大嘴刨幺 彩票微信群 事后分析 江西多乐彩任五遗漏 河南十一选五走势图带连线 浙江十一选五预测号码 河北11选5