CN103091570B - 一种人工电磁材料的结构单元参数仿真实验方法及装置 - Google Patents
一种人工电磁材料的结构单元参数仿真实验方法及装置 Download PDFInfo
- Publication number
- CN103091570B CN103091570B CN201110338295.XA CN201110338295A CN103091570B CN 103091570 B CN103091570 B CN 103091570B CN 201110338295 A CN201110338295 A CN 201110338295A CN 103091570 B CN103091570 B CN 103091570B
- Authority
- CN
- China
- Prior art keywords
- geometric parameter
- structural unit
- geometric
- parameter
- value
- 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
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提供了一种人工电磁材料的结构单元参数仿真实验方法和装置,该方法包括步骤:获取结构单元的第一几何参数;对多个水平的多个第一几何参数进行均匀设计,获取第一集合;获取结构单元的第二几何参数值;对获取的多个水平的多个第二几何参数与第一集合进行正交设计,获取第二集合;获取结构单元的第三几何参数的值;对获取的多个水平的多个第三几何参数与第二集合进行正交设计,获取第三集合。本发明的人工电磁材料的结构单元参数仿真实验方法及装置在保证实验点选取的空间“代表性”的同时,能较快、较准确地找到符合初始实验假设的人工电磁材料的电磁特性所对应的结构单元的参数,减少实验次数,提高了工作效率并节约了实验成本。
Description
技术领域
本发明涉及人工电磁材料领域,特别涉及一种人工电磁材料的结构单元参数仿真实验方法及装置。
背景技术
人工电磁材料(Metamaterials,或称超材料)是将具有特定几何形状的亚波长宏观结构单元周期性或非周期性地排列所构成的人工材料,其媒质特性取决于其结构单元的结构以及空间分布。
测量人工电磁材料结构单元的电磁特性,是人工电磁材料设计过程中的一个重要环节。而如何选取一定尺寸的人工电磁材料结构单元用于电磁仿真测量又是整个设计过程中必须要解决的问题。单个微结构单元的几何形状信息由一组参数描述,参数的取值代表其形状的大小。在合理的参数取值范围内,每个参数又对应多个取值水平。为了获得特定类型微结构单元的电磁响应特征规律,需要进行人工电磁材料结构单元电磁特性仿真实验,由于某些结构单元的复杂拓扑结构对应的几何参数较多,参数之间又相互约束,欲探明其电磁特性,常需要做极其大量的实验。
如将各参数各水平一一对应测量衍生材料的电磁特性,需要做出的拓扑结构模板将以百万计,如此大的实验数量给具体的实施带来了很大的难度。为了优化实验步骤,目前常用的方法仅为依靠之前积累起来的实验经验及对电磁原理方面的基本认知,对结构单元每个参数所允许的范围内,对其进行微调。该方法全凭经验,且不可预测整个探索过程所需时间,很难直接实现实验步骤的优化。
因此,需要探明一种人工电磁材料的结构单元参数仿真实验方法,用尽量少的实验次数获得对特定类型结构单元的电磁响应规律,对了解如何拓展对该微结构的应用有着十分重要的意义,另一方面也对促进人工电磁材料的产业化进程有着十分积极的作用。
发明内容
本发明主要解决的技术问题是提供一种人工电磁材料的结构单元参数仿真实验方法及装置,该方法及装置能较快、较准确地找到符合初始实验假设的人工电磁材料的电磁特性所对应的结构单元的参数。
本发明提供了一种人工电磁材料的结构单元参数仿真实验方法,该人工电磁材料的结构单元参数仿真实验方法包括步骤:获取结构单元的第一几何参数;对多个水平的多个第一几何参数进行均匀设计,获取第一集合;获取结构单元的第二几何参数,第二几何参数的取值范围受第一几何参数的取值约束,根据第一集合以及第一几何参数与第二几何参数之间的约束关系获取多个第二几何参数的值;对获取的多个水平的多个第二几何参数与第一集合进行正交设计,获取第二集合;获取结构单元的第三几何参数,第三几何参数的取值范围受第一几何参数与第二几何参数的取值约束,根据第二集合以及第三几何参数与第一几何参数及第二几何参数之间的约束关系获取第三几何参数的值;对获取的多个水平的多个第三几何参数与第二集合进行正交设计,获取第三集合。
根据本发明的一优选实施例,人工电磁材料的结构单元参数仿真实验方法包括步骤:根据第三集合进行人工电磁材料的电磁特性的仿真实验,并获取仿真实验结果以及仿真实验结果极值所对应的几何参数;对实验结果采用t检验进行分析,获取t检验结果大于预设显著性阈值的几何参数。
根据本发明的一优选实施例,根据第一集合以及第一几何参数与第二几何参数之间的约束关系获取第二几何参数的值的步骤包括如下步骤:根据第一集合的每一行中的第一几何参数的值以及第一几何参数与第二几何参数之间的约束关系获取第二几何参数的取值范围,并在每个第二几何参数的取值范围内对每个第二几何参数均匀取点。
根据本发明的一优选实施例,根据第二集合以及第三几何参数与第一几何参数及第二几何参数之间的约束关系获取第三几何参数的值的步骤包括:根据第二集合的每一行中的第一几何参数及第二几何参数的值以及第三几何参数与第一几何参数及第二几何参数之间的约束关系获取第三几何参数的取值范围,并在每个第三几何参数的取值范围内对每个第三几何参数进行均匀取点。
根据本发明的一优选实施例,在对多个水平的多个第一几何参数进行均匀设计,获取第一集合的步骤中,均匀设计采用好格子点法。
根据本发明的一优选实施例,对于n次电磁特性的仿真实验,好格子点法包括如下步骤:采用同余运算法生成均匀设计表的第j列,其公式如下表示:
uij=ihj[modn]
其中,i为均匀设计表的行数,n为各几何参数对应的采样样本点个数的最小公倍数,h为比n小的且与n互质的正整数组成的向量,即h=(h1,h2,...,hm),m由欧拉函数φ(n)决定,
当ihj大于n时,将ihj减去n的一个倍数使得差值落在[1,n]之中,均匀设计表的其他列可由如下公式递推生成:
u1j=hj
其中,i=1,…,n-1;
运用中心化L2偏差测度均匀设计表的均匀性,取偏差最小的均匀设计表作为第一集合,记做Un(ns)。
本发明还提供了一种人工电磁材料的结构单元参数仿真实验装置,该人工电磁材料的结构单元参数仿真实验装置包括:第一几何参数获取单元,获取结构单元的第一几何参数;第一集合获取单元,对多个水平的多个第一几何参数进行均匀设计,获取第一集合;第二几何参数获取单元,获取结构单元的第二几何参数,第二几何参数的取值范围受第一几何参数的取值约束,根据第一集合以及第一几何参数与第二几何参数之间的约束关系获取多个第二几何参数的值;第二集合获取单元,对获取的多个水平的多个第二几何参数与第一集合进行正交设计,获取第二集合;第三几何参数获取单元,获取结构单元的第三几何参数,第三几何参数的取值范围受第一几何参数与第二几何参数的取值约束,根据第二集合以及第三几何参数与第一几何参数及第二几何参数之间的约束关系获取第三几何参数的值;第三集合获取单元,对获取的多个水平的多个第三几何参数与第二集合进行正交设计,获取第三集合。
根据本发明的一优选实施例,人工电磁材料的结构单元参数仿真实验装置还包括:仿真实验单元,根据第三集合进行人工电磁材料的电磁特性的仿真实验,并获取仿真实验结果以及仿真实验结果极值所对应的几何参数;显著性参数获取单元,对实验结果采用t检验进行分析,获取t检验结果大于预设显著性阈值的几何参数。
根据本发明的一优选实施例,第二几何参数获取单元用于:根据第一集合的每一行中的第一几何参数的值以及第一几何参数与第二几何参数之间的约束关系获取第二几何参数的取值范围,并在每个第二几何参数的取值范围内对每个第二几何参数均匀取点。
根据本发明的一优选实施例,第三几何参数获取单元用于:根据第二集合的每一行中的第一几何参数及第二几何参数的值以及第三几何参数与第一几何参数及第二几何参数之间的约束关系获取第三几何参数的取值范围,并在每个第三几何参数的取值范围内对每个第三几何参数进行均匀取点。
根据本发明的一优选实施例,第一集合获取单元采用的均匀设计为好格子点法。
根据本发明的一优选实施例,对于n次电磁特性的仿真实验,好格子点法包括如下步骤:采用同余运算法生成均匀设计表的第j列,其公式如下表示:
uij=ihj[modn]
其中,i为均匀设计表的行数,n为各几何参数对应的采样样本点个数的最小公倍数,h为比n小的且与n互质的正整数组成的向量,即h=(h1,h2,...,hm),m由欧拉函数φ(n)决定,
当ihj大于n时,将ihj减去n的一个倍数使得差值落在[1,n]之中,均匀设计表的其他列可由如下公式递推生成:
u1j=hj
其中,i=1,…,n-1;运用中心化L2偏差测度均匀设计表的均匀性,取偏差最小的均匀设计表作为第一集合,记做Un(ns)。
本发明的有益效果是:区别于现有技术的情况,本发明的人工电磁材料的结构单元参数仿真实验方法及装置在保证实验点选取的空间“代表性”的同时,通过采用前次选取的实验方案的对应参数选取下次实验方案的对应参数的方式能较快、较准确地找到符合初始实验假设的人工电磁材料的电磁特性所对应的结构单元的参数,减少实验次数,提高了工作效率并节约了实验成本。
附图说明
图1是具有“电场耦合谐振器”型拓扑结构的人工电磁材料结构单元的立体结构示意图;
图2是本发明实施例的人工电磁材料的结构单元参数仿真实验方法的流程示意图;
图3是采用本发明实施例的人工电磁材料的结构单元参数仿真实验方法进行仿真获得的具有“电场耦合谐振器”型拓扑结构的人工电磁材料Str0的频段变化与损耗系数关系图;
图4是采用本发明实施例的人工电磁材料的结构单元参数仿真实验方法进行仿真获得的具有“电场耦合谐振器”型拓扑结构的人工电磁材料Str1的频段变化与损耗系数关系图;以及
图5是本发明实施例的人工电磁材料的结构单元参数仿真实验装置的结构示意图。
具体实施方式
下面结合附图和实施例对本发明进行详细说明。
人工电磁材料的单个微结构单元的几何形状信息由一组参数描述,复杂拓扑结构对应的几何参数较多,参数之间常常相互约束。以图1所示出的具有“电场耦合谐振器”型拓扑结构的人工电磁材料结构单元为例。如图1所示,该结构单元的拓扑结构中有7个对应的几何参数{L,L1,ax,ay,cellx,celly,hs}。
其中,cellx=celly,cellx、celly及hs之间的约束关系不明显,例如,可以定义参数cellx与celly的范围均为[1.5,5],厚度hs的范围为[1,4]。线宽W固定为0.1。而L,L1,ax及ay则分别满足如下关系式:
max(2L1+0.5,2L+0.5)≤ax≤min(cellx-0.1,celly-0.1)
max(2L1+0.5,2L+0.5)≤ay≤min(cellx-0.1,celly-0.1)
针对这种结构,在本发明中,cellx、celly及hs等参数的取值范围不受明显约束,因此,这类参数将采取均匀设计原理,使用Good lattice point法构造均匀设计表,设计相关实验安排方案,同时将该类参数定义为第一几何参数。L、L1、ax及ay等参数的取值范围受明显约束,L及L1取值范围受cellx、celly的取值约束,ax及ay取值范围受cellx、celly、L及L1的取值约束,因此,这类约束性参数将采取正交试验设计和分区设计定义实验安排方案,同时将L、L1类参数定义为第二几何参数,ax及ay等参数定义为第三几何参数,第二几何参数的取值范围受第一几何参数的取值约束,第三几何参数的取值范围受第一几何参数与第二几何参数的取值约束。
图2是本发明实施例的人工电磁材料的结构单元参数仿真实验方法的流程示意图。如图2所示,一种人工电磁材料的结构单元参数仿真实验方法包括步骤S1、步骤S2、步骤S3、步骤S4、步骤S5、步骤S6、步骤S7及步骤S8,具体如下:
步骤S1:获取结构单元的第一几何参数。
步骤S2:对多个水平的多个第一几何参数进行均匀设计,获取第一集合。
在本实施例中,步骤S2中的均匀设计采用好格子点法进行。对于n次电磁特性的仿真实验,好格子点法包括如下步骤:
采用同余运算法生成均匀设计表的第j列,其公式如下表示:
uij=ihj[modn]
其中,i为均匀设计表的行数,n为各几何参数对应的采样样本点个数的最小公倍数,h为比n小的且与n互质的正整数组成的向量,即h=(h1,h2,...,hm),m由欧拉函数φ(n)决定,
当ihj大于n时,将ihj减去n的一个倍数使得差值落在[1,n]之中,均匀设计表的其他列可由如下公式递推生成:
u1j=hj
其中,i=1,…,n-1;
运用中心化L2偏差测度均匀设计表的均匀性,取偏差最小的均匀设计表作为第一集合,记做Un(ns),并将Un(ns)定义为“代表性”前提条件。
步骤S3:获取结构单元的第二几何参数,第二几何参数的取值范围受第一几何参数的取值约束,根据第一集合以及第一几何参数与第二几何参数之间的约束关系获取多个第二几何参数的值。
在本实施例中,步骤S3中的根据第一集合以及第一几何参数与第二几何参数之间的约束关系获取第二几何参数的值的步骤包括如下步骤:根据第一集合的每一行中的第一几何参数的值以及第一几何参数与第二几何参数之间的约束关系获取第二几何参数的取值范围,并在每个第二几何参数的取值范围内对每个第二几何参数均匀取点。具体的,将Un(ns)的每一行定义成一个“代表性”前提条件,在每一个“代表性”前提下,分别获得L,L1的取值范围[Lbegin,Lend]&[L1 begin,L1 end],在这个范围内,采用均匀设计进行取点。
步骤S4:对获取的多个水平的多个第二几何参数与第一集合进行正交设计,获取第二集合。
具体的,依照传统“Lo.t0正交设计”实验构造思想(任两因素的水平组合做相同数量的实验),将该步骤所得内容与步骤S2中第一集合内容对应形成一个正交设计的排列,记做第二集合U2。
步骤S5:获取结构单元的第三几何参数,第三几何参数的取值范围受第一几何参数与第二几何参数的取值约束,根据第二集合以及第三几何参数与第一几何参数及第二几何参数之间的约束关系获取第三几何参数的值;
其中,在步骤S5中,根据第二集合以及第三几何参数与第一几何参数及第二几何参数之间的约束关系获取第三几何参数的值的步骤包括:根据第二集合的每一行中的第一几何参数及第二几何参数的值以及第三几何参数与第一几何参数及第二几何参数之间的约束关系获取第三几何参数的取值范围,并在每个第三几何参数的取值范围内对每个第三几何参数进行均匀取点。具体的,重新定义第二集合U2的每一行为一组新的“代表性”前提条件,在每一个“代表性”前提下,分别获得ax及ay的取值范围在这个范围内,采用分区设计进行取点,即在和区域内分别均匀采α个试验点。
步骤S6:对获取的多个水平的多个第三几何参数与所述第二集合进行正交设计,获取第三集合;
具体的,采取“Lo.t0正交设计”,将该步骤所得内容与第二集合U2中内容对应形成一个正交设计的排列,记做第三集合U3。
步骤S7:根据第三集合进行人工电磁材料的电磁特性的仿真实验,并获取仿真实验结果以及仿真实验结果极值所对应的几何参数。
步骤S7具体为:依照U3进行电磁仿真,并将代表每个拓朴结构模板损耗的结果S2,1读取出来,建立频段变化与S2,1关系的坐标系。
然后将每个坐标系中大于平均值的部分除去,取出余下部分的平均值A0,找出最小A0建立与A0对应的结构参数关系Str0(该结构相关频段变化与S2,1关系如图3的S1,1线所示,希望看到S1,1线整体趋势接近于0‘顶部区域’)。
步骤S8:对实验结果采用t检验进行分析,获取t检验结果大于预设显著性阈值的几何参数。
步骤S8中采用“t检验”测出对结果影响最大的几个参数,在参数允许的范围进行微调,以此获得最佳结果对应的参数。具体为:通过“t检验”可得知,各参数变化对S2,1值的影响如下表所示(‘1’表示最大,‘5’表示最小,‘+’正比关系,‘-’反比关系)。
参数 | cellx | L | L1 | ax | ay |
影响 | 1 | 4 | 5 | 2 | 3 |
关系 | + | - | + | - | - |
然后根据上表信息,对Str0中的结构参数值进行微调,并找出最优结构Str1(该结构相关频段变化与S2,1关系如图4的S1,1线所示)。
图5是本发明实施例的人工电磁材料的结构单元参数仿真实验装置的结构示意图。如图5所示,人工电磁材料的结构单元参数仿真实验装置100包括第一几何参数获取单元110、第一集合获取单元120、第二几何参数获取单元130、第二集合获取单元140、第三几何参数获取单元150、第三集合获取单元160、仿真实验单元170及显著性参数获取单元180。
其中,第一几何参数获取单元110,获取结构单元的第一几何参数。
第一集合获取单元120,对多个水平的多个第一几何参数进行均匀设计,获取第一集合。采用的均匀设计为好格子点法。好格子点法的具体技术特征请参照步骤S2中所述,此处不再赘述。
第二几何参数获取单元130,获取结构单元的第二几何参数,第二几何参数的取值范围受第一几何参数的取值约束,根据第一集合以及第一几何参数与第二几何参数之间的约束关系获取多个第二几何参数的值。其中,根据第一集合以及第一几何参数与第二几何参数之间的约束关系获取多个第二几何参数的值具体为:根据第一集合的每一行中的第一几何参数的值以及第一几何参数与第二几何参数之间的约束关系获取第二几何参数的取值范围,并在每个第二几何参数的取值范围内对每个第二几何参数均匀取点。
第二集合获取单元140,对获取的多个水平的多个第二几何参数与第一集合进行正交设计,获取第二集合。
第三几何参数获取单元150,获取结构单元的第三几何参数,第三几何参数的取值范围受第一几何参数与第二几何参数的取值约束,根据第二集合以及第三几何参数与第一几何参数及第二几何参数之间的约束关系获取第三几何参数的值。其中,根据第二集合以及第三几何参数与第一几何参数及第二几何参数之间的约束关系获取第三几何参数的值具体为:根据第二集合的每一行中的第一几何参数及第二几何参数的值以及第三几何参数与第一几何参数及第二几何参数之间的约束关系获取第三几何参数的取值范围,并在每个第三几何参数的取值范围内对每个第三几何参数进行均匀取点。
第三集合获取单元160,对获取的多个水平的多个第三几何参数与第二集合进行正交设计,获取第三集合。
仿真实验单元170,根据第三集合进行人工电磁材料的电磁特性的仿真实验,并获取仿真实验结果以及仿真实验结果极值所对应的几何参数。
显著性参数获取单元180,对实验结果采用t检验进行分析,获取t检验结果大于预设显著性阈值的几何参数。
本发明的有益效果是:区别于现有技术的情况,本发明的人工电磁材料的结构单元参数仿真实验方法及装置在保证实验点选取的空间“代表性”的同时,通过采用前次选取的实验方案的对应参数选取下次实验方案的对应参数的方式能较快、较准确地找到符合初始实验假设的人工电磁材料的电磁特性所对应的结构单元的参数,减少实验次数,提高了工作效率并节约了实验成本。
以上所述仅为本发明的实施例,并非因此限制本发明的专利范围,凡是利用本发明说明书及附图内容所作的等效结构或等效流程变换,或直接或间接运用在其他相关的技术领域,均同理包括在本发明的专利保护范围内。
Claims (12)
1.一种人工电磁材料的结构单元参数仿真实验方法,其特征在于,所述人工电磁材料的结构单元参数仿真实验方法包括步骤:
获取所述结构单元的第一几何参数;
对多个水平的多个第一几何参数进行均匀设计,获取第一集合;
获取所述结构单元的第二几何参数,所述第二几何参数的取值范围受所述第一几何参数的取值约束,根据所述第一集合以及所述第一几何参数与所述第二几何参数之间的约束关系获取多个第二几何参数的值;
对获取的多个水平的多个第二几何参数与所述第一集合进行正交设计,获取第二集合;
获取所述结构单元的第三几何参数,所述第三几何参数的取值范围受所述第一几何参数与所述第二几何参数的取值约束,根据所述第二集合以及所述第三几何参数与所述第一几何参数及所述第二几何参数之间的约束关系获取所述第三几何参数的值;
对获取的多个水平的多个第三几何参数与所述第二集合进行正交设计,获取第三集合。
2.根据权利要求1所述的人工电磁材料的结构单元参数仿真实验方法,其特征在于,所述人工电磁材料的结构单元参数仿真实验方法包括步骤:
根据所述第三集合进行所述人工电磁材料的电磁特性的仿真实验,并获取仿真实验结果以及所述仿真实验结果极值所对应的几何参数;
对所述实验结果采用t检验进行分析,获取t检验结果大于预设显著性阈值的几何参数。
3.根据权利要求1所述的人工电磁材料的结构单元参数仿真实验方法,其特征在于,所述根据所述第一集合以及所述第一几何参数与所述第二几何参数之间的约束关系获取第二几何参数的值的步骤包括如下步骤:
根据所述第一集合的每一行中的第一几何参数的值以及所述第一几何参数与所述第二几何参数之间的约束关系获取第二几何参数的取值范围,并在每个第二几何参数的取值范围内对每个第二几何参数均匀取点。
4.根据权利要求1所述的人工电磁材料的结构单元参数仿真实验方法,其特征在于,根据所述第二集合以及所述第三几何参数与所述第一几何参数及所述第二几何参数之间的约束关系获取所述第三几何参数的值的步骤包括:
根据所述第二集合的每一行中的第一几何参数及第二几何参数的值以及所述第三几何参数与所述第一几何参数及所述第二几何参数之间的约束关系获取第三几何参数的取值范围,并在每个第三几何参数的取值范围内对每个第三几何参数进行均匀取点。
5.根据权利要求1所述的人工电磁材料的结构单元参数仿真实验方法,其特征在于,在所述对多个水平的多个第一几何参数进行均匀设计,获取第一集合的步骤中,所述均匀设计采用好格子点法。
6.根据权利要求5所述的人工电磁材料的结构单元参数仿真实验方法,其特征在于,对于n次电磁特性的仿真实验,所述好格子点法包括如下步骤:
采用同余运算法生成均匀设计表的第j列,其公式如下表示:
uij=ihj[modn]
其中,i为均匀设计表的行数,n为各几何参数对应的采样样本点个数的最小公倍数,h为比n小的且与n互质的正整数组成的向量,即h=(h1,h2,...,hm),m由欧拉函数φ(n)决定,当ihj大于n时,将ihj减去n的一个倍数使得差值落在[1,n]之中,均匀设计表的其他列可由如下公式递推生成:
u1j=hj
其中,i=1,…,n-1;
运用中心化L2偏差测度均匀设计表的均匀性,取偏差最小的均匀设计表作为所述第一集合,记做Un(ns)。
7.一种人工电磁材料的结构单元参数仿真实验装置,其特征在于,所述人工电磁材料的结构单元参数仿真实验装置包括:
第一几何参数获取单元,获取所述结构单元的第一几何参数;
第一集合获取单元,对多个水平的多个第一几何参数进行均匀设计,获取第一集合;
第二几何参数获取单元,获取所述结构单元的第二几何参数,所述第二几何参数的取值范围受所述第一几何参数的取值约束,根据所述第一集合以及所述第一几何参数与所述第二几何参数之间的约束关系获取多个第二几何参数的值;
第二集合获取单元,对获取的多个水平的多个第二几何参数与所述第一集合进行正交设计,获取第二集合;
第三几何参数获取单元,获取所述结构单元的第三几何参数,所述第三几何参数的取值范围受所述第一几何参数与所述第二几何参数的取值约束,根据所述第二集合以及所述第三几何参数与所述第一几何参数及所述第二几何参数之间的约束关系获取所述第三几何参数的值;
第三集合获取单元,对获取的多个水平的多个第三几何参数与所述第二集合进行正交设计,获取第三集合。
8.根据权利要求7所述的人工电磁材料的结构单元参数仿真实验装置,其特征在于,所述人工电磁材料的结构单元参数仿真实验装置还包括:
仿真实验单元,根据所述第三集合进行所述人工电磁材料的电磁特性的仿真实验,并获取仿真实验结果以及所述仿真实验结果极值所对应的几何参数;显著性参数获取单元,对所述实验结果采用t检验进行分析,获取t检验结果大于预设显著性阈值的几何参数。
9.根据权利要求7所述的人工电磁材料的结构单元参数仿真实验装置,其特征在于,所述第二几何参数获取单元用于:
根据所述第一集合的每一行中的第一几何参数的值以及所述第一几何参数与所述第二几何参数之间的约束关系获取第二几何参数的取值范围,并在每个第二几何参数的取值范围内对每个第二几何参数均匀取点。
10.根据权利要求7所述的人工电磁材料的结构单元参数仿真实验装置,其特征在于,所述第三几何参数获取单元用于:
根据所述第二集合的每一行中的第一几何参数及第二几何参数的值以及所述第三几何参数与所述第一几何参数及所述第二几何参数之间的约束关系获取第三几何参数的取值范围,并在每个第三几何参数的取值范围内对每个第三几何参数进行均匀取点。
11.根据权利要求7所述的人工电磁材料的结构单元参数仿真实验装置,其特征在于,所述第一集合获取单元采用的所述均匀设计为好格子点法。
12.根据权利要求11所述的人工电磁材料的结构单元参数仿真实验装置,其特征在于,对于n次电磁特性的仿真实验,所述好格子点法包括如下步骤:
采用同余运算法生成均匀设计表的第j列,其公式如下表示:
uij=ihj[modn]
其中,i为均匀设计表的行数,n为各几何参数对应的采样样本点个数的最小公倍数,h为比n小的且与n互质的正整数组成的向量,即h=(h1,h2,...,hm),m由欧拉函数φ(n)决定,当ihj大于n时,将ihj减去n的一个倍数使得差值落在[1,n]之中,均匀设计表的其他列可由如下公式递推生成:
u1j=hj
其中,i=1,…,n-1;
运用中心化L2偏差测度均匀设计表的均匀性,取偏差最小的均匀设计表作为所述第一集合,记做Un(ns)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110338295.XA CN103091570B (zh) | 2011-10-31 | 2011-10-31 | 一种人工电磁材料的结构单元参数仿真实验方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110338295.XA CN103091570B (zh) | 2011-10-31 | 2011-10-31 | 一种人工电磁材料的结构单元参数仿真实验方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103091570A CN103091570A (zh) | 2013-05-08 |
CN103091570B true CN103091570B (zh) | 2015-09-16 |
Family
ID=48204380
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201110338295.XA Active CN103091570B (zh) | 2011-10-31 | 2011-10-31 | 一种人工电磁材料的结构单元参数仿真实验方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103091570B (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6157900A (en) * | 1994-10-19 | 2000-12-05 | Intellisense Corp. | Knowledge based system and method for determining material properties from fabrication and operating parameters |
CN101105841A (zh) * | 2007-02-12 | 2008-01-16 | 浙江大学 | 由大规模基因芯片表达谱数据构建基因调控亚网络的方法 |
CN101655525A (zh) * | 2009-07-27 | 2010-02-24 | 肖怀宝 | 基于支持向量机(svm)的人工电磁材料电磁参数提取方法 |
CN101923590A (zh) * | 2010-08-16 | 2010-12-22 | 北京理工大学 | 一种高效的拉丁超立方试验设计方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8638505B2 (en) * | 2008-05-30 | 2014-01-28 | The Invention Science Fund 1 Llc | Negatively-refractive focusing and sensing apparatus, methods, and systems |
-
2011
- 2011-10-31 CN CN201110338295.XA patent/CN103091570B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6157900A (en) * | 1994-10-19 | 2000-12-05 | Intellisense Corp. | Knowledge based system and method for determining material properties from fabrication and operating parameters |
CN101105841A (zh) * | 2007-02-12 | 2008-01-16 | 浙江大学 | 由大规模基因芯片表达谱数据构建基因调控亚网络的方法 |
CN101655525A (zh) * | 2009-07-27 | 2010-02-24 | 肖怀宝 | 基于支持向量机(svm)的人工电磁材料电磁参数提取方法 |
CN101923590A (zh) * | 2010-08-16 | 2010-12-22 | 北京理工大学 | 一种高效的拉丁超立方试验设计方法 |
Non-Patent Citations (1)
Title |
---|
一种新型负折射率微波媒质的设计与仿真模拟;曹云建 等;《科学通报》;20061130;第51卷(第22期);第2612页至2617页 * |
Also Published As
Publication number | Publication date |
---|---|
CN103091570A (zh) | 2013-05-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104200087B (zh) | 用于机器学习的参数寻优及特征调优的方法及系统 | |
CN103675629B (zh) | 基于电场特征量的电极起晕电压预测方法 | |
Gu et al. | Back analysis of mechanical parameters of roller compacted concrete dam | |
CN109033705A (zh) | 一种空间褶皱薄膜弹性应变能和主应力有限元分析方法 | |
CN101788425A (zh) | 一种结构件复合型裂纹前缘应力强度因子分离和分布的确定方法 | |
CN103091570B (zh) | 一种人工电磁材料的结构单元参数仿真实验方法及装置 | |
Barndorff-Nielsen et al. | Bipower variation for Gaussian processes with stationary increments | |
CN109918829A (zh) | 一种针对结构的轴向模态识别方法 | |
CN104077479B (zh) | 一种基于守恒迎风格式获取参量阵声场空间分布的方法 | |
CN108984914A (zh) | 一种求解复杂时域电磁问题的快速多尺度计算方法 | |
CN103091569B (zh) | 一种人工电磁材料的结构单元参数仿真实验方法及装置 | |
CN104297706A (zh) | 一种基于谐振的磁场发生装置及其设计方法 | |
CN104581734A (zh) | 一种网络社区划分方法及装置 | |
CN103136401B (zh) | 一种人造微结构的试验设计表获取方法及装置 | |
CN114912302A (zh) | 一种基于模态空间采样算法的加筋板模态密度获取方法 | |
CN103136396B (zh) | 一种超材料结构单元参数仿真选取方法及装置 | |
CN103135004B (zh) | 一种试验表构造的方法、一种测量电磁特性的方法及其装置 | |
CN110896218B (zh) | 一种用于建立集合性居民负荷的谐波建模方法及系统 | |
CN102768312B (zh) | 一种人工电磁材料单元的试验点选取方法及装置 | |
Beirão da Veiga et al. | A-priori and a-posteriori error analysis for a family of Reissner–Mindlin plate elements | |
CN107526856B (zh) | 可并行的显隐式混合不连续伽辽金时域有限元法 | |
CN102768688B (zh) | 一种超材料单元结构的推演方法和装置 | |
CN102890202B (zh) | 一种人工电磁材料单元的试验点选取方法及装置 | |
CN105868456B (zh) | 基于过滤器技术与细分矩形算法的飞行器约束优化方法 | |
Ayache et al. | Stochastic volatility and multifractional Brownian motion |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |