CN101135726A - 一种星载sar内定标信号处理平台***及实现方法 - Google Patents

一种星载sar内定标信号处理平台***及实现方法 Download PDF

Info

Publication number
CN101135726A
CN101135726A CNA2007101221069A CN200710122106A CN101135726A CN 101135726 A CN101135726 A CN 101135726A CN A2007101221069 A CNA2007101221069 A CN A2007101221069A CN 200710122106 A CN200710122106 A CN 200710122106A CN 101135726 A CN101135726 A CN 101135726A
Authority
CN
China
Prior art keywords
phase
module
unit
signal
error
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.)
Granted
Application number
CNA2007101221069A
Other languages
English (en)
Other versions
CN101135726B (zh
Inventor
陈杰
杨威
王鹏波
周荫清
李春升
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beihang University
Original Assignee
Beihang University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Beihang University filed Critical Beihang University
Priority to CN2007101221069A priority Critical patent/CN101135726B/zh
Publication of CN101135726A publication Critical patent/CN101135726A/zh
Application granted granted Critical
Publication of CN101135726B publication Critical patent/CN101135726B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明公开了一种星载SAR内定标信号处理平台***及实现方法,它具有高效、精确、完备、可靠性高的处理特点。该平台由六个子模块组成,分别是有效数据抽取模块、稳定性定标模块、相位误差定标模块、TR组件频域定标模块、可靠性检验模块和分析结果判读模块。该星载SAR内定标信号处理实现方法的处理对象为星载SAR内定标信号,采用协同并行处理的工作方式处理内定标信号,不仅可以得到更准确的处理结果,而且还能对雷达***工作状态进行定量的描述。

Description

一种星载SAR内定标信号处理平台***及实现方法
技术领域
本发明属于信号处理领域,涉及一种信号处理平台,特别涉及一种星载SAR内定标信号处理平台。
背景技术
合成孔径雷达(SAR)卫星近些年来发展迅速,由于SAR卫星不受天气、地理、时间等因素的限制,能够对地进行全天时的观测,且具有一定的穿透力,因而被广泛的应用于军事侦察、资源探测、海洋观测、生态监测等方面。
在雷达内部传输过程中,为了扩大SAR的应用领域,获取清晰度高、可视性高的SAR遥感图像,需要补偿雷达的回波信号。但由于雷达内部器件的非理想性会引起信号幅度和相位的畸变,而要补偿这种畸变所带来的误差,先要测出整个雷达***内部的传递函数,则必须首先保证雷达处于正常的工作状态,否则测得的传递函数就会存在较大的误差,造成卫星对畸变产生的误差的补偿能力也将很有限。这正是目前内定标信号处理过程中的一个缺陷,即忽略了雷达处于正常工作状态这个前提。
目前,对内定标信号的研究还较少,主要是针对补偿雷达回波信号在雷达内部传输时所引入的信号幅度和相位的畸变,以及后续的一些研究,而对雷达工作状态是否正常进行判断这方面基本属于空白。目前国内基本上都采用一种三回路信号定标的方法。三个定标回路是指发射定标回路,接收定标回路和参考定标回路。发射定标回路完成对发射定标回路信号的采样,分析发射回路信号的波形、测定发射通道的增益、稳定性,并获取发射定标回路的传递函数;接收定标回路完成对接收定标回路信号的采样,分析接收回路信号的波形,测定接收定标回路前端增益及其稳定性,获取接收定标回路的传递函数;参考定标回路完成对参考定标信号的复制,获取参考定标回路的传递函数。
这些处理方法的一个共同点就是要获取各个回路的传递函数,但是在获取这些回路的传递函数处理之前,对数据没有做任何的预处理,也没有对雷达的工作状态进行分析,因而获取的传递函数可靠性、精度都不高。
发明内容
本发明的目的是提供一种星载SAR内定标信号处理平台***及实现方法,它的处理对象为星载SAR内定标信号,通过输入信号时有效数据抽取模块对有效内定标信号的抽取,以及可靠性检验模块对稳定性定标模块、相位误差定标模块的输出结果对比映证,采用协同并行处理的工作方式处理内定标信号,不仅可以得到更准确的处理结果,而且还能对雷达***工作状态进行定量的描述,具有高效、精确、完备、可靠性高的处理特点。
本发明是一种星载SAR内定标信号处理平台***,包括有效数据抽取模块、分析结果判读模块,其特征在于,还包括稳定性定标模块、相位误差定标模块、TR组件频域定标模块和可靠性检验模块。有效数据抽取模块用自适应相邻对比法将输入平台的原始内定标信号中的有效数据抽取出来,并把这组有效数据分别输出到稳定性定标模块、相位误差定标模块、TR组件频域定标模块和可靠性检验模块;稳定性定标模块测定匹配滤波后的信号功率峰值处的幅度和相位,反映雷达***的工作状态,并分别将处理后的相位和幅度信息输出到分析结果判读模块中,幅度变化曲线输出到可靠性验证模块中;相位误差定标模块采用高精度正交域误差相位分解方法对误差相位进行分解,得到正交的各次相位误差的峰值,并分别将处理后的各高次相位误差峰峰值输出到分析结果判读模块,各高次相位误差拟合曲线输出到可靠性验证模块;TR组件频域定标模块主要检测TR组件的工作状态,并将处理后的结果输出到分析结果判读模块;可靠性验证模块功能为验证整个平台分析结果正确与否,并将结果输出到分析结果判读模块;最终由分析结果判读模块分析所有的处理结果,从而判断雷达的工作状态。
所述有效数据抽取模块利用自适应相邻对比法抽取有效数据。由输入接口单元、匹配滤波器和输出接口单元组成。原始内定标信号按帧的形式由输入接口单元输入,并成对进入匹配滤波器中进行匹配滤波;取其匹配滤波后功率的最大值之间差的绝对值。如果该值很小,可直接去除,不进行后续操作;如果该值很大,判定这两组数据为有效数据,通过输出接口单元,将信号输出到稳定性定标模块、相位误差定标模块、TR组件频域定标模块、可靠性检验模块;处理完毕后进行下一对数据帧的处理,并重复操作,直到处理完全部的原始数据。
所述稳定性定标模块将内定标数据进行匹配滤波,测定滤波后的信号功率峰值处的幅度和相位,分析其随时间的变化趋势。由输入接口单元、幅度特性保存单元、匹配滤波器、插值单元和输出接口单元组成。通过输入接口单元的信号帧分两路,分别进入幅度特性保存单元和匹配滤波器,进入幅度特性保存单元可将信号帧的信号幅度变换为曲线,并保存在该单元内;进入匹配滤波器的信号滤波后输出到插值单元,并取功率峰值处附近的数据进行插值;在插值序列中重新选取峰值的位置,得到该位置信号的幅度和相位并保存,然后重复插值操作,直到匹配滤波后的最大值处信号的相位变化在误差允许的微小范围内。输出接口单元将幅度变化曲线传递给可靠性检验模块,并把功率峰值处的相位和幅度传递给分析结果判读模块。
所述相位误差定标模块由输入接口单元、相位去除单元、相位拟合单元、相位误差预处理单元、相位误差分析单元和输出接口单元组成;从输入接口单元按帧读入有效的内定标信号到相位去除单元,在相位去除单元中,标准二次相位被除去,剩下的为误差相位;相位去除单元的输出信号输入到相位拟合单元,获取相位拟合曲线;相位拟合曲线通过相位误差预处理单元,得到相位误差曲线,并输入到相位误差分析单元,获取各高次相位误差的拟合曲线和峰峰值,并把拟合曲线通过输出接口单元传递给可靠性检验模块,把峰峰值传递给分析结果判读模块;重复上述操作,直到有效的内定标信号处理完毕。
所述可靠性检验模块通过借鉴点目标评估中的指标来对内定标信号匹配滤波后的结果进行分析,由输入接口单元、匹配滤波器、误差注入单元、减法器单元和输出接口单元组成;稳定性定标模块和相位误差定标模块的输出结果输入到误差注入单元注入到理想线性调频信号中,得到的合成信号在匹配滤波器中进行匹配滤波,获取其分辨率、积分旁瓣比、峰值旁瓣比与内定标信号帧经过匹配滤波器得到的分辨率、积分旁瓣比和峰值旁瓣比比较,得到相应差的绝对值,并通过输出接口单元输出到分析结果判读模块;再接收下一帧内定标数据,重复上述操作,直到处理完输入数据。
一种星载SAR内定标信号处理实现方法,其特征在于,包括以下步骤:
步骤1:原始内定标信号按帧的形式由有效数据抽取模块的输入接口单元输入,并成对进入匹配滤波器中进行匹配滤波,取其匹配滤波后功率的最大值之间差的绝对值;根据该值的大小判断是否有效;通过输出接口单元,将信号输出到稳定性定标模块、相位误差定标模块、TR组件频域定标模块、可靠性检验模块。
步骤2:通过稳定性定标模块的输入接口单元的信号帧分别进入幅度特性保存单元和匹配滤波器,进入幅度特性保存单元可将信号帧的信号幅度变换为曲线,并保存在该单元内;进入匹配滤波器的信号滤波后输出到插值单元,并取功率峰值处附近的数据进行插值;幅度变化曲线传递给可靠性检验模块,功率峰值处的相位和幅度传递给分析结果判读模块。
步骤3:读入有效的内定标信号到相位误差定标模块的相位去除单元,除去标准二次相位,剩下的为误差相位并输入到相位拟合单元,获取相位拟合曲线;相位拟合曲线通过相位误差预处理单元,得到相位误差曲线,并输入到相位误差分析单元,获取各高次相位误差的拟合曲线和峰峰值,并把拟合曲线传递给可靠性检验模块,峰峰值传递给分析结果判读模块。
步骤4:相关输出结果输入到可靠性检验模块的误差注入单元,并注入到理想线性调频信号中,得到的合成信号在匹配滤波器中进行匹配滤波,获取其分辨率、积分旁瓣比、峰值旁瓣比与内定标信号按帧经过匹配滤波器得到相应指标进行比较,得到相应差的绝对值,并通过输出接口单元输出到分析结果判读模块。
步骤5:有效的内定标信号输入到TR组件频域定标模块,用来检测TR组件的工作状态,并将处理后的结果和可靠性验证模块验证得出的分析结果一起输出到分析结果判读模块;最终由分析结果判读模块分析所有的输入的处理结果,从而判断雷达的工作状态。
步骤6:接收下一帧内定标数据,重复步骤1至步骤5,直到处理完输入数据。
所述步骤3中相位误差定标模块的相位去除单元中,产生的理想线性调频信号在经过取共轭装置取共轭后,在乘法器中与有效的内定标信号相乘,此时输出信号的标准二次相位被除去,剩下的为误差相位。
所述步骤3中相位误差定标模块中的相位拟合单元的处理步骤如下:
A、首先取相位去除单元的输出信号的第一个采样点,作为相位曲线的第一个点。
B、在取共轭装置中取这第一个采样点的共轭,然后再乘法器中和第二个采样点相乘,取其乘积结果的相位。
C、把这个相位和第一个采样点的相位通过加法器相加,作为相位曲线的第二个点。
D、重复操作,直至所有采样点处理完毕,获得相位拟合曲线。
所述步骤3中相位误差定标模块的相位误差预处理单元中,含有处理芯片可将输入的数据进行勒让德展开,相位拟合曲线通过输入该单元后即进行勒让德展开,获取一次相位误差和常数项;并利用减法器减去一次相位误差和常数项,如果一次相位误差峰值小于一个阈值,则该单元处理结束,进入相位误差分析单元;否则继续进行该拟合曲线的勒让德展开,重复以上操作。
所述步骤4中可靠性检验模块误差注入单元中,由一信号发生装置和一处理芯片组成,处理芯片把稳定性定标模块输出的幅度误差和相位误差定标模块输出的相位误差调制到理想线性调频信号中,得到合成信号。
本发明一种星载SAR内定标信号处理平台***及实现方法的优点在于:
(1)加入了有效数据抽取模块,同时利用自适应相邻对比法技术,大大的提高了对内定标信号的处理效率,效率是传统处理方法的几倍。
(2)在稳定性定标模块中,对数据进行匹配滤波后,采用回路延时抖动自适应处理技术,防止由于通道延时的抖动造成最大值的采样不准;在相位误差定标模块中,采用高精度正交域误差相位分解技术对误差相位进行分解,大大提高了处理的精度。
(3)TR组件频域定标模块中,将数据变化到频率之后,再反映和监测TR组件的性能,对TR组件的监测更加有效、可靠。
(4)稳定性定标模块、相位误差定标模块和可靠性检验这三个模块中的结果相互对比映证,提高了处理结果的可靠性,使整个平台构造更加合理、完善。
附图说明
图1是本发明一种星载SAR内定标信号处理平台***的总体结构图;
图2是本发明一种星载SAR内定标信号处理平台***的有效数据抽取模块结构图;
图3是本发明一种星载SAR内定标信号处理平台***的稳定性定标模块结构图;
图4是本发明一种星载SAR内定标信号处理平台***的相位误差定标模块结构图;
图5是本发明一种星载SAR内定标信号处理平台***的可靠性检验模块结构图;
图6是本发明一种星载SAR内定标信号处理平台***的分析结果判读模块处理流程图;
图7是本发明一种星载SAR内定标信号处理实现方法的流程图;
图8是本发明一种星载SAR内定标信号处理实现方法步骤3中相位误差定标模块的相位拟合单元的处理流程图;
图9a是本发明一种星载SAR内定标信号处理平台***和实现方法的稳定性定标模块功率峰值处的相位变化曲线;
图9b是本发明一种星载SAR内定标信号处理平台***和实现方法的稳定性定标模块功率峰值处的幅度变化曲线;
图10是本发明一种星载SAR内定标信号处理平台***和实现方法的相位误差定标模块输出的各高次相位误差峰峰值变化曲线;
图11a是本发明一种星载SAR内定标信号处理平台***和实现方法的可靠性检验模块的匹配滤波器中得到的分辨率比较图;
图11b是本发明一种星载SAR内定标信号处理平台***和实现方法的可靠性检验模块的匹配滤波器中得到的峰值旁瓣比比较图;
图11c是本发明一种星载SAR内定标信号处理平台***和实现方法的可靠性检验模块的匹配滤波器中得到的积分旁瓣比比较图。
图中:1.有效数据抽取模块   101.输入接口单元
102.匹配滤波器             103.输出接口单元     2.相位误差定标模块
201.输入接口单元           202.相位去除单元     203.相位拟合单元
204.相位误差预处理单元     205.相位误差分析单元 206.输出接口单元
3.稳定性定标模块           301.输入接口单元     302.匹配滤波器      303.插值单元
304.幅度特征保存单元       305.输出接口单元     4.可靠性检验模块
401.输入接口单元           402.匹配滤波器       403.减法滤波器
404.输出接口单元           405.误差注入单元     5.TR 组件频域定标模块
6.分析结果判读模块
具体实施方式
下面将结合附图对本发明作进一步的详细说明。
如图1所示,本发明是一种星载SAR内定标信号处理平台。包括有效数据抽取模块1、分析结果判读模块6、稳定性定标模块3、相位误差定标模块2、TR组件频域定标模块5和可靠性检验模块4。有效数据抽取模块1用自适应相邻对比法将输入平台的原始内定标信号中的有效数据抽取出来,并把这组有效数据分别输出到稳定性定标模块3、相位误差定标模块2、TR组件频域定标模块5和可靠性检验模块4;稳定性定标模块3将内定标数据进行匹配滤波,测定滤波后的信号功率峰值处的幅度和相位,并分别将处理后的相位和幅度信息输出到分析结果判读模块6中,幅度变化曲线输出到可靠性验证模块4中;相位误差定标模块2采用高精度正交域误差相位分解方法对误差相位进行分解,得到正交的各次相位误差的峰值,并分别将处理后的各高次相位误差峰峰值输出到分析结果判读模块6,各高次相位误差和拟合曲线输出到可靠性验证模块4;TR组件频域定标模块5主要检测TR组件的工作状态,并将处理后的结果输出到分析结果判读模块6;可靠性验证模块4功能为验证整个平台分析结果正确与否,并将结果输出到分析结果判读模块6;最终由分析结果判读模块6分析所有的处理结果,从而判断雷达的工作状态。
如图2所示,整个平台的入口是有效数据抽取模块1,数据通过有效数据抽取模块1后,无效内定标数据被滤除,只剩下有效内定标数据,
有效数据抽取模块1由输入接口单元101、匹配滤波器102和输出接口单元103组成。原始内定标信号按帧的形式由输入接口单元101输入,输入后的内定标帧成对进入匹配滤波器102中进行匹配滤波,即第一帧内定标信号和第二帧内定标信号先后由输入接口单元101输出至匹配滤波器102,取其匹配滤波后功率的最大值之间差的绝对值,如果绝对值之差小于0.5dB,说明这两帧数据不处于回路中放大器的线性放大区域,进而判断这两帧为无效信号;如果之差绝对值大于0.5dB,说明这两帧数据处于回路中放大器的线性放大区域,判断这两帧信号为有效信号。如果得出的是有效数据,则通过输出接口单元103,将信号输出到后面的稳定性定标模块3、相位误差定标模块2、TR组件频域定标模块5、可靠性检验模块4;反之,则抛弃这两帧数据。处理完毕后输入接口单元103再将第三帧内定标信号和第四帧内定标信号先后输出至匹配滤波器102,进入下一次处理,并重复操作,直到处理完全部的原始数据。
稳定性定标模块3中,为了提高数据处理的鲁棒性,首先将内定标数据进行匹配滤波,测定滤波后的信号功率峰值处的幅度和相位,分析其随时间的变化趋势。当变化很小时判定该回路工作正常。随着时间、温度的变化,回路的延时可能会有微小变化,致使功率峰值位置抖动,造成位置选取不准确,这里使用了回路延时抖动自适应处理技术。
如图3所示,稳定性定标模块3由输入接口单元301、幅度特性保存单元304、匹配滤波器302、插值单元303和输出接口单元305组成。从稳定性定标模块3的输入接口单元301按帧读入有效的内定标数据。通过输入接口单元301的第一帧信号分两路,第一路进入幅度特性保存单元304,由于理想情况下,幅度为一恒定值,相当于信号的幅度为矩形,则幅度特性保存单元304可将信号帧的幅度变换为曲线,并保存在该单元内;第二路进入匹配滤波器302,进行匹配滤波,滤波后的输出输入到插值单元303进行处理,即取功率峰值处附近的数据进行插值。在插值序列中重新选取功率峰值的位置,得到该位置处的相位,如果相位变化小于0.005rad,则结束这一步的处理,保存功率峰值处的相位和幅度,否则,重复插值操作。输出接口单元305将幅度变化曲线传递给可靠性检验模块4,并把功率峰值处的相位和幅度传递给分析结果判读模块6。
相位误差定标模块2中,使用了高精度正交域误差相位分解技术对误差相位进行分解。
如图4所示,相位误差定标模块2由输入接口单元201、相位去除单元202、相位拟合单元203、相位误差预处理单元204、相位误差分析单元205和输出接口单元206组成。从输入接口单元201按帧读入有效的内定标信号,并将输出按帧输入到相位去除单元202,相位去除单元202的输出信号输入到相位拟合单元203,获取相位拟合曲线。
相位拟合单元203输出的相位拟合曲线输入至相位误差预处理单元204,把预处理后相位误差曲线输入到相位误差分析单元205,利用勒让德展开式,获取各高次相位误差(一般为二次相位误差到五次相位误差)的拟合曲线和各高次相位误差的峰峰值,并把各高次相位误差拟合曲线通过输出接口单元206传递给可靠性检验模块4,把各高次相位误差峰峰值传递给分析结果判读模块6。重复上述操作,直到有效的内定标信号处理完毕。
TR组件频域定标模块5从频域的角度来分析内定标信号,通过把信号转换到频域,分析对比相同的一段微秒级时间内,不同TR组件之间的幅频变化,以及同一TR组件在不同时间内的幅频变化,来检测TR组件的性能。在比较幅频变化时,通常取通带内一小段频谱幅值的平均值进行比较,通过这种方式,在处理速度和精度上都能得到满意的结果。
可靠性检验模块4的功能主要是验证整个平台分析结果正确与否。它通过借鉴点目标评估中的一些指标来对内定标信号匹配滤波后的结果进行分析,获取其分辨率、积分旁瓣比、峰值旁瓣比然后和前面得到的幅度误差和相位误差进行对比验证。
如图5所示,可靠性检验模块4由输入接口单元401、匹配滤波器402、误差注入单元405、减法器单元403和输出接口单元404。通过输入接口单元401,按帧读入有效的内定标数据及稳定性定标模块3和相位误差定标模块2的输出结果,内定标信号帧输入匹配滤波器402,获取其分辨率、积分旁瓣比和峰值旁瓣比;稳定性定标模块3和相位误差定标模块2的输出结果输入到误差注入单元405,vv误差注入单元405由一信号发生装置和一处理芯片组成,处理芯片把稳定性定标模块3输出中的幅度误差和误差相位定标模块2输出的相位误差调制到理想线性调频信号中,并将得到的误差注入模块405输出的合成信号在匹配滤波器402中进行匹配滤波,也同样获取其分辨率、积分旁瓣比、峰值旁瓣比。将这两组指标通过减法器单元403,分别获取分辨率、积分旁瓣比和峰值旁瓣比相应差的绝对值,并将结果通过输出接口单元404,将输出到分析结果判读模块6。接收下一帧内定标数据,重复上述操作,直到处理完所有数据。
分析结果判读模块6接收来自稳定性定标模块3、相位误差定标模块2、TR组件频域定标模块5、可靠性验证模块4的数据流。该模块的功能主要是对处理的结果进行分析,判断雷达的工作状态。
如图6所示,对于稳定性定标模块3检验单元的输出结果,分析各帧数据功率峰值处的相位变化和幅度变化,如果幅度变化波动小于0.5dB则属于正常,相位变化波动小于3°时为正常状态;
对于相位误差定标检验单元2的输出结果,如果二次相位误差小于10°,其他高次相位误差小于6°时为正常状态;
对于TR组件频域定标模块5的输出结果,如果频谱幅度变化波动小于0.5dB,则为正常状态;
根据可靠性检测模块4中得到的两组指标的对比结果,包括:分辨率、积分旁瓣比和峰值旁瓣比,如果对应指标之差在5%以下,说明分析结果正确;否则,检查稳定性定标模块3和相位误差定标模块2。
如果上述指标都合格,则雷达工作状态正常,可进行***传递函数的测量,否则,返回相应的模块进行检查。
一种星载SAR内定标信号处理实现方法,如图7所示,包括以下步骤:
步骤(1):原始内定标信号按帧的形式由有效数据抽取模块1的输入接口单元101输入,并成对进入匹配滤波器102中进行匹配滤波,取其匹配滤波后功率的最大值之间差的绝对值;如果该值很小,可直接去除,不进行后续操作;如果该值很大,判定这两组数据为有效数据,通过输出接口单元103,将信号输出到稳定性定标模块3、相位误差定标模块2、TR组件频域定标模块5、可靠性检验模块4。
步骤(2):通过稳定性定标模块3的输入接口单元301的信号帧分两路,分别进入幅度特性保存单元304和匹配滤波器302,进入幅度特性保存单元304可将信号帧的信号幅度变换为曲线,并保存在该单元内;进入匹配滤波器302的信号滤波后输出到插值单元303,并取功率峰值处附近的数据进行插值;在插值序列中重新选取峰值的位置,得到该位置信号的幅度和相位并保存,然后重复插值操作,直到匹配滤波后的最大值处信号的相位变化在误差允许的微小范围内;输出接口单元305将幅度变化曲线传递给可靠性检验模块4,并把功率峰值处的相位和幅度传递给分析结果判读模块6。
步骤(3):从相位误差定标模块2的输入接口单元201按帧读入有效的内定标信号到相位去除单元202,在相位去除单元202中,标准二次相位被除去,剩下的为误差相位;相位去除单元202的输出信号输入到相位拟合单元203,获取相位拟合曲线;相位拟合曲线通过相位误差预处理单元204,得到相位误差曲线,并输入到相位误差分析单元205,获取各高次相位误差的拟合曲线和峰峰值,并把拟合曲线通过输出接口单元206传递给可靠性检验模块4,把峰峰值传递给分析结果判读模块6。
步骤(4):稳定性定标模块3和相位误差定标模块2的输出结果输入到可靠性检验模块4的误差注入单元405注入到理想线性调频信号中,得到的合成信号在匹配滤波器402中进行匹配滤波,获取其分辨率、积分旁瓣比、峰值旁瓣比与内定标信号按帧经过匹配滤波器得到的分辨率、积分旁瓣比和峰值旁瓣比比较,得到相应差的绝对值,并通过输出接口单元404输出到分析结果判读模块6。
步骤(5):有效的内定标信号输入到TR组件频域定标模块5,用来检测TR组件的工作状态,并将处理后的结果和可靠性验证模块4验证得出的分析结果一起输出到分析结果判读模块6;最终由分析结果判读模块6分析所有的输入的处理结果,从而判断雷达的工作状态。
步骤(6):接收下一帧内定标数据,重复步骤(1)至步骤(5),直到处理完输入数据。
步骤(3)中相位误差定标模块2的相位去除单元202中,产生的理想线性调频信号在经过取共轭装置取共轭后,在乘法器中与有效的内定标信号相乘,此时输出信号的标准二次相位被除去,剩下的为误差相位。
如图8所示,步骤(3)中相位拟合单元203的处理步骤如下:
A、首先取相位去除单元203的输出信号的第一个采样点,并取其相位,作为相位曲线的第一个点。
B、在取共轭装置中取这第一个采样点的共轭,然后再乘法器中和第二个采样点相乘,取其乘积结果的相位。
C、把这个相位和第一个采样点的相位通过加法器相加,作为相位曲线的第二个点。
D、重复操作,直至所有采样点处理完毕,获得相位拟合曲线。
步骤(3)中相位误差定标模块2的相位误差预处理单元204中含有处理芯片可将输入的数据进行勒让德展开,相位拟合曲线通过输入该单元后即进行勒让德展开,获取一次相位误差和常数项;并利用减法器减去一次相位误差和常数项,如果一次相位误差峰值小于一个阈值(一般为1度),则该单元处理结束,进入相位误差分析单元205;否则继续进行该拟合曲线的勒让德展开,重复以上操作。
相位误差预处理单元204中,勒让德展开式可按下式操作:
勒让德展开式为
其中,pn(t)为勒让德多项式的归一化形式,其中T为每帧信号持续时间,t为采样间隔时间,前五项为:
P 0 ( T ) = 1 T P 1 ( T ) = 3 T ( 2 t T ) P 2 ( T ) = 5 T ( 6 t 2 T 2 - 1 2 )
P 3 ( T ) = 7 T ( 20 t 3 T 3 - 3 t T ) P 4 ( T ) = 9 T ( 70 t 4 T 4 - 15 t 2 T 2 + 3 8 )
P 5 ( T ) = 11 T ( 252 t 5 T 5 - 15 t 3 T 3 + 15 t 4 T ) ,
因此,要求出误差曲线,只要求出系数bn即可。为了求出bn,可以利用展开式的正交性,即
∫ - T 2 T 2 P i ( t ) P j ( t ) dt = 0 i = j 1 i ≠ j , 则:
Figure A20071012210600169
这样就可以求出误差曲线,及相应误差的峰峰值。
各次相位误差峰峰值的表达式如下,其中Wj代表峰峰值,单位为弧度。
W 1 = b 1 × 3 T W 2 = b 2 × 5 T × 1.5 W 3 = b 3 × 7 T × 2
W 4 = b 4 × 1 T × 30 7 W 5 = b 5 × 11 T × 2 .
步骤(4)中可靠性检验模块4的误差注入单元405由一信号发生装置和一处理芯片组成,处理芯片把稳定性定标模块3输出中的幅度误差和误差相位定标模块2输出的相位误差调制到理想线性调频信号中,得到合成信号。
本实施例中,将一组内定标信号输入本发明中的平台,相位误差定标模块2中的相位拟合单元203中处理芯片型号为TMS320F28015,可靠性检验模块4的误差注入单元405中的处理芯片型号为TMS320LC549-80,随机取其中取相邻10帧内定标数据进行分析。
如图9a、图9b所示,为稳定性定标模块3功率峰值处的相位变化曲线和功率峰值处的幅度变化曲线。从图中可以看出,相位的变化范围为108.8°~109.3°之间,小于0.05°;幅度的变化范围为41.2dB~41.25dB之间,小于0.5dB。说明稳定性定标模块3输出的功率峰值处的幅度和相位的变化很小,即雷达***的工作状态十分的稳定。
从这10帧内定标数据中随机取一帧的进行相位分解,如图10所示,相位误差定标模块2输出的各高次相位误差峰峰值变化曲线,为内定标信号相位误差正交分解后得到的一次相位误差到五次相位误差。从图中看出,它们的峰峰值的变化范围为都在-1°~1°之间,都小于2°,说明各高次相位误差峰峰值相位误差变化很小,即雷达***工作正常。
如图11a、图11b、图11c所示,分别为这10帧内定标数据在可靠性检验模块4的匹配滤波器402中得到的分辨率、峰值旁瓣比和积分旁瓣比比较图。虚线分别表示稳定性定标模块3和相位误差定标模块2的输出结果输入到误差注入单元405注入到理想线性调频信号后,得到的合成信号经匹配滤波得到的分辨率、积分旁瓣比、峰值旁瓣比,实线分别表示内定标信号帧经匹配滤波得到的分辨率、积分旁瓣比和峰值旁瓣比。从图中可以看出,两者都基本重合。说明稳定性定标模块3和相位误差定标模块2分析得到的结果是正确的。附图中的曲线分别存在着一定的差异是由于:在相位分解中,本发明只分解处理了一次相位误差到五次相位误差的结果,但实际上还存在着更高次,即六次以上的相位误差,由于高次相位误差一般很小,故一般忽略不予以考虑,或者当随机噪声处理。

Claims (10)

1.一种星载SAR内定标信号处理平台***,包括TR组件频域定标模块、分析结果判读模块,其特征在于,还包括稳定性定标模块、相位误差定标模块、有效数据抽取模块和可靠性检验模块;
有效数据抽取模块用自适应相邻对比法将输入平台的原始内定标信号中的有效数据抽取出来,并把这组有效数据分别输出到稳定性定标模块、相位误差定标模块、TR组件频域定标模块和可靠性检验模块;稳定性定标模块测定匹配滤波后的信号功率峰值处的幅度和相位,反映雷达***的工作状态,并分别将处理后的相位和幅度信息输出到分析结果判读模块中,幅度变化曲线输出到可靠性验证模块中;相位误差定标模块采用高精度正交域误差相位分解方法对误差相位进行分解,得到正交的各次相位误差的峰值,并分别将处理后的各高次相位误差峰峰值输出到分析结果判读模块,各高次相位误差拟合曲线输出到可靠性验证模块;TR组件频域定标模块主要检测TR组件的工作状态,并将处理后的结果输出到分析结果判读模块;可靠性验证模块功能为验证整个平台分析结果正确与否,并将结果输出到分析结果判读模块;最终由分析结果判读模块分析所有的处理结果,从而判断雷达的工作状态。
2.根据权利要求1所述的一种星载SAR内定标信号处理平台***,其特征在于:所述有效数据抽取模块由输入接口单元、匹配滤波器和输出接口单元组成;
原始内定标信号按帧的形式由输入接口单元输入,并成对进入匹配滤波器中进行匹配滤波,取其匹配滤波后功率的最大值之间差的绝对值;如果该值很小,可直接去除,不进行后续操作;如果该值很大,判定这两组数据为有效数据,通过输出接口单元,将信号输出到稳定性定标模块、相位误差定标模块、TR组件频域定标模块、可靠性检验模块;处理完毕后进行下一对数据帧的处理,并重复操作,直到处理完全部的原始数据。
3.根据权利要求1所述的一种星载SAR内定标信号处理平台***,其特征在于:所述稳定性定标模块由输入接口单元、幅度特性保存单元、匹配滤波器、插值单元和输出接口单元组成;
通过输入接口单元的信号帧分两路,分别进入幅度特性保存单元和匹配滤波器,进入幅度特性保存单元可将信号帧的信号幅度变换为曲线,并保存在该单元内;进入匹配滤波器的信号滤波后输出到插值单元,并取功率峰值处附近的数据进行插值;在插值序列中重新选取峰值的位置,得到该位置信号的幅度和相位并保存,然后重复插值操作,直到匹配滤波后的最大值处信号的相位变化在误差允许的微小范围内;输出接口单元将幅度变化曲线传递给可靠性检验模块,并把功率峰值处的相位和幅度传递给分析结果判读模块。
4.根据权利要求1所述的一种星载SAR内定标信号处理平台***,其特征在于:所述相位误差定标模块由输入接口单元、相位去除单元、相位拟合单元、相位误差预处理单元、相位误差分析单元和输出接口单元组成;
从输入接口单元按帧读入有效的内定标信号到相位去除单元,在相位去除单元中,标准二次相位被除去,剩下的为误差相位;相位去除单元的输出信号输入到相位拟合单元,获取相位拟合曲线;相位拟合曲线通过相位误差预处理单元,得到相位误差曲线,并输入到相位误差分析单元,获取各高次相位误差的拟合曲线和峰峰值,并把拟合曲线通过输出接口单元传递给可靠性检验模块,把峰峰值传递给分析结果判读模块;重复上述操作,直到有效的内定标信号处理完毕。
5.根据权利要求1所述的一种星载SAR内定标信号处理平台***,其特征在于:所述可靠性检验模块通过借鉴点目标评估中的指标来对内定标信号匹配滤波后的结果进行分析,由输入接口单元、匹配滤波器、误差注入单元、减法器单元和输出接口单元组成;
稳定性定标模块和相位误差定标模块的输出结果输入到误差注入单元注入到理想线性调频信号中,得到的合成信号在匹配滤波器中进行匹配滤波,获取其分辨率、积分旁瓣比、峰值旁瓣比与内定标信号按帧经过匹配滤波器得到的分辨率、积分旁瓣比和峰值旁瓣比比较,得到相应差的绝对值,并通过输出接口单元输出到分析结果判读模块;再接收下一帧内定标数据,重复上述操作,直到处理完输入数据。
6.一种星载SAR内定标信号处理实现方法,其特征在于,包括以下步骤:
步骤一:原始内定标信号按帧的形式由有效数据抽取模块的输入接口单元输入,并成对进入匹配滤波器中进行匹配滤波,取其匹配滤波后功率的最大值之间差的绝对值;根据该值的大小判断是否有效;通过输出接口单元,将信号输出到稳定性定标模块、相位误差定标模块、TR组件频域定标模块、可靠性检验模块;
步骤二:通过稳定性定标模块的输入接口单元的信号帧分别进入幅度特性保存单元和匹配滤波器,进入幅度特性保存单元可将信号帧的信号幅度变换为曲线,并保存在该单元内;进入匹配滤波器的信号滤波后输出到插值单元,并取功率峰值处附近的数据进行插值;幅度变化曲线传递给可靠性检验模块,功率峰值处的相位和幅度传递给分析结果判读模块;
步骤三:读入有效的内定标信号到相位误差定标模块的相位去除单元,除去标准二次相位,剩下的为误差相位并输入到相位拟合单元,获取相位拟合曲线;相位拟合曲线通过相位误差预处理单元,得到相位误差曲线,并输入到相位误差分析单元,获取各高次相位误差的拟合曲线和峰峰值,并把拟合曲线传递给可靠性检验模块,峰峰值传递给分析结果判读模块;
步骤四:相关输出结果输入到可靠性检验模块的误差注入单元,并注入到理想线性调频信号中,得到的合成信号在匹配滤波器中进行匹配滤波,获取其分辨率、积分旁瓣比、峰值旁瓣比与内定标信号按帧经过匹配滤波器得到相应指标进行比较,得到相应差的绝对值,并通过输出接口单元输出到分析结果判读模块;
步骤五:有效的内定标信号输入到TR组件频域定标模块,用来检测TR组件的工作状态,并将处理后的结果和可靠性验证模块验证得出的分析结果一起输出到分析结果判读模块;最终由分析结果判读模块分析所有的输入的处理结果,从而判断雷达的工作状态;
步骤六:接收下一帧内定标数据,重复步骤一至步骤五,直到处理完输入数据。
7.根据权利要求6所述的一种星载SAR内定标信号处理实现方法,其特征在于:
所述步骤三中相位误差定标模块的相位去除单元中,产生的理想线性调频信号在经过取共轭装置取共轭后,在乘法器中与有效的内定标信号相乘,此时输出信号的标准二次相位被除去,剩下的为误差相位。
8.根据权利要求6所述的一种星载SAR内定标信号处理实现方法,其特征在于:所述步骤三中相位误差定标模块中的相位拟合单元的处理步骤如下:
A、首先取相位去除单元的输出信号的第一个采样点,作为相位曲线的第一个点;
B、在取共轭装置中取这第一个采样点的共轭,然后再乘法器中和第二个采样点相乘,取其乘积结果的相位;
C、把这个相位和第一个采样点的相位通过加法器相加,作为相位曲线的第二个点;
D、重复操作,直至所有采样点处理完毕,获得相位拟合曲线。
9.根据权利要求6所述的一种星载SAR内定标信号处理实现方法,其特征在于:所述步骤三中相位误差定标模块的相位误差预处理单元中,含有处理芯片可将输入的数据进行勒让德展开,相位拟合曲线通过输入该单元后即进行勒让德展开,获取一次相位误差和常数项;并利用减法器减去一次相位误差和常数项,如果一次相位误差峰值小于一个阈值,则该单元处理结束,进入相位误差分析单元;否则继续进行该拟合曲线的勒让德展开,重复以上操作。
10.根据权利要求6所述的一种星载SAR内定标信号处理实现方法,其特征在于:所述步骤四中可靠性检验模块的误差注入单元中,由一信号发生装置和一处理芯片组成,处理芯片把稳定性定标模块输出的幅度误差和误差相位定标模块输出的相位误差调制到理想线性调频信号中,得到合成信号。
CN2007101221069A 2007-09-21 2007-09-21 一种星载sar内定标信号处理平台***及实现方法 Expired - Fee Related CN101135726B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2007101221069A CN101135726B (zh) 2007-09-21 2007-09-21 一种星载sar内定标信号处理平台***及实现方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2007101221069A CN101135726B (zh) 2007-09-21 2007-09-21 一种星载sar内定标信号处理平台***及实现方法

Publications (2)

Publication Number Publication Date
CN101135726A true CN101135726A (zh) 2008-03-05
CN101135726B CN101135726B (zh) 2010-11-03

Family

ID=39159907

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2007101221069A Expired - Fee Related CN101135726B (zh) 2007-09-21 2007-09-21 一种星载sar内定标信号处理平台***及实现方法

Country Status (1)

Country Link
CN (1) CN101135726B (zh)

Cited By (20)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101236247B (zh) * 2008-03-07 2010-06-16 北京航空航天大学 一种星载多通道天线sar数据通道幅相误差校正平台
CN101915908A (zh) * 2010-06-29 2010-12-15 西安空间无线电技术研究所 一种星载微波散射计多耦合回路内定标方法
CN102183743A (zh) * 2011-03-04 2011-09-14 北京航空航天大学 一种长波长星载ctlr模式简缩极化sar的定标方法
CN102253374A (zh) * 2011-04-20 2011-11-23 北京航空航天大学 一种长波长星载π/4模式简缩极化SAR的定标方法
CN101788663B (zh) * 2009-01-23 2012-07-04 中国科学院电子学研究所 一种星载合成孔径雷达***预失真方法
CN101571593B (zh) * 2008-04-30 2012-11-21 北京航空航天大学 一种星载sar图像的严密共线方程模型
CN103091666A (zh) * 2011-11-07 2013-05-08 中国科学院电子学研究所 非理想定标器条件下的机载p波段极化sar定标方法
CN103760534A (zh) * 2014-01-19 2014-04-30 中国人民解放军国防科学技术大学 一种星载sar数据的电离层色散效应校正方法
CN105892304A (zh) * 2014-11-05 2016-08-24 中国空间技术研究院 曲线数据自动判读方法
CN106526588A (zh) * 2016-12-02 2017-03-22 上海无线电设备研究所 一种太赫兹主动成像雷达***相位补偿方法
CN107271994A (zh) * 2017-05-16 2017-10-20 上海卫星工程研究所 基于内定标数据的星载合成孔径雷达参考函数获取方法
CN108562880A (zh) * 2018-04-26 2018-09-21 西安空间无线电技术研究所 一种反射面星载sar***内定标网络单元和内定标方法
CN109884631A (zh) * 2019-02-21 2019-06-14 中国科学院电子学研究所 星载合成孔径雷达内定标数据处理的方法及***
CN110018460A (zh) * 2019-04-11 2019-07-16 上海卫星工程研究所 星载合成孔径雷达整星阶段板间相位差的远场测量方法
CN110146856A (zh) * 2019-04-11 2019-08-20 上海卫星工程研究所 星载sar内定标数据处理流程设计方法
CN111398921A (zh) * 2020-04-22 2020-07-10 中国人民解放军63816部队 一种减小雷达光电差标定不确定度的方法
CN112698283A (zh) * 2020-12-11 2021-04-23 中国科学院空天信息创新研究院 雷达测试***、方法、信号发生设备及信号反馈设备
CN113723225A (zh) * 2021-08-13 2021-11-30 陕西航天技术应用研究院有限公司 一种新型的定标信号处理方法、***、电子设备及存储介质
CN115629363A (zh) * 2022-12-16 2023-01-20 西安空间无线电技术研究所 星载探测雷达实时内定标方法及装置
CN116184340A (zh) * 2023-04-27 2023-05-30 中国科学院空天信息创新研究院 一种分布式合成孔径雷达验证***及方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1148067C (zh) * 2001-06-01 2004-04-28 中国科学院计算技术研究所 合成孔径雷达复数图像数据压缩方法

Cited By (29)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101236247B (zh) * 2008-03-07 2010-06-16 北京航空航天大学 一种星载多通道天线sar数据通道幅相误差校正平台
CN101571593B (zh) * 2008-04-30 2012-11-21 北京航空航天大学 一种星载sar图像的严密共线方程模型
CN101788663B (zh) * 2009-01-23 2012-07-04 中国科学院电子学研究所 一种星载合成孔径雷达***预失真方法
CN101915908A (zh) * 2010-06-29 2010-12-15 西安空间无线电技术研究所 一种星载微波散射计多耦合回路内定标方法
CN101915908B (zh) * 2010-06-29 2012-08-22 西安空间无线电技术研究所 一种星载微波散射计多耦合回路内定标方法
CN102183743A (zh) * 2011-03-04 2011-09-14 北京航空航天大学 一种长波长星载ctlr模式简缩极化sar的定标方法
CN102253374A (zh) * 2011-04-20 2011-11-23 北京航空航天大学 一种长波长星载π/4模式简缩极化SAR的定标方法
CN103091666A (zh) * 2011-11-07 2013-05-08 中国科学院电子学研究所 非理想定标器条件下的机载p波段极化sar定标方法
CN103091666B (zh) * 2011-11-07 2015-02-04 中国科学院电子学研究所 非理想定标器条件下的机载p波段极化sar定标方法
CN103760534A (zh) * 2014-01-19 2014-04-30 中国人民解放军国防科学技术大学 一种星载sar数据的电离层色散效应校正方法
CN105892304A (zh) * 2014-11-05 2016-08-24 中国空间技术研究院 曲线数据自动判读方法
CN105892304B (zh) * 2014-11-05 2019-05-31 中国空间技术研究院 曲线数据自动判读方法
CN106526588B (zh) * 2016-12-02 2019-04-26 上海无线电设备研究所 一种太赫兹主动成像雷达***相位补偿方法
CN106526588A (zh) * 2016-12-02 2017-03-22 上海无线电设备研究所 一种太赫兹主动成像雷达***相位补偿方法
CN107271994A (zh) * 2017-05-16 2017-10-20 上海卫星工程研究所 基于内定标数据的星载合成孔径雷达参考函数获取方法
CN107271994B (zh) * 2017-05-16 2019-08-02 上海卫星工程研究所 基于内定标数据的星载合成孔径雷达参考函数获取方法
CN108562880A (zh) * 2018-04-26 2018-09-21 西安空间无线电技术研究所 一种反射面星载sar***内定标网络单元和内定标方法
CN108562880B (zh) * 2018-04-26 2021-03-26 西安空间无线电技术研究所 一种反射面星载sar***内定标网络单元和内定标方法
CN109884631B (zh) * 2019-02-21 2020-11-10 中国科学院电子学研究所 星载合成孔径雷达内定标数据处理的方法及***
CN109884631A (zh) * 2019-02-21 2019-06-14 中国科学院电子学研究所 星载合成孔径雷达内定标数据处理的方法及***
CN110018460A (zh) * 2019-04-11 2019-07-16 上海卫星工程研究所 星载合成孔径雷达整星阶段板间相位差的远场测量方法
CN110146856A (zh) * 2019-04-11 2019-08-20 上海卫星工程研究所 星载sar内定标数据处理流程设计方法
CN110146856B (zh) * 2019-04-11 2023-04-14 上海卫星工程研究所 星载sar内定标数据处理流程设计方法
CN111398921A (zh) * 2020-04-22 2020-07-10 中国人民解放军63816部队 一种减小雷达光电差标定不确定度的方法
CN112698283A (zh) * 2020-12-11 2021-04-23 中国科学院空天信息创新研究院 雷达测试***、方法、信号发生设备及信号反馈设备
CN113723225A (zh) * 2021-08-13 2021-11-30 陕西航天技术应用研究院有限公司 一种新型的定标信号处理方法、***、电子设备及存储介质
CN115629363A (zh) * 2022-12-16 2023-01-20 西安空间无线电技术研究所 星载探测雷达实时内定标方法及装置
CN116184340A (zh) * 2023-04-27 2023-05-30 中国科学院空天信息创新研究院 一种分布式合成孔径雷达验证***及方法
CN116184340B (zh) * 2023-04-27 2023-07-11 中国科学院空天信息创新研究院 一种分布式合成孔径雷达验证***及方法

Also Published As

Publication number Publication date
CN101135726B (zh) 2010-11-03

Similar Documents

Publication Publication Date Title
CN101135726B (zh) 一种星载sar内定标信号处理平台***及实现方法
CN106646546A (zh) 一种卫星信号多维快速捕获方法及***
CN105259410B (zh) 一种强噪声干扰下的欠采样波形的频率估计方法及其装置
CN103941254A (zh) 一种基于地质雷达的土壤物性类别识别方法和装置
CN106443203A (zh) 一种脉冲信号检测***及方法
CN104391278A (zh) 一种利用极化对消的雷达抗干扰方法
EP2207044B1 (en) A signal phase difference measurement method, apparatus and system
CN103983617A (zh) 一种基于小波变换改进激光探针定量分析的方法
US20120098924A1 (en) Signal spectra detection system
CN101436879B (zh) 环境噪声下干扰信号幅度的提取方法
CN107561357A (zh) 一种基于信道化的高精度瞬时测频方法和装置
CN104375132A (zh) 数字阵列雷达多路模拟通道相对延时测量装置及方法
CN103117973A (zh) 基于正交检波的多路并行解调***
CN104483668A (zh) 一种高精度雷达信号检测和跟踪***及其方法
CN104392086A (zh) 一种基于皮尔逊秩次变量相关系数的信号检测电路及方法
Huang et al. Efficient angular chirp-Fourier transform and its application to high-speed target detection
CN107192689B (zh) 一种基于多尺度太赫兹光谱的原包装奶粉无损检测方法
CN204008440U (zh) 用于差分气体检测的信号处理装置及其气体检测***
CN107688188A (zh) 一种gps l1c/a与l1c信号联合捕获方法和卫星导航接收机
CN109495195A (zh) 无线电通信pcma信号幅度的联合估计方法及装置
CN104333423A (zh) 一种基于基尼相关的信号检测电路及方法
CN110083865B (zh) 一种抗脉冲噪声的已知信号检测方法及检测电路
CN101572564A (zh) 弱相关非高斯环境下基于局部最佳检测器的伪码捕获方法
US9800317B2 (en) Method for the detection of an electromagnetic signal by an antenna array with pattern diversity, and device implementing said method
CN103248368A (zh) 一种判断随机解调器压缩采样重构成败的方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
C17 Cessation of patent right
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20101103

Termination date: 20130921