0 引言

二次雷达系统（Second Surveillance Radar，SSR）发展到今天，已经成为空中交通管制系统的基础 ［1］，其系统中的信号模式也随着发展更新换代。最早的二次雷达系统包括模式1、2、3，随后系统功能不断扩充、不断改进后增加了模式A、C、4，其中模式1、2、3/A、C称为选择性特征识别（Selectivity Identification Function，SIF）模式。模式 1、2、3/A 用于身份认证，模式C提供飞机的飞行高度，在询问时根据需要可只发出SIF模式中的某一种模式信号，亦可几种模式信号交错发送［2］。而最新的二次雷达系统向上兼容，在原系统基础上增加了模式S。模式S大量装备用于民用航管二次雷达，包括点名询问和数据链通信等功能。随着日益加剧的空中流量压力，空间信号环境日益复杂，二次雷达系统在接收应答信号时，易产生多个应答机的应答信号相互混扰现象。而传统的解码方法难以分辨这些在时域和频域相互重叠的信号，严重影响了信号的模式识别和解码，降低了SSR系统的安全性和可靠性［3］。

基于曼彻斯特编码特性的分选算法利用了模式S应答信号的编码特性，可在多个模式S应答信号混扰时进行有效分离。然而，在实际环境中，混扰的二次雷达系统应答信号中常常不仅仅包含模式S应答信号，还会出现SIF模式信号，在这种情况下，如果直接使用基于曼彻斯特编码特性的分选算法对混扰信号进行分选，则SIF模式的应答信号会被忽略。因为SIF模式未采用曼彻斯特编码，即不满足该算法假设的曼彻斯特编码特性［6］。

从中国刑警学院回来后，王敬凯继续深入梳理案情，以取得更多的证据。 在那封匿名信及3张纸条中，多次出现“瓶装”这个常人颇感生疏的特殊词组。

本文针对曼彻斯特算法在混扰信号中掺杂SIF模式信号时不能有效分选信号的问题，提出一种优化方法。首先根据信号模式特征之间的不同，对混扰信号进行针对性预处理。然后通过求解得出分离矩阵，以实现对混扰信号中的SIF模式部分进行有效分选。然后对非SIF模式类信号进行分选。最后通过仿真验证了本文算法的有效性。

1 应答信号模型

1.1 应答信号格式

SIF模式的应答信号如图1所示，包括2个脉冲间隔为20.3 μs的框架脉冲F1和F2、12个数据脉冲、一个恒为0的备用位脉冲X以及一个特殊位置编码脉冲SPI，SPI一般不用，当两目标互相接近或者应答码相同时，可要求其中一个目标增加一个SPI脉冲以便识别。脉冲宽度为0.45 μs，相邻脉冲的间隔为1.45 μs，脉冲F2与SPI的间隔为4.35 μs，当SPI被激活时，该脉冲持续时间可达到15 s～30 s。信号经PAM调制后可以写作：

假设天线阵元数为M，进入接收机的应答信号个数为K且K＜M，采样点数为N，则观测到的混扰矩阵Y可以表示为：

 

式中，Ri（t）代表各个模式的应答信号；i=1，2，3/A，C代表信号模式；ai，1=1，ai，15=1 代表框架脉冲；ai，8=0，代表X 脉冲；T0=1.45 μs，代表相邻脉冲间隔；pR（t）代表宽度为0.45μs的矩形脉冲；PSPI（t）表示SPI脉冲；TSPI=24.65 μs。

模式S的应答信号包括4个脉宽为0.5 μs的前导脉冲和一个长度为56 μs（包含56 bit数据）或112 μs（包含112 bit数据）的应答数据块，如图2所示。数据采用曼彻斯特编码，即前0.5 μs为高电平、后0.5 μs为低电平时表示数据1，反之则表示数据0。经PAM调制后可表示为

 

  

图1 SIF模式应答信号

  

图2 模式S应答信号

式中，pS，R（t） 代表宽度为 0.5 μs 的矩形脉冲；TS，k（k=1，2，3，4） 分别取 0 μs、1 μs、3.5 μs、4.5 μs，TS=8 μs，TS，R=0.5 μs。

1.2 应答信号混扰模型

 

式中，为混扰信号矩阵；为引导矩阵；为应答信号矩阵；为噪声矩阵。

应答信号矩阵S可表示为：

 

式中表示矩阵Hadamard积。若询问机的接收天线采用均匀线性阵列天线，阵元间距为半波长，则阵列引导矢量ak为：

 

式中，θk为应答信号sk相对于阵列法线的入射角度。然而，在现实中引导矩阵A会受到多方面因素影响，如不精确的阵列天线刻度、天线耦合、阵列混扰等。

宋代的中央集权，大体上在“专务以矫失为得”的原则上建立起来的。政权之主要注意力集中于如何防范文武重臣的篡夺之祸；如何防止人事、财政、军政等大权旁落；如何禁制百官间凭藉种种因缘相互朋比，以致构成为专制政权的离心力量；……任人格局中“且要异论相搅，即各不敢为非”的良苦用心，诸般设施中所贯彻的文武相制、内外相维、上下相轧原则，都分明体现出“事为之防，曲为之制”的精神。［3］520

最后在固定信噪比为15 dB、两个信号的DOA为-3°和5°，变化两应答信号的相对延时，变化区间取0 μs至20 μs。4种算法分选的误码率如图11所示。可见，由于PA算法核心点在于利用两信号非重叠部分提取特征向量构造分离矩阵，故在相对时延较小时产生较高误码率；ESPRIT算法和ICA算法的性能因其算法特性几乎不受相对延时影响；本文算法始终保持较低误码率，但随着相对延时增大误码率稍有增加。

 

式中，分离矩阵WH近似为引导矩阵A的伪逆，即。

2 改进算法简介与流程

2.1 针对性预处理

在针对信号的预处理中，算法需要初步估计信源个数以及噪声模型建立。信源个数K可采用最小描述长度（Minimum Description Length，MDL）准则进行估计，该算法从信息论的角度出发，可对信源数进行有效估计，即

 

式中，λi为样本协方差矩阵的第i个大的特征值。同时，噪声矩阵，各向量服从零均值，方差为σ2的白复高斯分布且相互独立，且噪声与信源相互独立。

对观测到的混扰信号矩阵Y作奇异值分解（SVD）分解有

 

式中，均为酉矩阵；，为混扰矩阵Y的奇异值，且排序为。

令

取wk为分离矩阵WSIF的列向量，则，可有

 

 

式中，UK、VK分别为 U、V 的前 K 列，ΛK为 K×K 对角阵，为混扰信号矩阵Y的行空间Row（Y）的一组基。

在我国跨境电商物流的发展过程中，基础设施的不完善及体系不完整等，物流的成本较高，根据相关资料得出结论，从成本来看，我国跨境电商物流的成本远远高于其他国家甚至是国际标准。跨境电商物流需要考虑的问题多而复杂，涉及国内物流、国内海关、 国际运输、 国外海关、 国外物流众多环节，在实际操作中风险很大，因此，在一定程度上大大增加了我国跨境电商的物流成本。

定义

 

预处理后，维数等于应答信号个数K。

2.2 核心方法

算法的主要思想是，先从混扰信号中分选出SIF模式应答信号，将其滤除，再将基于曼彻斯特编码特性的分选算法应用于剩余混扰信号进行分选。

他朝着地上“呸呸呸”吐了三口，嘴里骂了几句，来给自己壮胆，然而未等骂完，忽然闻到了空气中似乎有什么味道。

分析SIF模式应答信号，可以得到如下特性

 

式中，TSIF=1.45 μs为SIF模式应答信号脉冲间隔。

农民通过对其土地的占有，在土地上进行相关的生产和生活，也通过生产取得的利益满足他们衣食住行以及娱乐生活的需要，这是农民对占有土地的主要需求。农民的土地收入权也可以分为增值收益的权利和获得孽息的权利。

对接收到的混扰信号R（t）进行采样，得到混扰信号矩阵YSIF

 

第二，在实施机制方面，管理机制对于学校导向的产学研协同育人模式促进最大。产学研合作主体中的一方学校没有建立起严格而科学的管理机制，没有在前期计划、过程组织、人员激励及结果控制中发挥管理职能。高校产学研协同育人大多停留在企业、科研单位签订产学研共建合作协议，企业、学校和科研单位都没有从产学研协同育人模式中获益。

 

所以

 

由kronecker积的性质，可得

 

定义矩阵PSIF有如下形式

 

所以有

 

若令，则有

 

显然在P的零空间内，即。对P作奇异值分解分解可得

 

求得分离矩阵WSIF后，可令矩阵QSIF为矩阵WSIF的列向量张成子空间的正交投影矩阵。那么矩阵为只包含模式S应答信号的混扰信号矩阵，可以通过曼彻斯特编码分选算法，对混扰信号矩阵进行分选。

 

可见为矩阵E的特征向量，可通过特征值分解求得，进而得到分离矩阵WSIF。则应答信号S的估计值为

令，，则

 

今天，我们主要来了解一下中国古代最具代表性的十大乐器，分别是 ：琵琶、二胡、编钟、箫、笛、瑟、琴、埙（）、笙、鼓。

评价是合作学习中不可缺少的重要环节。在每次小组活动后，教师要及时进行总结评价，一方面通过小组汇报的形式使学生对知识形成统一的认识，得出结论；另一方面，对各个学习小组的活动情况进行评价，如哪一组最团结协作，哪一组最有创新精神等，对小组成员也要进行适当合理的评价，比如谁在合作学习中最积极，谁能用正确的方式向其他学生提出自己不同的意见等，让学生尝到合作学习的甜头，享受学习的快乐。我们还可以在教室后面的黑板报上，营造精神家园，让小组之间的竞争看得见。小组之间合作好便给他们贴上一颗星星，每个月评比一次，这样更能激发学生的积极性。

显然，对于B0的每一列有

 

将 bi由 K2×1 的向量变形为 K×K 的矩阵˙i，即，则有

 

式中为对角矩阵。

可以定义矩阵˙0为

 

对矩阵作分解，可分为K个K×K的方阵，每个方阵可表示为

 

令

余荫山房中植物景观与园内的建筑关系紧密，其植物景观轴线的设计起着强化和补充园中建筑轴线的作用。余荫山房主要的一条植物景观轴线由“深柳堂”前的百年炮仗花与百年榆树的对植构成，两侧加盆栽以强化对称感，此为南北景观轴线，与园中自南向北的“临池别馆”、方形水池与“深柳堂”形成的建筑轴线相呼应[12]。余荫山房的另一条植物景观轴线由“玲珑水榭”左右两侧9株对列而植的桂花构成，此为东西景观轴线，与园中“玲珑水榭”、八角形水池及“烷红跨绿桥”垂直方向构成的建筑轴线相合。

 

令，易得

式中，，R 为 K维可逆矩阵。

 

式中，对角矩阵。

式中，VP（i）表示VP的第i列。

 

在对分选得到的应答信号S˜解码前，首先要对其取模并将其幅值归一化，以使其幅值相对稳定，然后采用6 dB阈值法等传统方法进行解码。

取VP的后K个列向量，则构成P零空间的一组基，即

宁波市江北区存量土地资源盘活实践与探索（何建军） .................................................................................9-17

2.3 分选步骤

曼彻斯特改进算法分选步骤如下：

Step1：通过MDL准则估计应答机个数K；

Step2：针对SIF模式信号特性进行预处理得到Y'；

Step3：由Y'构建矩阵P；

Step4：对矩阵P做奇异值分解，求解其零空间的一组基B0；

Step5：求解非对称非正交联合对角化问题，得到分离矩阵W；

Step6：求解应答信号S的估计值；

Step7：利用6 dB阈值法等方法解码出SIF模式信号；

Step8：通过曼彻斯特编码分选算法求解Y'中剩余模式S信号。

3 仿真分析

与曼彻斯特算法相似，影响该改进算法分选效果的因素主要有信噪比、不同应答信号的DOA差值与信号延时，现分别从这几个影响因素出发将本文算法分选性能与ESPRIT算法、ICA算法和PA算法对比。现假设测试中阵列天线阵元数为M=4，以2 MHz的采样率对下变频后的信号进行采样，仿真结果均取1 000次独立蒙特卡罗实验的平均值并进行比对。

3.1 分选过程仿真

假设混扰信号中包含1个模式S信号与1个SIF模式信号，应用本算法进行分选，效果如下。固定两个信号的DOA为-3°和5°，相对延时为30 μs，信噪比为15 dB。图3为混扰前的原始应答信号，图4为噪声环境下的混扰信号。可见两个信号在时域上重叠，如果直接对其解码会产生大量误码。图5为本文算法分选得到的两个应答信号。由图5对比图4可见，两信号成功分离，解码后与原信号码元对比，计算出误码率为1.79%，通过计算得到DOA估计值为-2.934 6°和5.041 7°。

  

图3 原始应答信号

  

图4 混扰信号

  

图5 算法信号分选结果

若设定成2个模式S信号与1个SIF模式信号进行混扰，本算法依然有效果。固定两个信号的DOA 为 -3°、1°和 5°，延时分别为 10 μs、20 μs与50 μs，信噪比为15 dB。图6为混扰前的原始应答信号，图7为噪声环境下的混扰信号。图8为本文算法分选得到的应答信号。由图8对比图7可见，信号成功分离，解码后与原信号码元对比，计算出误码率为1.93%，通过计算得到DOA估计值为-3.021 6°、1.070 9°和 4.998 3°。

  

图6 原始应答信号

  

图7 混扰信号

  

图8 算法信号分选结果

  

图9 不同信噪比下的误码率比较

3.2 性能分析

性能分析时将采用1个模式S信号与1个SIF模式信号混扰模型进行分析。首先，固定两个信号的DOA为-3°和5°、相对延时为10 μs时，变化信噪比从0 dB～30 dB。从图9中可以直观看出，随着信噪比增加，信号分选误码率明显降低。信噪比低于15 dB时，本文算法误码率明显低于另外3种算法，在低信噪比情况下保持较好分选效果；信噪比大于15 dB后，4种算法误码率均趋于零，均能对混扰信号实现完全分离。不难得出结论，在DOA估计值方面改进算法与PA算法、ESPRIT算法性能相似，均可实现有效DOA估计，而ICA算法则因其原理不同而表现最差。

在数据收集后，本研究采用方差膨胀因子(VIF)检测是否存在多重共线性问题。结果显示，所有参数估计的VIF均小于2(Max=1．88)，说明在高校创新能力与高校绩效的回归模型中，不存在多重共线性问题。

这种理论在比较法上的支持是日本的实践。日本高等法院在1995年BBS Wheels案中，曾认为日本传统的权利用尽原则的合理性依据即在于专利法的平衡原则，当他售出商品时即已经获得了补偿，因而不应二次获利。这同样适用在专利权人在国外销售产品的情况[11]。

接下来固定信噪比为15 dB、两信号相对延时为10 μs，且固定其中一个应答信号的DOA为0°时，另一个应答信号的DOA值相对变化，变化区间取0°～12°。4种算法分选的误码率如图10所示。易见，在DOA差值较小时，4种算法的分选性能均相对较低，但本文算法的分选性能受DOA差值影响相对最小。

  

图10 不同DOA差值下的误码率比较

为了实现信号的分选，需要找到一个分离矩阵W，得到应答信号S的估计值

  

图11 不同相对延时下的误码率比较

4 结论

本文针对曼彻斯特算法在混扰信号中掺杂SIF模式信号时不能有效分选信号的问题，提出一种优化方法。首先根据信号模式特征之间的不同，对混扰信号进行针对性预处理。然后通过求解得出分离矩阵，以实现对混扰信号中的SIF模式部分进行有效分选。最后对非SIF模式信号进行分选。仿真结果表明，该算法在低信噪比情况下保持较好的分选性能。并且相较其他算法而言，由于本算法内容的针对性，分选性能具有效果更佳。但与曼彻斯特算法类似，当应答信号相对延时较大时，本文算法的误码率略有增加，此问题将在下一步的工作中研究解决。

参考文献：

［1］张尉.二次雷达原理［M］.北京：国防工业出版社，2009：9-12.

［2］王洪，刘昌忠，汪学刚.二次雷达S模式综述［J］.电讯技术，2008，48（7）：113-118.

［3］GALATI G，BARTOLINI S，MENè L.Analysis of SSR signals by super resolution algorithms［C］//Proceedings of the Fourth IEEE International Symposium on Signal Processing and Information Technology，Roma，2004：166-170.

［4］CHAUMETTE E，COMON P，MULLER D.ICA-based technique for radiating sources estimation：application to airport surveillance［J］.IEE Proc.-F，1993，140（6）：395-401.

由于机会主义的存在，实际投资A10小于A1＊、A20小于A2＊，为了使利润最大化，中间组织应选择价格P0和W0，使买方和卖方在直接交易和通过中间组织交易之间至少没有感到差别：

［5］PETROCHILOS N，COMON P.A zero-cumulant random variable and its applications［J］.Signal Processing，2006，86（11）：3334-3338.

［6］PETROCHILOS N，GALATI G，PIRACCI E.Separation of SSR signals by array processing in multilateration systems［J］.IEEE Transactions on Aerospace and Electronic Systems，2009，45（3）：965-982.

［7］CAEDOSO J F，SOULOUMIAC A.Blind beamforming for non-Gaussian signals［C］//IEE Proceedings F-Radar and Signal Processing.IET，1993，140（6）：362-370.

［8］PETROCHILOS N，VANDER V A J.Algorithms to separate overlapping secondary surveillance radar replies［C］//IEEE International Conference on Acoustics Speech and Signal Processing（ICASSP），2004：II-49-52.

［9］GALATI G，LEONARDI M，PIRACCI E G，et al.The transponder data recorder：implementation and first results［J］.IEEE Aerospace and Electronic Systems magazine，2014，29（2）：405-412.

［10］COMON P.Tensor decompositions，state of the art and applications［J］.Mathematics in Signal Processing，Warwick，UK，2000（12）：18-20.

［11］YEREDOR A.Non-orthogonal joint diagonalization in the least-squares sense with application in blind source separation［J］.IEEE Transactions on Signal Processing，2002，50（7）：1545-1553.