基于过零点斜率加窗解调的ACARS接收方法研究

专利名称：基于调频斜率精搜索的低信噪比信号参数提取方法专利类型：发明专利
发明人：张磊,孙希平,张曼
申请号：CN201810144257.2
申请日：20180212
公开号：CN108415013B
公开日：
20220503
专利内容由知识产权出版社提供
摘要：本发明公开了一种基于调频斜率精搜索的低信噪比信号参数提取方法，属于雷达技术领域，其主要思路为：获取雷达脉冲信号，根据所述雷达脉冲信号得到数字信号，并根据数字信号确定参考信号，然后基于调频斜率精搜索得到最终的调频斜率估计值；根据参考信号和最终的调频斜率估计值，得到参考信号的时间间隔，并将所述参考信号的时间间隔作为雷达脉冲信号的脉冲重复时间估计值；根据所述数字信号和所述雷达脉冲信号的脉冲重复时间估计值，得到高信噪比的积累信号；根据高信噪比的积累信号，得到雷达脉冲信号的脉宽估计值、雷达脉冲信号的带宽估计值和雷达脉冲信号的载频，记为基于调频斜率精搜索的低信噪比信号参数提取结果。

申请人：西安电子科技大学,西安中电科西电科大雷达技术协同创新研究院有限公司
地址：710071 陕西省西安市太白南路2号
国籍：CN
代理机构：西安睿通知识产权代理事务所(特殊普通合伙)
代理人：惠文轩
更多信息请下载全文后查看。

基于ACARS数据链的MSK解调方法抗干扰性能研究作者：王敏杰文荣辉陈昱均来源：《现代信息科技》2020年第06期摘要：根據敌我识别系统中对ACARS数据链识别解调的需求，分析了ACARS的基带信号特点。

针对ACARS基带信号的MSK解调易受噪声和频偏影响的问题，对3种常用的MSK 解调方法进行了比较分析与仿真测试。

在结构分析中可以看出数字差分解调的构成易于工程实现，适合于FPGA上实现高速MSK解调。

解调正确率的测试结果表明，数字差分解调方法相较于几种常规MSK解调方法，其误码率更低，抗频偏能力更强。

关键词：ACARS;MSK;数字差分解调;频偏;误码率中图分类号：TN971.+1;TN919.2 文献标识码：A 文章编号：2096-4706（2020）06-0061-05Abstract：According to the requirement of ACARS data link identification and demodulation in IFF system，the baseband signal characteristics of ACARS are analyzed. Aiming at the problem that the MSK demodulation of ACARS baseband signal is easily affected by noise and frequency offset，three kinds of commonly used MSK demodulation methods are compared，analyzed and simulated. In the structural analysis，it can be seen that the structure of DDM is easy to implement in engineering，and it is suitable for high-speed MSK demodulation on FPGA. The test results of demodulation accuracy show that compared with several conventional MSK demodulation methods，the digital differential demodulation method has lower bit error rate and stronger anti frequency offset ability.Keywords：ACARS;MSK;digital differential demodulation;frequency offset;bit error rate0 引言ACARS数据链信号是敌我识别系统中一项重要的检测信号，它是由美国航空无线电公司研发的一种用于飞机与地面站双向通信的甚高频数据链，用以保障飞机的飞行安全和全球范围的空中交通管制[1]。

第34卷第2期指挥控制与仿真V ol.34 No.2 2012年4月 Command Control & Simulation Apr. 2012 文章编号:1673-3819(2012)02-0109-05基于频谱分析的雷达信号分选新方法孙鹏1,2,潘谊春1,郁春来1,吴彩华1(1.空军雷达学院,湖北 武汉 430019;2.95865部队,北京 100021)摘 要:相参性是现代雷达的普遍性质,利用雷达脉冲之间的相参性,提取相参特征参数,并运用于雷达信号分选,对降低当前信号分选中存在的“漏批”率具有重要的意义。

因此在基于频谱分析的相参特征模型基础上,研究了基于最小二乘模板匹配误差的相参特征,并提出了具体应用中的改进措施。

仿真结果表明,利用频谱主峰图像最小二乘模板匹配误差这一相参特征参数进行信号分选,可以有效地减少“漏批”率。

关键词:相参;信号分选;频谱分析;最小二乘模板匹配误差中图分类号:TP391.9 文献标识码:A DOI: 10.3969/j.issn.1673-3819.2012.02.026A Novel Sorting Method of Radar Signal based on Spectrum AnalysisSUN Peng1,2,PAN Yic-hun1,YU Chun-lai1,WU Cai-hua1(1. AFRA, Wuhan 430019; 2. the Unit 95865 of PLA, Beijing 100021,China)Abstract: Coherent is an important characteristics of the modern radar, by use of the coherent specific character between the Radar pulses, abstract coherent specific character parameter, and apply the parameter to the Radar signal sorting, is extremely valuable to reduce the probability of leak-alarm. this paper extracts the spectrum peak image of the reconstruct signal by the way of spectrum, and introduce a coherent characteristic parameter based on Least-Square Template Matching Error(LSTME). Theoretical analysis and simulation both indicate that the spectrum peak image of the reconstruct signal is one of the manifestations of coherent characteristics, use the coherent characteristic parameter to signal Sorting can reduce the probability of lead-alarm.Key words: coherent; signal sorting; spectrum analysis; LSTME信号分选是侦察系统中的关键处理进程,其分选识别水平是衡量信息处理技术先进程度的重要标志。

第17卷　第4期 地　球　物　理　学　进　展 V ol.17　N o.4 2002年12月(605～613) PROG RESS　I N　GE OPHY SICS Dec.2002接收函数方法及研究进展许卫卫　郑天愉(中国科学院地质与地球物理研究所,北京100101)[摘　要]　远震P波波形数据中包含了大量在台站下方地壳上地幔速度间断面所产生的P2S转换波及其多次反射波的信息,是研究台站下方局部区域S波速度分布理想的震相,由此产生的接收函数方法是反演台站下方S波速度结构的有效手段.接收函数方法可以通过波形反演拟合接收函数的径向分量,对观测台站下方地球介质的S波速度结构进行估计,也可以通过偏移叠加获得的接收函数道集(地震剖面图)追踪速度间断面.这种方法避免了对天然地震震源及其附近结构混响效应等复杂因素的影响,对S波速度的垂向分布敏感,垂向分辨率高.由于宽频带流动地震台阵的发展,用此方法还可获得研究区域速度结构的横向变化,横向分辨能力主要取决于台站的间距.本文回顾20年来接收函数研究的进展,探讨了方法研究的发展趋势,介绍了对地壳—上地幔结构的部分研究结果.[关键词]　接收函数;S波速度结构;速度间断面;台阵[中图分类号]　P315 [文献标识码]　A [文章编号]　100422903(2002)04206052090　引　言用地震波资料来探测地球内部结构是现代地震学的主要任务之一.近10年来,利用宽频带地震台阵记录到的远震资料获得地壳和地幔各间断面结构的接收函数方法发展迅速,获得一系列重要成果.早在1979年Langston就研究了这个问题[1].他给出了源等效假定,从长周期远震体波中分离出接收函数.Owens等(1984)把这一方法扩展到了宽频带记录,并发展了接收函数的线性反演方法,研究了坎帕兰郡高原的构造[2].K ind等(1995)研究了G DS N (全球数字化地震台网)台站下的地壳结构,并对接收函数方法做了综和阐述[3].为了提高机构成像结果的分辨率,Y uan等(1997)发展了接收函数的偏移叠加方法,并用以研究上地幔间断面的横向变化[4],他们用此方法研究青藏高原和安第斯地区[5]的地壳上地幔结构获得了成功.在国内,刘启元在宽频带地震台阵观测和接收函数方法研究方面做了许多工作,如接收函数的非线性反演和合成三维横向非均匀介质远震体波接收函数的方法等[6—8].吴庆举等[9]发展了时间域的反褶积方法,用接收函数研究了青藏高原的地壳结构.尽管接收函数方法早在上个世纪70年代就已经提了出来,但由于当时台站稀少、所用资料为人工数字化得到的长周期数据,很难得到研究地区下的精细结构,与其他方法相比优势并不明显[10—12].近年来,随着数字化观测技术及流动台站观测的迅速发展,接收函数方法开始被广泛采用,研究地震台阵下方地壳和上地幔结构的结果层出不穷[13—15],接收函数逐渐成为研究地壳上地幔结构的一种有效手段.[收稿日期]　2002202219;　[修回日期]　2002208210.[基金项目]　国家基础研究发展规划项目(G1999043301)资助.[作者简介]　许卫卫,男,1978年4月生,陕西眉县人,2000年毕业于北京航空航天大学,现为中国科学院地质与地球所博士研究生,主要从事流动台站的观测及地球内部结构的研究.(Email:w wxu@)1　接收函数方法接收函数方法主要包括接收函数的提取和接收函数的反演两大部分.接收函数的提取实质是一个反褶积的过程,得到所需的接收函数;接收函数的反演采用理论图与观测图均方误差最小原则,用依据理论介质模型得到的理论接收函数拟合观测得到的接收函数,反演得到台站下方的S 波速度结构.为了得到更高的分辨率,要对直接提取的接收函数进行偏移叠加.1.1　接收函数的提取在时间域,三分量远震P 波波形数据可表示为仪器响应、震源时间函数及介质结构响应的卷积D V (t )=I (t )3S (t )3E V (t ),D R (t )=I (t )3S (t )3E R (t ),D T (t )=I (t )3S (t )3E T (t ).(1)D (t )为远震P 波波形记录,I (t )为仪器的脉冲响应,S (t )为入射平面波的有效震源时间函数,E (t )为介质结构响应,下标V ,R ,T 分别表示垂直,径向,切向分量.理论计算与实际观测表明,远震P 波波形的垂直分量主要由直达波构成,尾随波列能量较弱,可忽略不记,故可令:E V (t )≈δ(t ).即D V (t )=I (t )3S (t ).将此式代入(1)式并变换到频率域可以得到E R (ω)=D R (ω)D V (ω),E T (ω)=D T (ω)D V (ω).(2) 若将E R (ω)、E T (ω)再变换到时间域得到的E R (t )、E T (t ),即被称为接收函数.在实际计算中,由于(2)式中的分母有可能为零或趋于零,导致频率域的除法不稳定,特做如下处理:E R (ω)∧=D R (ω)D 3V (ω)Φ(ω)G (ω).(3)其中:Φ(ω)=max {D 3V (ω)D V (ω),c max [D 3V (ω)D V (ω)]},G (ω)=exp (-ω2/4α2).0<c <1.Φ(ω)称为“水准值”,由经验选择,与数据噪声有关,并对结果的分辨率有重要影响.G (ω):高斯滤波器,压制高频分量,α由经验选择.噪声的不同频率分量可能并不相同,而上式中的c 不随频率分量变化,这样会降低接收函数的分辨率.为了处理这一问题,刘启元等提出了接收函数复谱比的最大或然性估计方法[6],吴庆举等在时间域应用了Wiener 滤波反褶积和最大熵谱反褶积算法[9],经证实都是行之有效的.如果射线路径相同,相应的接收函数应该也是一致的.因此,有选择的进行接收函数的叠加可减弱随机干扰的影响,增强P 2S 转换波信号.另一方面,不同震中距的地震,直达P 波与P 2S 转换波之间走时差不同,需要做偏移处理.这是个对地震图拉伸或者压缩的处理过程,会造成地震波形一些的畸变,但实践证明,与其增强信号的能力相比,是完全可略的.・606・ 地　球　物　理　学　进　展 17卷K ind 等在1996年对西藏的研究中采用了单台叠加的方法[16].对某一台站,选取一个参考震中距进行偏移处理之后,以直达P 波为参照,对接收函数进行叠加,得到此台站下方的平均接收函数.依次类推,可得到所有台站下方的平均接收函数.然而上面的方法产生了一个新的问题.在单台叠加中,信号产生的噪音是相关的,叠加会生成一些虚假界面.为此,Y uan 等在1997年对接收函数方法做了一些改进[4].他们对每道记录独立进行标准化后,按照台站、方位排列起来,形成一个接收函数道集,之后对所有台站的全部接收函数叠加,形成研究区域下结果的一个粗略面貌.结合接收函数道集剖面图和叠加曲线对结果进行解释.为了研究该区域下的横向变化,用相邻台站移动平均的办法来提高台站间的相关性.具体做法如下:假设有A1、A2、…A10共10个台站,对相邻台站A1、A2、A3的接收函数进行平均得到假想台站B1的接收函数,对A2、A3、A4进行平均得到B2,依次类推得到各假想台站下的接收函数图.实践证明这种方法比单台叠加更有优越性,可以更清楚、全面地看到研究区域下的情况.同年,Dueker 和shechan 提出了接收函数共中心点叠加方法(亦称共面元叠加)[17],这也是目前运用较多的一种方法.根据资料分辨率和研究目标选择一大小适中的面元,在待研究的某一深度将研究区域划分成相邻(可重叠)的许多面元,对观测到的地震追踪穿过该区域的P 2S 转换波射线,对转换点落在同一面元上的接收函数做偏移叠加,得到该面元在该深度叠加的接收函数.对于各种深度均可做上述处理,从而得到整个区域的接收函数.为了简化处理,我们可以采用沿深度分段的做法,对一定的深度范围,采用一个平均深度来进行偏移叠加.这种方法实际上是一种多台叠加方法,能够把单台叠加中的相关噪音通过多台不相关来消除.1.2　接收函数反演从偏移叠加的接收函数道集,通过追踪大振幅的震相,在时间域或深度域获得速度间断面.对于研究区域下的速度结构,需要进行接收函数反演来建立.反演过程中,先根据对研究区域已有的认识给出S 波速度结构初始模型,以理论地震图与观测地震图之间的均方根误差最小为准则,对S 波速度进行修改,P 波速度根据S 波速度推算得到.拟合远震P 波波形的方法多种多样.Phinney (1964)首次给出了远震P 波振幅谱的谱幅比拟合的结果[18],之后Burdick 和Langston (1977)指出了在频率域拟合的固有缺点(第一,从时间域到频率域转换过程中的截断误差不可避免;第二,特定震相的信息不够明显),并把波形拟合引入到时间域[19].Owens 等(1984,1987)利用了宽频带P 波波形展示的更多细节,进一步发展了时间域波形反演[2,20].在此基础上,Amm on 等(1990)[21]引入了Randall (1989)计算微分地震图的有效算法[22]以及Shaw 和Orcutt (1985)提出的“跳动”反演技术[23],加快了反演的速度,并且对反演得出的速度模型进行了光滑度的约束,这就是现在被广泛使用的时间域线性化反演技术.由于计算效率的提高,Amm on 对反演结果的非唯一性也作了研究,揭示了高、低速层以及它们的厚度对反演结果的影响.他指出,要减轻反演结果的非唯一性,除了使用更高分辨率的数据外,与其他手段(比如反射、折射方法,地球化学、地质学方法)所得出的结果进行比较也是非常必要的.上面提到的反演方法属于线性化反演,它要求充分接近“真实”的初始模型,这使得已有的接收函数反演技术有赖于其他方法所获得的先验模型,为克服上述困难,刘启元等在・706・4期 许卫卫:接收函数方法及研究进展 1996年提出了频率域的非线性反演方法[6].他们利用Sham way 和Der 的多道最大或然性反褶积原理[24],讨论了从单台三分量远震P 波波形中分离接收函数径向与垂向分量复谱比的方法,根据T arantola 的波形反演理论[25]发展了接收函数复谱比的非线性反演方法.在反演中除了拟合接收函数的复谱比之外,还需拟合时间域中接收函数垂向和径向分量的初至振幅比.用理论合成的波形数据进行数字实验的结果表明,这种方法的反演结果不依赖于初始模型的选择.接收函数反演基本按以下步骤进行:(1)数据准备.包括数据格式转换、滤波等工作.(2)坐标系的转换.把地震记录由地理直角坐标系变换到射线坐标系.(3)用地震记录的垂直分量对径向和切向分量进行反褶积,得到相应的接收函数.(4)选定一个参考震中距,对接收函数做偏移校正.(5)对接收函数进行叠加.(6)反演.(7)对结果进行解释.2　接收函数的应用接收函数方法经过20多年的发展,已逐渐成熟,在地壳上地幔速度结构研究中占有着重要的地位.近年来,应用接收函数方法研究地震台阵观测资料获得的主要结果,表现在以下几个方面.2.1　板块构造研究板块活动在地球内部留下的构造信息,是地壳上地幔速度结构研究的重要目标.为探寻青藏高原隆升和巨厚地壳形成的机制,1992年开始启动I NDEPTH 计划(西藏与喜马拉雅地区深部探测国际计划),开始了宽频带流动地震台阵观测.此期后第二和第三期探测也分别在1994年和1997年完成.从1992年持续至今的流动地震台阵观测,为青藏高原研究积累了大量的资料.K ind 等(1996)首先利用I NDEPTH 2Ⅱ的数据,用接收函数方法对藏布(西藏南部印度板块与欧亚板块的碰撞区)缝合带做了研究[16].结果表明,沿藏布缝合带莫霍面连续分布,在其北部发现了地壳低速带,这和宽角反射研究得到的强烈P 2S 转换波、地磁研究中的高导层这些所谓的“亮点”是一致的.Y uan 等(1997)[4]采用了包括I NDEPTH 2Ⅱ/GE DEPTH (西藏与喜马拉雅地区德国测深计划)和G S N (全球地震台网)拉萨台的数据资料,同时在接收函数方法上做了较大的改进.Y uan 等的研究结果表明,沿整个测线莫霍面深度在70km 到80km 之间,在藏布缝合带以北50km 处地下10到20km 深度发现显著的低速区,这被解释为由部分熔融引起.另有一个重要的结果就是在整个地区之下50～60km 深度存在“下地壳”,这个间断面在Owens 等1997年对该地区所做的研究中也可看到.吴庆举等(1998)根据PASSC A L (地壳与岩石层的地震台阵研究计划)流动地震台阵观测实验所得的资料对青藏高原的巨厚地壳做了研究[9].他们利用时间域广义线性反演的跳动算法,引入模型光滑度约束,并将合成地震图的K ennet 算法[26]及微分地震图的Randall 快速算法[22]用于接收函数的正演计算,由台站接收函数获得了各台站下方的一维S 波速度分布.反演结果表明,青藏高原的莫霍界面在班公2怒江缝合带附近存在明显的深度错断,可能的解释是班公2怒江缝合带是印度地壳向欧亚下地壳挤入的前沿;在日喀则、拉萨、桑雄和安东多等地的地壳内部,可以连续观测到三个显著的速度界面,存在壳内低速层.K osarev 等(1999)[14]集中利用上述多期观测获得的759道接收函数,进一步对该地区做了更深入的研究,得到一个十分引人注意的结果.结果显示了在从藏布缝合线以北50km 处・806・ 地　球　物　理　学　进　展 17卷的莫霍面(80km 深)有一向北倾斜到200km 的界面(藏布转换边界);而在高原北部,有一分段南倾的结构.他们认为藏布转换边界是印度岩石圈顶部与亚洲岩石圈底部的界面.中安第斯山脉是海洋板块与大陆板块相互作用形成的最高山脉.这种抬升运动发生于20Ma 年前,主要是大陆地壳加厚的结果,且有明显的构造削减作用.Y uan 等(2000)利用中安第斯山脉当地地震台网记录到的远震P 2S 转换波,得到了与削减过程紧密相关的深层结构图[27].他们发现莫霍界面在阿蒂普拉诺高原深75km ,在海拔高达4000m 的普纳高原(秘鲁的干冷高原)却只有50km.高海拔却又相对薄的地壳厚度说明普纳高原的上抬引起岩石圈变薄.纳斯卡海洋板块向下俯冲到120km 的深度,他们的研究在这儿却没有观测到远震转换波,分析认为可能是因为辉长岩向榴辉岩过渡已经完成,这是俯冲带中存在动力变质反应的直接证据,也是对大多数中等深度地震波在120km 深度俯冲板块中消失的可能解释.研究获得的存在于整个阿蒂普拉诺高原和普纳高原下的10～20km 厚的壳内低速带,被解释为下地壳向上地壳过渡的连续变质和部分熔融带.2.2　上地幔间断面和热点的研究热点和上地幔410km 和660km 间断面,是人们探寻地表构造运动的深部热动力学机制的重要观测信息.现有的重力观测资料和地球化学资料已经证实了热柱的存在,但是地震学在这方面的资料相对而言比较缺乏.尤其对于热柱的起源问题,一直是争论的焦点.由于410km 和660km 界面岩石组分的不同,导致两个界面的热学性质截然不同,410km 界面深度随温度升高而增加,660km 界面则正好相反.于是便可通过该地区410km 和610km 界面的研究来推测热柱的起源.全球的热点以冰岛、夏威夷为典型代表.体波走时层析成像发现,在冰岛中部低速分布延伸到了至少400km 的深度,但对于更深的深度则显得无能为力[28,29].Shen 等(1997)运用接收函数方法进行了研究[30].他们在冰岛布设了NW —EW 、W —E 、N —S 三条测线,通过对该地区的410km 界面和660km 界面横向变化的研究来分析热柱起源,主要用到的震相为P 410S 和P 660S.他们的研究发现,冰岛中部和南部下方410km 界面与660km 界面之间的过渡区比全球平均值薄20km 左右,温度高出约150K,但是在它的周围区域为正常值.这些结果和熔融物产生深度的热异常估计也非常吻合,此结果有力地证明了热柱来源于下地幔.Li 等(2000)对夏威夷热柱问题用转换地震波也做了研究[31].结果表明,夏威夷群岛下部中心地带130～140km 处存在一个剪切波低速带(<4km ),并且一直延伸到整个上地幔.在夏威夷南2西南部,410km 界面与660km 界面之间的过渡区厚度比平均值小40～50km ,其间的温度高出近300℃.热柱起源类似于冰岛的情形,不过温度要高出冰岛100℃左右.2.3　典型区域构造研究地壳和壳幔边界的细结构,是获取深部信息、开展综合性的天然实验室研究的重要组成部分.随着地震观测能力的提高,人们开始利用密集分布的台阵观测,提高探测地下结构的分辨能力,以实现对典型构造区域的综合研究.地球上最古老的大陆壳壳核一般都形成于太古代.津巴布韦和喀布瓦尔(一北一南)克拉通共同形成了非州南部的陆核,这些地区未受到后来一些构造运动的影响,为太古代的原始地壳.喀布瓦尔克拉通北部有世界闻名、全球最大的成层铁镁质侵入体,岩石学研究认为它是经历元古代几次大规模事件后形成的,在陆核的周围还有一些小的区域(如奥克瓦区,・906・4期 许卫卫:接收函数方法及研究进展 玛高恩迪区,喀易司区,纳马夸兰2纳塔尔区)均是在后太古代时期形成,夹在津巴布韦克拉通和喀布瓦尔克拉通之间的是太古代晚期这两个克拉通固化时形成的碰撞岩层———林波波河带.Nguuri 等(2001)利用1997年至1999间在该地区布设的81个宽频带地震仪,以及3个全球数字台站的资料,对该地区地下1～101km 深度范围的结构进行了研究[32].结果表明,太古代时期形成的津巴布韦和喀布瓦尔克拉通莫霍面非常清晰,界面速度差最少在1km/s ,喀布瓦尔克拉通的莫霍面在34～37km 范围内;津巴布韦克拉通稍微复杂一些,莫霍面深度范围在33～48km 内,平均38km ,在与南非灌木丛区、奥克瓦/玛高恩迪区相邻的地方深度偏大.后太古代时期形成的地壳和未经干扰的克拉通相比厚度要大,并且模糊不清,上地幔速度偏低,一个可能的解释是这和太古代地壳受到的古生代岩浆上覆有关.这些都与层析成像、地质学得到的结果比较相符,进一步加深了人们对该地区的认识.利用南加利福尼亚地区大量宽频带三分量地震台站积累的资料,Zhu 等(2000)用接收函数方法对该地区的莫霍面深度变化进行了研究[33].研究中发现,接收函数可以很好地实现台站下方莫霍面深度的“点”测量,结果对地壳P 波速度不太敏感,主要依赖于S 波速度,或者说是依赖于V P /V S .由P 2S 转换波测定的莫霍面深度与V P /V S 值关联性强,他们利用较晚到达的转换震相P p P s 及P p S s +P s P s 做为约束来解决这个问题.以地壳厚度H 和P 、S 波速度比k =V P /V S 为自变量构造函数s (H ,k )=w 1r (t 1)+w 2r (t 2)-w 3r (t 3),其中r (t )是径向接收函数,t 1、t 2、t 3分别为由H 和V P /V S 决定的转换震相P t S 、P p P s 及P p S s +P s P s 的到时(由初始模型预先估计),w i 为加权系数且∑w i =1,当H 和V P /V S 正确时,这三个震相相关叠加,s (H ,k )可达到最大值.此方法省去了对转换震相到时的拾取,可以方便地处理大量资料,结果的可信度还可以通过对s (H ,k )最大值的方差估计得知.研究结果表明,整个区域莫霍面深度在21km 到37km 之间,平均29km.在东横断山、半岛山脉、内华达山脉莫霍面较深,加利福尼亚内部和萨儿顿山谷较浅.在莫哈韦沙漠西部和中部莫霍面比较平坦,接近于平均值,向东则变浅.Zhu 等(2000)还利用这些观测的远震P 2S 转换波接收函数,获得了穿过圣安德烈斯断层的地壳结构[34].结果表明,在安德烈斯断层两侧莫霍面连续分布,平均深度30km ;安德烈斯断层之下莫霍面错开6～8km.在东加州剪切带下,莫霍面也有间断,这被解释为北美板块与太平洋板块边界向东扩展形成的新的边界.这些都表明圣安德烈斯断层和东加州剪切带贯穿了整个地壳.同时还发现,圣加百利山下的莫霍面上隆,这和重力观测以及地质学方面的研究成果都较为相符,据分析是水平挤压的结果.2.4　地震构造研究延怀盆地位于汾渭断陷盆地带的北段,由雁行式排列的涿鹿、怀来、延庆三个新生代裂谷组成,是一个地质结构复杂的交汇区.为研究延怀盆地及其临近区域内的地震构造,刘启元等(1997)[7]在沙城以东以延庆盆地为中心大约46000km 2范围内布设了15个测点组成的G DS 21000宽频带流动地震台阵,台站间距5～10km ,跨越了观测区域各主要构造边界,大体构成东西向和南北向交汇的测线.他们对方位不同、震中距变化不大的接收函数进行了平滑,压制了横向结构不均匀造成的散射效应,并且进一步降低了背景噪音.研究结果表明,延怀盆地临近区上地壳普遍存在低速层,这就为进一步研究该区域内地震孕育的机理提供了依据.・016・ 地　球　物　理　学　进　展 17卷1997年新疆伽师地区30km 范围内连续发生了7次6级以上强烈地震,引起了地震学家的广泛关注.刘启元等(2000)运用与延怀盆地研究中类似的方法对该地区进行了研究[8],得出了以下结论:(1)伽师强震群区的莫霍界面深度在塔里木盆地一侧为40～52km ,而在天山一侧的褶皱变形区为60～76km ,南天山山前褶皱变形区的地壳厚度明显增大,地壳内部产生明显变形;(2)塔里木盆地北缘存在明显的不均匀性,台阵下方壳幔界面的三维地形与地表的地貌特征有较好的相关性,塔里木盆地与天山之间的接触变形关系显示了塔里木盆地向北北西方向的挤压作用,山前弧形褶皱带是在特定的地下深部及上述挤压作用下形成的结果,块体深部构造边界的不均匀性在碰撞造山带可能具有某种普遍性.(3)伽师强震群分布在壳幔界面梯度带的上方,该震群的成因可能与震源区附近的附伏断裂活动有关,这可能是塔里木盆地北缘的非均匀变形及现金构造运动造成的结果.(4)由于震源区处于褶皱变形区的沉积层向结晶地壳的过渡区,相应的介质剪切模量较小,这可能是伽师震群应力降明显偏低的主要原因.3　讨论随着流动台阵观测资料的增加,新的研究结果将不断涌现.同时,如何从这些地震记录中更有效地获取更多来自地球内部的信息,亦是一个研究主题.目前接收函数方法主要集中于对P 2S 转换波震相的分析来研究地壳上地幔结构,即P 波接收函数方法.选用P 2S 震相,主要是因为该震相初动尖锐,能量较强.但因地壳多次反射波与P 2S 到时相当,降低了接收函数的分辨率,人们也开展了对其它类型震相的接收函数研究[35],这其中又以S 2P 为典型代表.S 2P 震相走时比P 波的地壳多次反射波晚,但比S 波的地壳多次反射波要早,可以克服P 2S 在这方面的缺点.但它的初动比较微弱,需要研究相应的信号提取技术.P 2S 震相反演的得到的是S 波速度结构,而S 2P 波反演得到的是P 波速度结构,同时利用P 2S 及S 2P 转换波,对P 波和S 波速度进行联合反演,可以更准确地反映地球内部介质的物性特征,是有发展潜力的研究工作.在反演方法上,由于频率域固有的缺点,主要还是时间域的线性反演,其结果对初始模型有依赖性.台阵观测可以获得空间分布的大量接收函数道集,有待发展相应的非线性反演方法.在叠加手段上,可以借鉴已发展较为成熟的地震勘探的许多方法.目前接收函数叠加中所用的偏移实质上是走时的动校正,叠加只是在共转换面元上的叠加,代表一种平均效应,与勘探地球物理中所用的方法还存在许多不同点,这主要受地震数量与台站间距的制约.但随着观测周期的增加与台站间距的减少,将勘探地球物理中的方法技术(如抽道集、速度分析、叠前偏移等)引入到接收函数的偏移叠加中是一种新的发展方向[36],它会提高处理结果的精度和信噪比.观测地震图中含有大量理论地震图中所没有的信息,在直接从接收函数道集中追踪界面信息时,可结合对观测地震图中这些信息的判断,选择合适的方法和资料加以连续追踪可以得到许多有意义的结果[37].另外,利用其他资料或方法得到的结果,约束接收函数反演存在的解答非唯一性问题,是正在探寻的手段,目前提出的区域P 波接收函数方法是一个新的思路.借助多锥形谱分析方法,利用区域P 波可以得到许多关于地壳结构的信息,和远震接收函数方法结合起来对其可起到很好的约束和指导作用[38].・116・4期 许卫卫:接收函数方法及研究进展 。

一种改进的基于串行干扰抵消的Rake接收机杨朝建;高勇【摘要】在基于串行干扰抵消的Rake接收机的基础上作了改进,提出了一种新的Rake接收机结构,对其进行了原理和性能分析,给出了Matlab仿真结果,并与基于串行干扰抵消的Rake接收机性能进行了比较.仿真结果表明,改进的基于串行干扰抵消的Rake接收机性能明显得到提高,同时减小了误码率.【期刊名称】《现代电子技术》【年(卷),期】2007(030)017【总页数】3页(P26-28)【关键词】多径信号;多径干扰;Rake接收机;串行干扰抵消;最大比值合并【作者】杨朝建;高勇【作者单位】四川大学,电子信息学院,四川,成都,610065;四川大学,电子信息学院,四川,成都,610065【正文语种】中文【中图分类】TN851 引言信号经过无线信道传播，不可避免地会产生多径效应，使通信系统的性能明显下降。

Rake接收机是一种有效的抗多径技术，他充分利用了各条多径信号能量以获得分集增益。

对于Rake接收机的原理和性能，文献[1-5]都有详细的分析和推导。

从中可以看出各种传统的Rake接收机对多径干扰的抑制效果没有达到期望值。

因此，研究人员提出了一些改进的Rake接收机，如：RLS-Rake接收机[3] 、Pre-Rake[5] 、采用串行干扰抵消的Rake接收机结构[6] 。

相比于传统的Rake接收机，他们的性能有了很大的提高。

文献[6]中的接收机结构，假定各条多径信号的功率与时延之间有确定的关系，也就是说时延越小信号的功率就越大。

但实际上，多径传播的损耗包含了路径上的损耗、物体的吸收损耗，因此，多径时延越小，多径的功率就越大这种对应关系并不严密。

在Cost-231模型中，根据物体的材质不同，损耗是不同的。

例如：经过加钢筋的混凝土墙的损耗是单纯的混凝土、砖墙的3倍，是有大量窗户的墙的6倍[7]。

因此文献[6]中的Rake接收机结构并不能保证消除来自功率最大的一条多径信号的干扰，而对于各条多径信号而言，受到的最大的干扰往往来自功率最大的那条多径信号。

一种保留相位谱的DS信号接收方法
王永民;郭建新
【期刊名称】《空军工程大学学报（自然科学版）》
【年(卷),期】2009(010)004
【摘要】为抑制接收DS信号中的窄带干扰,研究了基于DFT只保留相位谱的DS 信号接收方法,给出了该方法的一种迭代实现算法,有效降低了相位谱计算的运算量.讨论了该算法在DS系统抑制窄带干扰中的应用原理,分析了该算法在加性高斯白噪声信道中的性能损失,给出了一种计算性能损失的方法.对采用相位谱保留算法接收DS信号进行了数值仿真,给出了在窄带干扰和多址干扰情况下,该算法与常规DS 信号接收方法的对比数值仿真结果.仿真结果表明:采用文中的算法处理后,可有效抑制窄带干扰,在多址干扰情况下,保留相位谱处理后没有引起明显的性能损失.
【总页数】6页(P48-52,57)
【作者】王永民;郭建新
【作者单位】空军工程大学,电讯工程学院,陕西,西安,710077;空军工程大学,电讯工程学院,陕西,西安,710077
【正文语种】中文
【中图分类】TN973.2
【相关文献】
1.一种提高色谱指纹谱保留时间重现性的新方法 [J], 王龙星;肖红斌;梁鑫淼
2.一种基于频差补偿的相位谱时延估计方法 [J], 邱天爽;尤国红;沙岚;赵小平;高阳
3.DDS相位截断杂散谱精确分析方法的改进 [J], 袁福;郑林华;袁继兵
4.一种互谱相位在时延估计中的应用方法 [J], 丁超;陈喆;张宗堂;程玉胜
5.强地震动场相位谱的一种工程预测方法 [J], 丁海平;廖振鹏
因版权原因，仅展示原文概要，查看原文内容请购买。