CN113282553B - 一种电磁环境远场数据的快速无损压缩存储方法 - Google Patents

一种电磁环境远场数据的快速无损压缩存储方法 Download PDF

Info

Publication number
CN113282553B
CN113282553B CN202110566228.7A CN202110566228A CN113282553B CN 113282553 B CN113282553 B CN 113282553B CN 202110566228 A CN202110566228 A CN 202110566228A CN 113282553 B CN113282553 B CN 113282553B
Authority
CN
China
Prior art keywords
data
far
field
nbase
field data
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
Application number
CN202110566228.7A
Other languages
English (en)
Other versions
CN113282553A (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 CN202110566228.7A priority Critical patent/CN113282553B/zh
Publication of CN113282553A publication Critical patent/CN113282553A/zh
Application granted granted Critical
Publication of CN113282553B publication Critical patent/CN113282553B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
    • G06F16/10—File systems; File servers
    • G06F16/17—Details of further file system functions
    • G06F16/174—Redundancy elimination performed by the file system
    • G06F16/1744—Redundancy elimination performed by the file system using compression, e.g. sparse files
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Evolutionary Computation (AREA)
  • Computer Hardware Design (AREA)
  • Data Mining & Analysis (AREA)
  • Databases & Information Systems (AREA)
  • Computer Graphics (AREA)
  • Software Systems (AREA)
  • Mram Or Spin Memory Techniques (AREA)

Abstract

本发明公开了一种电磁环境远场数据的快速无损压缩存储方法，包括以下步骤：S1.加载CAD网格模型；S2.转换形成系统矩阵和右侧激励项；S3.构建电磁环境的远场模型，通过步骤S2中得到的数据得到每个方向上垂直极化方向
Figure DDA0003081071360000011
水平极化方向
Figure DDA0003081071360000012
的远场；S4.将步骤S2得到面元导出成文件；并将每个网格顶点的远场电场和磁场数据导出成文件；S5.将导出远场数据和导出的网格数据读取进内；S6.对内存中的远场数据和网格数据进行压缩；S7.借助于rar压缩工具将得到的远场数据和网格数据进行二次压缩，得到最终的压缩存储文件。本发明能够有效降低远场数据的存储空间，且数据格式能够适用于海量数据的存储。

Description

一种电磁环境远场数据的快速无损压缩存储方法
技术领域
本发明涉及电磁环境数据存储，特别是涉及一种电磁环境远场数据的快速无损压缩存储方法。
背景技术
电磁环境数据库在进行设备测试、场景模拟等情况下具有着非常重要的作用，电磁环境数据库中的数据主要包含天线的近场数据和远场数据；
对于远场数据而言，一般通过传统矩量法得到的天线远场数据，但是这些数据占用硬盘空间很大，具有很大的冗余，且数据格式不适用于海量数据存储，也不适用于对存储长度要求非常严格的场合。
发明内容
本发明的目的在于克服现有技术的不足，提供一种电磁环境远场数据的快速无损压缩存储方法，能够有效降低远场数据的存储空间，且数据格式能够适用于海量数据的存储。
本发明的目的是通过以下技术方案来实现的：一种电磁环境远场数据的快速无损压缩存储方法，包括以下步骤：
S1.加载CAD网格模型：
从nastran格式的目标对象网格模型文件中提取CAD面元网格模型的点列表和点连接列表，设置材料电磁参数信息和模型激励源参数，完成电磁模型创建工作；
S2.将提取的CAD面元网格模型、设置的材料电磁参数信息，以及设置的激励源参数转换成系统矩阵[A]nbase×nbase和右侧激励项[rhs]nbase；
S3.构建电磁环境的远场模型，通过步骤S2中获得的数据得到每个方向上垂直极化方向
Figure BDA0003081071340000011
水平极化方向
Figure BDA0003081071340000012
的远场；
S4.将步骤S1得到CAD面元网格模型的点列表和点连接列表导出网格数据文件；步骤S3得到的每个方向上垂直极化方向
Figure BDA0003081071340000013
水平极化方向
Figure BDA0003081071340000014
的远场导出成远场数据文件；
S5.将导出的网格数据文件和远场数据文件读取进内存；
S6.对内存中的远场数据和网格数据进行压缩；
S7.借助于rar压缩工具将得到的远场数据和网格数据进行二次压缩，得到最终的压缩存储文件。
进一步地，所述步骤S2包括：
S201.将CAD面元网格模型中的点列表和点连接列表通过RWG基函数的格式转换成矩量法所需的基函数，设nbase为电磁模型中所有基函数未知量总数，材料的电磁参数信息包括磁导率μl和介电常数εl；激励源参数包括源点所在基函数区域的位置矢量和入射场分布；
S202.转换得到系统矩阵[A]nbase×nbase：设amn为[A]nbase×nbase的第m行第n列元素，则amn填充公式如下：
Figure BDA0003081071340000021
其中，G为三维格林函数；
Figure BDA0003081071340000022
ω＝2πf0为角频率；f0为频率；
Figure BDA0003081071340000023
为哈密尔顿算子；fm和fn为第m个和第n个基函数；r′为源点所在区域即fn区域内的位置矢量；r为场点所在区域即fm区域内的位置矢量；
S203.转换得到右侧激励项[rhs]nbase：设rhsm为右侧激励项[rhs]nbase中的第m个元素，则填充公式如下：
Figure BDA0003081071340000024
式中，
Figure BDA0003081071340000025
为第m个基函数域内的入射场分布；
S204.通过矩阵求解计算求解[A]nbase×nbasex＝[rhs]nbase得到x，作为每个基函数ibase未知量上的电流量x[ibase]，ibase为1到nbase的整数下标。
进一步地，所述步骤S3中，所述步骤S3中，通过步骤S2中得到的模型中每个基函数上的电流的积分，得到每个方向上垂直极化方向
Figure BDA0003081071340000026
水平极化方向
Figure BDA0003081071340000027
的远场，针对计算域S′的面积分的远场计算公式为：
Figure BDA0003081071340000028
Figure BDA0003081071340000029
其中，
Figure BDA00030810713400000210
为电磁模型远场的垂直极化的单位方向矢量，
Figure BDA00030810713400000211
为电磁模型远场的水平极化的单位方向矢量，
Figure BDA00030810713400000212
电磁模型远场的径向单位矢量；S′为电磁模型基函数所组成的计算域；r′为S′上任何一点矢量；R为远场的径向距离；k为自由空间波数；n为面元法向单位矢量；k为自由空间波数；
Figure BDA00030810713400000213
进一步地，所述步骤S4中，远场数据文件采用维度为
Figure BDA00030810713400000214
的三维浮点数组
Figure BDA00030810713400000215
其中，mθ为θ角总数，
Figure BDA00030810713400000216
为
Figure BDA00030810713400000217
角总数；设θ在0～180°内取值，取值间隔为
Figure BDA00030810713400000218
在0～360°内取值，取值间隔为
Figure BDA00030810713400000219
则：
θ角的取值θi为：
Figure BDA0003081071340000031
Figure BDA0003081071340000032
角的取值
Figure BDA0003081071340000033
为：
Figure BDA0003081071340000034
每个
Figure BDA0003081071340000035
位置的远场存储格式
Figure BDA0003081071340000036
均包含4列数据：
第一列：Mag(E⊥)；第二列Deg(E⊥)；第三列Mag(E||)；第四列Deg(E||)；
其中，E⊥为垂直极化电场，E||为水平极化电场，Mag(x)为幅度函数，Deg(x)为相位函数；
由于θi有mθ个取值，
Figure BDA0003081071340000037
有
Figure BDA0003081071340000038
个取值，故实际上总共有
Figure BDA0003081071340000039
个
Figure BDA00030810713400000310
又由于
Figure BDA00030810713400000311
均包含4列数据，故形成了维度为
Figure BDA00030810713400000312
的三维浮点数组
Figure BDA00030810713400000313
进一步地，所述步骤S6包括：
S601.对内存中加载的远场数据格式作如下归一化处理：
按照用户需求设定如下参数：空间分辨率:Sr、动态范围:Dr、幅度分辨率:Mr、相位分辨率:Pr；
其中，20≤Sr≤300；20≤Dr≤300；Dr≤Sr；0.01≤Mr≤1；0.1≤Pe≤5
在计算前对Mr和Pr的值进行处理：
Nmr＝ceil(log2Sr/Me)
Npr＝ceil(log2360/Pr)
Figure BDA00030810713400000314
Figure BDA00030810713400000315
其中，Nmr为幅度存储位数，Npr为相位存储位数；
S602.将远场数据离散压缩处理:
第一步、从文件中读取维度为
Figure BDA00030810713400000316
的三维浮点数组
Figure BDA00030810713400000317
该数组由共有
Figure BDA00030810713400000318
个
Figure BDA00030810713400000319
所对应的远场数据
Figure BDA00030810713400000320
构成；每个
Figure BDA00030810713400000321
包括四列数据，分别为：
Mag(E⊥)、Deg(E⊥)、Mag(E||)、Deg(E||)；
并从中恢复出E⊥、E||；
对于每一个
Figure BDA00030810713400000322
所对应的远场数据
Figure BDA00030810713400000323
进行总场计算，总场计算公式：
Etot＝10*log10(|E⊥|2+|E|||2+1e-30)
其中，1e-30的作用是为了保证log10取值的有效性；
第二步、取Etot最大值，记为Emax；
第三步、对任一个远场数据
Figure BDA0003081071340000041
进行处理；
A01.判断是否满足：Etot＜Emax-Dr,若满足，则将E⊥、E||置零，若不满足，则E⊥、E||不变；
A02.对步骤A01得到的数据再次进行处理：
判断是否满足20*log10(|E⊥|+1e-30)<Etot-Sr,若满足，E⊥置0，若不满足，E⊥不变；
判断是否满足20*log10(|E|||+1e-30)<Etot-Sr,若满足，E||置0，若不满足，E||不变；
A03、对于步骤A02得到的数据采用的场值离散化流程进行处理，公式为：
Figure BDA0003081071340000042
Figure BDA0003081071340000043
Figure BDA0003081071340000044
Figure BDA0003081071340000045
进行上述处理后，将每个远场数据
Figure BDA0003081071340000046
作为一行用4*(Nmr+Npr)个二进制数位表示；
第四步、将远场数据
Figure BDA0003081071340000047
的处理结果写入二进制文件：
二进制文件中，数据头的数据格式
数位(byte) 4 4 4 4 4
数值 S<sub>r</sub> D<sub>r</sub> M<sub>e</sub> P<sub>e</sub> E<sub>max</sub>
类型 浮点 浮点 浮点 浮点 浮点
数据行的每行数据格式：
Figure BDA0003081071340000048
第五步、对于不同
Figure BDA0003081071340000049
对应的远场数据
Figure BDA00030810713400000410
重复执行第三步到第四步，得到所有远场数据的压缩文件。
本发明的有益效果是：本发明能够有效降低远场数据的存储空间，且数据格式能够适用于海量数据的存储和一些对存储长度要求非常严格的场合。
附图说明
图1为本发明的方法流程图；
图2为电磁环境的远场模型示意图。
具体实施方式
下面结合附图进一步详细描述本发明的技术方案，但本发明的保护范围不局限于以下所述。
如图1所示，一种电磁环境远场数据的快速无损压缩存储方法，包括以下步骤：
S1.加载CAD网格模型：
从nastran格式的目标对象网格模型文件中提取CAD面元网格模型的点列表和点连接列表，设置材料电磁参数信息和模型激励源参数，完成电磁模型创建工作；在本申请的实施例中，所述目标对象包括天线、电路，以及天线与电路形成的复杂组合体中的一种。
S2.将提取的CAD面元网格模型、设置的材料电磁参数信息，以及设置的激励源参数转换成系统矩阵[A]nbase×nbase和右侧激励项[rhs]nbase；
所述步骤S2包括：
S201.将CAD面元网格模型中的点列表和点连接列表通过RWG基函数的格式转换成矩量法所需的基函数，设nbase为电磁模型中所有基函数未知量总数，材料的电磁参数信息包括磁导率μl和介电常数εl；激励源参数包括源点所在基函数区域的位置矢量和入射场分布；
S202.转换得到系统矩阵[A]nbase×nbase：设amn为[A]nbase×nbase的第m行第n列元素，则amn填充公式如下：
Figure BDA0003081071340000051
其中，G为三维格林函数；
Figure BDA0003081071340000052
ω＝2πf0为角频率；f0为频率；
Figure BDA0003081071340000053
为哈密尔顿算子；fm和fn为第m个和第n个基函数；r′为源点所在区域即fn区域内的位置矢量；r为场点所在区域即fm区域内的位置矢量；
S203.转换得到右侧激励项[rhs]nbase：设rhsm为右侧激励项[rhs]nbase中的第m个元素，则填充公式如下：
Figure BDA0003081071340000054
式中，
Figure BDA0003081071340000061
为第m个基函数域内的入射场分布；
S204.通过矩阵求解计算求解[A]nbase×nbasex＝[rhs]nbase得到x，作为每个基函数ibase未知量上的电流量x[ibase]，ibase为1到nbase的整数下标。
S3.构建电磁环境的远场模型，通过步骤S2中获得的数据得到每个方向上垂直极化方向
Figure BDA0003081071340000062
水平极化方向
Figure BDA0003081071340000063
的远场；其中，电磁环境的远场模型如图2所示；
通过步骤S2中得到的模型中每个基函数上的电流的积分，得到每个方向上垂直极化方向
Figure BDA0003081071340000064
水平极化方向
Figure BDA0003081071340000065
的远场，针对计算域S′的面积分的远场计算公式为：
Figure BDA0003081071340000066
Figure BDA0003081071340000067
其中，
Figure BDA0003081071340000068
为电磁模型远场的垂直极化的单位方向矢量，
Figure BDA0003081071340000069
为电磁模型远场的水平极化的单位方向矢量，
Figure BDA00030810713400000610
电磁模型远场的径向单位矢量；S′为电磁模型基函数所组成的计算域；r′为S′上任何一点矢量；R为远场的径向距离；k为自由空间波数；n为面元法向单位矢量；k为自由空间波数；
Figure BDA00030810713400000611
S4.将步骤S1得到CAD面元网格模型的点列表和点连接列表导出网格数据文件；步骤S3得到的每个方向上垂直极化方向
Figure BDA00030810713400000612
水平极化方向
Figure BDA00030810713400000613
的远场导出成远场数据文件；
以常见的三角形面单元网格为例，一般的存储格式包括如下四个部分：
1)点总数nPoint
2)点坐标列表Points
3)单元总数nElement顶点数nVertex
4)点连接列表Connects
其中，nPoint、nElement和nVertex为整型数据，点坐标列表存储的是维度为nPoint*3的二维浮点数组，Connects存储的维度为nElement*(nVertex+1)的二维整型数组，前nVertex列存储的是每个单元的各个顶点所在Points中的序号(序号以1开始计数)，最后一列存储的是每个单元的材料编号(一般0表示完全电导体材料)；
在本申请的实施例中，远场数据文件采用维度为
Figure BDA00030810713400000614
的三维浮点数组
Figure BDA00030810713400000615
其中，mθ为θ角总数，
Figure BDA00030810713400000616
为
Figure BDA00030810713400000617
角总数；设θ在0～180°内取值，取值间隔为
Figure BDA00030810713400000618
在0～360°内取值，取值间隔为
Figure BDA00030810713400000619
则：
θ角的取值θi为：
Figure BDA0003081071340000071
Figure BDA0003081071340000072
角的取值
Figure BDA0003081071340000073
为：
Figure BDA0003081071340000074
每个
Figure BDA0003081071340000075
位置的远场存储格式
Figure BDA0003081071340000076
均包含4列数据：
数据列 第一列 第二列 第三列 第四列
数据 Mag(E<sup>⊥</sup> Deg(E<sup>⊥</sup>) Mag(E<sup>||</sup>) Deg(E<sup>||</sup>)
其中，E⊥为垂直极化电场，E||为水平极化电场，，Mag(x):幅度函数(dB)，Deg(x):相位函数(度)；
由于θi有mθ个取值，
Figure BDA0003081071340000077
有
Figure BDA0003081071340000078
个取值，故实际上总共有
Figure BDA0003081071340000079
个
Figure BDA00030810713400000710
又由于
Figure BDA00030810713400000711
均包含4列数据，故形成了维度为
Figure BDA00030810713400000712
的三维浮点数组
Figure BDA00030810713400000713
S5.将导出的网格数据文件和远场数据文件读取进内存，分别建立文件与内存数据的映射，该映射即每个
Figure BDA00030810713400000714
与对应远场
Figure BDA00030810713400000715
的映射；
S6.对内存中的远场数据和网格数据进行压缩；
S7.借助于rar压缩工具将得到的远场数据和网格数据进行二次压缩，得到最终的压缩存储文件。
其中，步骤S6中远场数据的处理过程如下：
S601.对内存中加载的远场数据格式做如下归一化处理：
按照用户需求设定如下参数：空间分辨率(dB):Sr、动态范围(dB):Dr、幅度分辨率(dB):Mr、相位分辨率(度):Pr；
参数约定：(1)20≤Sr≤300；(2)20≤Dr≤300；(3)Dr≤Sr；(4)0.01≤Mr≤1；(5)0.1≤Pe≤5；
Sr默认取120dB；Dr默认取60dB；Me默认取0.3；Pe默认取0.36；
由下列公式对Mr和Pr的值在计算前进行处理：
Nmr＝ceil(log2Sr/Me)
Npr＝ceil(log2360/Pr)
Figure BDA00030810713400000716
Figure BDA00030810713400000717
其中，Nmr为幅度存储位数，Npr为相位存储位数；
S602.将远场数据离散压缩处理:
第一步、从文件中读取维度为
Figure BDA0003081071340000081
的三维浮点数组
Figure BDA0003081071340000082
该数组由共有
Figure BDA0003081071340000083
个
Figure BDA0003081071340000084
所对应的远场数据
Figure BDA0003081071340000085
构成；每个
Figure BDA0003081071340000086
包括四列数据，分别为：
Mag(E⊥)、Deg(E⊥)、Mag(E||)、Deg(E||)；
并从中恢复出E⊥、E||；
对于每一个
Figure BDA0003081071340000087
所对应的远场数据
Figure BDA0003081071340000088
进行总场计算，总场计算公式：
Etot＝10*log10(|E⊥|2+|E|||2+1e-30)
其中，1e-30的作用是为了保证log10取值的有效性；
第二步、取Etot最大值，记为Emax；
第三步、对任一个远场数据
Figure BDA0003081071340000089
进行处理；
A01.判断是否满足：Etot＜Emax-Dr,若满足，则将E⊥、E||置零，若不满足，则E⊥、E||不变；
A02.对步骤A01得到的数据再次进行处理：
判断是否满足20*log10(|E⊥|+1e-30)<Etot-Sr,若满足，E⊥置0，若不满足，E⊥不变；
判断是否满足20*log10(|E|||+1e-30)<Etot-Sr,若满足，E||置0，若不满足，E||不变；
A03、对于步骤A02得到的数据采用的场值离散化流程进行处理，公式为：
Figure BDA00030810713400000810
Figure BDA00030810713400000811
Figure BDA00030810713400000812
Figure BDA00030810713400000813
进行上述处理后，将每个远场数据
Figure BDA00030810713400000814
作为一行用4*(Nmr+Npr)个二进制数位表示；第四步、将远场数据
Figure BDA00030810713400000815
的处理结果写入二进制文件：
二进制文件中，数据头的数据格式
数位(byte) 4 4 4 4 4
数值 S<sub>r</sub> D<sub>r</sub> M<sub>e</sub> P<sub>e</sub> E<sub>max</sub>
类型 浮点 浮点 浮点 浮点 浮点
数据行的每行数据格式：
Figure BDA0003081071340000091
第五步、对于不同
Figure BDA0003081071340000092
对应的远场数据
Figure BDA0003081071340000093
重复执行第三步到第四步，得到所有远场数据的压缩文件。
上述说明示出并描述了本发明的一个优选实施例，但如前所述，应当理解本发明并非局限于本文所披露的形式，不应看作是对其他实施例的排除，而可用于各种其他组合、修改和环境，并能够在本文所述发明构想范围内，通过上述教导或相关领域的技术或知识进行改动。而本领域人员所进行的改动和变化不脱离本发明的精神和范围，则都应在本发明所附权利要求的保护范围内。

Claims (5)

1.一种电磁环境远场数据的快速无损压缩存储方法，其特征在于：包括以下步骤：
S1.加载CAD网格模型：
从nastran格式的目标对象网格模型文件中提取CAD面元网格模型的点列表和点连接列表，设置材料电磁参数信息和模型激励源参数，完成电磁模型创建工作；
S2.将提取的CAD面元网格模型、设置的材料电磁参数信息，以及设置的激励源参数转换成系统矩阵[A]nbase×nbase和右侧激励项[rhs]nbase；
在进行转换时，将CAD面元网格模型中的点列表和点连接列表通过RWG基函数的格式转换成矩量法所需的基函数，nbase为电磁模型中所有基函数未知量总数，右侧激励项[rhs]nbase中的第m个元素为rhsm，填充公式如下：
Figure FDA0003268283150000011
fm为第m个基函数，r为场点所在区域即fm区域内的位置矢量，
Figure FDA0003268283150000012
为第m个基函数域内的入射场分布；
S3.构建电磁环境的远场模型，通过步骤S2中获得的数据得到每个方向上垂直极化方向
Figure FDA0003268283150000013
水平极化方向
Figure FDA0003268283150000014
的远场；
S4.将步骤S1得到CAD面元网格模型的点列表和点连接列表导出网格数据文件；将步骤S3得到的每个方向上垂直极化方向
Figure FDA0003268283150000015
水平极化方向
Figure FDA0003268283150000016
的远场导出成远场数据文件；
S5.将导出的网格数据文件和远场数据文件读取进内存；
S6.对内存中的远场数据和网格数据进行压缩；
S7.借助于rar压缩工具将得到的远场数据和网格数据进行二次压缩，得到最终的压缩存储文件。
2.根据权利要求1所述的一种电磁环境远场数据的快速无损压缩存储方法，其特征在于：所述步骤S2包括：
S201.将CAD面元网格模型中的点列表和点连接列表通过RWG基函数的格式转换成矩量法所需的基函数，设nbase为电磁模型中所有基函数未知量总数，材料的电磁参数信息包括磁导率μl和介电常数εl；激励源参数包括源点所在基函数区域的位置矢量和入射场分布；
S202.转换得到系统矩阵[A]nbase×nbase：设amn为[A]nbase×nbase的第m行第n列元素，则amn填充公式如下：
Figure FDA0003268283150000017
其中，G为三维格林函数；
Figure FDA0003268283150000021
ω＝2πf0为角频率；f0为频率；
Figure FDA0003268283150000022
为哈密尔顿算子；fm和fn为第m个和第n个基函数；r′为源点所在区域即fn区域内的位置矢量；r为场点所在区域即fm区域内的位置矢量；
S203.转换得到右侧激励项[rhs]nbase：设rhsm为右侧激励项[rhs]nbase中的第m个元素，则填充公式如下：
Figure FDA0003268283150000023
式中，
Figure FDA0003268283150000024
为第m个基函数域内的入射场分布；
S204.通过矩阵求解计算求解[A]nbase×nbasex＝[rhs]nbase得到x，作为每个基函数ibase未知量上的电流量x[ibase]，ibase为1到nbase的整数下标。
3.根据权利要求1所述的一种电磁环境远场数据的快速无损压缩存储方法，其特征在于：所述步骤S3中，通过步骤S2中得到的模型中每个基函数上的电流的积分x[ibase]，得到每个方向上垂直极化方向
Figure FDA0003268283150000025
水平极化方向
Figure FDA0003268283150000026
的远场，针对计算域S′的面积分的远场计算公式为：
Figure FDA0003268283150000027
Figure FDA0003268283150000028
其中，
Figure FDA0003268283150000029
为电磁模型远场的垂直极化的单位方向矢量，
Figure FDA00032682831500000210
为电磁模型远场的水平极化的单位方向矢量，
Figure FDA00032682831500000211
电磁模型远场的径向单位矢量；S′为电磁模型基函数所组成的计算域；r′为S′上任何一点矢量；R为远场的径向距离；k为自由空间波数；n为面元法向单位矢量；k为自由空间波数；
Figure FDA00032682831500000212
4.根据权利要求1所述的一种电磁环境远场数据的快速无损压缩存储方法，其特征在于：所述步骤S4中，远场数据文件采用维度为
Figure FDA00032682831500000213
的三维浮点数组
Figure FDA00032682831500000214
其中，mθ为θ角总数，
Figure FDA00032682831500000215
为
Figure FDA00032682831500000216
角总数；设θ在0～180°内取值，取值间隔为
Figure FDA00032682831500000217
Figure FDA00032682831500000218
在0～360°内取值，取值间隔为
Figure FDA00032682831500000219
则：
θ角的取值θi为：
Figure FDA00032682831500000220
Figure FDA00032682831500000221
角的取值
Figure FDA00032682831500000222
为：
Figure FDA0003268283150000031
每个
Figure FDA0003268283150000032
位置的远场存储格式
Figure FDA0003268283150000033
均包含4列数据：
第一列：Mag(E⊥)；第二列Deg(E⊥)；第三列Mag(E||)；第四列Deg(E||)；
其中，E⊥为垂直极化电场，E||为水平极化电场，Mag(x)为幅度函数，Deg(x)为相位函数；
由于θi有mθ个取值，
Figure FDA0003268283150000034
有
Figure FDA0003268283150000035
个取值，故实际上总共有
Figure FDA0003268283150000036
个
Figure FDA0003268283150000037
又由于
Figure FDA0003268283150000038
均包含4列数据，故形成了维度为
Figure FDA0003268283150000039
的三维浮点数组
Figure FDA00032682831500000310
5.根据权利要求1所述的一种电磁环境远场数据的快速无损压缩存储方法，其特征在于：所述步骤S6包括：
S601.对内存中加载的远场数据格式作如下归一化处理：
按照用户需求设定如下参数：空间分辨率:Sr、动态范围:Dr、幅度分辨率:Mr、相位分辨率:Pr；
其中，20≤Sr≤300；20≤Dr≤300；Dr≤Sr；0.01≤Mr≤1；0.1≤Pe≤5
在计算前对Mr和Pr的值进行处理：
Nmr＝ceil(log2Sr/Me)
Npr＝ceil(log2360/Pr)
Figure FDA00032682831500000311
Figure FDA00032682831500000312
其中，Nmr为幅度存储位数，Npr为相位存储位数；
S602.将远场数据离散压缩处理:
第一步、从文件中读取维度为
Figure FDA00032682831500000313
的三维浮点数组
Figure FDA00032682831500000314
该数组由共有
Figure FDA00032682831500000315
个
Figure FDA00032682831500000316
所对应的远场数据
Figure FDA00032682831500000317
构成；每个
Figure FDA00032682831500000318
包括四列数据，分别为：
Mag(E⊥)、Deg(E⊥)、Mag(E||)、Deg(E||)；
并从中恢复出E⊥、E||；
对于每一个
Figure FDA00032682831500000319
所对应的远场数据
Figure FDA00032682831500000320
进行总场计算，总场计算公式：
Etot＝10*log10(|E⊥|2+|E|||2+1e-30)
其中，1e-30的作用是为了保证log10取值的有效性；
第二步、取Etot最大值，记为Emax；
第三步、对任一个远场数据
Figure FDA0003268283150000041
进行处理；
A01.判断是否满足：Etot<Emax-Dr,若满足，则将E⊥、E||置零，若不满足，则E⊥、E||不变；
A02.对步骤A01得到的数据再次进行处理：
判断是否满足20*log10(|E⊥|+1e-30)<Etot-Sr,若满足，E⊥置0，若不满足，E⊥不变；
判断是否满足20*log10(|E|||+1e-30)<Etot-Sr,若满足，E||置0，若不满足，E||不变；
A03、对于步骤A02得到的数据采用的场值离散化流程进行处理，公式为：
Figure FDA0003268283150000042
Figure FDA0003268283150000043
Figure FDA0003268283150000044
Figure FDA0003268283150000045
进行上述处理后，将每个远场数据
Figure FDA0003268283150000046
作为一行用4*(Nmr+Npr)个二进制数位表示；
第四步、将远场数据
Figure FDA0003268283150000047
的处理结果写入二进制文件；
第五步、对于不同
Figure FDA0003268283150000048
对应的远场数据
Figure FDA0003268283150000049
重复执行第三步到第四步，得到所有远场数据的压缩文件。
CN202110566228.7A 2021-05-24 2021-05-24 一种电磁环境远场数据的快速无损压缩存储方法 Active CN113282553B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110566228.7A CN113282553B (zh) 2021-05-24 2021-05-24 一种电磁环境远场数据的快速无损压缩存储方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110566228.7A CN113282553B (zh) 2021-05-24 2021-05-24 一种电磁环境远场数据的快速无损压缩存储方法

Publications (2)

Publication Number Publication Date
CN113282553A CN113282553A (zh) 2021-08-20
CN113282553B true CN113282553B (zh) 2021-11-02

Family

ID=77281156

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110566228.7A Active CN113282553B (zh) 2021-05-24 2021-05-24 一种电磁环境远场数据的快速无损压缩存储方法

Country Status (1)

Country Link
CN (1) CN113282553B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113987826B (zh) * 2021-11-09 2022-04-26 北京航空航天大学 基于波导端口的天线隔离度全波计算方法
CN114067079B (zh) * 2021-11-19 2022-05-13 北京航空航天大学 复杂曲面电磁波矢量动态可视化方法
CN115358282A (zh) * 2022-10-21 2022-11-18 北京智芯微电子科技有限公司 电磁环境数据压缩方法及装置

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1012778A1 (en) * 1997-07-14 2000-06-28 Sarnoff Corporation Apparatus and method for macroblock based rate control in a coding system
CN103279589A (zh) * 2013-04-18 2013-09-04 南京理工大学 基于矩阵嵌套压缩的旋转对称体电磁散射特性仿真方法
CN105786765A (zh) * 2016-02-25 2016-07-20 南京航空航天大学 一种快速自适应地生成激励无关特征基函数的方法
CN106199524A (zh) * 2016-06-23 2016-12-07 中国科学院电子学研究所 基于基追踪去噪的远场宽带rcs数据采集与压缩方法
CN109376485A (zh) * 2018-12-03 2019-02-22 上海无线电设备研究所 基于aca-mlfma加速的区域分解非共形网格的快速仿真建模方法
US10282500B1 (en) * 2016-10-17 2019-05-07 Northrop Grumman Systems Corporation Efficient modeling of electromagnetic properties of an assembly
CN110990909A (zh) * 2019-11-29 2020-04-10 北京仿真中心 一种三维cad模型数据处理方法

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2387004B1 (en) * 2010-05-11 2016-12-14 Dassault Systèmes Lossless compression of a structured set of floating point numbers, particularly for CAD systems
CN102156764B (zh) * 2010-12-27 2013-06-26 南京理工大学 一种分析天线辐射和电磁散射的多分辨预条件方法
US10275547B2 (en) * 2015-06-03 2019-04-30 The Mathworks, Inc. Method and system for assessing performance of arbitrarily large arrays
CN108038313A (zh) * 2017-12-13 2018-05-15 南京理工大学 一种剖分不均匀的目标电磁散射特性的分析方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1012778A1 (en) * 1997-07-14 2000-06-28 Sarnoff Corporation Apparatus and method for macroblock based rate control in a coding system
CN103279589A (zh) * 2013-04-18 2013-09-04 南京理工大学 基于矩阵嵌套压缩的旋转对称体电磁散射特性仿真方法
CN105786765A (zh) * 2016-02-25 2016-07-20 南京航空航天大学 一种快速自适应地生成激励无关特征基函数的方法
CN106199524A (zh) * 2016-06-23 2016-12-07 中国科学院电子学研究所 基于基追踪去噪的远场宽带rcs数据采集与压缩方法
US10282500B1 (en) * 2016-10-17 2019-05-07 Northrop Grumman Systems Corporation Efficient modeling of electromagnetic properties of an assembly
CN109376485A (zh) * 2018-12-03 2019-02-22 上海无线电设备研究所 基于aca-mlfma加速的区域分解非共形网格的快速仿真建模方法
CN110990909A (zh) * 2019-11-29 2020-04-10 北京仿真中心 一种三维cad模型数据处理方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
A Novel Multilevel Matrix Compression Method for Analysis of Electromagnetic Scattering from PEC Targets;Mengmeng Li等;《IEEE Transactions on Antennas and Propagation》;20111120;第60卷(第3期);1390-1399页 *
Analyzing electromagnetic scattering using characteristic basis function method with compressed sensing;QuanQuan Wang等;《2013 6th international Congress on image and signal processing》;20140220;1633-1637页 *
电磁散射矩阵压缩算法研究;朱伟光;《中国优秀硕士学位论文全文数据库 基础科学辑》;20210515;A005-266 *
非均匀各向异性阻抗目标电磁散射建模方法研究;张凡;《中国博士学位论文全文数据库 基础科学辑》;20180115;A005-37 *

Also Published As

Publication number Publication date
CN113282553A (zh) 2021-08-20

Similar Documents

Publication Publication Date Title
CN113282553B (zh) 一种电磁环境远场数据的快速无损压缩存储方法
CN113297763B (zh) 一种适用于矩量法的近场数据快速无损压缩存储方法
KR102541743B1 (ko) 학습 장치, 추론 장치 및 학습 완료 모델
CN106991222B (zh) 一种基于叠层矩阵分解的低频电磁特性仿真方法
CN115618479A (zh) 一种基于bim的幕墙高封闭性安装方法
CN106909539A (zh) 图像检索系统、服务器、数据库及相关的方法
CN116030923B (zh) 获取材料动态本构关系的方法、装置、设备及存储介质
CN114330042A (zh) 基于sn曲线的程序载荷谱编制方法、系统和存储介质
CN112017219A (zh) 一种激光点云配准方法
CN112825199A (zh) 碰撞检测方法、装置、设备及存储介质
CN105786760B (zh) 一种基于稀疏块状矩阵压缩存储结构的预条件共轭梯度区域网平差方法
CN107194994B (zh) 一种无标定曲面点云数据重建圆柱面的方法及装置
CN109544557A (zh) 基于区块的主成分分析转换方法及其装置
US9898834B2 (en) Method and apparatus for generating a bitstream of repetitive structure discovery based 3D model compression
CN115457202B (zh) 一种三维模型更新的方法、装置及存储介质
CN116957955A (zh) 晶圆图像校正方法、装置、计算机设备和存储介质
CN113776632B (zh) 称重偏载补偿方法、装置、计算机设备和存储介质
CN116662463A (zh) Ndt地图数据的压缩方法、装置及计算机程序
CN108319717A (zh) 一种矢量点转换为栅格点时最佳栅格尺寸范围的计算方法
CN115049813A (zh) 一种基于一阶球谐的粗配准方法、装置及系统
CN116385571B (zh) 一种基于多维化动态可变分辨率的点云压缩方法及系统
CN116168300B (zh) 一种高分辨率遥感影像等gsd的核线重采样方法和装置
CN114034938B (zh) 一种介质空间中天线性能评估方法、装置及电子设备
CN115063432B (zh) 对准标记图形的识别方法及装置、计算机设备和存储介质
CN114579926B (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