• 四川郎酒股份有限公司获第十二届人民企业社会责任奖年度环保奖 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
    • / 12
    • 下载费用:30 金币  

    重庆时时彩总和路珠: GPS测站坐标时间序列周期性探测方法及系统.pdf

    摘要
    申请专利号:

    重庆时时彩单双窍门 www.4mum.com.cn CN201510173563.5

    申请日:

    2015.04.14

    公开号:

    CN104765055A

    公开日:

    2015.07.08

    当前法律状态:

    撤回

    有效性:

    无权

    法律详情: 发明专利申请公布后的视为撤回IPC(主分类):G01S 19/37申请公布日:20150708|||实质审查的生效IPC(主分类):G01S 19/37申请日:20150414|||公开
    IPC分类号: G01S19/37(2010.01)I 主分类号: G01S19/37
    申请人: 武汉大学
    发明人: 周晓慧; 姜卫平; 马一方
    地址: 430072湖北省武汉市武昌区珞珈山武汉大学
    优先权:
    专利代理机构: 武汉科皓知识产权代理事务所(特殊普通合伙)42222 代理人: 胡艳
    PDF完整版下载: PDF下载
    法律状态
    申请(专利)号:

    CN201510173563.5

    授权公告号:

    ||||||

    法律状态公告日:

    2017.09.01|||2015.08.05|||2015.07.08

    法律状态类型:

    发明专利申请公布后的视为撤回|||实质审查的生效|||公开

    摘要

    本发明提供了一种GPS测站坐标时间序列周期性探测方法及系统,包括步骤:步骤1,获取GPS测站坐标时间序列观测值;步骤2,剔除GPS测站坐标时间序列观测值粗差并修正天线相位中心偏差;步骤3,采用谐函数描述GPS测站坐标时间序列获得GPS测站坐标时间序列的谐函数模型,并设计谐函数模型矩阵;步骤4,采用最小二乘法解算谐函数模型矩阵,探测GPS测站坐标时间序列周期性。本发明可能处理含有时变周期性信号的GPS测站坐标时间序列,并保证探测准确性。

    权利要求书

    权利要求书
    1.  GPS测站坐标时间序列周期性探测方法,其特征是,包括步骤:
    步骤1,获取GPS测站坐标时间序列观测值;
    步骤2,剔除GPS测站坐标时间序列观测值粗差并修正天线相位中心偏差;
    步骤3,采用谐函数描述GPS测站坐标时间序列获得GPS测站坐标时间序列的谐函数模型,并设计谐函数模型矩阵其中,y表示GPS测站坐标时间序列观测值;A=1t11t2......1tm,]]>Ak=cosωkt1sinωkt1cosωkt2sinωkt2......cosωktmsinωktm,]]>ti为观测历元,i=1,2,...,m,m为观测历元数量,ωk表示可变周期性信号的频率;xk=akbk,]]>ak、bk分别为正弦函数分量、余弦函数分量;x=y(0)a0,]]>y(0)为GPS测站参考历元坐标,a0为GPS测站的线性速度;k表示谐函数编号,q为谐函数数量;
    步骤4,采用最小二乘法解算谐函数模型矩阵,获得可变周期性信号的频率估计值,并采用假设检验方法验证频率估计值;
    步骤5,基于假设检验验证后的频率估计值,采用最小二乘方法解算谐函数模型矩阵获得正弦函数分量ak、余弦函数分量bk和线性速度a0。

    2.  如权利要求1所述的GPS测站坐标时间序列周期性探测方法,其特征是:
    步骤1中通过双差定位软件工具、精密单点定位软件工具或IGS分析中心获取GPS测站坐标时间序列观测值。

    3.  如权利要求1所述的GPS测站坐标时间序列周期性探测方法,其特征是:
    所述的GPS测站坐标时间序列的谐函数模型为:
    y(ti)=y(0)+a0ti+Σk=1q[akcos(ωkti)+bksin(ωkti)]]]>
    其中,y(ti)为观测历元ti对应的GPS测站坐标观测值,y(0)为GPS测站参考历元坐标;a0为GPS测站的线性速度,ti为观测历元;ak、bk分别为正弦函数分量、余弦函数分量,ωk表示可变周期性信号的频率;k表示谐波函数编号,q为谐波函数数量。

    4.  如权利要求1所述的GPS测站坐标时间序列周期性探测方法,其特征是:
    步骤4进一步包括子步骤:

    4.  1解算以下最大化问题以探测频率ωj及其对应的Aj:
    ωj=argmaxωl||PA‾ly||Qy-12,A‾l=PA‾Al]]>
    PA‾=I-A‾(A‾TQy-1A‾)-1A‾TQy-1]]>
    PA‾l=A‾l(A‾lTQy-1A‾l)-1A‾lTQy-1]]>
    其中,I为单位矩阵,是L2范数运算符,Qy是GPS测站坐标时间序列观测值的方差协方差矩阵;A‾=[A,A1,...,Aj-1],]]>Al=cosωlt1sinωlt1cosωlt2sinωlt2......cosωltmsinωltm,]]>l=1,2,...j-1;

    4.  2采用假设H0:y=Ax+Σk=1j-1(Akxk)Ha:y=Ax+Σk=1j(Akxk)]]>判断是否向谐函数模型增加新频率,即根据步骤4.1获得的探测频率ωj及其对应的Aj判定是否接受零假设,若接受,则增加新频率,然后执行步骤4.1;否则,执行步骤4.3;H0为零假设,Ha为备选假设;

    4.  3采用Qy=σ2I和统计量评估假设检验,其中,为零假设下的最小二乘残差;为备选假设下的验后方差,为备选假设下的最小二乘残差,统计量T2≈F(2,m-n-2s),s为假设检验过程中确定的备选谐函数数量。

    5.  如权利要求1所述的GPS测站坐标时间序列周期性探测方法,其特征是:
    步骤5具体为:
    采用最小二乘法,通过解算以下矩阵获得正弦函数分量ak和余弦函数分量bk:
    y=Bx+v
    x^=(BTQy-1B)-1BTQy-1y]]>
    Qx^=(BTQY-1B)-1]]>
    其中,矩阵B=[A Ak],k=1,...,j;为x的最小二乘估值;Qy是GPS测站坐标时间序列观测值的方差协方差矩阵;y是GPS测站坐标时间序列观测值,v是随机误差矢量,采用最小二乘法解算获得。

    6.  GPS测站坐标时间序列周期性探测系统,其特征是,包括:
    观测值获得???,用来获取GPS测站坐标时间序列观测值;
    修正???,用来剔除GPS测站坐标时间序列观测值粗差并修正天线相位中心偏差;
    谐函数模型构建???,用来采用谐函数描述GPS测站坐标时间序列获得GPS测站坐标时间序列的谐函数模型,并设计谐函数模型矩阵其中,y表示GPS测站坐标时间序列观测值;A=1t11t2......1tm,]]>Ak=cosωkt1sinωkt1cosωkt2sinωkt2......cosωktmsinωktm,]]>ti为观测历元,i=1,2,...,m,m为观测历元数量,ωk表示可变周期性信号的频率;xk=akbk,]]>ak、bk分别为正弦函数分量、余弦函数分量;x=y(0)a0,]]>y(0)为GPS测站参考历元坐标,a0为GPS测站的线性速度;k表示谐函数编号,q为谐函数数量;
    频率估计???,用来采用最小二乘法解算谐函数模型矩阵,获得可变周期性信号的频率估计值,并采用假设检验方法验证频率估计值;
    周期性探测???,用来基于假设检验验证后的频率估计值,采用最小二乘方法解算谐函数模型矩阵获得正弦函数分量ak、余弦函数分量bk和线性速度a0。

    说明书

    说明书GPS测站坐标时间序列周期性探测方法及系统
    技术领域
    本发明属于GPS数据精密处理技术领域,具体涉及一种GPS测站坐标时间序列周期性探测方法及系统。
    背景技术
    GPS测站坐标时间序列主要用于噪声模型和速度模型的建立。现有的GPS测站坐标时间序列预测估计模型通常包含季节性信号,并且将周期性信号如周年信号和半周年等作为常数参数广泛应用(Blewitt and Lavallée,2002),一般采用周年周期及其谐波和常数振幅的正弦曲线建模(Nicolaidis,2002)。
    GPS测站坐标时间序列中存在的季节性信号是对环境变化的响应,而隐藏在季节性信号中的环境噪声是不稳定的。研究成果表明,环境噪声的表现形式具有周期性可变特征(Davis,Wernicke et al.,2012)。
    然而,数据处理分析策略与GPS相关的技术因素、固体地球质量负荷、大气压负荷、水储量变化、岩石热膨胀、降雨、多路径效应等因素的综合影响,导致GPS测站坐标时间序列中的季节性信号呈现周期性可变特征(Davis,Wernicke et al.,2012)。田云锋(2011)对青藏高原和喜马拉雅地区的GPS垂向分量分析发现存在明显的相位变化,推断是受到地表水体因子控制;而在CMONOC及周边IGS台站的位置时间序列中均发现了周期约为351/n(n=1,…,6)天的“异?!敝芷谙?,并认为地表质量负荷并非该异常周期项的来源。
    目前GPS测站坐标时间序列周期性探测方法存在不足:将周期性时变信号当作时不变信号进行描述,对速率估值及其不确定性都会产生明显影响。
    发明内容
    针对现有技术存在的问题,本发明提供了一种利用谐函数和最小二乘准则、顾及周期性可变信号的GPS测站坐标时间序列周期性探测方法及系统。
    本发明采用谐函数对GPS测站坐标时间序列建模,通过最小二乘准则估计GPS测站周期性可变信号,从而实现GPS测站坐标时间序列周期性探测。
    为解决上述技术问题,本发明采用如下的技术方案:
    一、GPS测站坐标时间序列周期性探测方法,包括步骤:
    步骤1,获取GPS测站坐标时间序列观测值;
    步骤2,剔除GPS测站坐标时间序列观测值粗差并修正天线相位中心偏差;
    步骤3,采用谐函数描述GPS测站坐标时间序列获得GPS测站坐标时间序列的谐函数模型,并设计谐函数模型矩阵其中,y表示GPS测站坐标时间序列观测值;A=1t11t2......1tm,]]>Ak=cosωkt1sinωkt1cosωkt2sinωkt2......cosωktmsinωktm,]]>ti为观测历元,i=1,2,...,m,m为观测历元数量,ωk表示可变周期性信号的频率;xk=akbk,]]>ak、bk分别为正弦函数分量、余弦函数分量;x=y(0)a0,]]>y(0)为GPS测站参考历元坐标,a0为GPS测站的线性速度;k表示谐函数编号,q为谐函数数量;
    步骤4,采用最小二乘法解算谐函数模型矩阵,获得可变周期性信号的频率估计值,并采用假设检验方法验证频率估计值;
    步骤5,基于假设检验验证后的频率估计值,采用最小二乘方法解算谐函数模型矩阵获得正弦函数分量ak、余弦函数分量bk和线性速度a0。
    步骤1中通过双差定位软件工具、精密单点定位软件工具或IGS分析中心获取GPS测站坐标时间序列观测值。
    步骤3中,所述的GPS测站坐标时间序列的谐函数模型为:
    y(ti)=y(0)+a0ti+Σk=1q[akcos(ωkti)+bksin(ωkti)]]]>
    其中,y(ti)为观测历元ti对应的GPS测站坐标观测值,y(0)为GPS测站参考历元坐标;a0为GPS测站的线性速度,ti为观测历元;ak、bk分别为正弦函数分量、余弦函数分量,ωk表示可变周期性信号的频率;k表示谐波函数编号,q为谐波函数数量。
    步骤4进一步包括子步骤:
    4.1解算以下最大化问题以探测频率ωj及其对应的Aj:
    ωj=argmaxωl||PA‾ly||Qy-12,A‾l=PA‾Al]]>
    PA‾=I-A‾(A‾TQy-1A‾)-1A‾TQy-1]]>
    PA‾l=A‾l(A‾lTQy-1A‾l)-1A‾lTQy-1]]>
    其中,I为单位矩阵,是L2范数运算符,Qy是GPS测站坐标时间序列观测值的方差协方差矩阵;A‾=[A,A1,...,Aj-1],]]>Al=cosωlt1sinωlt1cosωlt2sinωlt2......cosωltmsinωltm,]]>l=1,2,...j-1;
    4.2采用假设H0:y=Ax+Σk=1j-1(Akxk)Ha:y=Ax+Σk=1j(Akxk)]]>判断是否向谐函数模型增加新频率,即根据步骤4.1获得的探测频率ωj及其对应的Aj判定是否接受零假设,若接受,则增加新频率,然后执行步骤4.1;否则,执行步骤4.3;H0为零假设,Ha为备选假设;
    4.3采用Qy=σ2I和统计量评估假设检验,其中,为零假设下的最小二乘残差;为备选假设下的验后方差,为备选假设下的最小二乘残差,统计量T2≈F(2,m-n-2s),s为假设检验过程中确定的备选谐函数数量。
    步骤5具体为:
    采用最小二乘法,通过解算以下矩阵获得正弦函数分量ak和余弦函数分量bk:
    y=Bx+v
    x^=(BTQy-1B)-1BTQy-1y]]>
    Qx^=(BTQY-1B)-1]]>
    其中,矩阵B=[A Ak],k=1,...,j;为x的最小二乘估值;Qy是GPS测站坐标时间序列观测值的方差协方差矩阵;y是GPS测站坐标时间序列观测值,v是随机误差矢量,采用最小二乘法解算获得。
    二、GPS测站坐标时间序列周期性探测系统,包括:
    观测值获得???,用来获取GPS测站坐标时间序列观测值;
    修正???,用来剔除GPS测站坐标时间序列观测值粗差并修正天线相位中心偏差;
    谐函数模型构建???,用来采用谐函数描述GPS测站坐标时间序列获得GPS测站坐标时间序列的谐函数模型,并设计谐函数模型矩阵其中,y表示GPS测站坐标时间序列观测值;A=1t11t2......1tm,]]>Ak=cosωkt1sinωkt1cosωkt2sinωkt2......cosωktmsinωktm,]]>ti为观测历元,i=1,2,...,m,m为观测历元数量,ωk表示可变周期性信号的频率;xk=akbk,]]>ak、bk分别为正弦函数分量、余弦函数分量;x=y(0)a0,]]>y(0)为GPS测站参考历元坐标,a0为GPS测站的线性速度;k表示谐函数编号,q为谐函数数量;
    频率估计???,用来采用最小二乘法解算谐函数模型矩阵,获得可变周期性信号的频率估计值,并采用假设检验方法验证频率估计值;
    周期性探测???,用来基于假设检验验证后的频率估计值,采用最小二乘方法解算谐函数模型矩阵获得正弦函数分量ak、余弦函数分量bk和线性速度a0。
    与现有技术相比,本发明具有特点:
    本发明采用谐函数建立GPS测站坐标时间序列预估模型,突破了现有技术中仅采用周期性时不变信号描述坐标时间序列的现状;建立顾及周期性可变特征信号模型,一方面有助于进一步精化已有的GPS测站坐标时间序列的速度模型,提高速度估值的可靠性;另一方面,顾及周期性可变特征信号模型有助于真实反映与当地环境效应相关的(如降雨、温度、地表负荷及含水层抽水等)季节性地球物理信号,以进一步开展地球物理过程等方面的建模及解释。
    现有GPS测站坐标时间序列周期性探测中,将周期性信号当作时不变信号进行估计,这样导致GPS测站坐标时间序列速度模型的不完善、速率估值及其不确定性的较大偏差,本发明可解决上述问题,可处理含有时变周期性信号的GPS测站坐标时间序列,并保证周期性探测的准确性。
    附图说明
    图1为本发明方法的具体流程示意图。
    具体实施方式
    为了使本发明目的、技术方案及有益效果更加清楚明白,下面将结合附图及具体实施方式,进一步说明本发明。应当理解,以下描述的具体实施方式仅用以解释本发明,并不用于限定本发明。
    本发明具体步骤如下:
    步骤1,获取GPS测站坐标时间序列观测值。
    本步骤属于本技术领域的常规技术,具体可通过双差定位软件(如成熟的GAMIT/GLOBK数据处理软件)、精密单点定位软件(如成熟的GIPSY/OASIS数据处理软件)、或IGS分析中心获取GPS测站坐标时间序列观测值。
    步骤2,剔除GPS测站坐标时间序列观测值中粗差,修正与更换硬件相关的天线相位中心偏差。
    本步骤属于本技术领域内常规的数据预处理方式,在此不作赘述。
    步骤3,采用谐函数描述GPS测站坐标时间序列,获得GPS测站坐标时间序列谐函数模型。
    GPS测站坐标时间序列谐函数模型如下:
    y(ti)=y(0)+a0ti+Σk=1q[akcos(ωkti)+bksin(ωkti)]---(1)]]>
    式(1)中:
    y(ti)为观测历元ti对应的GPS测站坐标观测值;
    y(0)为GPS测站参考历元坐标;
    a0为GPS测站的线性速度;
    系数ak、bk分别为正弦函数分量、余弦函数分量,用来描述GPS测站周期性运动的振幅,对应于频率ωk和观测历元ti,其中,i=1,2,...,m,m为观测历元数量;
    q为谐函数数量,k为谐函数编号,这里的谐函数指akcos(ωkti)+bksin(ωkti)。
    采用矩阵表达公式(1),获得GPS测站坐标时间序列谐函数模型矩阵:
    y=Ax+Σk=1q(Akxk)---(2)]]>
    式(2)中:
    y表示GPS测站坐标时间序列观测值;
    A为GPS测站坐标时间序列谐函数模型的线性回归系数构成的2列矩阵,即A=1t11t2......1tm;]]>
    Ak为cos(ωkt)和sin(ωkt)构成的2列矩阵,即Ak=cosωkt1sinωkt1cosωkt2sinωkt2......cosωktmsinωktm;]]>
    xk=akbk;]]>
    x=y(0)a0,]]>y(0)为GPS测站参考历元坐标。
    GPS测站坐标时间序列谐函数模型中待估参数包括频率ωk及系数ak、bk。
    本具体实施方式中,计算GPS测站坐标时间序列时不仅考虑测站的线性速 度和周年/半周年等周期性信号影响,更考虑了可变周期性信号。
    与以往GPS测站坐标时间序列建模方法不同,本发明充分考虑了不同因素对GPS测站坐标时间序列造成的可变周期性,同时利用利用三角函数系数及其对应的频率将这些可变周期性纳入计算体系,并通过假设检验方法对待估频率ωk(即可变周期性信号的频率)进行数值计算、验证。
    步骤4,解算GPS测站坐标时间序列谐函数模型矩阵获得待估频率ωk,并采用统计检验方法验证待估频率ωk。
    采用最小二乘法确定待估频率ωk,具体如下:
    H0:y=Ax+Σk=1j-1(Akxk)Ha:y=Ax+Σk=1j(Akxk)---(3)]]>
    式(3)中:
    H0表示零假设,Ha表示备选假设,j=1、2、…、q。
    该过程迭代完成,j从1到q取值;每次运行,如果接受零假设,即向函数模型增加一个新的频率。
    4.1解算以下最大化问题以探测待估频率ωj(Teunissen,2000)及其对应的Aj:
    ωj=argmaxωl||PA‾ly||Qy-12,A‾l=PA‾Al---(4)]]>
    PA‾=I-A‾(A‾TQy-1A‾)-1A‾TQy-1---(5)]]>
    PA‾l=A‾l(A‾lTQy-1A‾l)-1A‾lTQy-1---(6)]]>
    公式(4)~(6)中:
    I为单位矩阵;
    运算符表示L2范数;
    Qy是GPS测站坐标时间序列观测值y的方差协方差矩阵;
    Al为的子矩阵,l=1,2,...j-1,Al结构同Ak相同, 即Al=cosωlt1sinωlt1cosωlt2sinωlt2......cosωltmsinωltm,]]>用于假设检验待估频率ωl对应的系数al和bl。满足公式(4)-(6)中取最大值时的ωl即待估频率ωj,ωl对应的Al被设置为系数矩阵Aj。
    采用公式(4)~(6)所示的最大化问题解算待估频率ωj,采用公式(3)的假设检验判定是否增加新的待估频率。
    采用数值方法解算上述问题,GPS测站坐标时间序列的功率谱通过应用不同频率ωj的谱值具有最大谱值的连续频率用来构建矩阵Aj。
    4.2采用Qy=σ2I评估假设检验,先验单位权方差σ2未知。采用下述统计量进行评估:

    A‾i=PA‾Ai---(8)]]>
    式(7)~(8)中:
    为零假设下的最小二乘残差;
    为备选假设下的验后方差,采用计算,其中,df为自由度,为备选假设下的最小二乘残差;
    公式(7)中统计量T2的分布为自由度为2和m-n-2s的F-分布,即T2≈F(2,m-n-2s),m为观测历元数量,n为未知参数数量,本发明中n为3,m-n表示模型冗余度;s为假设检验过程中确定的备选谐函数数量。
    本步骤为假设检验统计标准过程,为便于理解,特说明本步骤过程:T2为统计量,根据统计量T2对应的T2≈F(2,m-n-2s)分布可以得到一个具体分布值; 另根据公式(7)也可计算T2的分布值,将两个分布值比较,判断假设检验中是接受零假设、还是接受备选假设。
    步骤5,基于假设检验验证后的频率ωk,采用最小二乘方法解算系数ak、bk。
    采用最小二乘法,解算以下矩阵:
    y=Bx+v        (9)
    x^=(BTQy-1B)-1BTQy-1y---(10)]]>
    Qx^=(BTQY-1B)-1---(11)]]>
    式(9)~(11)中:
    设计矩阵B包括矩阵A和Ak,k=1,...,j,即B=[A Ak];
    为x的最小二乘估值;
    v是随机误差矢量,采用最小二乘法解算获得。
    采用上述方法估计任一GPS测站可变周期性信号频率及其对应振幅、线性速度等参数,从而实现更完善的速度模型的建立,提高速度估值可靠性。
    本文中所描述的具体实施仅仅是对本发明精神作举例说明。本发明所属技术领域的技术人员可以对所描述的具体实施做各种的修改或补充或采用类似的方式替代,但并不会偏离本发明的精神或者超越所附权利要求书所定义的范围。

    关 键 词:
    GPS 坐标 时间 序列 周期性 探测 方法 系统
      专利查询网所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
    关于本文
    本文标题:GPS测站坐标时间序列周期性探测方法及系统.pdf
    链接地址://www.4mum.com.cn/p-5894810.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