CN109709509A - 一种信号方位趋势的判断方法 - Google Patents
一种信号方位趋势的判断方法 Download PDFInfo
- Publication number
- CN109709509A CN109709509A CN201811585147.6A CN201811585147A CN109709509A CN 109709509 A CN109709509 A CN 109709509A CN 201811585147 A CN201811585147 A CN 201811585147A CN 109709509 A CN109709509 A CN 109709509A
- Authority
- CN
- China
- Prior art keywords
- data
- group
- follow
- trend
- observation
- 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
Links
Landscapes
- Mobile Radio Communication Systems (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种信号方位趋势的判断方法，接收一组有效序列数据，命名为A组，作为趋势判断的基准数据；接收下一个观察数据DATA Now，将观察数据DATA Now加入A组，形成一个观察数据队列，将观察数据队列命名为B组；继续接收下一个后续数据DATA Next，将后续数据DATA Next加入A组，形成一个后续数据队列，将后续数据队列命名为C组。本发明通过趋势判断处理，滤除偶然跳动的干扰数据，实际工程应用表明该方法简单可实现性强，滤波效果良好。
Description
技术领域
本发明涉及无线电测向领域，尤其涉及一种信号方位趋势的判断方法。
背景技术
随着无线电磁波在生产生活中的广泛运用，同时无线电信号辐射是有方向性特点，可通过无线电技术手段确定无线电信号辐射源方向，诞生了无线电测向技术。在实际应用中通过测向设备对对一种信号进行方位(范围0到360度)测量，即可判断信号的方位。
但是由于测向设备以及环境等因素的干扰，导致实际的测量时存在偶然跳动，造成测量结果可能会受到影响。
发明内容
为了解决上述问题，本发明提出一种信号方位趋势的判断方法。
一种信号方位趋势的判断方法，包括以下趋势判断步骤：
S1. 接收一组有效序列数据，命名为A组，作为趋势判断的基准数据；接收下一个观察数据DATA Now，将观察数据DATA Now加入A组，形成一个观察数据队列，将观察数据队列命名为B组；继续接收下一个后续数据DATA Next，将后续数据DATA Next加入A组，形成一个后续数据队列，将后续数据队列命名为C组；
S2. 计算A组、B组和C组的方差值，分别记为sd、sd1和sd2；并根据所述的sd、sd1和sd2计算观察数据在基准数据上引入的波动Fluct1和后续数据在观察数据基准上引入的波动Fluct2；
S3. 根据步骤S2所述的Fluct1和所述的Fluct2，计算Fluct1的判定门限值Threshold1和计算Fluct2的判定门限值Threshold2；
S4. 判断Fluct1与Threshold1之间的大小关系和判断Fluct2与Threshold2之间的大小关系；
S5. 若Fluct1> Threshold1且Fluct2> Threshold2，表示观察数据和后续数据的变化趋势均超出允许范围，则剔除本次观察数据，本次输出上一次输出值，同时保存后续数据作为下一次判断的观察值；
若Fluct1>Threshold1且Fluct2<=Threshold2，表示观察数据变化趋势在基准数据变化趋势范围以内，则输出本次观察数据，将本次观察数据加入A组；同时保存后续数据作为下一次判断的观察值；
若Fluct1<=Threshold1且Fluct2>Threshold2，表示观察数据和后续数据变化趋势一致，实际的数据变化趋势已经偏离基准数据的变化趋势，此变化趋势是允许的，则输出本次观察数据，将本次观察数据加入A组；同时保存后续数据作为下一次判断的观察值；
若Fluct1<Threshold1且Fluct2<Threshold2，表示观察数据和后续数据的变化趋势均未超出允许范围，则输出本次观察数据，将本次观察数据加入A组；同时保存后续数据作为下一次判断的观察值。
信号方位趋势的判断方法还包括以下初始化步骤：
S21. 设置基准队列中的起始数据个数和基准队列中的数据个数上限值，同时设置门限判定百分比per；
S22. 将测向设备接收的方位数据进行保存，形成一个新的队列，并判断新的队列中的数据个数与起始数据个数之间的关系；
S23. 若新的队列中的数据个数小于起始数据个数，则返回至步骤S22，同时输出本次接收方位数据；
若新的队列中的数据个数等于起始数据个数，则保存本次接受数据作为观察值，然后返回至步骤S22，接收后续数据进行判断；
若新的队列中的数据个数大于起始数据个数，则判断新的队列中的数据个数是否超过数据个数上限值；
S33. 若否，则执行所述的趋势判断步骤，然后返回至步骤S22；
若超过数据个数上限值，则将新的队列中的首个数据移除，确保新的队列的数据个数不超过数据个数上限值，然后执行所述的趋势判断步骤，然后返回至步骤S22。
所述步骤S2中的Fluct1计算公式为Fluct1=abs（sd1-sd），其中abs（）为取绝对值函数。
所述步骤S2中的Fluct2计算公式为Fluct2=abs（sd2-sd1），其中abs（）为取绝对值函数。
所述步骤S3中的Threshold1计算公式为Threshold1=min（sd1，sd）*per，其中min（）表示取sd1与sd中的最小值，per为门限判定百分比。
所述步骤S3中的Threshold2计算公式为Threshold2=min（sd2，sd1）*per，其中min（）表示取sd2与sd1中的最小值，per为门限判定百分比。
本发明的有益效果在于：本发明能够很好的滤除偶然出现的干扰数据，滤波以后输出数据符合实际变化趋势，实际工程应用效果良好。
附图说明
图1是趋势判断方法流程图；
图2是初始化流程图
图3是趋势判断滤波前后输出数据曲线对比图。
具体实施方式
为了使本领域的技术人员更好地理解本发明的技术方案，下面结合附图和具体实施例对本发明作进一步的详细说明。
实施方案一：一种信号方位趋势的判断方法，包括以下趋势判断步骤：
S1. 接收一组有效序列数据，命名为A组，作为趋势判断的基准数据；接收下一个观察数据DATA Now，将观察数据DATA Now加入A组，形成一个观察数据队列，将观察数据队列命名为B组；继续接收下一个后续数据DATA Next，将后续数据DATA Next加入A组，形成一个后续数据队列，将后续数据队列命名为C组；
S2. 计算A组、B组和C组的方差值，分别记为sd、sd1和sd2；并根据所述的sd、sd1和sd2计算观察数据在基准数据上引入的波动Fluct1和后续数据在观察数据基准上引入的波动Fluct2；
S3. 根据步骤S2所述的Fluct1和所述的Fluct2，计算Fluct1的判定门限值Threshold1和计算Fluct2的判定门限值Threshold2；
S4. 判断Fluct1与Threshold1之间的大小关系和判断Fluct2与Threshold2之间的大小关系；
S5. 若Fluct1> Threshold1且Fluct2> Threshold2，表示观察数据和后续数据的变化趋势均超出允许范围，则剔除本次观察数据，本次输出上一次输出值，同时保存后续数据作为下一次判断的观察值；
若Fluct1>Threshold1且Fluct2<=Threshold2，表示观察数据变化趋势在基准数据变化趋势范围以内，则输出本次观察数据，将本次观察数据加入A组；同时保存后续数据作为下一次判断的观察值；
若Fluct1<=Threshold1且Fluct2>Threshold2，表示观察数据和后续数据变化趋势一致，实际的数据变化趋势已经偏离基准数据的变化趋势，此变化趋势是允许的，则输出本次观察数据，将本次观察数据加入A组；同时保存后续数据作为下一次判断的观察值；
若Fluct1<Threshold1且Fluct2<Threshold2，表示观察数据和后续数据的变化趋势均未超出允许范围，则输出本次观察数据，将本次观察数据加入A组；同时保存后续数据作为下一次判断的观察值。
所述A组方差sd表示基准数据的波动程度，所述B组方差sd1表示DATA Now引入的波动程度，所述C组方差sd2表示DATA Next引入的波动程度。
所述步骤S2中的Fluct1计算公式为Fluct1=abs（sd1-sd），其中abs（）为取绝对值函数。
所述步骤S2中的Fluct2计算公式为Fluct2=abs（sd2-sd1），其中abs（）为取绝对值函数。
所述步骤S2中的Threshold1计算公式为Threshold1=min（sd1，sd）*per，其中min（）表示取sd1与sd中的最小值，per为门限判定百分比。
所述步骤S2中的Threshold2计算公式为Threshold2=min（sd2，sd1）*per，其中min（）表示取sd2与sd1中的最小值，per为门限判定百分比。
如图2所示，信号方位趋势的判断方法还包括以下初始化步骤：
S21. 设置基准队列中的起始数据个数和基准队列中的数据个数上限值，同时设置门限判定百分比per；
S22. 将测向设备接收的方位数据进行保存，形成一个新的队列，并判断新的队列中的数据个数与起始数据个数之间的关系；
S23. 若新的队列中的数据个数小于起始数据个数，则返回至步骤S22，同时输出本次接收方位数据；
若新的队列中的数据个数等于起始数据个数，则保存本次数据作为观察值，然后返回至步骤S22，接收后续数据进行判断；
若新的队列中的数据个数大于起始数据个数，则判断新的队列中的数据个数是否超过数据个数上限值；
S33. 若否，则执行所述的趋势判断步骤，然后返回至步骤S22；
若超过数据个数上限值，则将新的队列中的首个数据移除，确保新的队列的数据个数不超过数据个数上限值，然后执行所述的趋势判断步骤，然后返回至步骤S22。
实施方案二：如队列起始数据个数设置为3，那么数据队列数据这时是123，继续接收数据4进行保存（观察值），输出数据3（输出上一次数据），再接收数据5来判断决定4这个数据是不是偶然跳动的值，如果是就剔除，若数据4在允许范围内，保存数据4到数据队列，此时数据队列为1234。若基准队列数据个数上限值为10，当第11个数据11在允许范围内，保存数据11到数据队列，同时将首元素数据1移除，此时数据队列为234…11。
方差是衡量一组数据的波动大小，方差越小，说明数据越稳定。在这里判断数据的趋势就是计算3组数据（123、1234、1235）的方差分别是sd、sd1、sd2。通过对sd和sd1进行比较，对sd1和sd2进行比较，如果都在允许范围（设置值）之内，直接输出观察值同时保存到数据队列作为采样数据，作为后续数据判断的基准。否则剔除观察值，将当前值作为观察值，等接收下一次数据进行判断，依次循环。
如图3所示，右图表示采用趋势判断滤波输出数据曲线，左图表示原始数据曲线，10_30%表示趋势判断设定参数，其中10表示基准队列数据个数上限值，30%表示门限判定时最小值设定的百分比。从对比曲线可以看出本发明能够很好的滤除偶然出现的干扰数据，滤波以后输出数据符合实际变化趋势。
本发明的有益效果在于：本发明能够很好的滤除偶然出现的干扰数据，滤波以后输出数据符合实际变化趋势，实际工程应用效果良好。
需要说明的是，对于前述的各个方法实施例，为了简单描述，故将其都表述为一系列的动作组合，但是本领域技术人员应该知悉，本申请并不受所描述的动作顺序的限制，因为依据本申请，某一些步骤可以采用其他顺序或者同时进行。其次，本领域技术人员也应该知悉，说明书中所描述的实施例均属于优选实施例，所涉及的动作和单元并不一定是本申请所必须的。
在上述实施例中，对各个实施例的描述都各有侧重，某个实施例中没有详细描述的部分，可以参见其他实施例的相关描述。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成，所述的程序可存储于计算机可读取存储介质中，该程序在执行时，可包括如上述各方法的实施例的流程。其中，所述的存储介质可为磁碟、光盘、ROM、RAM等。
以上所揭露的仅为本发明较佳实施例而已，当然不能以此来限定本发明之权利范围，因此依本发明权利要求所作的等同变化，仍属本发明所涵盖的范围。
Claims (6)
1.一种信号方位趋势的判断方法，其特征在于：包括以下趋势判断步骤：
S1. 接收一组有效序列数据，命名为A组，作为趋势判断的基准数据；接收下一个观察数据DATA Now，将观察数据DATA Now加入A组，形成一个观察数据队列，将观察数据队列命名为B组；继续接收下一个后续数据DATA Next，将后续数据DATA Next加入A组，形成一个后续数据队列，将后续数据队列命名为C组；
S2. 计算A组、B组和C组的方差值，分别记为sd、sd1和sd2；并根据所述的sd、sd1和sd2计算观察数据在基准数据上引入的波动Fluct1和后续数据在观察数据基准上引入的波动Fluct2；
S3. 根据步骤S2所述的Fluct1和所述的Fluct2，计算Fluct1的判定门限值Threshold1和计算Fluct2的判定门限值Threshold2；
S4. 判断Fluct1与Threshold1之间的大小关系和判断Fluct2与Threshold2之间的大小关系；
S5. 若Fluct1> Threshold1且Fluct2> Threshold2，表示观察数据和后续数据的变化趋势均超出允许范围，则剔除本次观察数据，本次输出上一次输出值，同时保存后续数据作为下一次判断的观察值；
若Fluct1>Threshold1且Fluct2<=Threshold2，表示观察数据变化趋势在基准数据变化趋势范围以内，则输出本次观察数据，将本次观察数据加入A组；同时保存后续数据作为下一次判断的观察值；
若Fluct1<=Threshold1且Fluct2>Threshold2，表示观察数据和后续数据变化趋势一致，实际的数据变化趋势已经偏离基准数据的变化趋势，此变化趋势是允许的，则输出本次观察数据，将本次观察数据加入A组；同时保存后续数据作为下一次判断的观察值；
若Fluct1<Threshold1且Fluct2<Threshold2，表示观察数据和后续数据的变化趋势均未超出允许范围，则输出本次观察数据，将本次观察数据加入A组；同时保存后续数据作为下一次判断的观察值。
2.如权利要求1所述的一种信号方位趋势的判断方法，其特征在于：还包括以下初始化步骤：
S21. 设置基准队列中的起始数据个数和基准队列中的数据个数上限值，同时设置门限判定百分比per；
S22. 将测向设备接收的方位数据进行保存，形成一个新的队列，并判断新的队列中的数据个数与起始数据个数之间的关系；
S23. 若新的队列中的数据个数小于起始数据个数，则返回至步骤S22，同时输出本次接收方位数据；
若新的队列中的数据个数等于起始数据个数，则保存本次接收数据作为观察值，同时返回至步骤S22，接收后续数据进行判断；
若新的队列中的数据个数大于起始数据个数，则判断新的队列中的数据个数是否超过数据个数上限值；
S33. 若否，则执行所述的趋势判断步骤，返回至步骤S22；
若超过数据个数上限值，则将新的队列中的首个数据移除，确保新的队列的数据个数不超过数据个数上限值，然后执行所述的趋势判断步骤，返回至步骤S22。
3.如权利要求1所述的一种信号方位趋势的判断方法，其特征在于：所述步骤S2中的Fluct1计算公式为Fluct1=abs（sd1-sd），其中abs（）为取绝对值函数。
4.如权利要求1所述的一种信号方位趋势的判断方法，其特征在于：所述步骤S2中的Fluct2计算公式为Fluct2=abs（sd2-sd1），其中abs（）为取绝对值函数。
5.如权利要求1所述的一种信号方位趋势的判断方法，其特征在于：所述步骤S2中的Threshold1计算公式为Threshold1=min（sd1，sd）*per，其中min（）表示取sd1与sd中的最小值，per为门限判定百分比。
6.如权利要求1所述的一种信号方位趋势的判断方法，其特征在于：所述步骤S2中的Threshold2计算公式为Threshold2=min（sd2，sd1）*per，其中min（）表示取sd2与sd1中的最小值，per为门限判定百分比。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811585147.6A CN109709509B (zh) | 2018-12-24 | 2018-12-24 | 一种信号方位趋势的判断方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811585147.6A CN109709509B (zh) | 2018-12-24 | 2018-12-24 | 一种信号方位趋势的判断方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109709509A true CN109709509A (zh) | 2019-05-03 |
| CN109709509B CN109709509B (zh) | 2022-01-11 |
Family
ID=66257285
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811585147.6A Active CN109709509B (zh) | 2018-12-24 | 2018-12-24 | 一种信号方位趋势的判断方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109709509B (zh) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20060009909A1 (en) * | 2004-07-12 | 2006-01-12 | Aviation Communication & Surveillance Systems Llc | Systems and methods for determining bearing |
| JP2014010483A (ja) * | 2012-06-27 | 2014-01-20 | Yanmar Co Ltd | 予測装置、予測方法及びコンピュータプログラム |
| CN103746750A (zh) * | 2013-08-23 | 2014-04-23 | 西华大学 | 无线电监测电磁态势预测系统 |
| CN104220896A (zh) * | 2012-04-13 | 2014-12-17 | 高通股份有限公司 | 用于估计到达方向的系统、方法和设备 |
| CN108344895A (zh) * | 2018-02-22 | 2018-07-31 | 国网辽宁省电力有限公司电力科学研究院 | 信号检测方法及装置 |
-
2018
- 2018-12-24 CN CN201811585147.6A patent/CN109709509B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20060009909A1 (en) * | 2004-07-12 | 2006-01-12 | Aviation Communication & Surveillance Systems Llc | Systems and methods for determining bearing |
| CN104220896A (zh) * | 2012-04-13 | 2014-12-17 | 高通股份有限公司 | 用于估计到达方向的系统、方法和设备 |
| JP2014010483A (ja) * | 2012-06-27 | 2014-01-20 | Yanmar Co Ltd | 予測装置、予測方法及びコンピュータプログラム |
| CN103746750A (zh) * | 2013-08-23 | 2014-04-23 | 西华大学 | 无线电监测电磁态势预测系统 |
| CN108344895A (zh) * | 2018-02-22 | 2018-07-31 | 国网辽宁省电力有限公司电力科学研究院 | 信号检测方法及装置 |
Non-Patent Citations (3)
| Title |
|---|
| DON TORRIERI: "Direction Finding of a Compromised Node in a Wireless Network", 《IEEE TRANSACTIONS ON MOBILE COMPUTING》 * |
| 杨军佳 等: "相干干扰下无源测向系统方位测量误差模型", 《探测与控制学报》 * |
| 章辉 等: "一种无线信号变化趋势的预测方法", 《2006通信理论与技术新进展——第十一届全国青年通信学术会议论文集》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109709509B (zh) | 2022-01-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US10742027B2 (en) | Method and device for online limit early-warning to shunt capacitor bank | |
| CN109740687A (zh) | 一种基于dlae的发酵过程故障监测方法 | |
| CN105424190A (zh) | 产品外观的灰度检测方法 | |
| CN103942444B (zh) | 一种基于dtw的平权1‑范数多传感器不等长序列相似度度量方法 | |
| CN107784314A (zh) | 一种多变量报警系统的正常异常数据划分方法及系统 | |
| CN112043258A (zh) | 动态心率预测方法、装置、设备及可读存储介质 | |
| CN109709509A (zh) | 一种信号方位趋势的判断方法 | |
| CN104199417A (zh) | 一种半导体镀膜工艺的统计过程监控方法 | |
| CN103995180A (zh) | 一种考虑不等式约束的电力系统频率估计方法 | |
| CN102034022B (zh) | 一种基于倍频分析的信号处理方法及系统 | |
| CN111524036B (zh) | 一种暂态电能质量扰动分类方法 | |
| CN113255541B (zh) | 一种基于本征模态函数重组信号相对熵的自适应流程工业过程的工艺参数去噪方法 | |
| Paul et al. | How price signals in pulses are transmitted across regions and value chain? Examining horizontal and vertical market price integration for major pulses in India | |
| Gupta et al. | Extreme Learning Machines with frequency based noise filtering for prediction of critical digressions in a noisy industrial process | |
| CN107256045A (zh) | 大体积混凝土智能温控的方法及装置 | |
| CN104535033B (zh) | 准确测量不同流程产品阻焊厚度的方法和系统 | |
| CN106595628B (zh) | 数字光纤陀螺信号的野值实时剔除方法 | |
| Mocci et al. | SOM analysis of IMPOC data in relation to process parameters in strip steel production | |
| Fritsche et al. | Bayesian Bhattacharyya Bound for discrete-time filtering revisited | |
| Ramezani et al. | Closed Loop Identification of Multi-Rate System by Expectation-Maximization | |
| TWI736999B (zh) | 資料分析裝置、資料分析方法以及與其相關的品質提升系統 | |
| Aguero et al. | Virtual closed loop identification: A Subspace Approach | |
| CN116247734B (zh) | 面向边缘侧弱通信环境的分布式一致性功率控制方法 | |
| CN108594676A (zh) | 一类工业受控过程特征参数快速辨识方法 | |
| CN102181322B (zh) | 一种药用茶油制备的在线监测方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |