[go: up one dir, main page]
More Web Proxy on the site http://driver.im/

KR102198714B1 - System and method for estimating property of noise control material using genetic algorithm - Google Patents

System and method for estimating property of noise control material using genetic algorithm Download PDF

Info

Publication number
KR102198714B1
KR102198714B1 KR1020190033519A KR20190033519A KR102198714B1 KR 102198714 B1 KR102198714 B1 KR 102198714B1 KR 1020190033519 A KR1020190033519 A KR 1020190033519A KR 20190033519 A KR20190033519 A KR 20190033519A KR 102198714 B1 KR102198714 B1 KR 102198714B1
Authority
KR
South Korea
Prior art keywords
solution
fitness
solutions
parent
fit
Prior art date
Application number
KR1020190033519A
Other languages
Korean (ko)
Other versions
KR20200113392A (en
Inventor
박영호
이인석
Original Assignee
창원대학교 산학협력단
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 창원대학교 산학협력단 filed Critical 창원대학교 산학협력단
Priority to KR1020190033519A priority Critical patent/KR102198714B1/en
Publication of KR20200113392A publication Critical patent/KR20200113392A/en
Application granted granted Critical
Publication of KR102198714B1 publication Critical patent/KR102198714B1/en

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Data Mining & Analysis (AREA)
  • General Engineering & Computer Science (AREA)
  • Software Systems (AREA)
  • Computational Mathematics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Databases & Information Systems (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Algebra (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Abstract

본 발명은 유전알고리즘을 이용한 소음제어물질의 물성치(예컨대, 유동저항, 공극률, 뒤틀림도, 열특성길이, 점성특성길이 등) 추정 시스템 및 방법에 관한 것으로, 특히 소음제어물질의 물성치에 대해서 해집단을 구성하고, 해집단의 적합도를 계산하며, 적합도 계산 결과와 해집단을 병합하여 새로운 배열을 구성하여 부모해를 선택하며, 부모해를 2진수로 부호화하여 교배하여서 자식해를 생성하며, 자식해의 다양성을 증가시키기 위한 돌연변이를 발생시키며, 부모해와 자식해를 10진수로 복호화하여 가정을 구성하고 그 가정에 대한 적합도를 계산하고, 해집단과 가정을 비교하여 해를 대치해주며, 가장 높은 적합도를 갖는 해를 추정치로 결정하는, 유전알고리즘을 이용한 소음제어물질의 물성치 추정 시스템 및 방법에 관한 것이다.The present invention relates to a system and method for estimating the physical properties of a noise control material (e.g., flow resistance, porosity, warpage, thermal length, viscous length, etc.) using a dielectric algorithm. Is constructed, the fitness of the solution group is calculated, the result of the fitness calculation and the solution group are merged to form a new array to select the parent solution, and the parent solution is coded in binary and crossed to create a child solution. It generates a mutation to increase the diversity of the parent and child solutions, and constructs an assumption by decoding the parent and child solutions into decimal numbers, calculates the fit for the assumption, compares the solution group and the assumption, and replaces the solution. The present invention relates to a system and method for estimating the property values of noise control materials using a genetic algorithm that determines a solution with goodness of fit as an estimate.

Description

유전알고리즘을 이용한 소음제어물질의 물성치 추정 시스템 및 방법{SYSTEM AND METHOD FOR ESTIMATING PROPERTY OF NOISE CONTROL MATERIAL USING GENETIC ALGORITHM}A system and method for estimating the properties of noise control materials using genetic algorithms {SYSTEM AND METHOD FOR ESTIMATING PROPERTY OF NOISE CONTROL MATERIAL USING GENETIC ALGORITHM}

본 발명은 유전알고리즘을 이용한 소음제어물질의 물성치(예컨대, 유동저항, 공극률, 뒤틀림도, 열특성길이, 점성특성길이 등) 추정 시스템 및 방법에 관한 것으로, 특히 소음제어물질의 물성치에 대해서 해집단을 구성하고, 해집단의 적합도를 계산하며, 적합도 계산 결과와 해집단을 병합하여 새로운 배열을 구성하여 부모해를 선택하며, 부모해를 2진수로 부호화하여 교배하여서 자식해를 생성하며, 자식해의 다양성을 증가시키기 위한 돌연변이를 발생시키며, 부모해와 자식해를 10진수로 복호화하여 가정을 구성하고 그 가정에 대한 적합도를 계산하고, 해집단과 가정을 비교하여 해를 대치해주며, 가장 높은 적합도를 갖는 해를 추정치로 결정하는, 유전알고리즘을 이용한 소음제어물질의 물성치 추정 시스템 및 방법에 관한 것이다.The present invention relates to a system and method for estimating the physical properties of a noise control material (e.g., flow resistance, porosity, warpage, thermal length, viscous length, etc.) using a dielectric algorithm. Is constructed, the fitness of the solution group is calculated, the result of the fitness calculation and the solution group are merged to form a new array to select the parent solution, and the parent solution is coded in binary and crossed to create a child solution. It generates a mutation to increase the diversity of the parent and child solutions, and constructs an assumption by decoding the parent and child solutions into decimal numbers, calculates the fit for the assumption, compares the solution group and the assumption, and replaces the solution. The present invention relates to a system and method for estimating the property values of noise control materials using a genetic algorithm that determines a solution with goodness of fit as an estimate.

일반적으로, 공학적인 정적인 방음 기재로서 크게 흡음재와 차음재로 나뉜다. 능동적인 소음 제거를 꾀하는 액티브 노이즈 캔슬링과는 차이가 있으며 방음막 설치를 통해서 소리를 차단하여 방음 효과를 꾀하는 것이다. 흡음재는 소음과 진동을 흡수하는 효과를 가지며, 차음재의 경우 소음을 차단하여 투과되는 소음을 감소시키는 특징이 있다. 시공하는 장소에 따라 다르지만 보통 차음재와 흡음재를 함께 시공해야만 최대의 방음 효과를 얻을 수 있다. In general, as an engineering static soundproofing substrate, it is largely divided into a sound absorbing material and a sound insulating material. It is different from active noise canceling, which aims to remove active noise, and it is intended to block sound through the installation of a sound insulation film to achieve a sound insulation effect. The sound-absorbing material has an effect of absorbing noise and vibration, and in the case of the sound-absorbing material, it has a characteristic of reducing transmitted noise by blocking noise. Although it depends on the place of construction, the maximum sound insulation effect can be obtained only when a sound insulation material and a sound absorption material are usually installed together.

한편, 이와 같은 흡음재 및 차음재를 적재적소에 사용하기 위해서는, 그것에 대한 물성을 정확히 파악할 필요가 있으며, 더욱이 이와 같은 소음제어물질의 물성에 대해서 수치로 추정할 필요성이 있었다.On the other hand, in order to use the sound-absorbing material and the sound-insulating material in the right place, it is necessary to accurately grasp the physical properties thereof, and furthermore, it is necessary to estimate the physical properties of the noise control material by numerical values.

국제 특허 공개 WO 2009/119344호 공보(공개일:2009.10.01)International Patent Publication WO 2009/119344 (published date: 2009.10.01)

따라서 본 발명은 상기와 같은 종래의 문제점을 해결하기 위해 이루어진 것으로서, 본 발명의 목적은 소음제어물질의 물성치를 수치로 정확히 추정할 수 있는, 유전알고리즘을 이용한 소음제어물질의 물성치 추정 시스템 및 방법을 제공하는 데에 있다.Accordingly, the present invention has been made to solve the conventional problems as described above, and an object of the present invention is to provide a system and method for estimating the properties of a noise control material using a genetic algorithm that can accurately estimate the properties of the noise control material in numerical values. To provide.

상기의 목적을 달성하기 위해, 본 발명의 실시형태에 의한 유전알고리즘을 이용한 소음제어물질의 물성치 추정 시스템은 소음제어물질의 소정개수(n개)의 물성치를 그 물성치가 존재 가능한 영역에서 중복을 허용하지 않고 무작위로 설정 개수(m개)만큼 선택한 뒤 열이 물성치가 되도록 하여 m행 n열의 배열 형태로 해 집단을 구성하는 해 집단 구성부; 상기 해 집단의 모든 해에 대한 적합도를 계산하고 m행 1열의 배열 형태로서 결과를 저장하는 해 집단의 적합도 계산부; 상기 적합도의 계산 결과와 상기 해 집단을 병합하여 m행 n+1 열의 새로운 배열을 구성하고 1열(적합도 열)을 기준으로 오름차순으로 정렬하고, 모든 해 집단에서 2개의 해(부모해)를 중복을 허용하여 선택하는 부모해 선택부; 교배를 위해 선택된 상기 부모해를 2진수로 부호화하는 2진 부호화부; 2진수로 부호화된 상기 부모해를 교배하여 자식해를 생성하는 자식해 생성부; 상기 자식해의 다양성을 증가시키기 위한 돌연변이를 발생시키는 돌연변이 발생부; 상기 부모해 및 자식해를 10진수로 복호화하고 하나의 가정을 구성하는 가정 구성부; 상기 가정의 적합도를 상기 해 집단의 적합도 계산부에서의 적합도 계산 방식과 동일한 방법으로 계산하는 가정의 적합도 계산부; 상기 해 집단에서 가장 높은 적합도(열성)를 갖는 해(해 1)와 상기 가정에서 가장 낮은 적합도(우성)를 갖는 해(해 2)를 비교하여, 열성보다 우성이 낮은 적합도를 가질 경우 상기 해 집단에서 해 1을 해 2로 대치하는 해 대치부; 및 가장 높은 적합도를 갖는 해를 추정치로 결정하는 추정치 결정부;를 포함하는 것을 특징으로 한다.In order to achieve the above object, the system for estimating the physical properties of the noise control material using the genetic algorithm according to the embodiment of the present invention allows the overlapping of the property values of a predetermined number (n) of the noise control material in an area where the property can exist. A solution group constructing unit configured to form a solution group in an array of m rows and n columns by randomly selecting as many as a set number (m) and making the columns become material properties; A goodness-of-fit calculation unit for calculating the fitness for all the solutions of the solution group and storing the result as an array of m rows and 1 column; By merging the calculation result of the goodness-of-fit and the solution group, a new array of m rows and n+1 columns is formed and sorted in ascending order based on column 1 (goodness-of-fit column), and two solutions (parent solutions) are duplicated in all solution groups. Parental selection unit to select by allowing; A binary encoder for encoding the parent solution selected for mating into binary numbers; A child solution generator for generating a child solution by crossing the parent solution encoded in binary numbers; A mutation generating unit that generates a mutation for increasing the diversity of the offspring; A home construction unit that decodes the parent and child solutions into decimal numbers and configures one assumption; A fitness of the assumption that calculates the fitness of the assumption by the same method as the fitness of the solution group; Compare the solution with the highest goodness of fit (recessive) in the solution group (Solution 1) and the solution with the lowest goodness of fit (dominant) in the assumption (Solution 2). The year replacement unit replacing year 1 with year 2 in; And an estimate determination unit that determines a solution having the highest fitness as an estimate.

상기 실시형태에 의한 유전알고리즘을 이용한 소음제어물질의 물성치 추정 시스템에 있어서, 상기 물성치는 유동저항, 공극률, 열 특성 길이(thermal characteristic length) 및 점성 특성 길이를 포함할 수 있다.In the system for estimating the properties of the noise control material using the dielectric algorithm according to the above embodiment, the property values may include flow resistance, porosity, thermal characteristic length, and viscous characteristic length.

상기 실시형태에 의한 유전알고리즘을 이용한 소음제어물질의 물성치 추정 시스템에 있어서, 상기 해 집단의 적합도 계산부 및 가정의 적합도 계산부는 다음의 [수학식 1]에 의해 적합도를 계산할 수 있다.In the system for estimating the property values of the noise control material using the genetic algorithm according to the above embodiment, the fitness of the solution group and the fitness of assumption calculation unit may calculate the fitness by the following [Equation 1].

[수학식 1][Equation 1]

Figure 112019030304587-pat00001
Figure 112019030304587-pat00001

[여기서,

Figure 112019030304587-pat00002
는 적합도를 나타내고,
Figure 112019030304587-pat00003
는 I번째 해를 나타내고,
Figure 112019030304587-pat00004
는 x축 최상단 주파수를 나타내고,
Figure 112019030304587-pat00005
는 상기 x축 최 하단 주파수를 나타내고,
Figure 112019030304587-pat00006
비오(Biot) 이론을 활용하여 추정한 주파수에서의 흡음계수를 나타내며,
Figure 112019030304587-pat00007
는 시험을 통해 계측한 주파수에서의 흡음계수를 나타낸다][here,
Figure 112019030304587-pat00002
Represents the fitness,
Figure 112019030304587-pat00003
Represents the I-th solution,
Figure 112019030304587-pat00004
Represents the uppermost frequency of the x-axis,
Figure 112019030304587-pat00005
Represents the lowermost frequency of the x-axis,
Figure 112019030304587-pat00006
It represents the sound absorption coefficient at the frequency estimated using Biot theory,
Figure 112019030304587-pat00007
Represents the sound absorption coefficient at the frequency measured through the test]

상기 실시형태에 의한 유전알고리즘을 이용한 소음제어물질의 물성치 추정 시스템에 있어서, 상기 부모해 선택부는 다음의 [수학식 2]를 이용하여 모든 해 집단에서 2개의 해(부모해)를 중복을 허용하여 선택할 수 있다.In the system for estimating the properties of the noise control material using the genetic algorithm according to the above embodiment, the parental solution selection unit allows overlapping of two solutions (parent solutions) in all solution groups using the following [Equation 2] You can choose.

[수학식 2][Equation 2]

Figure 112019030304587-pat00008
Figure 112019030304587-pat00008

[여기서,

Figure 112019030304587-pat00009
Figure 112019030304587-pat00010
순위의 선택확률을 나타내고, sp는 선택압을 나타내며, n은 총 해의 개수를 나타내며, i는 순위를 나타낸다][here,
Figure 112019030304587-pat00009
Is
Figure 112019030304587-pat00010
It represents the probability of selection of the ranking, sp represents the selection pressure, n represents the total number of solutions, and i represents the rank]

상기의 목적을 달성하기 위해, 본 발명의 다른 실시형태에 의한 유전알고리즘을 이용한 소음제어물질의 물성치 추정 방법은, 해 집단 구성부가 소음제어물질의 소정개수(n개)의 물성치를 그 물성치가 존재 가능한 영역에서 중복을 허용하지 않고 무작위로 설정 개수(m개)만큼 선택한 뒤 열이 물성치가 되도록 하여 m행 n열의 배열 형태로 해 집단을 구성하는 제 1 단계; 해 집단의 적합도 계산부가 상기 해 집단의 모든 해에 대한 적합도를 계산하고 m행 1열의 배열 형태로서 결과를 저장하는 제 2 단계; 부모해 선택부가 상기 적합도의 계산 결과와 상기 해 집단을 병합하여 m행 n+1 열의 새로운 배열을 구성하고 1열(적합도 열)을 기준으로 오름차순으로 정렬하고, 모든 해 집단에서 2개의 해(부모해)를 중복을 허용하여 선택하는 제 3 단계; 2진 부호화부가 교배를 위해 선택된 상기 부모해를 2진수로 부호화하는 제 4 단계; 자식해 생성부가 2진수로 부호화된 상기 부모해를 교배하여 자식해를 생성하는 제 5 단계; 돌연변이 발생부가 상기 자식해의 다양성을 증가시키기 위한 돌연변이를 발생시키는 제 6 단계; 가정 구성부가 상기 부모해 및 자식해를 10진수로 복호화하고 하나의 가정을 구성하는 제 7 단계; 가정의 적합도 계산부가 상기 가정의 적합도를 상기 제 2 단계에서의 적합도 계산 방식과 동일한 방법으로 계산하는 제 8 단계; 해 대치부가 상기 해 집단에서 가장 높은 적합도(열성)를 갖는 해(해 1)와 상기 가정에서 가장 낮은 적합도(우성)를 갖는 해(해 2)를 비교하여, 열성보다 우성이 낮은 적합도를 가질 경우 상기 해 집단에서 해 1을 해 2로 대치하는 제 9 단계; 추정치 결정부가 상기 제 2 단계 내지 제 9 단계가 설정 횟수 반복되었는지의 여부를 결정하는 제 10 단계; 및 상기 제 10 단계에서 상기 제 2 단계 내지 제 9 단계가 설정 횟수 반복되었다면, 상기 추정치 결정부에 의해 가장 높은 적합도를 갖는 해가 추정치로 결정되는 제 11 단계;를 포함하는 것을 특징으로 한다.In order to achieve the above object, in the method for estimating the physical properties of the noise control material using a genetic algorithm according to another embodiment of the present invention, the solution group component has a predetermined number (n) of the material properties of the noise control material. A first step of forming a group in an array of m rows and n columns by randomly selecting as many as a set number (m) of possible areas without allowing overlapping, and making the columns become physical properties; A second step of calculating, by a fitness calculation unit of a solution group, the fitness for all solutions of the solution group and storing the result as an array of m rows and 1 column; The parental solution selection unit merges the calculation result of the goodness-of-fit and the solution group to form a new array of m rows and n+1 columns, and sorts them in ascending order based on column 1 (goodness-of-fit column), and in all solution groups, two solutions (parent A third step of selecting a solution by allowing redundancy; A fourth step of encoding, by a binary encoding unit, the parent solution selected for mating into binary numbers; A fifth step of generating a child solution by crossing the parent solution encoded as a binary number by a child solution generator; A sixth step of generating a mutation for increasing the diversity of the offspring by the mutation generating unit; A seventh step of decoding the parent solution and the child solution to a decimal number by a family construction unit and constructing one assumption; An eighth step of calculating the fitness of the assumption by the fitness calculation unit in the same manner as the fitness calculation method in the second step; If the solution replacement part has a goodness-of-fit that is less dominant than recessive by comparing the solution with the highest goodness-of-fit (recessive) in the solution group (Solution 1) and the solution with the lowest goodness-of-fit (dominant) in the assumption (Solution 2) A ninth step of replacing solution 1 with solution 2 in the solution group; A tenth step of determining whether the second to ninth steps have been repeated a set number of times, by an estimate determination unit; And an eleventh step of determining, by the estimate determination unit, a solution having the highest degree of fitness as an estimate, if the second to ninth steps are repeated a set number of times in the tenth step.

상기 다른 실시형태에 의한 유전알고리즘을 이용한 소음제어물질의 물성치 추정 방법은, 상기 제 10 단계에서 상기 제 2 단계 내지 제 9 단계가 설정 횟수 반복되지 않았다면 상기 제 2 단계로 진행될 수 있다.The method for estimating the property values of the noise control material using the genetic algorithm according to the other embodiment may proceed to the second step if the second to ninth steps are not repeated a set number of times in the tenth step.

본 발명의 실시형태들에 의한 유전알고리즘을 이용한 소음제어물질의 물성치 추정 시스템 및 방법에 의하면, 소음제어물질의 소정개수(n개)의 물성치를 그 물성치가 존재 가능한 영역에서 중복을 허용하지 않고 무작위로 설정 개수(m개)만큼 선택한 뒤 열이 물성치가 되도록 하여 m행 n열의 배열 형태로 해 집단을 구성하고, 상기 해 집단의 모든 해에 대한 적합도를 계산하고 m행 1열의 배열 형태로서 결과를 저장하며, 상기 적합도의 계산 결과와 상기 해 집단을 병합하여 m행 n+1 열의 새로운 배열을 구성하고 1열(적합도 열)을 기준으로 오름차순으로 정렬하고, 모든 해 집단에서 2개의 해(부모해)를 중복을 허용하여 선택하며, 교배를 위해 선택된 상기 부모해를 2진수로 부호화하며, 2진수로 부호화된 상기 부모해를 교배하여 자식해를 생성하며, 상기 자식해의 다양성을 증가시키기 위한 돌연변이를 발생시키며, 상기 부모해 및 자식해를 10진수로 복호화하고 하나의 가정을 구성하며, 상기 가정의 적합도를 상기 해 집단의 적합도 계산 방식과 동일한 방법으로 계산하며, 상기 해 집단에서 가장 높은 적합도(열성)를 갖는 해(해 1)와 상기 가정에서 가장 낮은 적합도(우성)를 갖는 해(해 2)를 비교하여, 열성보다 우성이 낮은 적합도를 가질 경우 상기 해 집단에서 해 1을 해 2로 대치하며, 가장 높은 적합도를 갖는 해를 추정치로 결정하도록 구성함으로써, 소음제어물질의 물성치를 수치로 정확히 추정할 수 있다는 뛰어난 효과가 있다. According to the system and method for estimating the property values of noise control materials using the genetic algorithm according to the embodiments of the present invention, the property values of a predetermined number (n) of the noise control materials are randomly applied without allowing overlapping in an area where the property values exist. After selecting as many as the set number (m) and making the columns become material properties, the solution group is formed in an array of m rows and n columns, and the goodness of fit for all solutions of the solution group is calculated, and the result is calculated as an array of m rows and 1 column. Save, and merge the calculation result of the goodness-of-fit and the solution group to form a new array of m rows and n+1 columns, and sort them in ascending order based on column 1 (goodness-of-fit column), and in all solution groups, two solutions (parent solution ) Is selected by allowing duplicates, the parent solution selected for crossing is encoded in binary, the parent solution encoded in binary is crossed to generate a child solution, and mutation to increase the diversity of the child solution Is generated, the parent and child solutions are decoded into decimal numbers to form an assumption, and the fitness of the assumption is calculated in the same way as the fitness of the solution group, and the highest fitness ( Comparing the solution with recessive) (Solution 1) and the solution with the lowest goodness of fit (dominant) in the above assumption (Solution 2), and replacing the solution 1 with the solution 2 in the above solution group if the fit has a lower dominance than recessive And, by configuring the solution with the highest fit to be determined as an estimate, there is an excellent effect that the physical properties of the noise control material can be accurately estimated numerically.

도 1은 본 발명의 실시예에 의한 유전알고리즘을 이용한 소음제어물질의 물성치 추정 시스템에 대한 블록 구성도이다.
도 2는 도 1의 유전알고리즘을 이용한 소음제어물질의 물성치 추정 시스템을 이용한 소음제어물질의 물성치 추정 방법을 설명하기 위한 플로우챠트이다.
도 3은 도 2의 스텝(S500)에서 아버지 및 어머니 염색체가 열별로 2개의 요소로 분할되는 것을 나타낸 도면이다.
도 4는 도 2의 스텝(S500)에서 실험계획법의 직교배열표를 이용하여 1024개의 경우의 수를 16개의 경우의 수로 축약한 것을 나타낸 도면이다.
도 5는 도 2의 스텝(S500)에서 효과적이고 빠른 계산을 위해 4개의 경우의 수를 추가하여 최종적으로 자식해를 구성하는 것을 나타낸 도면이다.
도 6은 도 2의 스텝(S600)에서 돌연변이가 자식해의 염색체 중 하나의 유전자 값을 반대되는 것으로 임의로 변경해주는 형태로 구현되는 것을 나타낸 도면이다.
1 is a block diagram of a system for estimating a property value of a noise control material using a genetic algorithm according to an embodiment of the present invention.
FIG. 2 is a flowchart illustrating a method of estimating a property value of a noise control material using the system for estimating a property value of a noise control material using the genetic algorithm of FIG. 1.
3 is a diagram showing that the father and mother chromosomes are divided into two elements for each column in step S500 of FIG. 2.
FIG. 4 is a diagram showing that the number of 1024 cases is reduced to the number of 16 cases using the orthogonal arrangement table of the experimental design method in step S500 of FIG. 2.
FIG. 5 is a diagram showing that in step S500 of FIG. 2, a child solution is finally constructed by adding the number of four cases for effective and fast calculation.
6 is a view showing that the mutation is implemented in a form in which the gene value of one of the chromosomes of the offspring is changed to the opposite one in step S600 of FIG. 2.

이하, 본 발명의 실시예를 도면을 참조하여 상세히 설명하기로 한다.Hereinafter, embodiments of the present invention will be described in detail with reference to the drawings.

도 1은 본 발명의 실시예에 의한 유전알고리즘을 이용한 소음제어물질의 물성치 추정 시스템에 대한 블록 구성도이다.1 is a block diagram of a system for estimating a property value of a noise control material using a genetic algorithm according to an embodiment of the present invention.

본 발명의 실시예에 의한 유전알고리즘을 이용한 소음제어물질의 물성치 추정 시스템은, 도 1에 도시된 바와 같이, 해 집단 구성부(100), 해 집단의 적합도 계산부(200), 부모해 선택부(300), 2진 부호화부(400), 자식해 생성부(500), 돌연변이 발생부(600), 가정 구성부(700), 가정의 적합도 계산부(800), 해 대치부(900) 및 추정치 결정부(1000)를 포함한다. The system for estimating the property values of noise control materials using a genetic algorithm according to an embodiment of the present invention, as shown in FIG. 1, includes a solution group configuration unit 100, a suitability calculation unit 200, and a parent solution selection unit. (300), binary encoding unit 400, child solution generation unit 500, mutation generation unit 600, assumption construction unit 700, family suitability calculation unit 800, solution replacement unit 900, and It includes an estimated value determining unit 1000.

해 집단 구성부(100)는 소음제어물질의 소정개수(n개)(예컨대, 5열)의 물성치[예컨대, 유동저항, 공극률, 열 특성 길이(thermal characteristic length), 점성 특성 길이 등]를 그 물성치가 존재 가능한 영역에서 중복을 허용하지 않고 무작위로 설정 개수(m개)(예컨대, 150개)만큼 선택한 뒤 열(column)이 물성치가 되도록 하여 m행 n열의 배열 형태로 해 집단을 구성하는 역할을 한다.The solution group component 100 draws a predetermined number (n) of noise control materials (e.g., 5 columns) of physical properties (e.g., flow resistance, porosity, thermal characteristic length, viscous characteristic length, etc.) The role of forming a group in the form of an array of m rows and n columns after selecting as many as the set number (m) (e.g., 150) randomly in the area where material properties can exist, and making the columns become material properties. Do it.

상기 물성치가 존재 가능한 영역은 예컨대, 다음과 같다.Areas in which the physical properties can exist are as follows, for example.

Figure 112019030304587-pat00011
유동저항: 1000 ~ 10000000 MKS Rayls/m
Figure 112019030304587-pat00011
Flow Resistance: 1000 ~ 10000000 MKS Rayls/m

Figure 112019030304587-pat00012
공극률: 0 ~ 1
Figure 112019030304587-pat00012
Porosity: 0 to 1

Figure 112019030304587-pat00013
뒤틀림도: 1 ~ 5
Figure 112019030304587-pat00013
Warp degree: 1 to 5

Figure 112019030304587-pat00014
열특성길이: (1~400)*e-6 m
Figure 112019030304587-pat00014
Thermal characteristic length: (1~400)*e-6 m

Figure 112019030304587-pat00015
점성특성길이: (1~400)*e-6 m
Figure 112019030304587-pat00015
Viscous character length: (1~400)*e-6 m

해 집단의 적합도 계산부(200)는 해 집단 구성부(100)에서 구성된 해 집단의 모든 해에 대한 적합도를 계산하고, m행 1열의 배열 형태로서 결과를 저장하는 역할을 한다.The fitness of the solution group 200 calculates the fitness for all the solutions of the solution group configured in the solution group configuration unit 100 and stores the result as an arrangement of m rows and 1 column.

해 집단의 적합도 계산부(200)는 다음의 [수학식 1]에 의해 적합도를 계산한다.The fitness of the solution group 200 calculates the fitness according to the following [Equation 1].

[수학식 1][Equation 1]

Figure 112019030304587-pat00016
Figure 112019030304587-pat00016

[여기서,

Figure 112019030304587-pat00017
는 적합도를 나타내고,
Figure 112019030304587-pat00018
는 I번째 해를 나타내고,
Figure 112019030304587-pat00019
는 x축 최상단 주파수를 나타내고,
Figure 112019030304587-pat00020
는 상기 x축 최 하단 주파수를 나타내고,
Figure 112019030304587-pat00021
비오(Biot) 이론을 활용하여 추정한 주파수에서의 흡음계수를 나타내며,
Figure 112019030304587-pat00022
는 시험을 통해 계측한 주파수에서의 흡음계수를 나타낸다][here,
Figure 112019030304587-pat00017
Represents the fitness,
Figure 112019030304587-pat00018
Represents the I-th solution,
Figure 112019030304587-pat00019
Represents the uppermost frequency of the x-axis,
Figure 112019030304587-pat00020
Represents the lowermost frequency of the x-axis,
Figure 112019030304587-pat00021
It represents the sound absorption coefficient at the frequency estimated using Biot theory,
Figure 112019030304587-pat00022
Represents the sound absorption coefficient at the frequency measured through the test]

해 집단의 적합도 계산부(200)에서 이루어지는 적합도 계산에 대해 예를 들어 설명하면 다음과 같다.The fitness calculation performed by the fitness calculation unit 200 of the solution group will be described as an example as follows.

총 150개의 해에 대해, 비오(Biot) 이론을 활용하여 흡음계수 곡선을 추정한다. 흡음계수 곡선의 y 축은 흡음계수, x 축은 주파수로써 주어진다. 이때, y 축은 0~1 사이의 값을 가지고, x축은 16 Hz ~ 16000 Hz 사이의 값을 가지며, 주파수별 간격은 1/3 옥타브 밴드를 따른다. 어떤 해를 기준으로 추정한 흡음계수 곡선과 시험결과 얻어진 흡음계수 곡선 간 주파수별 오차 절댓값들의 x축 전체에 대한 합으로서 적합도를 정의한다(수학식 1 참조). 적합도가 "작아진다"는 것은 정답(시험결과)에 근접해간다라는 의미와 상통한다.For a total of 150 solutions, the sound absorption coefficient curve is estimated using Biot theory. The y-axis of the sound absorption coefficient curve is given as the sound absorption coefficient and the x-axis is the frequency. At this time, the y-axis has a value between 0 and 1, the x-axis has a value between 16 Hz and 16000 Hz, and the interval for each frequency follows a 1/3 octave band. The goodness of fit is defined as the sum of the absolute values of errors for each frequency between the sound-absorbing coefficient curve estimated based on a certain solution and the sound-absorbing coefficient curve obtained from the test result on the entire x-axis (see Equation 1). The "smaller" of fit is in common with the meaning of getting closer to the correct answer (test result).

모든 해에 대해 상기와 같은 방법으로 적합도를 계산하고 150행 1열의 배열 형태로서 결과를 저장한다.For all solutions, the fit is calculated in the same way as above, and the result is stored as an array of 150 rows and 1 column.

부모해 선택부(300)는 해 집단의 적합도 계산부(200)에서의 적합도의 계산 결과와 해 집단 구성부(100)에 의해서 구성된 해 집단을 병합하여 m행 n+1 열의 새로운 배열을 구성하여 1열(적합도 열)을 기준으로 오름차순으로 정렬하고, 모든 해 집단에서 2개의 해(부모해)를 중복을 허용하여 선택하는 역할을 한다.The parental solution selection unit 300 merges the calculation result of the fitness of the solution group and the solution group configured by the solution group construction unit 100 to form a new array of m rows and n+1 columns. It sorts in ascending order based on column 1 (conformity column), and selects two solutions (parent solutions) from all solution groups by allowing overlapping.

부모해 선택부(300)는 다음의 [수학식 2]를 이용하여 모든 해 집단에서 2개의 해(부모해)를 중복을 허용하여 선택한다.The parental solution selection unit 300 selects two solutions (parent solutions) by allowing overlapping in all solution groups using the following [Equation 2].

[수학식 2][Equation 2]

Figure 112019030304587-pat00023
Figure 112019030304587-pat00023

[여기서,

Figure 112019030304587-pat00024
Figure 112019030304587-pat00025
순위의 선택확률을 나타내고, sp는 선택압을 나타내며, n은 총 해의 개수를 나타내며, i는 순위를 나타낸다][here,
Figure 112019030304587-pat00024
Is
Figure 112019030304587-pat00025
It represents the probability of selection of the ranking, sp represents the selection pressure, n represents the total number of solutions, and i represents the rank]

부모해 선택부(300)에 의해 2개의 해(부모해)를 중복을 허용하여 선택하는 과정에 대해 예를 들어 설명하면 다음과 같다.A process of selecting two solutions (parent solutions) by allowing overlapping by the parental solution selection unit 300 will be described as follows.

해 집단의 적합도 계산 결과(150행 1열 데이터)와 해 집단(150행 5열)을 서로 병합하여 150행 6열(1열이 적합도)의 새로운 배열을 구성하고, 이를 1열(적합도)을 기준으로 오름차순 정렬한다. 이 배열에서, 1행이 가장 낮은 적합도를 갖는 "우성"의 해가 되며, 150행이 가장 높은 적합도를 갖는 "열성"의 해가 된다. 확률적으로 우성의 해가 열성의 해 보다 높은 확률로 선택되도록 하기 위해 룰렛-휠 연산을 활용하되, 각 해의 선택 확률은 각각의 적합도 순위를 기반으로 정해준다. 이때, 적합도 순위는 행의 번호가 된다. 즉, 1행이 1순위로 가장 높으며 150행이 150순위로 가장 낮다. 이와 같이 순위에 따라 선택 확률을 부여하되, 1순위 해의 선택 확률이 150순위 해의 선택 확률의 4배를 넘어가지 않도록 선택에 압력을 가한다. 여기서 선택압(selection pressure,sp)은 4라고 표현한다. [수학식 2]에 따라 150개의 해 집단에서 2개의 해를 중복을 허용하여 선택한다.The result of calculating the fitness of the solution group (data in 150 rows and 1 column) and the solution group (data in 150 rows and 5 columns) are merged to form a new array of 150 rows and 6 columns (column 1 is goodness-of-fit). Sort by ascending order. In this arrangement, row 1 is the "dominant" solution with the lowest goodness of fit, and row 150 is the "recessive" solution with the highest goodness of fit. The roulette-wheel operation is used to ensure that the dominant solution is selected with a higher probability than the recessive solution probabilistically, but the selection probability of each solution is determined based on each fitness ranking. At this time, the order of fitness is the number of the row. In other words, row 1 is the highest with first priority and row 150 is the lowest with 150th. In this way, the selection probability is given according to the ranking, but pressure is applied to the selection so that the selection probability of the first-ranked solution does not exceed 4 times the selection probability of the 150-ranked solution. Here, the selection pressure (sp) is expressed as 4. According to [Equation 2], two solutions are selected by allowing duplicates from 150 solution groups.

2진 부호화부(400)는 부모해 선택부(300)에서 교배를 위해 선택된 부모해를 2진수로 부호화하는 역할을 한다.The binary encoding unit 400 serves to encode a parent solution selected for mating by the parent solution selection unit 300 into binary numbers.

2진 부호화부(400)에 의한 부호화 과정을 예를 들어 상세히 설명하면 다음과 같다.The coding process by the binary encoder 400 will be described in detail as follows.

부모해로써 선택된 두 해(2행 5열)의 각 요소를 10진수에서 2진수로 부호화한다. 여기서, 1행에 해당하는 유동저항은 10 bit, 2행에 해당하는 공극률은 10 bit, 3행에 해당하는 뒤틀림도는 9 bit, 4행에 해당하는 열특성길이는 16 bit, 그리고 5행에 해당하는 점성특성길이는 16 bit로 부호화 한다. 이렇게 2진수로 구성된 행을 염색체(Chromosome)라 하고 염색체를 구성하는 한 요소를 유전자(Gene)라 한다.Each element of the two solutions selected as parent solutions (row 2, column 5) is encoded from decimal to binary. Here, the flow resistance corresponding to the first row is 10 bits, the porosity corresponding to the second row is 10 bits, the twisting degree corresponding to the third row is 9 bits, the column characteristic length corresponding to the fourth row is 16 bits, and the fifth row is The corresponding viscous characteristic length is encoded as 16 bits. The row composed of binary numbers in this way is called a chromosome, and one element constituting the chromosome is called a gene.

예를 들어 이진수 1011에 대해,For example, for the binary number 1011,

1011은 4 bit 이진수이고,1011 is a 4 bit binary number,

1011은 염색체이며,1011 is a chromosome,

염색체를 구성하는 1, 0, 1, 1 은 각각 유전자이다.Each of 1, 0, 1, and 1 constituting a chromosome is a gene.

자식해 생성부(500)는 2진 부호화부(400)에 의해 2진수로 부호화된 부모해를 교배하여 자식해를 생성하는 역할을 한다. The child solution generator 500 serves to generate a child solution by crossing parent solutions encoded in binary numbers by the binary encoding unit 400.

자식해 생성부(500)에 의해 자식해를 생성하는 과정에 대해 예를 들어 설명하기로 한다.A process of generating a child solution by the child solution generation unit 500 will be described as an example.

부모해의 교배를 위해 각 염색체 내부에 임의로 한 지점을 선택한다. 아버지(1행) 및 어머니(2행)의 염색체는, 도 3에 도시된 바와 같이, 열 별로 2개의 요소로 분할된다.For mating of the parental solution, a point is randomly selected inside each chromosome. The chromosomes of father (row 1) and mother (row 2) are divided into two elements for each column, as shown in FIG. 3.

도 3을 기준으로 설명하면, 유동저항에 대한 자식해의 존재 경우의 수는 CF1CF2, CF1CM2, CM1CF2, 그리고 CM1CM2 총 4개가 있다. 이는 유동저항 등 5개의 물성(제어변수)에 대해, 각각이 4개의 변수(제어수준)를 가짐을 의미한다. 5개의 제어변수들이 모두 4의 제어수준을 갖는다면, 조합 가능한 자식 염색체 경우의 수는 45(=1024)개가 된다. 모든 경우의 수를 전부 고려하기 어렵고 비효율적이므로, 모든 경우의 수를 고려한 경우와 거의 동일한 효과를 갖는 실험계획법(Design Of Experiment, DOE)의 직교배열표(Orthogonal array)를 이용하여 1024개의 경우의 수를 단 16개의 경우로 축약한다(다음의 [수학식 3] 및 도 4 참조)Referring to FIG. 3, the number of cases in which child solutions for the flow resistance exist are four in total: CF1CF2, CF1CM2, CM1CF2, and CM1CM2. This means that for 5 properties (control variables) such as flow resistance, each has 4 variables (control levels). If all 5 control variables have a control level of 4, the number of possible child chromosome cases is 4 5 (=1024). Since it is difficult and inefficient to consider the number of all cases, the number of 1024 cases using an orthogonal array of Design Of Experiment (DOE) that has almost the same effect as the case of considering all cases. Is reduced to only 16 cases (see [Equation 3] and Fig. 4 below)

[수학식 3][Equation 3]

Figure 112019030304587-pat00026
Figure 112019030304587-pat00026

[여기서, L은 직교배열표 임을 의미하고, V는 제어수준을 나타내고, m은 제어변수의 개수를 나타내고, n은 시험 개수로서

Figure 112019030304587-pat00027
로 정의된다][Here, L means an orthogonal array table, V represents the control level, m represents the number of control variables, and n represents the number of tests.
Figure 112019030304587-pat00027
Is defined as]

이때, 3열의 뒤틀림도와 5열의 점성특성길이가 가장 마지막에 수렴하는 점에 착안하여 효과적이고 빠른 계산을 위해, 도 5에 도시된 바와 같이, 4개의 경우의 수를 추가하여 최종적으로 자식해를 구성한다.At this time, for an effective and fast calculation, focusing on the point at which the distortion of the 3rd column and the viscous characteristic length of the 5th column converge at the last, as shown in Fig. 5, the number of 4 cases is added to finally construct a child solution. do.

돌연변이 발생부(600)는 자식해 생성부(500)에 의해 구성된 자식해의 다양성을 증가시키기 위한 돌연변이를 발생시키는 역할을 한다.The mutation generating unit 600 serves to generate a mutation to increase the diversity of the child solutions constituted by the child solution generating unit 500.

돌연변이를 발생시키는 과정에 대해서 예를 들어 상세히 설명하기로 한다.The process of generating the mutation will be described in detail with an example.

초기 해 집단이 정답과 동떨어진 영역에서 모두 선택되었다면, 이들을 아무리 교배해도 정답에 수렴할 수 없는 국소수렴(Local converging)의 문제가 있으므로 국소해로의 수렴 탈출을 위해 자식해의 일부에 대해 돌연변이를 일으킨다. 정답으로 어느 정도 수렴하는 1세대부터 150세대까지는 부모의 유전정보를 가능한 보존하기 위해 2%의 돌연변이 확률을 사용하고, 정답에 어느 정도 수렴한 150세대 이후에는 기본적으로 60%의 높은 돌연변이 확률을 주어 해의 다양성을 최대한 확보한다. 돌연변이는, 도 6에 도시된 바와 같이, 자식해의 염색체 중 하나의 유전자 값을 반대되는 것으로 임의 변경해주는 형태로 구현한다If the initial solution group is all selected in a region that is far from the correct answer, there is a problem of local converging that cannot converge to the correct answer no matter how much they are crossed, so mutations are generated for some of the child solutions to escape convergence to the local solution. . From the 1st generation to the 150th generation, which converges to some extent as the correct answer, a 2% mutation probability is used to preserve the parent's genetic information as much as possible, and after the 150th generation when the correct answer converges to some extent, a high mutation probability of 60% is basically given. Maximize the diversity of the year. Mutation, as shown in Fig. 6, is implemented in the form of randomly changing the value of one of the chromosomes of the offspring to the opposite one.

가정 구성부(700)는 부모해 선택부(300) 및 자식해 생성부(500)에 의해 생성된 부모해 및 자식해를 10진수로 복호화하고 하나의 가정을 구성하는 역할을 한다.The home configuration unit 700 decodes the parent solution and the child solution generated by the parent solution selection unit 300 and the child solution generation unit 500 into decimal numbers, and configures one home.

예컨대, 2개의 부모해와 20개의 자식해를 10진수로 복호화하여 하나의 가정을 구성할 수 있다(22행 5열).For example, two parent solutions and 20 child solutions can be decoded into a decimal number to form an assumption (row 22, column 5).

가정의 적합도 계산부(800)는 가정의 적합도를 해 집단의 적합도 계산부(200)에서의 적합도 계산 방식과 동일한 방법으로 계산하는 역할을 한다.The fitness of the assumption calculation unit 800 serves to calculate the fitness of the assumption in the same manner as the fitness calculation method of the fitness calculation unit 200 of the group.

해 대치부(900)는 해 집단의 적합도 계산부(200)에 의해 계산된 해 집단에서 가장 높은 적합도(열성)를 갖는 해(해 1)와 가정의 적합도 계산부(800)에 의해 계산된 가정에서 가장 낮은 적합도(우성)를 갖는 해(해 2)를 비교하여, 열성보다 우성이 낮은 적합도를 가질 경우 해 집단에서 해 1을 해 2로 대치하는 역할을 한다.The solution replacement unit 900 is a solution (solution 1) having the highest fitness (recessiveness) in the solution group calculated by the fitness calculation unit 200 of the solution group and an assumption calculated by the fitness calculation unit 800 The solution (Solution 2) with the lowest fit (dominant) in is compared, and if the dominant fit has a lower fit than the recessive, it plays a role of replacing Solution 1 with Solution 2.

추정치 결정부(1000)는 가장 높은 적합도를 갖는 해를 추정치로 결정하는 역할을 한다.The estimate determination unit 1000 serves to determine a solution having the highest fitness as an estimate.

이하, 상기한 바와 같이 구성된 본 발명의 실시예에 의한 유전알고리즘을 이용한 소음제어물질의 물성치 추정 시스템을 이용한 소음제어물질의 물성치 추정 방법을 설명하기로 한다.Hereinafter, a method for estimating the properties of the noise control material using the system for estimating the properties of the noise control material using the genetic algorithm according to an embodiment of the present invention constructed as described above will be described.

도 2는 도 1의 유전알고리즘을 이용한 소음제어물질의 물성치 추정 시스템을 이용한 소음제어물질의 물성치 추정 방법을 설명하기 위한 플로우챠트로서, 여기서 S는 스텝(step)을 나타낸다.FIG. 2 is a flowchart illustrating a method of estimating a property value of a noise control material using the system for estimating property value of a noise control material using the genetic algorithm of FIG. 1, where S represents a step.

먼저, 해 집단 구성부(100)가 소음제어물질의 소정개수(n개)의 물성치를 그 물성치가 존재 가능한 영역에서 중복을 허용하지 않고 무작위로 설정 개수(m개)만큼 선택한 뒤 열이 물성치가 되도록 하여 m행 n열의 배열 형태로 해 집단을 구성한다(S100)First, the solution group component 100 randomly selects a predetermined number (n) of material properties of the noise control material as much as a set number (m) in an area where the material properties can exist without allowing overlapping, So that the group is formed in an array of m rows and n columns (S100)

이어서, 해 집단의 적합도 계산부(200)가 상기 스텝(S100)에서 구성된 해 집단의 모든 해에 대한 적합도를 계산하고, m행 1열의 배열 형태로서 결과를 저장한다(S200).Subsequently, the fitness calculation unit 200 of the solution group calculates the fitness for all the solutions of the solution group configured in step S100, and stores the result as an arrangement of m rows and 1 column (S200).

이후, 부모해 선택부(300)가 상기 스텝(S200)에서의 적합도의 계산 결과와 상기 스텝(S100)에서 구성된 해 집단을 병합하여 m행 n+1 열의 새로운 배열을 구성하고 1열(적합도 열)을 기준으로 오름차순으로 정렬하고, 모든 해 집단에서 2개의 해(부모해)를 중복을 허용하여 선택한다(S300).Thereafter, the parental solution selection unit 300 merges the calculation result of the fitness at step S200 and the solution group configured at step S100 to form a new array of m rows and n+1 columns, ), and selects by allowing overlapping of two solutions (parent solutions) in all solution groups (S300).

그리고 2진 부호화부(400)가 상기 스텝(S300)에서 교배를 위해 선택된 부모해를 2진수로 부호화하고(S400), 자식해 생성부(500)가 2진수로 부호화된 부모해를 교배하여 자식해를 생성한다(S500).Then, the binary encoding unit 400 encodes the parent solution selected for crossing in the step (S300) into binary numbers (S400), and the child solution generation unit 500 crosses the parent solution encoded in binary numbers, A solution is generated (S500).

이어서, 돌연변이 발생부(600)가 상기 스텝(S500)에서 생성된 자식해의 다양성을 증가시키기 위한 돌연변이를 발생시키고(S600), 가정 구성부(700)가 부모해 및 자식해를 10진수로 복호화하고 하나의 가정을 구성한다(S700).Subsequently, the mutation generator 600 generates a mutation to increase the diversity of the child solutions generated in the step S500 (S600), and the home construction unit 700 decodes the parent and child solutions into decimal numbers. And constitutes one assumption (S700).

이후, 가정의 적합도 계산부(800)가 상기 스텝(S700)에서 구성된 가정의 적합도를 상기 스텝(S200)에서의 적합도 계산 방식과 동일한 방법으로 계산한다(S800).Thereafter, the fitness of the assumption calculation unit 800 calculates the fitness of the assumption configured in the step S700 in the same manner as the calculation method of the fitness in the step S200 (S800).

이어서, 해 대치부(900)가 상기 스텝(S200)에서 계산된 해 집단에서 가장 높은 적합도(열성)를 갖는 해(해 1)와 상기 스텝(S800)에서 계산된 가정에서 가장 낮은 적합도(우성)를 갖는 해(해 2)를 비교하여, 열성보다 우성이 낮은 적합도를 가지는 지의 여부를 결정한다(S900).Subsequently, the solution replacement unit 900 has the highest fitness (recessive) in the solution group calculated in step S200 (solution 1) and the lowest fitness (dominant) in the assumption calculated in step S800 By comparing the solution with (Solution 2), it is determined whether the dominance has a lower fit than the recessive (S900).

상기 스텝(S900)에서 열성보다 우성이 낮은 적합도를 가지면(YES), 해 집단에서 해 1을 해 2로 대치한다(S1000). In the step S900, if the dominance is lower than the recessive fit (YES), the solution 1 is replaced by the solution 2 in the solution group (S1000).

한편, 상기 스텝(S900)에서 우성이 열성보다 낮은 적합도를 가지지 않으면(NO) 본 스텝(S900)을 반복수행한다. On the other hand, in the step S900, if the dominance does not have a degree of fitness lower than the recessive (NO), the step S900 is repeated.

이어서, 추정치 결정부(1000)가 상기 스텝(S200) 내지 스텝(S1000)이 설정 횟수 반복되는 지의 여부를 결정한다(S1100).Next, the estimate determination unit 1000 determines whether the steps S200 to S1000 are repeated a set number of times (S1100).

상기 스텝(S1100)에서 상기 스텝(S200) 내지 스텝(S1000)이 설정 횟수 반복되었다면(YES), 추정치 결정부(1000)에 의해 가장 높은 적합도를 갖는 해가 추정치로 결정된다(S1200). If the steps S200 to S1000 are repeated a set number of times (YES) in the step S1100, the solution having the highest degree of fit is determined by the estimate determining unit 1000 as an estimate value (S1200).

한편, 상기 스텝(S1100)에서 상기 스텝(S200) 내지 스텝(S1000)이 설정 횟수 반복되지 않는다면(NO), 상기 스텝(S200)으로 진행된다.On the other hand, if the steps S200 to S1000 are not repeated the set number of times (NO) in the step S1100, the process proceeds to the step S200.

본 발명의 실시예에 의한 유전알고리즘을 이용한 소음제어물질의 물성치 추정 시스템 및 방법에 의하면, 소음제어물질의 소정개수(n개)의 물성치를 그 물성치가 존재 가능한 영역에서 중복을 허용하지 않고 무작위로 설정 개수(m개)만큼 선택한 뒤 열이 물성치가 되도록 하여 m행 n열의 배열 형태로 해 집단을 구성하고, 상기 해 집단의 모든 해에 대한 적합도를 계산하고 m행 1열의 배열 형태로서 결과를 저장하며, 상기 적합도의 계산 결과와 상기 해 집단을 병합하여 m행 n+1 열의 새로운 배열을 구성하고 1열(적합도 열)을 기준으로 오름차순으로 정렬하고, 모든 해 집단에서 2개의 해(부모해)를 중복을 허용하여 선택하며, 교배를 위해 선택된 상기 부모해를 2진수로 부호화하며, 2진수로 부호화된 상기 부모해를 교배하여 자식해를 생성하며, 상기 자식해의 다양성을 증가시키기 위한 돌연변이를 발생시키며, 상기 부모해 및 자식해를 10진수로 복호화하고 하나의 가정을 구성하며, 상기 가정의 적합도를 상기 해 집단의 적합도 계산 방식과 동일한 방법으로 계산하며, 상기 해 집단에서 가장 높은 적합도(열성)를 갖는 해(해 1)와 상기 가정에서 가장 낮은 적합도(우성)를 갖는 해(해 2)를 비교하여, 열성보다 우성이 낮은 적합도를 가질 경우 상기 해 집단에서 해 1을 해 2로 대치하며, 가장 높은 적합도를 갖는 해를 추정치로 결정하도록 구성함으로써, 소음제어물질의 물성치를 수치로 정확히 추정할 수 있다. According to the system and method for estimating the property values of noise control materials using a genetic algorithm according to an embodiment of the present invention, the property values of a predetermined number (n) of the noise control materials are randomly obtained without allowing overlapping in an area where the property values exist. After selecting as many as the set number (m), configure the solution group in the form of an array of m rows and n columns by making the columns become material properties, calculate the fit for all solutions of the solution group, and store the result as an array of m rows and 1 column. And, by merging the calculation result of the goodness of fit and the solution group, a new array of m rows and n+1 columns is formed and sorted in ascending order based on column 1 (the fitness column), and two solutions (parent solutions) in all solution groups The parent solution selected for mating is selected by allowing duplication, encoding the parent solution selected for mating in binary, generating a child solution by crossing the parent solution coded in binary, and mutating mutations to increase the diversity of the child solution. The parent and child solutions are decoded into decimal numbers to form an assumption, and the fitness of the assumption is calculated in the same way as the fitness of the solution group, and the highest fitness (recessive) is calculated in the solution group. ) With the solution (Solution 1) and the solution with the lowest goodness of fit (dominant) in the assumption (Solution 2), and if the fit has a lower dominance than the recessive, the solution 1 is replaced by the solution 2 In addition, by configuring the solution with the highest fitness to be determined as an estimate, the properties of the noise control material can be accurately estimated numerically.

도면과 명세서에는 최적의 실시예가 개시되었으며, 특정한 용어들이 사용되었으나 이는 단지 본 발명의 실시형태를 설명하기 위한 목적으로 사용된 것이지 의미를 한정하거나 특허청구범위에 기재된 본 발명의 범위를 제한하기 위하여 사용된 것은 아니다. 그러므로 본 기술 분야의 통상의 지식을 가진 자라면 이로부터 다양한 변형 및 균등한 타 실시예가 가능하다는 점을 이해할 수 있을 것이다. 따라서 본 발명의 진정한 기술적 보호범위는 첨부된 특허청구범위의 기술적 사상에 의해 정해져야 할 것이다.In the drawings and specification, the best embodiments have been disclosed, and specific terms are used, but these are only used for the purpose of describing the embodiments of the present invention, and are used to limit the meaning or the scope of the present invention described in the claims. Was not done. Therefore, those of ordinary skill in the art will appreciate that various modifications and equivalent other embodiments are possible therefrom. Therefore, the true technical scope of the present invention should be determined by the technical spirit of the appended claims.

100: 해 집단 구성부
200: 해 집단의 적합도 계산부
300: 부모해 선택부
400: 2진 부호화부
500: 자식해 생성부
600: 돌연변이 발생부
700: 가정 구성부
800: 가정의 적합도 계산부
900: 해 대치부
1000: 추정치 결정부
100: sea group composition
200: fit calculation unit of solution group
300: parental choice
400: binary encoding unit
500: child year generator
600: mutation generator
700: home building block
800: family fitness calculation unit
900: Sun replacement
1000: estimate determination unit

Claims (9)

유전알고리즘을 이용하여 소음제어물질의 물성치에 대한 해들중 가장 높은 적합도를 갖는 해를 추정치로 결정하는 유전알고리즘을 이용한 소음제어물질의 물성치 추정 시스템으로서,
소음제어물질의 소정개수(n개)의 물성치를 그 물성치가 존재 가능한 영역에서 중복을 허용하지 않고 무작위로 설정 개수(m개)만큼 선택한 뒤 열이 물성치가 되도록 하여 m행 n열의 배열 형태로 해 집단을 구성하는 해 집단 구성부;
상기 해 집단의 모든 해에 대한 적합도를 계산하고 m행 1열의 배열 형태로서 결과를 저장하는 해 집단의 적합도 계산부;
상기 적합도의 계산 결과와 상기 해 집단을 병합하여 m행 n+1 열의 새로운 배열을 구성하고 1열(적합도 열)을 기준으로 오름차순으로 정렬하고, 모든 해 집단에서 2개의 해(부모해)를 중복을 허용하여 선택하는 부모해 선택부;
교배를 위해 선택된 상기 부모해를 2진수로 부호화하는 2진 부호화부;
2진수로 부호화된 상기 부모해를 교배하여 자식해를 생성하는 자식해 생성부;
상기 자식해의 다양성을 증가시키기 위한 돌연변이를 발생시키는 돌연변이 발생부;
상기 부모해 및 자식해를 10진수로 복호화하고 하나의 가정을 구성하는 가정 구성부;
상기 가정의 적합도를 상기 해 집단의 적합도 계산부에서의 적합도 계산 방식과 동일한 방법으로 계산하는 가정의 적합도 계산부;
상기 해 집단에서 가장 높은 적합도(열성)를 갖는 해(해 1)와 상기 가정에서 가장 낮은 적합도(우성)를 갖는 해(해 2)를 비교하여, 열성보다 우성이 낮은 적합도를 가질 경우 상기 해 집단에서 해 1을 해 2로 대치하는 해 대치부; 및
가장 높은 적합도를 갖는 해를 추정치로 결정하는 추정치 결정부;를 포함하고,
상기 해 집단의 적합도 계산부 및 가정의 적합도 계산부는 다음의 [수학식 1]에 의해 적합도를 계산하는, 유전알고리즘을 이용한 소음제어물질의 물성치 추정 시스템.

[수학식 1]
Figure 112020503455042-pat00055


[여기서,
Figure 112020503455042-pat00056
는 적합도를 나타내고,
Figure 112020503455042-pat00057
는 I번째 해를 나타내고,
Figure 112020503455042-pat00058
는 x축 최상단 주파수를 나타내고,
Figure 112020503455042-pat00059
는 x축 최하단 주파수를 나타내고,
Figure 112020503455042-pat00060
비오(Biot) 이론을 활용하여 추정한 주파수에서의 흡음계수를 나타내며,
Figure 112020503455042-pat00061
는 시험을 통해 계측한 주파수에서의 흡음계수를 나타낸다]
As a system for estimating the properties of the noise control material using a genetic algorithm that determines the solution with the highest fitness as an estimate among the solutions for the properties of the noise control material using the genetic algorithm,
The mass properties of a predetermined number (n) of the noise control material are randomly selected as many as the set number (m) without allowing duplication in the area where the properties can exist, and the columns are set to be the mass properties, so that they are arranged in m rows and n columns. A group composition unit that constitutes a group;
A goodness-of-fit calculation unit for calculating the fitness for all the solutions of the solution group and storing the result as an array of m rows and 1 column;
By merging the calculation result of the goodness-of-fit and the solution group, a new array of m rows and n+1 columns is formed and sorted in ascending order based on column 1 (goodness-of-fit column), and two solutions (parent solutions) are duplicated in all solution groups. Parental selection unit to select by allowing;
A binary encoder for encoding the parent solution selected for mating into binary numbers;
A child solution generator for generating a child solution by crossing the parent solution encoded in binary numbers;
A mutation generating unit that generates a mutation for increasing the diversity of the offspring;
A home construction unit that decodes the parent and child solutions into decimal numbers and configures one assumption;
A fitness of the assumption that calculates the fitness of the assumption by the same method as the fitness of the solution group;
Compare the solution with the highest goodness of fit (recessive) in the solution group (Solution 1) and the solution with the lowest goodness of fit (dominant) in the assumption (Solution 2). The year replacement unit replacing year 1 with year 2 in; And
Including; an estimate determination unit for determining a solution having the highest fitness as an estimate,
A system for estimating a property value of a noise control material using a genetic algorithm, which calculates the fitness according to the following [Equation 1].

[Equation 1]
Figure 112020503455042-pat00055


[here,
Figure 112020503455042-pat00056
Represents the fitness,
Figure 112020503455042-pat00057
Represents the I-th solution,
Figure 112020503455042-pat00058
Represents the uppermost frequency of the x-axis,
Figure 112020503455042-pat00059
Represents the lowest frequency of the x-axis,
Figure 112020503455042-pat00060
It represents the sound absorption coefficient at the frequency estimated using Biot theory,
Figure 112020503455042-pat00061
Represents the sound absorption coefficient at the frequency measured through the test]
제 1 항에 있어서,
상기 물성치는
유동저항, 공극률, 열 특성 길이(thermal characteristic length) 및 점성 특성 길이를 포함하는, 유전알고리즘을 이용한 소음제어물질의 물성치 추정 시스템.
The method of claim 1,
The above properties are
A system for estimating the properties of a noise control material using a dielectric algorithm, including flow resistance, porosity, thermal characteristic length, and viscous characteristic length.
삭제delete 제 1 항에 있어서,
상기 부모해 선택부는 다음의 [수학식 2]를 이용하여 모든 해 집단에서 2개의 해(부모해)를 중복을 허용하여 선택하는, 유전알고리즘을 이용한 소음제어물질의 물성치 추정 시스템.

[수학식 2]
Figure 112019030304587-pat00035


[여기서,
Figure 112019030304587-pat00036
Figure 112019030304587-pat00037
순위의 선택확률을 나타내고, sp는 선택압을 나타내며, n은 총 해의 개수를 나타내며, i는 순위를 나타낸다]
The method of claim 1,
The parental solution selection unit selects two solutions (parent solutions) in all solution groups using the following [Equation 2] by allowing overlapping, a system for estimating material properties of a noise control material using a genetic algorithm.

[Equation 2]
Figure 112019030304587-pat00035


[here,
Figure 112019030304587-pat00036
Is
Figure 112019030304587-pat00037
It represents the probability of selection of the ranking, sp represents the selection pressure, n represents the total number of solutions, and i represents the rank]
제 1 항에 기재된 유전알고리즘을 이용한 소음제어물질의 물성치 추정 시스템을 이용한 소음제어물질의 물성치 추정 방법으로서,
해 집단 구성부가 소음제어물질의 소정개수(n개)의 물성치를 그 물성치가 존재 가능한 영역에서 중복을 허용하지 않고 무작위로 설정 개수(m개)만큼 선택한 뒤 열이 물성치가 되도록 하여 m행 n열의 배열 형태로 해 집단을 구성하는 제 1 단계;
해 집단의 적합도 계산부가 상기 해 집단의 모든 해에 대한 적합도를 계산하고 m행 1열의 배열 형태로서 결과를 저장하는 제 2 단계;
부모해 선택부가 상기 적합도의 계산 결과와 상기 해 집단을 병합하여 m행 n+1 열의 새로운 배열을 구성하고 1열(적합도 열)을 기준으로 오름차순으로 정렬하고, 모든 해 집단에서 2개의 해(부모해)를 중복을 허용하여 선택하는 제 3 단계;
2진 부호화부가 교배를 위해 선택된 상기 부모해를 2진수로 부호화하는 제 4 단계;
자식해 생성부가 2진수로 부호화된 상기 부모해를 교배하여 자식해를 생성하는 제 5 단계;
돌연변이 발생부가 상기 자식해의 다양성을 증가시키기 위한 돌연변이를 발생시키는 제 6 단계;
가정 구성부가 상기 부모해 및 자식해를 10진수로 복호화하고 하나의 가정을 구성하는 제 7 단계;
가정의 적합도 계산부가 상기 가정의 적합도를 상기 제 2 단계에서의 적합도 계산 방식과 동일한 방법으로 계산하는 제 8 단계;
해 대치부가 상기 해 집단에서 가장 높은 적합도(열성)를 갖는 해(해 1)와 상기 가정에서 가장 낮은 적합도(우성)를 갖는 해(해 2)를 비교하여, 열성보다 우성이 낮은 적합도를 가질 경우 상기 해 집단에서 해 1을 해 2로 대치하는 제 9 단계;
추정치 결정부가 상기 제 2 단계 내지 제 9 단계가 설정 횟수 반복되었는지의 여부를 결정하는 제 10 단계; 및
상기 제 10 단계에서 상기 제 2 단계 내지 제 9 단계가 설정 횟수 반복되었다면, 상기 추정치 결정부에 의해 가장 높은 적합도를 갖는 해가 추정치로 결정되는 제 11 단계;를 포함하는, 소음제어물질의 물성치 추정 방법.
As a method for estimating the properties of a noise control material using the system for estimating properties of the noise control material using the genetic algorithm described in claim 1,
The solution group composition part randomly selects a predetermined number (n) of material properties of the noise control material as much as the set number (m) without allowing overlapping in the area where the material properties can exist. A first step of constructing a solution group in an array form;
A second step of calculating, by a fitness calculation unit of a solution group, the fitness for all solutions of the solution group and storing the result as an array of m rows and 1 column;
The parental solution selection unit merges the calculation result of the goodness-of-fit and the solution group to form a new array of m rows and n+1 columns, and sorts them in ascending order based on column 1 (goodness-of-fit column), and in all solution groups, two solutions (parent A third step of selecting a solution by allowing redundancy;
A fourth step of encoding, by a binary encoding unit, the parent solution selected for mating into binary numbers;
A fifth step of generating a child solution by crossing the parent solution encoded as a binary number by a child solution generator;
A sixth step of generating a mutation for increasing the diversity of the offspring by the mutation generating unit;
A seventh step of decoding the parent solution and the child solution to a decimal number by a family construction unit and constructing one assumption;
An eighth step of calculating the fitness of the assumption by the fitness calculation unit in the same manner as the fitness calculation method in the second step;
If the solution replacement part has a goodness-of-fit that is less dominant than recessive by comparing the solution with the highest goodness-of-fit (recessive) in the solution group (Solution 1) and the solution with the lowest goodness-of-fit (dominant) in the assumption (Solution 2) A ninth step of replacing solution 1 with solution 2 in the solution group;
A tenth step of determining whether the second to ninth steps have been repeated a set number of times, by an estimate determining unit; And
In the tenth step, if the second to ninth steps are repeated a set number of times, an eleventh step of determining, by the estimate determination unit, a solution having the highest fitness as an estimate value; estimating the physical properties of the noise control material, including: Way.
제 5 항에 있어서,
상기 제 10 단계에서 상기 제 2 단계 내지 제 9 단계가 설정 횟수 반복되지 않았다면 상기 제 2 단계로 진행되는, 소음제어물질의 물성치 추정 방법.
The method of claim 5,
In the tenth step, if the second to ninth steps are not repeated a set number of times, the method proceeds to the second step.
제 5 항에 있어서,
상기 물성치는
유동저항, 공극률, 열 특성 길이(thermal characteristic length) 및 점성 특성 길이를 포함하는, 소음제어물질의 물성치 추정 방법.
The method of claim 5,
The above properties are
A method for estimating properties of a noise control material, including flow resistance, porosity, thermal characteristic length, and viscous characteristic length.
제 5 항에 있어서,
상기 해 집단의 적합도 계산부 및 가정의 적합도 계산부는 다음의 [수학식 1]에 의해 적합도를 계산하는, 소음제어물질의 물성치 추정 방법.

[수학식 1]
Figure 112019030304587-pat00038


[여기서,
Figure 112019030304587-pat00039
는 적합도를 나타내고,
Figure 112019030304587-pat00040
는 I번째 해를 나타내고,
Figure 112019030304587-pat00041
는 x축 최상단 주파수를 나타내고,
Figure 112019030304587-pat00042
는 상기 x축 최 하단 주파수를 나타내고,
Figure 112019030304587-pat00043
비오(Biot) 이론을 활용하여 추정한 주파수에서의 흡음계수를 나타내며,
Figure 112019030304587-pat00044
는 시험을 통해 계측한 주파수에서의 흡음계수를 나타낸다]
The method of claim 5,
The goodness-of-fit calculation unit of the solution group and the goodness-of-fit calculation unit of the assumption calculate the fitness by the following [Equation 1].

[Equation 1]
Figure 112019030304587-pat00038


[here,
Figure 112019030304587-pat00039
Represents the fitness,
Figure 112019030304587-pat00040
Represents the I-th solution,
Figure 112019030304587-pat00041
Represents the uppermost frequency of the x-axis,
Figure 112019030304587-pat00042
Represents the lowermost frequency of the x-axis,
Figure 112019030304587-pat00043
It represents the sound absorption coefficient at the frequency estimated using Biot theory,
Figure 112019030304587-pat00044
Represents the sound absorption coefficient at the frequency measured through the test]
제 5 항에 있어서,
상기 부모해 선택부는 다음의 [수학식 2]를 이용하여 모든 해 집단에서 2개의 해(부모해)를 중복을 허용하여 선택하는, 소음제어물질의 물성치 추정 방법.

[수학식 2]
Figure 112019030304587-pat00045


[여기서,
Figure 112019030304587-pat00046
Figure 112019030304587-pat00047
순위의 선택확률을 나타내고, sp는 선택압을 나타내며, n은 총 해의 개수를 나타내며, i는 순위를 나타낸다]



The method of claim 5,
The parental solution selection unit selects two solutions (parent solutions) in all solution groups by allowing overlapping using the following [Equation 2], a method for estimating material properties of a noise control material.

[Equation 2]
Figure 112019030304587-pat00045


[here,
Figure 112019030304587-pat00046
Is
Figure 112019030304587-pat00047
It represents the probability of selection of the ranking, sp represents the selection pressure, n represents the total number of solutions, and i represents the rank]



KR1020190033519A 2019-03-25 2019-03-25 System and method for estimating property of noise control material using genetic algorithm KR102198714B1 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
KR1020190033519A KR102198714B1 (en) 2019-03-25 2019-03-25 System and method for estimating property of noise control material using genetic algorithm

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
KR1020190033519A KR102198714B1 (en) 2019-03-25 2019-03-25 System and method for estimating property of noise control material using genetic algorithm

Publications (2)

Publication Number Publication Date
KR20200113392A KR20200113392A (en) 2020-10-07
KR102198714B1 true KR102198714B1 (en) 2021-01-05

Family

ID=72884417

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1020190033519A KR102198714B1 (en) 2019-03-25 2019-03-25 System and method for estimating property of noise control material using genetic algorithm

Country Status (1)

Country Link
KR (1) KR102198714B1 (en)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR102480727B1 (en) * 2021-01-26 2022-12-26 주식회사 에이엘 Method for optimizing sound absorbing material structure for weight reduction of silencer using genetic algorithm
IT202100005120A1 (en) * 2021-03-04 2022-09-04 Milano Politecnico COMPUTER IMPLEMENTED METHOD AND SYSTEM TO MINIMIZE STRUCTURAL COSTS, MAXIMIZE FREE SPACE AND MINIMIZE ENVIRONMENTAL IMPACT IN A CONCEPTUAL BUILDING DESIGN
KR102560723B1 (en) * 2022-02-28 2023-07-27 한국전력기술 주식회사 Apparatus and method for operating cooling and heating facilities using genetic algorithm
KR20240074097A (en) 2022-11-18 2024-05-28 세종대학교산학협력단 Material and process optimization method by active learning

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4604477B2 (en) * 2003-11-18 2011-01-05 横浜ゴム株式会社 Method for estimating physical properties, method for optimal design of sound transmission body

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2009236630A (en) 2008-03-26 2009-10-15 Mitsubishi Heavy Ind Ltd Method of estimating physical property value of ceramic, method of estimating physical property value of heat-insulating coating material, method of estimating remaining life of heat-insulating coating material, and method of estimating remaining life of high-temperature member, and method of obtaining value of physical property
KR100986160B1 (en) * 2008-09-12 2010-10-07 인하대학교 산학협력단 System and method for manufacturing mask using genetic algorithm and dna computing

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4604477B2 (en) * 2003-11-18 2011-01-05 横浜ゴム株式会社 Method for estimating physical properties, method for optimal design of sound transmission body

Also Published As

Publication number Publication date
KR20200113392A (en) 2020-10-07

Similar Documents

Publication Publication Date Title
KR102198714B1 (en) System and method for estimating property of noise control material using genetic algorithm
Alba et al. The exploration/exploitation tradeoff in dynamic cellular genetic algorithms
Bungartz et al. Multivariate quadrature on adaptive sparse grids
Wright Genetic algorithms for real parameter optimization
Levitin The universal generating function in reliability analysis and optimization
Esfahanizadeh et al. Finite-length construction of high performance spatially-coupled codes via optimized partitioning and lifting
Ceylan Developing combined genetic algorithm—hill-climbing optimization method for area traffic control
CN110059405B (en) Construction method of high-quality Steiner minimum tree with differential evolution under X structure
CN109858698B (en) Vehicle supply and demand optimization method, device, equipment and storage medium for target area
CN110690907B (en) Method for estimating deletion mode of known branch information turbo code
Kaveh et al. Simultaneous topology and size optimization of structures by genetic algorithm using minimal length chromosome
EP1055196A2 (en) Code compaction by evolutionary algorithm
Darmawan et al. Low-depth random Clifford circuits for quantum coding against Pauli noise using a tensor-network decoder
CN109257050A (en) A kind of reparation binary code generator matrix building method and restorative procedure
CN111542069B (en) Method for realizing wireless AP deployment optimization based on rapid non-dominant genetic algorithm
Lang Space-partitioned ND-trees for the dynamic nondominance problem
Bang et al. Improved genetic algorithm for minimum latency problem
CN114329702B (en) High-stability method and device for calibrating and designing reaction spectrum based on improved differential evolution algorithm
Dhahri et al. Hierarchical learning algorithm for the beta basis function neural network
CN107623595B (en) Network server deployment method and system
Pitman et al. Growth of the Brownian forest
Boyer et al. Genetic algorithm with crossover based on confidence interval as an alternative to traditional nonlinear regression methods
Yen et al. A method for automatic generation of a fuzzy model
CN112686337A (en) Cross-platform credit information fusion method based on depth-first algorithm
Johanyák Clonal selection based parameter optimization for sparse fuzzy systems

Legal Events

Date Code Title Description
E701 Decision to grant or registration of patent right
GRNT Written decision to grant