CN103136401B - The EXPERIMENTAL DESIGN list acquirement method of a kind of man-made microstructure and device - Google Patents
The EXPERIMENTAL DESIGN list acquirement method of a kind of man-made microstructure and device Download PDFInfo
- Publication number
- CN103136401B CN103136401B CN201110396037.7A CN201110396037A CN103136401B CN 103136401 B CN103136401 B CN 103136401B CN 201110396037 A CN201110396037 A CN 201110396037A CN 103136401 B CN103136401 B CN 103136401B
- Authority
- CN
- China
- Prior art keywords
- value
- geometric parameter
- factor
- parameter
- test
- 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
The embodiment of the invention discloses EXPERIMENTAL DESIGN list acquirement method and the device of a kind of man-made microstructure, described method, comprise determining that the number of the geometric parameter factor of man-made microstructure to be measured, the value number of each geometric parameter factor and the parameter value of correspondence thereof, and obtain initial trial design table according to Latin Hypercube Sampling method simultaneously；Electromagnetic response delta data when testing according to described initial trial design table, obtains active geometric parameter scope in each geometric parameter factor；According to Latin Hypercube Sampling method and the geometric parameter scope of each geometric parameter factor of acquisition, obtain test parameters value；Test parameters value structure EXPERIMENTAL DESIGN table according to uniform Design principle and described acquisition.Use the present invention, the test procedure can be reduced, it is achieved in the test that quantity is few, draw distortionless result of the test in the case of not affecting result of the test.
Description
Technical field
The present invention relates to man-made microstructure test field, particularly relate to the EXPERIMENTAL DESIGN table of a kind of man-made microstructure
Acquisition methods and device.
Background technology
Artificial electromagnetic material is also referred to as Meta Materials, is a kind of synthetic material that electromagnetic wave can produce response
Material, is made up of substrate and the man-made microstructure being attached on substrate, owing to man-made microstructure is typically by conducting electricity
The structure with certain geometrical pattern that material arrangement becomes, therefore, it is possible to electromagnetic wave is produced response, so that
Meta Materials entirety embodies the electromagnetic property being different from substrate.Man-made microstructure is that the performance of electromagnetic material plays
Vital effect.
Measure the electromagnetic property of man-made microstructure, be an important step in artificial electromagnetic material design process.
In order to obtain the electromagnetic response characteristic rule to certain types of man-made microstructure, it is often necessary to do substantial amounts of reality
Test.At present, " orthogonal design " that this field generally uses obtains EXPERIMENTAL DESIGN table, then according to this examination
Test design table and carry out the electromagnetic response characteristic rule of the test certain types of man-made microstructure of acquisition.
As a example by the man-made microstructure of " snowflake " shape shown in Fig. 1, in prior art, use orthogonal design
Method generates the flow process of EXPERIMENTAL DESIGN table:
Determine the geometric parameter factor and the value number thereof of the man-made microstructure of " snowflake " shape, and determine each
Parameter value.It is determined here that the geometric parameter factor of " snowflake " shape man-made microstructure includes in Fig. 1: a1,
A2, b1, b2, calculate for convenience, can defined parameters factor s1 and s2 again, wherein, described s1=a1+b1,
Described s2=a2+b2.
Can first determine a numerical range, then in this numerical range, determine each factor
Corresponding each parameter value.
Herein, the span of certainty factor s1 is [1.60,3.00], and value 4, each parameter value
For s1=[1.6000,2.0667,2.5333,3.0000], the span of factor s2 is [0.05,1.45],
And value 4, each parameter value is s2=[0.0500,0.5167,0.9833,1.4500].
According to orthogonal design, available following orthogonal table:
Table 1:
| Sequence number | Parameter 1 | Parameter 2 |
| 1 | 1 | 1 |
| 2 | 1 | 2 |
| 3 | 1 | 3 |
| 4 | 1 | 4 |
| 5 | 2 | 1 |
| 6 | 2 | 2 |
| 7 | 2 | 3 |
| 8 | 2 | 4 |
| 9 | 3 | 1 |
| 10 | 3 | 2 |
| 11 | 3 | 3 |
| 12 | 3 | 4 |
| 13 | 4 | 1 |
| 14 | 4 | 2 |
| 15 | 4 | 3 |
| 16 | 4 | 4 |
Factor s1 is corresponding with described orthogonal table with the parameter value of s2, and combine factor w that numerical value is fixing,
I.e. can get the EXPERIMENTAL DESIGN table described in table 2.
Table 2:
After obtaining above-mentioned EXPERIMENTAL DESIGN table, according still further to the content of EXPERIMENTAL DESIGN table, make 16 corresponding ginsengs
Then these 16 man-made microstructure are carried out electromagnetism test by " snowflake " the shape man-made microstructure corresponding to number.
From the foregoing, it will be observed that existing orthogonal design needs to make substantial amounts of man-made microstructure, and need to carry out greatly
The Electromagnetic Simulation test of amount, such as, when choose " snowflake " shape man-made microstructure parameter s1 and s2
When value number is respectively 30,30, need 30*30=900 Electromagnetic Simulation test according to orthogonal design,
Need the man-made microstructure of 900 corresponding parametric values combinations.Arrange so many Electromagnetic Simulation test and people
Make micro structure, the most time-consuming, effort, and cost is the highest.
Summary of the invention
Embodiment of the present invention technical problem to be solved is, how the EXPERIMENTAL DESIGN table to man-made microstructure
Obtain, can make in the case of not affecting result of the test, reduce test number (TN).
In order to solve above-mentioned technical problem, embodiments provide the EXPERIMENTAL DESIGN of a kind of man-made microstructure
List acquirement method, including:
Determine the value of the number of the geometric parameter factor of man-made microstructure to be measured, each geometric parameter factor
Several and corresponding parameter value；
Number, the value number of each geometric parameter factor and correspondence thereof according to the described geometric parameter factor
Parameter value, and obtain initial trial design table according to Latin Hypercube Sampling method；
Electromagnetic response delta data when testing according to described initial trial design table, obtains each geometric parameters
Active geometric parameter scope in the number factor；
According to Latin Hypercube Sampling method and the geometric parameter scope of each geometric parameter factor of acquisition, obtain
Test parameters value；
Test parameters value structure EXPERIMENTAL DESIGN table according to uniform Design principle and described acquisition.
Wherein, described according to Latin Hypercube Sampling method with the geometric parameter of each geometric parameter factor of acquisition
Scope, obtains test parameters value, including:
According to the value number of each geometric parameter factor, use Latin Hypercube Sampling method to calculate acquisition and take out
Sample value；
According to calculated sample value and the geometric parameter scope of each geometric parameter factor of described acquisition, structure
Make and obtain test parameters value.
Wherein, the described test parameters value according to uniform Design principle and described acquisition constructs EXPERIMENTAL DESIGN table,
Including:
According to good grid point method structure uniform designs table in uniform Design principle；
Described test parameters value is substituted in described uniform designs table and obtain EXPERIMENTAL DESIGN table.
Wherein, described according in uniform Design principle good grid point method structure uniform designs table, including:
Arranging test number (TN) value n, described test number (TN) value n is the value number of each geometric parameter factor
Least common multiple；
According to described test number (TN) value n, extract less than described test number (TN) value n and with described test number (TN)
Value n becomes all numerical value of relatively prime relation；
According to test number (TN) value n and all numerical value of described extraction, generate uniform Design according to good grid point method
Each row of table, it is thus achieved that uniform designs table.
Wherein, the number of the described geometric parameter factor determining man-made microstructure to be measured, each geometric parameter because of
The value number of son and the parameter value of correspondence thereof include:
Determine the value of the number of the geometric parameter factor of man-made microstructure to be measured, each geometric parameter factor
Number；
Determine the span of each geometric parameter factor to determine each several according to described value number
The parameter value what parameter factors is corresponding.
Correspondingly, the embodiment of the present invention additionally provides the EXPERIMENTAL DESIGN table acquisition device of a kind of man-made microstructure,
Including:
Parameter determination module, for determining the number of the geometric parameter factor of man-made microstructure to be measured, Mei Geji
The value number of what parameter factors and the parameter value of correspondence thereof；
First constructing module, for according to the number of the described geometric parameter factor, each geometric parameter factor
Value number and the parameter value of correspondence thereof, and obtain initial trial design table according to Latin Hypercube Sampling method；
Parameter acquisition module, the electromagnetic response change in time testing according to described initial trial design table
Data, obtain active geometric parameter scope in each geometric parameter factor, and according to Latin hypercube
The geometric parameter scope of each geometric parameter factor of sampling approach and acquisition, obtains test parameters value；
Second constructing module, constructs test for the test parameters value according to uniform Design principle and described acquisition
Design table.
Wherein, described parameter acquisition module includes:
Parameter area acquiring unit, the electromagnetic response in time testing according to described initial trial design table
Delta data, obtains active geometric parameter scope in each geometric parameter factor；
Test parameters acquiring unit, for according to each geometric parameter of Latin Hypercube Sampling method and acquisition because of
The geometric parameter scope of son, obtains test parameters value.
Wherein, described second constructing module includes:
Structural unit, for according to good grid point method structure uniform designs table in uniform Design principle；
Acquiring unit, obtains EXPERIMENTAL DESIGN table for described test parameters value being substituted in described uniform designs table.
Wherein, described structural unit specifically includes:
Subelement is set, is used for arranging test number (TN) value n, and extracts little according to described test number (TN) value n
All numerical value of relatively prime relation are become in described test number (TN) value n and with described test number (TN) value n；
Constructor unit, for according to test number (TN) value n and all numerical value of described extraction, according to good grid
Point method generates each row of uniform designs table, it is thus achieved that uniform designs table.
Wherein, described parameter determination module is additionally operable to determine the span of each geometric parameter factor.
Implement the embodiment of the present invention, have the advantages that
During the EXPERIMENTAL DESIGN table of man-made microstructure is chosen, obtain according to Latin Hypercube Sampling method and treat
Survey parameter less initial trial design table, then after testing according to this initial trial design table, obtain
Take active geometric parameter scope, and carry out Latin hypercube again according to this geometric parameter scope and take out
After sample calculates, obtain final EXPERIMENTAL DESIGN table according to uniform Design principle, carry out according to this EXPERIMENTAL DESIGN table
Test contributes in the case of not affecting experimental result, reduces experimental procedure, it is achieved few in quantity
Experiment in draw distortionless experimental result.
Accompanying drawing explanation
In order to be illustrated more clearly that the embodiment of the present invention or technical scheme of the prior art, below will be to enforcement
In example or description of the prior art, the required accompanying drawing used is briefly described, it should be apparent that, describe below
In accompanying drawing be only some embodiments of the present invention, for those of ordinary skill in the art, do not paying
On the premise of going out creative work, it is also possible to obtain other accompanying drawing according to these accompanying drawings.
Fig. 1 is a kind of " snowflake " shape man-made microstructure to be measured schematic diagram；
Fig. 2 is the first embodiment flow process signal of the EXPERIMENTAL DESIGN list acquirement method of the man-made microstructure of the present invention
Figure；
Fig. 3 is the second embodiment flow process signal of the EXPERIMENTAL DESIGN list acquirement method of the man-made microstructure of the present invention
Figure；
Fig. 4 is the structure composition signal of the EXPERIMENTAL DESIGN table acquisition device of the man-made microstructure of the embodiment of the present invention
Figure.
Detailed description of the invention
Below in conjunction with the accompanying drawing in the embodiment of the present invention, the technical scheme in the embodiment of the present invention is carried out clearly
Chu, be fully described by, it is clear that described embodiment be only a part of embodiment of the present invention rather than
Whole embodiments.Based on the embodiment in the present invention, those of ordinary skill in the art are not making creation
The every other embodiment obtained under property work premise, broadly falls into the scope of protection of the invention.
Refer to Fig. 2, be the first embodiment of the EXPERIMENTAL DESIGN list acquirement method of the man-made microstructure of the present invention
Schematic flow sheet；The described method of the present embodiment includes:
S101: determine the number of the geometric parameter factor of man-made microstructure to be measured, each geometric parameter factor
Value number and the parameter value of correspondence thereof.
The concrete each parameter factors in described S101 can be configured as required by user, such as, right
In " snowflake " the shape man-made microstructure shown in Fig. 1, the geometric parameter factor that user can be arranged is: a1, a2
With b1, b2, for convenience of calculation, can further determine parameter factors s1 and s2, described s1=a1+b1,
Described s2=a2+b2, the acquisition so making test table is the simplest.
User also can be first according to the size of required Meta Materials, it is first determined go out the value model of each parameter factors
Enclose, in the range of this, extract some parameter values so that subsequent construction EXPERIMENTAL DESIGN table.
Described S101 specifically comprises the steps that the number of the geometric parameter factor determining man-made microstructure to be measured, every
The value number of the individual geometric parameter factor；Determine the span of each geometric parameter factor so that according to institute
State the parameter value that value number determines that each geometric parameter factor pair is answered.
S102: according to the number of the described geometric parameter factor, the value number of each geometric parameter factor and
Corresponding parameter value, and obtain initial trial design table according to Latin Hypercube Sampling method.
Concrete, Latin Hypercube Sampling method belongs to a kind of affined sampling method, Latin Hypercube Sampling
First [0,1] interval can be divided into the subinterval of n non-overlapping copies by method, then divides in each subinterval
Do not carry out independent sampling with equal probability, obtain n the random number i.e. sample value more than 0 less than 1.It is concrete
The formula of foundation is as follows:
AndI=0 therein, 1 ..., n.
uiValue be i-th subinterval and sample the sample value that obtains, owing to there is described relational expression so that every
One subinterval is only capable of producing a random number, taking of the geometric parameter factor needed for being i.e. only capable of according to user
Value number n obtains n random number.
N is the value number of aforesaid geometric parameter factor s1 and s2 herein, as s1 and s2 takes 4 values,
The most described n=4.
After obtaining sample value, to each geometric parameter factor, again according to formula: 1+ (n-1) * ui%,
I=1,2 ..., n.The most immediate integer is taken, using this integer as setting according to the principle rounded up
The value of the corresponding list item of the meter every string of table, wherein, the parameter corresponding to each parameter factors of this integer representation
The sequence number of value.Then corresponding parameter value is substituted in described design table and i.e. can get initial trial design table.
S103: electromagnetic response delta data when testing according to described initial trial design table, obtains each
Active geometric parameter scope in the geometric parameter factor.
User can design table according to the initial trial of described structure, makes the most artificial according to the parameter value of each row
The model of micro structure, and these models are carried out respectively electromagnetism test, obtain each artificial microstructure model
The situation of change of electromagnetic response data, described S103 determines parameter area according to changing violent situation, can
Further according to changing violent situation making change curve, then determine parameter area according to curve chart.
Such as, for " snowflake " shape man-made microstructure of Fig. 2, after a test, it may be determined that go out described geometry
Parameter factors s1 is in the range of [2.5333,3.000], and incorporating parametric factor s2 is in the range of [0.5167,0.9833]
Time, change is relatively violent, and changing value is above certain threshold value, and wherein, this threshold value can be according to user to super material
The test requirements document of material is configured.
S104: according to Latin Hypercube Sampling method and the geometric parameter scope of each geometric parameter factor of acquisition,
Obtain test parameters value.
Within the above range, again carry out Latin Hypercube Sampling, determine multiple sample value, and according to sampling
Value and described geometric parameter scope, determine the value number of each geometric parameter factor made new advances and new again
Parameter value, thus structure obtain EXPERIMENTAL DESIGN table.
S105: construct EXPERIMENTAL DESIGN table according to the test parameters value of uniform Design principle and described acquisition.
Uniform Design principle be based on testing site in whole trial stretch uniformly dispersing from uniformity angle
A kind of test design method set out, uses all according to the value number of described new each geometric parameter factor
Even design principle can get a design table, and on this table, the parameters factor is dispersed.Described S105
In design table according to uniform Design principle structure, the numerical value of each list item represents the corresponding geometric parameter factor
The sequence number of parameter value, the new parameter value correspondence of the parameters factor is entered in corresponding list item,
EXPERIMENTAL DESIGN table described in construction complete.
The present invention is during the EXPERIMENTAL DESIGN table of man-made microstructure is chosen, according to Latin Hypercube Sampling method
Obtain parameter to be measured less initial trial design table, then testing according to this initial trial design table
After, obtain active geometric parameter scope, and it is super to carry out Latin again according to this geometric parameter scope
After cube sample calculation, obtain final EXPERIMENTAL DESIGN table according to uniform Design principle, according to this EXPERIMENTAL DESIGN
Table carries out test and contributes in the case of not affecting experimental result, reduces experimental procedure, it is achieved in quantity
Few experiment draws distortionless experimental result.
Refer to Fig. 3 again, be the second enforcement of the EXPERIMENTAL DESIGN list acquirement method of the man-made microstructure of the present invention
Example schematic flow sheet；The described method of the present embodiment specifically includes:
S201: determine the number of the geometric parameter factor of man-made microstructure to be measured, each geometric parameter factor
Value number and the parameter value of correspondence thereof.
S202: according to the number of the described geometric parameter factor, the value number of each geometric parameter factor and
Corresponding parameter value, and obtain initial trial design table according to Latin Hypercube Sampling method.
S203: electromagnetic response delta data when testing according to described initial trial design table, obtains each
Active geometric parameter scope in the geometric parameter factor.
S204: use Latin Hypercube Sampling method to calculate and obtain sample value；According to calculated sample value
With the geometric parameter scope of each geometric parameter factor of described acquisition, obtain test parameters value.
I.e. it is sampled in the range of the geometric parameter determined again according to described S204, according to the examination of user
Test the new geometric parameter number needing to obtain each geometric parameter factor, and each new parameter value is i.e.
Test parameters value.
S205: according to good grid point method structure uniform designs table in uniform Design principle.
Described S205 specifically may include that and arranges test number (TN) value n, and described test number (TN) value n is each
The least common multiple of the value number of the geometric parameter factor；According to described test number (TN) value n, extract and be less than
Described test number (TN) value n and become all numerical value of relatively prime relation with described test number (TN) value n, according to test time
Numerical value n and all numerical value of described extraction, generate each row of uniform designs table according to good grid point method, it is thus achieved that
Uniform designs table.As user needs to carry out 8 tests, then n=8, described S205 can be set and obtain simultaneously
Less than 8 and the numerical value relatively prime with 8 is: 1,3,5,7.Specifically according to the equation below in good grid point method:
u1j=hj；
Wherein, i=1,2 ..., n-1, hjIt is designated as jth and n
Become the numerical value of relatively prime relation；
Every string of tectonic sieving table, obtains design table, then extracts the s row new matrix of composition from this design table,
And each new matrix is carried out the calculating of L2 centralization deviation, obtain the matrix that deviation is minimum, then according to being somebody's turn to do
The numerical value of each list item in matrix construction uniform designs table uniform designs table indicates each parameter factors institute
The sequence number of corresponding new parameter value.Wherein, the value of described s is the geometry that the parameter value of user setup is variable
The number of parameter factors.
S206: described test parameters value is substituted in described uniform designs table and obtain EXPERIMENTAL DESIGN table.
According to the numerical value of list item each in uniform designs table, by new parameter value i.e. test parameters value typing wherein,
I.e. can get EXPERIMENTAL DESIGN table.
Below with the instantiation EXPERIMENTAL DESIGN table choosing method to the man-made microstructure described in the present embodiment
Illustrate.
S1: arrange geometric parameter factor s1 and s2, all takes 4 values, each parameter value be s1=[1.6000,
2.0667,2.5333,3.0000], s2=[0.0500,0.5167,0.9833,1.4500].
S2: use Latin Hypercube Sampling method construct to obtain initial trial design table.
According to equation below:
AndI=1 therein, 2,3, n；The value of n herein is each
The value number of the geometric parameter factor, n=4 the most now.
It is calculated corresponding sample value u by above-mentioned formulai, then according to 1+ (n-1) * ui%, and therefrom
Take immediate integer, obtain the value of each list item of design table, thus it is super vertical to obtain the Latin shown in table 3
Side's EXPERIMENTAL DESIGN table.
Table 3:
| Sequence number | s1 | s2 |
| 1 | 3 | 2 |
| 2 | 4 | 3 |
| 3 | 2 | 1 |
| 4 | 1 | 4 |
Corresponding numerical value is substituted in table 3, obtains initial trial design table as shown in table 4.
Table 4:
| Sequence number | s1 | s2 |
| 1 | 2.5333 | 0.5167 |
| 2 | 3.0000 | 0.9833 |
| 3 | 2.0667 | 0.0500 |
| 4 | 1.6000 | 1.4500 |
S3: arrange man-made microstructure model to test according to above-mentioned table 4, determine man-made microstructure model
The span that middle data variation is the most violent is: s1 is [2.5333,3.000], and s2 is [0.5167,0.9833].
S4: use Latin Hypercube Sampling method to calculate and obtain sample value, and be calculated test parameters value.
Need to carry out 8 tests according to test needs, such as user, can be again according to Latin Hypercube Sampling side
Method obtain respectively s1 8 sample values sampling 1 (i=1,2 ..., 8) and 8 sample values of s2 sample 2
(i=1,2 ..., 8), then according to s1iValue=2.5333+ (3.0000-2.5333) × sampling 1 (i),
s2iValue=0.5167+ (0.9833-0.5167) × sampling 2 (i).
8 the new parameter values i.e. test parameters value that thus can get geometric parameter factor s1 is: s1=
[2.53330000000000,2.59997142857143,2.66664285714286,2.73331428571429,
2.79998571428571,2.86665714285714,2.93332857142857,3]；
8 the new parameter values i.e. test parameters value of available geometric parameter factor s2 is:
S2=[0.516700000000000,0.583357142857143,0.650014285714286,
0.716671428571429,0.783328571428571,0.849985714285714,
0.916642857142857,0.983300000000000].
S5: according to good grid point method structure uniform designs table in uniform Design principle.
Default needs carry out frequency n 2=8 (the most above-mentioned user needs to carry out 8 tests) tested, and obtain
To becoming the numerical value of relatively prime relation to include less than 8 and with 8: 1,3,5,7.
Equation below according to good grid point method:
u1j=hj；
Wherein, i=1,2 ..., n-1, hjIt is designated as jth and n
Become the numerical value of relatively prime relation；
Obtain uniform designs table as shown in table 5.
Table 5:
| 1 | 3 | 5 | 7 |
| 2 | 6 | 2 | 6 |
| 3 | 1 | 7 | 5 |
| 4 | 4 | 4 | 4 |
| 5 | 7 | 1 | 3 |
| 6 | 2 | 6 | 2 |
| 7 | 5 | 3 | 1 |
| 8 | 8 | 8 | 8 |
Obtaining s=2 therein row from table 5 and constitute multiple matrixes, s therein is variable parameter factors
Number, this example includes s1 and s2 two, therefore, s=2, and it is inclined that each matrix is carried out centralization L2
Difference calculates, it is thus achieved that the matrix of deviation minimum is as uniform designs table.Concrete, first row and the 3rd row are constituted
New matrix deviation value be: 0.0713 is minimum, therefore, using the first row in table 4 and the 3rd row as equal
Even design table, the most following table 6.
Table 6:
| Sequence number | Parameter 1 | Parameter 2 |
| 1 | 1 | 5 |
| 2 | 2 | 2 |
| 3 | 3 | 7 |
| 4 | 4 | 4 |
| 5 | 5 | 1 |
| 6 | 6 | 6 |
| 7 | 7 | 3 |
| 8 | 8 | 8 |
The test parameters value of corresponding s1 and the test parameters value of s2 are substituted in table 6, i.e. available such as table
EXPERIMENTAL DESIGN table shown in 7.
Table 7:
| Sequence number | s1 | s2 |
| 1 | 2.53330000000000 | 0.783328571428571 |
| 2 | 2.59997142857143 | 0.583357142857143 |
| 3 | 2.66664285714286 | 0.916642857142857 |
| 4 | 2.73331428571429 | 0.716671428571429 |
| 5 | 2.79998571428571 | 0.516700000000000 |
| 6 | 2.86665714285714 | 0.849985714285714 |
| 7 | 2.93332857142857 | 0.650014285714286 |
| 8 | 3.00000000000000 | 0.983300000000000 |
After obtaining table 7, user can design corresponding artificial micro-knot according to the parameter size provided in table 7
Structure model, carries out required test.
By above-described embodiment, according to orthogonal design mode, carrying out this test needs to carry out 4*4
I.e. 16 times, and use the present invention only to need according to user's needs, arrange to test for 8 times, greatly reduce
Test number (TN), and, the present invention is during the EXPERIMENTAL DESIGN table of man-made microstructure is chosen, based on uniformly
Design principle and Latin Hypercube Sampling method obtain EXPERIMENTAL DESIGN table, also do not interfere with experimental result.
Refer to Fig. 4, be the structure of the EXPERIMENTAL DESIGN table acquisition device of the man-made microstructure of the embodiment of the present invention
Composition schematic diagram.The described device of the present embodiment includes:
Parameter determination module 1, for determining the number, each of the geometric parameter factor of man-made microstructure to be measured
The value number of the geometric parameter factor and the parameter value of correspondence thereof；
Each parameter factors that concrete described parameter determination module 1 determines can be configured as required by user,
Such as, for " snowflake " the shape man-made microstructure shown in Fig. 1, the geometric parameter factor that user can be arranged is:
A1, a2 and b1, b2, for convenience of calculation, can further determine parameter factors s1 and s2, described
S1=a1+b1, described s2=a2+b2, the acquisition so making test table is the simplest.
User also can be first according to the size of required Meta Materials, it is first determined go out the value model of each parameter factors
Enclose, in the range of this, extract some parameter values so that subsequent construction EXPERIMENTAL DESIGN table.
First constructing module 2, for according to the number of the described geometric parameter factor, each geometric parameter factor
Value number and the parameter value of correspondence, and according to Latin Hypercube Sampling method obtain initial trial design
Table；
Concrete, Latin Hypercube Sampling method belongs to a kind of affined sampling method, Latin Hypercube Sampling
First [0,1] interval can be divided into the subinterval of n non-overlapping copies by method, then divides in each subinterval
Do not carry out independent sampling with equal probability, obtain n the random number i.e. sample value more than 0 less than 1.It is concrete
The formula of foundation is as follows:
AndI=0 therein, 1 ..., n.
uiValue be i-th subinterval and sample the sample value that obtains.N is the aforesaid geometric parameter factor herein
The value number of s1 and s2, as s1 and s2 takes 4 values, the most described n=4.
After obtaining sample value, to each geometric parameter factor, again according to formula: 1+ (n-1) * ui%,
I=1,2 ..., n.The most immediate integer is taken according to the principle rounded up, this integer is every as design table
The value of the corresponding list item of string, wherein, the sequence of the parameter value corresponding to each parameter factors of this integer representation
Number.Described first constructing module 2, will be corresponding after designing table well according to Latin Hypercube Sampling method construct
Parameter value substitute in described design table and i.e. can get initial trial design table.
Parameter acquisition module 3, the electromagnetic response in time testing according to described initial trial design table becomes
Change data, obtain active geometric parameter scope in each geometric parameter factor, and super vertical according to Latin
The geometric parameter scope of each geometric parameter factor of side's sampling approach and acquisition, obtains test parameters value；
User can design table according to the initial trial of described structure, makes the most artificial according to the parameter value of each row
The model of micro structure, and these models are carried out respectively electromagnetism test, obtain each artificial microstructure model
The situation of change of electromagnetic response data, described parameter acquisition module 3 determines parameter according to changing violent situation
Scope, then can determine according to curve chart further according to changing violent situation making change curve
Parameter area.Such as, for " snowflake " shape man-made microstructure of Fig. 2, after a test, it may be determined that
Go out described geometric parameter factor s1 in the range of [2.5333,3.000], incorporating parametric factor s2 in the range of
Time [0.5167,0.9833], change is relatively violent, and changing value is above certain threshold value.
After determining parameter area, described parameter acquisition module 3 carries out Latin Hypercube Sampling again, determines
Multiple sample values, and according to sample value and described geometric parameter scope, again determine each geometric parameters made new advances
The value number of the number factor and new parameter value thereof, thus structure obtains EXPERIMENTAL DESIGN table.
Second constructing module 4, constructs examination for the test parameters value according to uniform Design principle and described acquisition
Test design table.
Uniform Design principle be based on testing site in whole trial stretch uniformly dispersing from uniformity angle
A kind of test design method set out, uses all according to the value number of described new each geometric parameter factor
Even design principle can get a design table, and on this design table, the parameters factor is dispersed, and described
In this design table that two constructing modules 4 construct according to uniform Design principle, the numerical value of each list item represents phase
Answer the sequence number of the parameter value of the geometric parameter factor, the new parameter value correspondence of the parameters factor is entered into phase
Answer in list item, get final product EXPERIMENTAL DESIGN table described in construction complete.
Concrete, as shown in Figure 4, described parameter acquisition module 3 includes:
Parameter area acquiring unit 31, the electromagnetism in time testing according to described initial trial design table rings
Answer delta data, obtain active geometric parameter scope in each geometric parameter factor；
Test parameters acquiring unit 32, for according to Latin Hypercube Sampling method and each geometric parameter of acquisition
The geometric parameter scope of the factor, obtains test parameters value.
Concrete, as shown in Figure 4, described second constructing module 4 includes:
Structural unit 41, for according to good grid point method structure uniform designs table in uniform Design principle；
Acquiring unit 42, obtains EXPERIMENTAL DESIGN for described test parameters value being substituted in described uniform designs table
Table.
Further, the structural unit 41 in described second constructing module 4 may include that
Subelement is set, is used for arranging test number (TN) value n, and extracts little according to described test number (TN) value n
All numerical value of relatively prime relation are become in described test number (TN) value n and with described test number (TN) value n；
Constructor unit, for according to test number (TN) value n and all numerical value of described extraction, according to good grid
Point method generates each row of uniform designs table, it is thus achieved that uniform designs table.
Concrete, user can preset, by the described subelement that arranges, frequency n 2=8 needing to carry out testing, and
Yield less than 8 and become with 8 the numerical value of relatively prime relation to be: 1,3,5,7.
Equation below according to good grid point method:
u1j=hj；
Wherein, i=1,2 ..., n-1, hjIt is designated as jth and n
Become the value of relatively prime relation.
Every string of tectonic sieving table, obtains design table, then extracts the s row new matrix of composition from this design table,
And each new matrix is carried out the calculating of L2 centralization deviation, obtain the matrix that deviation is minimum, then according to being somebody's turn to do
The numerical value of each list item in matrix construction uniform designs table uniform designs table indicates each parameter factors institute
The sequence number of corresponding new parameter value.Wherein, the value of described s is the geometric parameters of changeable parameters of user setup
The number of the number factor.Finally described test parameters value is substituted into and described uniform designs table obtains EXPERIMENTAL DESIGN table.
The present invention is during the EXPERIMENTAL DESIGN table of man-made microstructure is chosen, according to Latin Hypercube Sampling method
Obtain parameter to be measured less initial trial design table, then testing according to this initial trial design table
After, obtain active geometric parameter scope, and it is super to carry out Latin again according to this geometric parameter scope
After cube sample calculation, obtain final EXPERIMENTAL DESIGN table according to uniform Design principle, according to this EXPERIMENTAL DESIGN
Table carries out test and contributes in the case of not affecting experimental result, reduces experimental procedure, it is achieved in quantity
Few experiment draws distortionless experimental result.
One of ordinary skill in the art will appreciate that all or part of flow process realizing in above-described embodiment method,
Can be by computer program and complete to instruct relevant hardware, described program can be stored in a calculating
In machine read/write memory medium, this program is upon execution, it may include such as the flow process of the embodiment of above-mentioned each method.
Wherein, described storage medium can be magnetic disc, CD, read-only store-memory body (Read-Only Memory,
Or random store-memory body (Random Access Memory, RAM) etc. ROM).
Above disclosed be only present pre-ferred embodiments, certainly can not with this limit the present invention it
Interest field, the equivalent variations therefore made according to the claims in the present invention, still belong to the scope that the present invention is contained.
Claims (10)
1. the EXPERIMENTAL DESIGN list acquirement method of a man-made microstructure, it is characterised in that including:
Determine the number of the geometric parameter factor of man-made microstructure to be measured, the value number of each geometric parameter factor
And the parameter value of correspondence；
Number, the value number of each geometric parameter factor and the ginseng of correspondence thereof according to the described geometric parameter factor
Numerical value, and obtain initial trial design table according to Latin Hypercube Sampling method；
Electromagnetic response delta data when testing according to described initial trial design table, obtains each geometric parameter
Active geometric parameter scope in the factor；
Active geometric parameters in each geometric parameter factor according to Latin Hypercube Sampling method and acquisition
Number scope, obtains test parameters value；
Test parameters value structure EXPERIMENTAL DESIGN table according to uniform Design principle and described acquisition；
The formula of the described concrete foundation of Latin Hypercube Sampling method is as follows:
AndI=0 therein, 1 ..., n；
uiThe sample value obtained for the sampling of i-th subinterval；
N is the value number of the described geometric parameter factor；
U is the random number more than 0 less than 1.
2. the method described in claim 1, it is characterised in that described according to Latin Hypercube Sampling method and
Active geometric parameter scope in each geometric parameter factor obtained, obtains test parameters value, including:
According to the value number of each geometric parameter factor, use Latin Hypercube Sampling method to calculate and obtain sampling
Value；
According to calculated sample value and each active geometric parameters of the geometric parameter factor of described acquisition
Number scope, structure obtains test parameters value.
3. method as claimed in claim 2, it is characterised in that described according to uniform Design principle with described
The test parameters value structure EXPERIMENTAL DESIGN table obtained, including:
According to good grid point method structure uniform designs table in uniform Design principle；
Described test parameters value is substituted in described uniform designs table and obtain EXPERIMENTAL DESIGN table.
4. method as claimed in claim 3, it is characterised in that described according to lattice good in uniform Design principle
Sub-some method structure uniform designs table, including:
Test number (TN) value n is set, described test number (TN) value n be the value number of each geometric parameter factor
Little common multiple；
According to described test number (TN) value n, extract less than described test number (TN) value n and with described test number (TN) value
N becomes all numerical value of relatively prime relation；
According to test number (TN) value n and all numerical value of described extraction, generate uniform designs table according to good grid point method
Each row, it is thus achieved that uniform designs table.
5. method as claimed in claim 4, it is characterised in that described determine the several of man-made microstructure to be measured
What number of parameter factors, the value number of each geometric parameter factor and the parameter value of correspondence thereof, including:
Determine the number of the geometric parameter factor of man-made microstructure to be measured, the value number of each geometric parameter factor；
Determine the span of each geometric parameter factor to determine each geometry according to described value number
The parameter value that parameter factors is corresponding.
6. the EXPERIMENTAL DESIGN table acquisition device of a man-made microstructure, it is characterised in that including:
Parameter determination module, for determining the number of the geometric parameter factor of man-made microstructure to be measured, each geometry
The value number of parameter factors and the parameter value of correspondence thereof；
First constructing module, for according to the number of the described geometric parameter factor, the taking of each geometric parameter factor
Value number and the parameter value of correspondence thereof, and obtain initial trial design table according to Latin Hypercube Sampling method；
Parameter acquisition module, the electromagnetic response change number in time testing according to described initial trial design table
According to, obtain active geometric parameter scope in each geometric parameter factor, and take out according to Latin hypercube
Active geometric parameter scope in each geometric parameter factor of quadrat method and acquisition, obtains test parameters
Value；
Second constructing module, constructs test for the test parameters value according to uniform Design principle and described acquisition and sets
Meter table；
The formula of the described concrete foundation of Latin Hypercube Sampling method is as follows:
AndI=0 therein, 1 ..., n；
uiThe sample value obtained for the sampling of i-th subinterval；
N is the value number of the described geometric parameter factor；
U is the random number more than 0 less than 1.
7. device as claimed in claim 6, it is characterised in that described parameter acquisition module includes:
Parameter area acquiring unit, the electromagnetic response in time testing according to described initial trial design table becomes
Change data, obtain active geometric parameter scope in each geometric parameter factor；
Test parameters acquiring unit, for according to Latin Hypercube Sampling method and each geometric parameter factor of acquisition
Active geometric parameter scope, obtains test parameters value.
8. device as claimed in claim 7, it is characterised in that described second constructing module includes:
Structural unit, for according to good grid point method structure uniform designs table in uniform Design principle；
Acquiring unit, obtains EXPERIMENTAL DESIGN table for described test parameters value being substituted in described uniform designs table.
9. device as claimed in claim 8, it is characterised in that described structural unit specifically includes:
Subelement is set, is used for arranging test number (TN) value n, and extract according to described test number (TN) value n and be less than
Described test number (TN) value n and become all numerical value of relatively prime relation with described test number (TN) value n；
Constructor unit, for according to test number (TN) value n and all numerical value of described extraction, according to good grid point
Method generates each row of uniform designs table, it is thus achieved that uniform designs table.
10. device as claimed in claim 9, it is characterised in that described parameter determination module is additionally operable to really
The span of each geometric parameter factor fixed.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201110396037.7A CN103136401B (en) | 2011-12-02 | 2011-12-02 | The EXPERIMENTAL DESIGN list acquirement method of a kind of man-made microstructure and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201110396037.7A CN103136401B (en) | 2011-12-02 | 2011-12-02 | The EXPERIMENTAL DESIGN list acquirement method of a kind of man-made microstructure and device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN103136401A CN103136401A (en) | 2013-06-05 |
| CN103136401B true CN103136401B (en) | 2016-08-17 |
Family
ID=48496224
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201110396037.7A Active CN103136401B (en) | 2011-12-02 | 2011-12-02 | The EXPERIMENTAL DESIGN list acquirement method of a kind of man-made microstructure and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103136401B (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106709082A (en) * | 2015-11-12 | 2017-05-24 | 北京仿真中心 | Universal interface-based simulation test design system and method |
| CN112507468B (en) * | 2021-02-01 | 2021-04-23 | 中国人民解放军国防科技大学 | Parameter optimization method and device for suspension system of magnetic-levitation train and electronic equipment |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101923590A (en) * | 2010-08-16 | 2010-12-22 | 北京理工大学 | High-efficiency Latin hypercube experimental design method |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE102006024097A1 (en) * | 2006-05-18 | 2007-11-22 | E.G.O. Elektro-Gerätebau GmbH | Use of left-handed metamaterials as a display, in particular on a cooktop, and display and display method |
-
2011
- 2011-12-02 CN CN201110396037.7A patent/CN103136401B/en active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101923590A (en) * | 2010-08-16 | 2010-12-22 | 北京理工大学 | High-efficiency Latin hypercube experimental design method |
Non-Patent Citations (2)
| Title |
|---|
| 同时实现介电常数和磁导率为负的H型结构单元左手材料;刘亚红 等;《物理学报》;20071031;第56卷(第10期);第5884页第2.1和3.1节和图1(a)、3 * |
| 均匀设计表的MATLAB实现;曹慧荣 等;《统计与决策》;20080331(第6期);正文第1页左栏第3行至第3页右栏第19行 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN103136401A (en) | 2013-06-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Peng et al. | Evaluating the effectiveness of landscape metrics in quantifying spatial patterns | |
| Xiao et al. | Emergence of symmetry in complex networks | |
| CN106338440B (en) | The true and false three axis material properties of rock test device of one kind and its method | |
| CN103136401B (en) | The EXPERIMENTAL DESIGN list acquirement method of a kind of man-made microstructure and device | |
| CN106610482A (en) | Method for analyzing impact on accuracy of capacitive-type voltage transformer from stray capacitance | |
| CN102968813A (en) | Surface sampling method of triangular patch mesh model | |
| CN105302992A (en) | Quick optimization and simulation design method for stirrer in reverberation chamber | |
| CN116449106A (en) | Impedance detection method, device, equipment and storage medium of network-structured converter | |
| CN103136402B (en) | The EXPERIMENTAL DESIGN list acquirement method of a kind of man-made microstructure and device | |
| CN102208027A (en) | Method for evaluating land utilization spatial pattern based on clearance degree dimension | |
| CN106951918A (en) | A kind of individual particle image clustering method analyzed for Ice mapping | |
| CN107219483A (en) | A kind of radial direction kurtosis anisotropic quantitative approach based on diffusion kurtosis imaging | |
| Rahman et al. | Applications of random finite element method in bearing capacity problems | |
| CN102914705B (en) | Method and method for selecting test points of artificial electromagnetic material unit | |
| Williams et al. | First investigation of void statistics in numerical relativity simulations | |
| CN107526856B (en) | Parallel explicit-implicit mixed discontinuous Galerkin time domain finite element method | |
| CN103136403B (en) | A kind of test design list acquirement method of man-made microstructure and device | |
| CN103135004A (en) | Method of constructing test table and method and device of measuring electromagnetic property | |
| Zhang et al. | Evaluation of role of Home Market Effects in China’s manufacturing industries | |
| CN102890202B (en) | Method and device for selecting test points of artificial electromagnetic material unit | |
| CN102768312A (en) | Test point selection method and device for artificial electromagnetic material unit | |
| Lefik et al. | Artificial neural network as a numerical form of effective constitutive law for composites with parametrized and hierarchical microstructure | |
| CN103091570B (en) | A kind of structural unit parameters simulation experimental technique of artificial electromagnetic material and device | |
| CN102890203B (en) | Method and device for selecting test point of artificial electromagnetic material unit | |
| CN103164557A (en) | Metamaterial design device and design method thereof |
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 |