CN110598289B - 一种不完备信息下的天线温度场测量方法 - Google Patents
一种不完备信息下的天线温度场测量方法 Download PDFInfo
- Publication number
- CN110598289B CN110598289B CN201910817061.XA CN201910817061A CN110598289B CN 110598289 B CN110598289 B CN 110598289B CN 201910817061 A CN201910817061 A CN 201910817061A CN 110598289 B CN110598289 B CN 110598289B
- Authority
- CN
- China
- Prior art keywords
- temperature field
- model
- antenna
- temperature
- calculating
- 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.)
- Active
Links
Images
Abstract
Description
技术领域
本发明属于雷达天线技术领域,涉及一种不完备信息下的天线温度场测量方法。
背景技术
天线在服役过程中,服役环境等不可抗因素会导致天线阵面温度产生不均匀变化;在天线内部,较小的空间内高密度集成了多个T/R组件,而单个T/R组件内多个收发通道上的多个芯片也会产生密集的热量,产生较大的热流密度,同样会导致温度的不均匀变化。温度的不均匀变化将会影响包括芯片S参数、收发通道驻波与增益乃至整个天线系统的电讯指标,因此热不均匀性给天线电性能带来的影响不可忽略。因此,为防止温度不均匀变化导致天线电性能恶化情况的发生,需要对天线阵面温度场进行监测。但在实际应用中,由于天线本身结构复杂,并且,天线阵面比较大,在每个位置布置传感器是不切实际的,所以只能依靠部分位置布置温度传感器获得的不完备信息进而计算出阵面的温度场,为后续天线电性能补偿做基础。
现有的基于不完备信息的温度场测量方法为基于本真正交分解的重构方法,基于POD分解应用于反射面天线的结构热类比方法,结构热类比方法是先确定反射面天线的结构方案,温度传感器的布局方案以及确定天线温度分布信息库,将天线划分环域,确定环域内的传感器的位置和数量,确定,提取当前工况下的实际数值,计算当前环域的初始温度场和初始温度值,计算环域类比系数,并修正初始温度场,重构环域的实际温度场,就目前技术而言,制约温度场重构精度的主要原因是对模型进行大量的实验与检测,过程中会产生较大的测量误差。例如申请公布号为CN108153954A,名称为“一种基于结构热类比的反射面天线温度快速重构方法”的专利申请,公开了一种基于结构热类比的反射面天线温度快速重构方法,通过划分环域,确定环域内温度传感器的位置和数量,实现了反射面天线的温度场快速重构。但其不足之处是需要对反射面天线的环境以及温度做大量的检测工作,对模型依赖度高,而且,测量过程中会带来较大的误差,无法满足天线温度场测量的精度。
发明内容
本发明的目的在于针对上述现有技术存在的不足,提供了一种不完备信息下的温度场测量方法,用于解决现有技术中存在的测量精度较低的技术问题。
为实现上述目的,本发明采取的技术方案包括如下步骤:
(1)建立待测天线模型及其对应的虚拟模型:
建立表面设置有p个温度传感器的待测天线模型,并建立与待测天线模型结构参数和材料属性相同的虚拟模型,在虚拟模型与待测天线模型中温度传感器对应的位置设置p个标记点,p≥2;
(2a)对虚拟模型施加热源和边界条件,并对虚拟模型进行网格划分,得到k个网格点,k>>p;
(3)计算y时刻待测天线模型的温度场T(y):
(3a)采用随机线性估计LSE方法,测量y-1时刻待测天线模型的温度,计算y-1时刻的POD模式系数αy-1,其中,y∈M;
(4)获取目标数据Ht和源数据Hs:
再对待测天线模型进行N次温度测量,得到p个位置处的温度,并将这些温度组成目标数据中标记点的温度Tnt,表示对待测天线模型的第j组测量数据,表示根据对待天线测模型第j组的测量数据计算出的待测天线模型的温度场,并将和组成目标数据表示虚拟模型中第i个时刻p个标记点处的温度值,Ti s表示虚拟模型第i个时刻的温度场,将Ti s和Ti ns组成源数据Hs=[Ti ns,Ti s],i=1,...,M;同时,组成待测天线模型在所有测量时刻的温度集合Tt,M>>N
(5)计算目标数据Ht的仿源数据Hs':
(6)获取训练数据H:
(7)建立温度计算伪预测模型pτ:
(7c)计算调整误差μτ:
计算训练误差系数eτ,并通过eτ计算调整误差μτ:
eτ=F/Dτ
(7d)判断μτ≥0.5或τ≥J是否成立,若成立,将步骤(7b)获取的pτ作为训练好的温度场计算伪预测模型pτ,否则,执行步骤(7e);
(7e)令τ=τ+1,同时对系数wτ进行更新,并执行步骤(7b),其中,wτ的更新公式为:
(8)计算天线的温度场:
(8a)将源温度信息Ti ns作为训练好的温度计算伪预测模型pτ(Tns)的输入,得到预测温度场Ts p,Ts p=pτ(Ti ns),并将预测温度场Ts p与目标温度场组合成伪目标温度场数据同时将源温度场Ti s与输出值组合成伪源温度场数据然后计算与之间的误差δ,
本发明与现有技术相比,具有如下优点:
1.本发明采用优化函数缩小了源数据Hs和目标数据Ht之间的差距,得到了逼近于目标温度值Tnt的伪源温度场值相当于把少量的目标数据Ht的性质迁移给了大量的源数据Hs,再用系数wτ进一步增强了目标数据Ht的影响力,从而,从虚拟模型中提取的源数据Hs更加接近目标数据Ht,最后得到了温度场重构方程来计算温度场,而现有技术依然通过基于插值的方法来重构温度场,并没有考虑构建的模型误差带来的影响,所以本发明相较于现有技术而言,考虑了模型误差和测量误差带来的影响,使重构精度更高;
2.本发明采用了递推最小二乘法计算POD模式系数,从开始测量温度信息的时刻起,每个时刻计算出该时刻对应的POD模式系数,将这些POD模式系数组成矩阵,采用递推最小二乘法根据POD模式系数矩阵计算下一时刻的POD模式系数,这样提前一个时刻计算温度场,能够提前知道下一时刻的温度信息,判断天线电性能的变化情况,对后续天线电性能补偿工作奠定基础。
附图说明
图1是本发明的实现流程图。
图2为本发明与验证现有技术温度场重构精度对比使用的模型和所处实验环境图。
具体实施方式
下面结合附图和具体实施例,对本发明作进一步的详细描述。
参照图1,本发明包括如下步骤:
步骤1)建立待测天线模型及其对应的虚拟模型:
建立表面设置有p个温度传感器的待测天线模型,并建立与待测天线模型结构参数和材料属性相同的虚拟模型,在虚拟模型与待测天线模型中温度传感器对应的位置设置p个标记点,p≥2;
步骤2a)对虚拟模型施加热源和边界条件,并对虚拟模型进行网格划分,得到k个网格点,k>>p;
步骤2b)对虚拟模型进行温度场瞬态仿真,得到M个时刻的瞬态温度场快照信息Ts和p个标记点处的M个时刻瞬态信息Tns,并对Ts进行本征正交基分解,得到待测天线模型的本征正交分解POD模式其中,M>>2;
(3)计算y时刻待测天线模型的温度场T(y):
步骤3a)采用随机线性估计LSE方法,测量y-1时刻待测天线模型的温度,计算y-1时刻的POD模式系数αy-1,其中,y∈M;
步骤4)获取目标数据Ht和源数据Hs:
再对待测天线模型进行N次温度测量,得到p个位置处的温度,并将这些温度组成目标数据中标记点的温度Tnt,表示对待测天线模型的第j组测量数据,表示根据对待天线测模型第j组的测量数据计算出的待测天线模型的温度场,并将和组成目标数据Ti ns表示虚拟模型中第i个时刻p个标记点处的温度值,Ti s表示虚拟模型第i个时刻的温度场,将Ti s和Ti ns组成源数据Hs=[Ti ns,Ti s],i=1,...,M;同时,组成待测天线模型在所有测量时刻的温度集合Tt,M>>N
步骤5)计算目标数据Ht的仿源数据Hs':
步骤6)获取训练数据H:
步骤7)建立温度计算伪预测模型pτ:
步骤7c)计算调整误差μτ:
计算训练误差系数eτ,并通过eτ计算调整误差μτ:
eτ=F/Dτ
步骤7d)判断μτ≥0.5或τ≥J是否成立,若成立,将步骤(7b)获取的pτ作为训练好的温度场计算伪预测模型pτ,否则,执行步骤(7e);
步骤7e)令τ=τ+1,同时对系数wτ进行更新,并执行步骤(7b),其中,wτ的更新公式为:
步骤8)计算天线的温度场:
步骤8a)将源温度信息Ti ns作为训练好的温度计算伪预测模型pτ(Tns)的输入,得到预测温度场Ts p,Ts p=pτ(Ti ns),并将预测温度场Ts p与目标温度场组合成伪目标温度场数据同时将源温度场Ti s与输出值组合成伪源温度场数据然后计算与之间的误差δ,
以下结合具体实验,对本发明的技术效果作进一步说明:
1、实验条件和内容:
在COMSOL Multiphysics 5.3进行热仿真分析,在MATLAB R2016a下运行算法程序,在图2所示的实验模型中进行实验。对本发明和现有的基于POD方法的温度场计算方法进行对比验证,其结果如表1所示。
表1
2、实验结果分析:
Claims (1)
1.一种不完备信息下的天线温度场测量方法,其特征在于,包括以下步骤:
(1)建立待测天线模型及其对应的虚拟模型:
建立表面设置有p个温度传感器的待测天线模型,并建立与待测天线模型结构参数和材料属性相同的虚拟模型,在虚拟模型与待测天线模型中温度传感器对应的位置设置p个标记点,p≥2;
(2a)对虚拟模型施加热源和边界条件,并对虚拟模型进行网格划分,得到k个网格点,k>>p;
(3)计算y时刻待测天线模型的温度场T(y):
(3a)采用随机线性估计LSE方法,测量y-1时刻待测天线模型的温度,计算y-1时刻的POD模式系数αy-1,其中,y∈M;
(4)获取目标数据Ht和源数据Hs:
再对待测天线模型进行N次温度测量,得到p个位置处的温度,并将这些温度组成目标数据中标记点的温度Tnt,表示对待测天线模型的第j组测量数据, 表示根据对待天线测模型第j组的测量数据计算出的待测天线模型的温度场,j=1,...,N,并将和组成目标数据j=1,...,N;Ti ns表示虚拟模型中第i个时刻p个标记点处的温度值,Ti s表示虚拟模型第i个时刻的温度场,将Ti s和Ti ns组成源数据Hs=[Ti ns,Ti s],i=1,...,M;同时,j=1,...,N,组成待测天线模型在所有测量时刻的温度集合Tt,M>>N;
(5)计算目标数据Ht的仿源数据Hs':
(6)获取训练数据H:
(7)建立温度计算伪预测模型pτ:
(7c)计算调整误差μτ:
计算训练误差系数eτ,并通过eτ计算调整误差μτ:
eτ=F/Dτ
(7d)判断μτ≥0.5或τ≥J是否成立,若成立,将步骤(7b)获取的pτ作为训练好的温度场计算伪预测模型pτ,否则,执行步骤(7e);
(7e)令τ=τ+1,同时对系数wτ进行更新,并执行步骤(7b),其中,wτ的更新公式为:
(8)计算天线的温度场:
(8a)将源温度信息Ti ns作为训练好的温度计算伪预测模型pτ(Tns)的输入,得到预测温度场Ts p,Ts p=pτ(Ti ns),并将预测温度场Ts p与目标温度场组合成伪目标温度场数据 同时将源温度场Ti s与输出值组合成伪源温度场数据 然后计算与之间的误差δ,
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910817061.XA CN110598289B (zh) | 2019-08-30 | 2019-08-30 | 一种不完备信息下的天线温度场测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910817061.XA CN110598289B (zh) | 2019-08-30 | 2019-08-30 | 一种不完备信息下的天线温度场测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110598289A CN110598289A (zh) | 2019-12-20 |
CN110598289B true CN110598289B (zh) | 2023-04-07 |
Family
ID=68857041
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910817061.XA Active CN110598289B (zh) | 2019-08-30 | 2019-08-30 | 一种不完备信息下的天线温度场测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110598289B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111461922B (zh) * | 2020-04-02 | 2023-04-21 | 国网冀北电力有限公司唐山供电公司 | 一种基于极限学习机的变压器热点温度实时预测方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108153954A (zh) * | 2017-12-19 | 2018-06-12 | 西安电子科技大学 | 一种基于结构热类比的反射面天线温度快速重构方法 |
CN109813967A (zh) * | 2017-11-21 | 2019-05-28 | 深圳市通用测试系统有限公司 | 一种阵列天线方向图的测量方法、设备、系统以及计算机可读存储介质 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101510229B (zh) * | 2009-03-20 | 2011-09-21 | 西安电子科技大学 | 基于机电热三场耦合的电子设备机箱结构优化设计方法 |
CN106650101A (zh) * | 2016-12-23 | 2017-05-10 | 西安电子科技大学 | 基于机电耦合模型的空间网状反射面天线温度载荷分析方法 |
-
2019
- 2019-08-30 CN CN201910817061.XA patent/CN110598289B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109813967A (zh) * | 2017-11-21 | 2019-05-28 | 深圳市通用测试系统有限公司 | 一种阵列天线方向图的测量方法、设备、系统以及计算机可读存储介质 |
CN108153954A (zh) * | 2017-12-19 | 2018-06-12 | 西安电子科技大学 | 一种基于结构热类比的反射面天线温度快速重构方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110598289A (zh) | 2019-12-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110940949B (zh) | 强冲击噪声环境下基于量子企鹅搜索机制的互质阵列doa估计方法 | |
CN101804581A (zh) | 一种机床热变形自动补偿的实现方法 | |
CN110749373A (zh) | 一种汽车噪声源检测方法 | |
CN107014339B (zh) | 一种用于大数据的高精度角度传感器误差补偿方法 | |
CN111929549B (zh) | 基于局部放电光学信号的gil局部放电源定位方法和系统 | |
CN112698112B (zh) | 电磁频谱地图构建方法、装置、计算机设备和存储介质 | |
CN110598289B (zh) | 一种不完备信息下的天线温度场测量方法 | |
CN104102836A (zh) | 一种电力系统快速抗差状态估计方法 | |
CN104614714B (zh) | 一种基于加权均方误差最小化的双重定标处理方法 | |
CN110470236B (zh) | 一种嵌入光纤光栅的柔性结构形变重构方法 | |
CN107422342A (zh) | Gnss卫星钟差实时估计质量控制方法 | |
CN107220450A (zh) | 一种非均匀材料连续分布力学参数场间接获取方法 | |
CN113835063B (zh) | 一种无人机阵列幅相误差与信号doa联合估计方法 | |
Su et al. | Parameter extraction of photovoltaic single-diode model using integrated current–voltage error criterion | |
CN107204616B (zh) | 基于自适应稀疏伪谱法的电力系统随机状态估计方法 | |
CN113503813A (zh) | 六自由度运动平台线位移定位精度测量与误差补偿方法 | |
CN116227389A (zh) | 气动热数据的预测方法及装置 | |
KR100316314B1 (ko) | 비정질화영역 결정방법 및 장치 | |
CN110175357B (zh) | 一种基于基准化分析的门级敏感性电路单元定位方法 | |
CN105699043B (zh) | 一种提高风洞传感器测量稳定性和精准度的方法 | |
CN110850366B (zh) | 混合视距和非视距环境下基于接收信号强度的定位方法 | |
CN115079090A (zh) | 基于esprit与加权降维搜索的多阵列非圆源快速定位方法 | |
CN114578197A (zh) | 基于平面uhf传感器的变压器套管局放三维定位方法及系统 | |
CN110657833B (zh) | 一种用于高精度源表一体化测量设备的新型校准方法 | |
CN113467590A (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 |