KR102102517B1 - 최적화 계산 장치 및 방법 - Google Patents

최적화 계산 장치 및 방법 Download PDF

Info

Publication number
KR102102517B1
KR102102517B1 KR1020180011562A KR20180011562A KR102102517B1 KR 102102517 B1 KR102102517 B1 KR 102102517B1 KR 1020180011562 A KR1020180011562 A KR 1020180011562A KR 20180011562 A KR20180011562 A KR 20180011562A KR 102102517 B1 KR102102517 B1 KR 102102517B1
Authority
KR
South Korea
Prior art keywords
solution
similarity
vector
similarity determination
optimization calculation
Prior art date
Application number
KR1020180011562A
Other languages
English (en)
Other versions
KR20190092135A (ko
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 KR1020180011562A priority Critical patent/KR102102517B1/ko
Priority to PCT/KR2018/011819 priority patent/WO2019151606A1/ko
Priority to US16/965,251 priority patent/US11915155B2/en
Publication of KR20190092135A publication Critical patent/KR20190092135A/ko
Application granted granted Critical
Publication of KR102102517B1 publication Critical patent/KR102102517B1/ko

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N5/00Computing arrangements using knowledge-based models
    • G06N5/04Inference or reasoning models
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/12Computing arrangements based on biological models using genetic models
    • G06N3/126Evolutionary algorithms, e.g. genetic algorithms or genetic programming
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N20/00Machine learning
    • 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
    • G06F17/18Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Biophysics (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Computing Systems (AREA)
  • Artificial Intelligence (AREA)
  • Mathematical Physics (AREA)
  • Data Mining & Analysis (AREA)
  • Evolutionary Computation (AREA)
  • General Physics & Mathematics (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Evolutionary Biology (AREA)
  • Computational Linguistics (AREA)
  • Molecular Biology (AREA)
  • General Health & Medical Sciences (AREA)
  • Physiology (AREA)
  • Biomedical Technology (AREA)
  • Genetics & Genomics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Medical Informatics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

최적화 계산 장치 및 방법이 제공된다. 최적화 계산 장치는 복수 개의 입력 데이터로부터 복수 개의 제1 해(S1)를 구하고, 상기 제1 해(S1)로부터 다시 제2 해(S2)를 구하며, 상기 과정을 반복 수행하여 최적 해(Sm)를 도출하는 알고리즘 모듈 및 상기 알고리즘 모듈에 연결되어 (i+1) 번째 해를 구하기 위해 i 번째 해의 유사도를 계산하는 유사도 판별 모듈을 포함할 수 있다.( 1 ≤ i ≤ m-1 )

Description

최적화 계산 장치 및 방법{APPARATUS AND METHOD OF OPTIMIZATION CALCULATION}
본 발명은 최적화 계산 장치 및 방법에 관한 것이다.
유전 알고리즘(Genetic Algorithm)은 자연계의 유전학에 바탕을 두며, 특히 다윈의 적자생존 이론을 기본 개념으로 한다. 유전자 프로그래밍에서는 문제에 대한 가능한 해(솔루션)들을 나열한 뒤, 유전자들을 변화시켜 정확도가 높은 해들을 만들어 낸다. 즉, 유전 알고리즘은 더 좋은 답을 찾아 가기 위해 진화를 모방한 탐색 알고리즘이라고 할 수 있다.
이처럼 유전 알고리즘은 견고함과 수렴성 때문에 다양한 최적화 문제에서 솔루션을 얻는데 널리 사용된다. 하지만, 휴리스틱 기반(heuristic-based) 알고리즘이기 때문에 특히 종료 단계에서 사용자의 임의적인 판단이 관련될 수 있는 가능성이 크고, 문제의 복잡성이 커짐에 따라 계산 시간이 급격히 늘어나는 단점을 가진다.
상기와 같은 문제점을 해결하기 위하여, 본 발명은 보다 효율적인 최적화 계산 장치를 제공한다.
상기와 같은 문제점을 해결하기 위하여, 본 발명은 상기 최적화 계산 장치를 이용한 최적화 계산 방법을 제공한다.
본 발명의 다른 목적들은 다음의 상세한 설명과 첨부한 도면으로부터 명확해 질 것이다.
본 발명의 실시예들에 따른 최적화 계산 장치는 복수 개의 입력 데이터로부터 복수 개의 제1 해(S1)를 구하고, 상기 제1 해(S1)로부터 다시 제2 해(S2)를 구하며, 상기 과정을 반복 수행하여 최적 해(Sm)를 도출하는 알고리즘 모듈 및 상기 알고리즘 모듈에 연결되어 (i+1) 번째 해를 구하기 위해 i 번째 해의 유사도를 계산하는 유사도 판별 모듈을 포함할 수 있다.( 1 ≤ i ≤ m-1 )
상기 알고리즘 모듈은 상기 유사도 판별 모듈에서 계산된 유사도에 따라 i 번째 해가 클러스터링되고, 상기 클러스터링된 i 번째 해에서 임의의 해가 선택되어 (i+1) 번째 해를 구할 수 있다.
상기 알고리즘 모듈은 상기 유사도의 평균, 표준편차, 신뢰도, 및 신뢰구간을 이용한 통계적 추정 방법으로 상기 클러스터링이 수행될 수 있다.
상기 알고리즘 모듈은 상기 신뢰도에 의해 상기 클러스터링이 제어될 수 있다.
상기 유사도 판별 모듈은 유사도 판별 함수를 이용하여 상기 유사도가 계산되고, 상기 유사도 판별 함수는 (i-1) 번째의 해 벡터의 무게 중심을 나타내는 기준 벡터(Oi), 상기 기준 벡터와 i 번째 해 벡터 사이의 반지름 (ri,jk, ri,jl), i 번째의 서로 다른 두 해 벡터 사이의 인접거리(di,jk,jl), 및 상기 두 해 벡터와 상기 기준 벡터를 기준으로 이루는 각(θi,jk,jl)을 이용할 수 있다.
상기 유사도 판별 함수는 상기 기준 벡터, 상기 반지름, 상기 인접거리, 상기 각을 이용한 다음 식
Figure 112018010716615-pat00001
일 수 있다.
본 발명의 실시예들에 따른 최적화 계산 방법은 복수 개의 입력 데이터로부터 복수 개의 제1 해(S1)를 구하는 제1 단계, 상기 제1 해(S1)를 이용하여 유사도를 계산하고 클러스터링 하는 제2 단계, 상기 클러스터링된 제1 해(S1)에서 임의의 해를 선택하여 제2 해(S2)를 구하는 제3 단계 및 상기 제2 단계 및 상기 제3 단계를 반복적으로 수행하여 최적 해(Sm)를 도출하는 제4 단계를 포함할 수 있다.
상기 제2 단계는 (i+1) 번째 해를 구하기 위해 i 번째 해의 유사도가 계산되고, 상기 유사도에 따라 i 번째 해가 클러스터링되며, 상기 클러스터링된 i 번째 해에서 임의의 해가 선택되어 (i+1) 번째 단계가 진행될 수 있다.( 1 ≤ i ≤ m-1 )
상기 제2 단계는 상기 유사도의 평균, 표준편차, 신뢰도, 및 신뢰구간을 이용한 통계적 추정 방법으로 상기 클러스터링이 수행될 수 있다.
상기 제2 단계는 상기 신뢰도에 의해 상기 클러스터링이 제어될 수 있다.
상기 제2 단계는 유사도 판별 함수를 이용하여 상기 유사도가 계산되고, 상기 유사도 판별 함수는 (i-1) 번째의 해 벡터의 무게 중심을 나타내는 기준 벡터(Oi), 상기 기준 벡터와 i 번째 해 벡터 사이의 반지름 (ri,jk, ri,jl), i 번째의 서로 다른 두 해 벡터 사이의 인접거리(di,jk,jl), 및 상기 두 해 벡터와 상기 기준 벡터를 기준으로 이루는 각(θi,jk,jl)을 이용할 수 있다.
상기 유사도 판별 함수는 상기 기준 벡터, 상기 반지름, 상기 인접거리, 상기 각을 이용한 다음 식
Figure 112018010716615-pat00002
일 수 있다.
본 발명의 실시예들에 따르면 계산 속도가 향상된 최적화 계산 장치를 제공할 수 있다.
본 발명의 실시예들에 따르면 수학적 신뢰도를 부여하여 결과적 신뢰도를 높힐 수 있는 최적화 계산 장치를 제공할 수 있다.
본 발명의 실시예들에 따르면 상기 장치를 이용한 최적화 계산 방법을 제공할 수 있다.
도 1은 본 발명의 일 실시예에 따른 흐름도이다.
도 2는 본 발명의 일 실시예에 따른 Z confidence 의 프로파일을 나타낸다.
도 3은 본 발명의 다른 실시예에 따른 알고리즘을 수행하기 위한 벡터값이다.
도 4는 도 3의 알고리즘을 수행한 결과를 나타낸다.
도 5는 본 발명의 일 실시예에 따른 실제 해 벡터와의 각을 나타낸다.
도 6은 본 발명의 일 실시예에 따른 실제 해 벡터와의 거리를 나타낸다.
도 7은 본 발명의 일 실시예에 따른 목적 함수 값을 나타낸다.
이하, 실시예들을 통하여 본 발명을 상세하게 설명한다. 본 발명의 목적, 특징, 장점은 이하의 실시예들을 통해 쉽게 이해될 것이다. 본 발명은 여기서 설명되는 실시예들에 한정되지 않고, 다른 형태로 구체화될 수도 있다. 여기서 소개되는 실시예들은 개시된 내용이 철저하고 완전해질 수 있도록 그리고 본 발명이 속하는 기술 분야에서 통상의 지식을 가진 자에게 본 발명의 사상이 충분히 전달될 수 있도록 하기 위해 제공되는 것이다. 따라서, 이하의 실시예들에 의하여 본 발명이 제한되어서는 안 된다.
본 명세서에서 제1, 제2 등의 용어가 다양한 요소들(elements)을 기술하기 위해서 사용되었지만, 상기 요소들이 이 같은 용어들에 의해서 한정되어서는 안 된다. 이러한 용어들은 단지 상기 요소들을 서로 구별시키기 위해서 사용되었을 뿐이다. 또, 어떤 요소가 다른 요소 위에 있다고 언급되는 경우에 그것은 다른 요소 위에 직접 형성될 수 있거나 또는 그들 사이에 제3의 요소가 개재될 수도 있다는 것을 의미한다.
도면들에서 요소의 크기, 또는 요소들 사이의 상대적인 크기는 본 발명에 대한 더욱 명확한 이해를 위해서 다소 과장되게 도시될 수 있다. 또, 도면들에 도시된 요소의 형상이 제조 공정상의 변이 등에 의해서 다소 변경될 수 있을 것이다. 따라서, 본 명세서에서 개시된 실시예들은 특별한 언급이 없는 한 도면에 도시된 형상으로 한정되어서는 안 되며, 어느 정도의 변형을 포함하는 것으로 이해되어야 한다.
본 발명에서 i는 입력 데이터로부터 해를 구하는 알고리즘 한 회분의 수행 과정으로서, 상기 알고리즘의 성격에 따라 단계 또는 세대로 표현될 수 있고, 상기 표현에 한정되지 않는다.
본 발명에서 사용되는 알고리즘과 관련된 실시에서는 네 가지 요건이 필요하다. 상기 요건은 유전 알고리즘 사용에 대한 기본 전제 조건에 해당될 수 있다.
첫째, 모든 해는 유전적 표현을 가진다. 상기 유전 알고리즘을 사용하기 위한 기본 요건이다.
둘째, 목적 함수는 해 도메인의 모든 값에 대해 발산하지 않는다. 상기 목적 함수가 상기 해 도메인의 일부 요소에 대해 발산하면 최적의 해를 얻는 것이 부적절하다.
셋째, 최적화 문제를 해결하는 데는 한정된 시간이 소요된다. 기존의 알고리즘과 본 발명에서 사용되는 유전 알고리즘의 효율성을 비교하기 위해 처리해야 할 최적화 문제는 한정된 시간 내에 해결되어야 한다.
넷째, 주어진 제약 조건 하에서 최적의 해가 존재한다. 상기 최적의 해가 주어진 범위 내에 있지 않으면, 상기 유전 알고리즘 사용의 필요성이 사라진다.
[최적화 계산 장치]
도 1은 본 발명의 일 실시예에 따른 흐름도이다.
도 1을 참조하면, 최적화 계산 장치는 알고리즘 모듈 및 유사도 판별 모듈을 포함할 수 있다.
상기 알고리즘 모듈은 복수 개의 입력 데이터로부터 복수 개의 제1 해(S1)를 구하고, 상기 제1 해(S1)로부터 다시 제2 해(S2)를 구하며, 상기 과정을 반복 수행하여 최적 해(Sm)를 도출할 수 있다. 상기 알고리즘 모듈은 복수의 개체를 진화시켜 세대를 거듭하며 해를 도출하는 유전 알고리즘이 사용될 수 있다. 상기 알고리즘 모듈은 모집단 초기화부, 평가부, 선택부, 교차 및 돌연변이부, 및 종결부를 포함할 수 있다.
모집단 초기화(Initialization of population)부는 모집단을 결정할 수 있다. 상기 모집단의 크기는 최적화할 문제의 성격을 고려하여 사용자에 의해 결정될 수 있다. 일반적으로 상기 모집단은 실현 가능한 해의 전체 범위에서 무작위로 생성될 수 있다.
평가(Evaluation)부는 생성된 상기 모집단에 대해 목적 함수의 값이 평가 된다. 상기 목적 함수는 다목적 함수가 될 수 있다.
선택(Selection)부는 상기 평가를 바탕으로, 엘리트 해 개체군이 새로운 세대를 낳기 위해 선정될 수 있다.
교차 및 돌연변이(Crossover and Mutation)부는 선택된 상기 해 개체군이 새로운 세대의 개체군을 생성하기 위해 교차 및 돌연변이법이 사용될 수 있다. 그 결과, 이전의 해 집단과 다른 새로운 해 집단이 생성될 수 있다.
종결(Termination)부는 최소 조건을 만족하는 해를 찾아 특정 수치의 생성에 도달하거나 사용자의 임의 조건을 만족하면, 상기 알고리즘이 종결될 수 있다.
상기 유사도 판별 모듈은 상기 알고리즘 모듈에 연결되어, (i+1) 번째 해를 구하기 위해 i 번째 해의 유사도를 계산할 수 있다. 상기 유사도 판별 모듈은 클러스터링(Clustering)을 위한 부목적 함수(Sub-objective function), 즉 유사도 판별 함수를 가질 수 있다. 상기 유사도 판별 함수는 각 세대에 대한 기준이 검토될 때 해를 획득하기 위해 요구되는 계산 시간이 쉽게 증가되기 때문에, 쉽게 계산될 수 있는 기본 함수를 사용할 수 있다.
상기 유사도 판별 함수를 만들기 위해서 기준 벡터, 반지름, 인접거리, 및 각 등의 수학 개념이 사용될 수 있다.
기준 벡터(Reference vector, O i )는 해 벡터 간의 유사성을 결정하는 벡터로 식 (1)과 같이 교차부와 돌연변이부의 각 (i-1) th 단계에서 얻은 해 벡터, s i-1,jk 의 평균값으로 정의될 수 있다. 즉, 세대가 갱신될 때, 이전 세대의 무게 중심(Center of gravity) 벡터로, i 번째 기준 벡터는 (i-1) 번째 세대의 해 벡터들의 무게 중심으로 표현될 수 있다. 초기 값은 사용자가 입력한 값으로 정의될 수 있다. 각 세대의 해 벡터의 개수 j k 의 값은 각 세대에 따라, 즉 i의 값에 따라 변할 수 있다.
Figure 112018010716615-pat00003
(1)
반지름(Radius, r i,jk )은 식 (2)와 같이 각 세대에서 생성된 모든 k 해 벡터와 i th 세대의 기준 벡터 사이의 유클리드 놈(Euclidean norm)으로 정의될 수 있다. 상기 반지름이 감소할수록 유사도는 주어진 기준 벡터에 대해 증가한다.
Figure 112018010716615-pat00004
(2)
인접 거리(Adjacent distance, d i,jk,jl )는 식 (3)과 같이 동일 세대의 서로 다른 두 해 벡터 s i,jk s i,jl 사이의 유클리드 놈으로 정의될 수 있다. 상기 인접 거리가 가까울수록 유사도가 높아진다.
Figure 112018010716615-pat00005
(3)
각(Angle)은 식 (4)와 같이 동일 세대의 서로 다른 두 해 벡터가 상기 기준 벡터 O i 를 기준으로 하여 이루는 각 θ i,jk,ji 로 정의될 수 있다. 두 벡터 사이의 각이 작을수록 유사도가 높아진다.
Figure 112018010716615-pat00006
(4)
상기 수학적 개념들을 이용하여 상기 유사도 판별 함수는 식 (5)와 같이 나타낼 수 있다.
Figure 112018010716615-pat00007
(5)
상기 유사도 판별 함수는 i th 세대에 속한 k 개의 해 벡터에 대하여 k C 2 번 계산될 수 있다. 만약 두 해 벡터가 같을 경우, 상기 함수 값은 0 이 된다.
상기 유사도 판별 모듈에는 드롭 아웃 정책(Drop out policy)이 포함될 수 있다. 상기 드롭 아웃 정책은 계산된 유사도 판별 함수를 기반으로 설정될 수 있다. J i,jk,jl ≤K 이면, s i,jk s i,jl 가 동일한 그룹에 속하는 것으로 취급되고, 그 중 하나는 드롭 아웃된다. K는 후술하는 통계적 추정에 의해 결정될 수 있다.
[통계적 추정(Statistical Estimation)]
통계적 추정은 표본 값으로부터 모집단 평균을 추정하는데 사용될 수 있다.표본 수가 증가함에 따라 신뢰도는 증가하지만 효율성은 감소할 수 있다.
식 (6),(7), 및 (8)은 평균 및 분산(Mean and Variance)을 나타낸다.
Figure 112018010716615-pat00008
(6)
Figure 112018010716615-pat00009
(7)
Figure 112018010716615-pat00010
(8)
상기 평균 및 분산을 기반으로 표준화 방법을 사용하면 간격을 설정(Interval generation)할 수 있다. 해가 특정한 수학적 신뢰도, Z confidence 를 가진 동일한 그룹에 속하는 것처럼 취급되는 가상 간격을 갖게 된다.
Figure 112018010716615-pat00011
(9)
여기서,
Figure 112018010716615-pat00012
이다. 예를 들어, 95% 신뢰도를 사용하면 Z confidence = 1.96이다. 여기서 n은 U variable 의 수이다.
K는 여러 가지 방법으로 정의될 수 있는데, 본 발명에서는 식 (10)과 같이 가장 간단한 형태로 정의될 수 있다.
Figure 112018010716615-pat00013
(10)
유사도 판별 함수의 통계적 추정을 기반으로 상기 해가 상기 가상 간격만큼 클러스터링될 수 있다. 상기 클러스터링된 해 그룹에서, 동일 그룹내에 속한 해들은 유사한 것으로 판단되어, 임의의 해 하나만을 추출하여 다음 세대로 계산을 진행시킬 수 있다. 상기 세대가 거듭될수록 상기 클러스터링되는 해의 그룹의 수는 줄어들 수 있고, 그에 따라 계산 속도도 빨라질 수 있다. 사용자가 설정한 일정 기준에 도달되면 상기 계산이 종결될 수 있다.
[Z confidence 의 프로파일]
도 2는 본 발명의 일 실시예에 따른 Z confidence 의 프로파일을 나타낸다. 새로운 세대가 만들어지면, Z confidence 값의 제어에 의해 효율적으로 추출된 해 벡터의 수를 줄일 수 있다. 시스템에 따라 다양한 Z confidence 의 프로파일이 가능해진다. 예를 들어, 단조 증가 함수 및 단조 감소 함수가 사용될 수 있다.
단조 증가 함수(Monotone increasing function)는 Z confidence 값이 계속 증가함에 따라 추출된 해의 수와 계산 시간을 줄일 수 있다. 세대에 따른 상기 Z confidence 의 감소 정도는 사용자의 필요에 따라 적절하게 조정될 수 있다. 상기 함수는 획득된 해 사이의 유사도가 높을 때 사용될 수 있다.
단조 감소 함수 (Monotone decreasing function)가 적용될 때, 추출된 해 벡터의 수는 세대에 걸쳐 증가될 수 있다. 모든 간격에 대해 상기 함수를 사용하면 관측 범위가 증가하는 동안 계산량이 증가한다는 단점이 있다. 상기 함수는 해 벡터의 수가 너무 작아 유전 알고리즘으로 진행할 수 없을 때 사용될 수 있다.
상기 획득된 해와 그로부터 계산된 목적 함수의 값 사이의 유사성에 기초하여 상기 Z confidence 의 프로파일이 유연하게 설정될 수 있다. 도 2를 참조하면, 본 발명에서는 지수 함수형의 단조 감소 함수가 사용될 수 있다.
[유사도 판별 함수의 적용 과정]
최적화 계산 장치를 이용한 계산 방법은 복수 개의 입력 데이터로부터 복수 개의 제1 해(S1)를 구하는 제1 단계, 상기 제1 해(S1)를 이용하여 유사도를 계산하고 클러스터링 하는 제2 단계, 상기 클러스터링된 제1 해(S1)에서 임의의 해를 선택하여 제2 해(S2)를 구하는 제3 단계 및 상기 제2 단계 및 상기 제3 단계를 반복적으로 수행하여 최적 해(Sm)를 도출하는 제4 단계를 포함할 수 있다.
상기 제2 단계는 (i+1) 번째 해를 구하기 위해 i 번째 해의 유사도가 계산되고 상기 유사도에 따라 i 번째 해가 클러스터링되며, 상기 클러스터링된 i 번째 해에서 임의의 해가 선택되어 (i+1) 번째 단계가 진행될 수 있다.( 1 ≤ i ≤ m-1 )
도 3은 본 발명의 다른 실시예에 따른 알고리즘을 수행하기 위한 벡터값들을 나타내고, 도 4는 도 3의 알고리즘을 수행한 결과를 나타낸다. 특정 세대에 속한 7차원 벡터 10개에 대해 본 발명의 알고리즘이 수행되었다.
도 3을 참조하면, 이전 세대의 무게 중심인 기준 벡터를 바탕으로 유사도 판별 함수를 계산하고, 상기 계산된 유사도 판별 함수값들의 평균과 표준 편차를 계산한다. 여기서 평균은 8.99, 표준 편차는 6.09로 계산된다.
통계적 추정에 사용할 신뢰도를 선택하고, 그에 따른 Z 값을 얻는다. 예를 들어, 95% 신뢰도일 경우, Z는 1.96, 99% 신뢰도일 경우, Z는 2.58 값을 가질 수 있다. 여기서는 95% 신뢰도를 선택하여 1.96의 Z 값을 가진다.
상기 평균, 상기 표준 편차, 상기 Z 값을 바탕으로 신뢰구간 L을 계산하고, 상기 유사도 판별 함수 값들 중 최대값과 최소값의 차이를 계산하여 범위를 구하고 R로 정의한다. 여기서 L은 3.77, R은 18.45 이다.
상기 L 값과 상기 R 값을 이용하여 식 (11)과 같이 그룹 수 N을 구할 수 있다.
Figure 112018010716615-pat00014
(11)
여기서 N은 5이다.
도 4를 참조하면, 상기 유사도 판별 함수 값들을 오름차순으로 정렬한 후 최소값으로부터 L을 더한 값 이하까지를 첫 번째 그룹으로 분할할 수 있다. 나머지 유사도 판별 함수값들에 대해서도 상기 분할 과정을 반복하면, 10개의 벡터가 5개의 그룹으로 분할될 수 있다.
상기 분할된 해 벡터를 바탕으로 유전 알고리즘의 다음 세대 계산을 수행할 수 있다. 상기 5개의 그룹에서 임의의 벡터 하나씩을 선택하여 앞선 과정들을 반복하여 수행할 수 있다. 상기 수행의 결과로 사용자가 원하는 신뢰도 하에서 사용자가 결정한 일정 그룹수 이하로 N값이 줄어들면 상기 유전 알고리즘은 종료될 수 있다.
[시뮬레이션]
기존의 유전 알고리즘과 본 발명에서 사용된 유전 알고리즘을 비교하기 위해 예시로 수동 선택적 촉매 환원 시스템(pSCR)에 상기 유전 알고리즘들을 적용했다. 상기 적용을 통해 매개 변수를 예측하여, 알려진 기존의 매개 변수 값과 본 발명의 유전 알고리즘을 이용한 결과 값을 비교할 수 있다. 실제 해 벡터와의 각, 거리, 및 목적 함수의 값 등 주요한 비교를 진행할 수 있다.
상기 목적 함수는 수동 선택적 촉매 환원 시스템에서 최소화되어야 하는 질소산화물(NOx)과 암모니아(NH3)의 누적 농도로 설정되고, 세대 수는 35로 고정되었다.
도 5는 본 발명의 일 실시예에 따른 실제 해 벡터와의 각을 나타내고, 도 6은 본 발명의 일 실시예에 따른 실제 해 벡터와의 거리를 나타내며, 도 7은 본 발명의 일 실시예에 따른 목적 함수 값을 나타낸다.
도 5 내지 도 7 에서, 적색은 기존의 유전 알고리즘을 이용한 것을 나타내고, 청색은 본 발명에서 사용된 유전 알고리즘을 이용한 것을 나타낸다. 그리고 본 발명의 유전 알고리즘은 27 세대 이후에 해를 생성하지 않는 것으로 나타난다.
도 5를 참조하면, 본 발명의 유전 알고리즘은 2 세대와 27 세대를 제외하고는 모든 부분에서 실제 해 벡터와의 각이 기존의 유전 알고리즘보다 훨씬 작은 값으로 나타난다. 즉, 본 발명의 유전 알고리즘이 기존 유전 알고리즘보다 기준 해 벡터와의 유사도가 더 높은 것으로 나타난다.
도 6을 참조하면, 본 발명의 유전 알고리즘이 1 세대부터 17 세대에서 기존의 유전 알고리즘보다 실제 해 벡터의 거리가 작게 나타나, 유사도가 높은 것으로 판단될 수 있고, 그 이후는 기존의 유전 알고리즘 보다 거리가 크게 나타나 유사도가 낮은 것으로 판단될 수 있다.
도 7을 참조하면, 상기 목적 함수는 기존 유전 알고리즘과 본 발명의 유전 알고리즘 방법 모두에서 감소하는 경향을 보이고, 15 세대 이후에는 유의한 변화를 보이지 않는 것으로 나타난다. 상기 15 세대 이후에는 기존의 유전 알고리즘으로부터 얻어진 해 벡터의 목적 함수값이 본 발명의 유전 알고리즘의 목적 함수의 값의 약 절반에 해당하는 것으로 나타난다.
본 발명에서 사용된 유전 알고리즘이 기존의 유전 알고리즘보다 각도와 거리에서 더 나은 성능을 보였음에도 불구하고 목적 함수의 계산된 값은 본 발명의 유전 알고리즘이 기존의 유전 알고리즘의 값보다 큰 것으로 나타났다. 그러나 이것은 현재 설정된 Z confidence 프로파일이 단조 증가 함수이므로, 상기 Z 프로파일이 제어되면 본 발명의 유전 알고리즘으로 27 세대 이후에 더 나은 결과를 얻을 수 있다.
게다가, 소요 시간을 비교하면, 기존 유전 알고리즘은 1.6시간/세대까지 소요되었지만, 본 발명의 유전 알고리즘은 0.57시간/세대가 소요되어, 기존 유전 알고리즘보다 3배 더 빠른 속도를 보였다. (계산은 Intel i5-4670 3.40 GHz 프로세서에서 수행되었다.) 따라서 요구되는 시간이 유전 알고리즘의 사용에 포함되면, 본 발명의 유전 알고리즘은 기존 유전 알고리즘에 비해 충분히 경쟁력이 있는 것으로 판단될 수 있다. 계산 시간이 긴 유전 알고리즘은 견고성을 확인하기가 어렵다. 본 발명의 유전 알고리즘은 동일한 시간 간격 동안 상기 유전 알고리즘을 여러 번 실행할 수 있으므로 견고성이 보장될 수 있다.
이제까지 본 발명에 대한 구체적인 실시예들을 살펴보았다. 본 발명이 속하는 기술 분야에서 통상의 지식을 가진 자는 본 발명이 본 발명의 본질적인 특성에서 벗어나지 않는 범위에서 변형된 형태로 구현될 수 있음을 이해할 수 있을 것이다. 그러므로 개시된 실시예들은 한정적인 관점이 아니라 설명적인 관점에서 고려되어야 한다. 본 발명의 범위는 전술한 설명이 아니라 특허청구범위에 나타나 있으며, 그와 동등한 범위 내에 있는 모든 차이점은 본 발명에 포함된 것으로 해석되어야 할 것이다.

Claims (12)

  1. 복수 개의 입력 데이터로부터 복수 개의 제1 해(S1)를 구하는 제1 과정과 상기 제1 해(S1)로부터 제2 해(S2)를 구하는 제2 과정을 수행하며, 상기 제2 과정을 반복 수행하여 최적 해(Sm)를 도출하는 알고리즘 모듈; 및
    상기 알고리즘 모듈에 연결되어, (i+1) 번째 해(1≤i≤m-1)를 구하기 위해 i 번째 해의 유사도를 계산하는 유사도 판별 모듈을 포함하고,
    상기 알고리즘 모듈은,
    상기 유사도 판별 모듈에서 계산된 유사도에 따라 i 번째 해를 클러스터링하고, 상기 클러스터링된 i 번째 해에서 임의의 해를 선택하여 (i+1) 번째 해를 구하는 것을 특징으로 하는 최적화 계산 장치.
  2. 삭제
  3. 제 1 항에 있어서,
    상기 알고리즘 모듈은,
    상기 유사도의 평균, 표준편차, 신뢰도, 및 신뢰구간을 이용한 통계적 추정 방법으로 상기 클러스터링을 수행하는 것을 특징으로 하는 최적화 계산 장치.
  4. 제 3 항에 있어서,
    상기 알고리즘 모듈은,
    상기 신뢰도에 의해 상기 클러스터링을 제어하는 것을 특징으로 하는 최적화 계산 장치.
  5. 복수 개의 입력 데이터로부터 복수 개의 제1 해(S1)를 구하는 제1 과정과 상기 제1 해(S1)로부터 제2 해(S2)를 구하는 제2 과정을 수행하며, 상기 제2 과정을 반복 수행하여 최적 해(Sm)를 도출하는 알고리즘 모듈; 및
    상기 알고리즘 모듈에 연결되어, (i+1) 번째 해(1≤i≤m-1)를 구하기 위해 i 번째 해의 유사도를 계산하는 유사도 판별 모듈을 포함하고,
    상기 유사도 판별 모듈은 유사도 판별 함수를 이용하여 상기 유사도를 계산하고,
    상기 유사도 판별 함수는,
    (i-1) 번째의 해 벡터의 무게 중심을 나타내는 기준 벡터(Oi), 상기 기준 벡터와 i 번째 해 벡터 사이의 반지름 (ri,jk, ri,jl), i 번째의 서로 다른 두 해 벡터 사이의 인접거리(di,jk,jl), 및 상기 두 해 벡터와 상기 기준 벡터를 기준으로 이루는 각(θi,jk,jl)을 이용하는 것을 특징으로 하는 최적화 계산 장치.
  6. 제 5 항에 있어서,
    상기 유사도 판별 함수는,
    상기 기준 벡터, 상기 반지름, 상기 인접거리, 상기 각을 이용한 다음 식
    Figure 112018010716615-pat00015
    인 것을 특징으로 하는 최적화 계산 장치.
  7. 제 1 항 또는 제 5 항의 최적화 계산 장치를 이용하여 최적화를 계산하는 방법으로서,
    복수 개의 입력 데이터로부터 복수 개의 제1 해(S1)를 구하는 제1 단계;
    상기 제1 해(S1)를 이용하여 유사도를 계산하고 클러스터링 하는 제2 단계;
    상기 클러스터링된 제1 해(S1)에서 임의의 해를 선택하여 제2 해(S2)를 구하는 제3 단계; 및
    상기 제2 단계 및 상기 제3 단계를 반복적으로 수행하여 최적 해(Sm)를 도출하는 제4 단계를 포함하고,
    상기 제2 단계에서,
    (i+1) 번째 해(1≤i≤m-1)를 구하기 위해 i 번째 해의 유사도가 계산되고,
    상기 유사도에 따라 i 번째 해가 클러스터링되며, 상기 클러스터링된 i 번째 해에서 임의의 해가 선택되어 (i+1) 번째 단계가 진행되는 것을 특징으로 하는 최적화 계산 방법.
  8. 삭제
  9. 제 7 항에 있어서,
    상기 제2 단계에서,
    상기 유사도의 평균, 표준편차, 신뢰도, 및 신뢰구간을 이용한 통계적 추정 방법으로 상기 클러스터링이 수행되는 것을 특징으로 하는 최적화 계산 방법.
  10. 제 9 항에 있어서,
    상기 제2 단계에서,
    상기 신뢰도에 의해 상기 클러스터링이 제어되는 것을 특징으로 하는 최적화 계산 방법.
  11. 제 7 항에 있어서,
    상기 제2 단계에서, 유사도 판별 함수를 이용하여 상기 유사도가 계산되고,
    상기 유사도 판별 함수는,
    (i-1) 번째의 해 벡터의 무게 중심을 나타내는 기준 벡터(Oi), 상기 기준 벡터와 i 번째 해 벡터 사이의 반지름 (ri,jk, ri,jl), i 번째의 서로 다른 두 해 벡터 사이의 인접거리(di,jk,jl), 및 상기 두 해 벡터와 상기 기준 벡터를 기준으로 이루는 각(θi,jk,jl)을 이용하는 것을 특징으로 하는 최적화 계산 방법.
  12. 제 11 항에 있어서,
    상기 유사도 판별 함수는,
    상기 기준 벡터, 상기 반지름, 상기 인접거리, 상기 각을 이용한 다음 식
    Figure 112018010716615-pat00016
    인 것을 특징으로 하는 최적화 계산 방법.
KR1020180011562A 2018-01-30 2018-01-30 최적화 계산 장치 및 방법 KR102102517B1 (ko)

Priority Applications (3)

Application Number Priority Date Filing Date Title
KR1020180011562A KR102102517B1 (ko) 2018-01-30 2018-01-30 최적화 계산 장치 및 방법
PCT/KR2018/011819 WO2019151606A1 (ko) 2018-01-30 2018-10-08 최적화 계산 장치 및 방법
US16/965,251 US11915155B2 (en) 2018-01-30 2018-10-08 Optimization calculation apparatus and optimization calculation method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
KR1020180011562A KR102102517B1 (ko) 2018-01-30 2018-01-30 최적화 계산 장치 및 방법

Publications (2)

Publication Number Publication Date
KR20190092135A KR20190092135A (ko) 2019-08-07
KR102102517B1 true KR102102517B1 (ko) 2020-04-21

Family

ID=67479374

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1020180011562A KR102102517B1 (ko) 2018-01-30 2018-01-30 최적화 계산 장치 및 방법

Country Status (3)

Country Link
US (1) US11915155B2 (ko)
KR (1) KR102102517B1 (ko)
WO (1) WO2019151606A1 (ko)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20230046245A (ko) 2021-09-29 2023-04-05 전북대학교산학협력단 세포 배양 최적 조건 검출 방법 및 이를 지원하는 전자 장치

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR102102517B1 (ko) * 2018-01-30 2020-04-21 서울대학교산학협력단 최적화 계산 장치 및 방법

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005202960A (ja) 2004-01-12 2005-07-28 Honda Research Inst Europe Gmbh 最適化方法及び最適化プログラム

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070220586A1 (en) 2006-03-01 2007-09-20 Norman Salazar Computing resource assignment method and apparatus using genetic algorithms
EP2207119A1 (en) 2009-01-06 2010-07-14 Koninklijke Philips Electronics N.V. Evolutionary clustering algorithm
KR101113903B1 (ko) * 2010-03-31 2012-03-06 (주)퀀텀바이오솔루션즈 양자역학 스코어링을 포함하는 단백질-리간드 도킹 구조의 예측방법
US8700548B2 (en) 2010-04-28 2014-04-15 Indian Statistical Institute Optimization technique using evolutionary algorithms
US9785472B2 (en) * 2010-06-11 2017-10-10 Hewlett Packard Enterprise Development Lp Computing cluster performance simulation using a genetic algorithm solution
KR101522306B1 (ko) * 2013-05-09 2015-05-26 서울대학교산학협력단 유사도 특성을 이용한 메타휴리스틱 알고리즘에 기반한 시스템 및 그 제어방법
KR102102517B1 (ko) * 2018-01-30 2020-04-21 서울대학교산학협력단 최적화 계산 장치 및 방법

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005202960A (ja) 2004-01-12 2005-07-28 Honda Research Inst Europe Gmbh 最適化方法及び最適化プログラム

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
Kalynmoy Deb 외 3명. A Fast and Elitist Multiobjective Genetic Algorithm: NSGA-II. 2002년

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20230046245A (ko) 2021-09-29 2023-04-05 전북대학교산학협력단 세포 배양 최적 조건 검출 방법 및 이를 지원하는 전자 장치

Also Published As

Publication number Publication date
WO2019151606A1 (ko) 2019-08-08
KR20190092135A (ko) 2019-08-07
US11915155B2 (en) 2024-02-27
US20210224667A1 (en) 2021-07-22

Similar Documents

Publication Publication Date Title
Prajwala A comparative study on decision tree and random forest using R tool
Ekárt et al. A metric for genetic programs and fitness sharing
Oldewage et al. The merits of velocity clamping particle swarm optimisation in high dimensional spaces
Neshat et al. A new cooperative algorithm based on PSO and k-means for data clustering
KR102102517B1 (ko) 최적화 계산 장치 및 방법
CN110991518A (zh) 一种基于进化多任务的两阶段特征选择方法及系统
KR101693405B1 (ko) Ldfgb 알고리즘을 이용한 비정상 침입 탐지 장치 및 방법
CN106156854A (zh) 一种基于dna编码的支持向量机参数预测方法
Sun et al. Hierarchical recognition of sparse patterns in large-scale simultaneous inference
Peignier et al. Data-driven gene regulatory networks inference based on classification algorithms
Windarto et al. A comparison of continuous genetic algorithm and particle swarm optimization in parameter estimation of Gompertz growth model
JP7440395B2 (ja) 最適解探索装置及び最適解探索プログラム
CN114386508A (zh) 一种基于样本点显著性的高斯过程回归序贯设计方法
Teymourifar et al. Dynamic priority rule selection for solving multi-objective job shop scheduling problems
Girsang et al. Robust adaptive genetic K-means algorithm using greedy selection for clustering
JP2021131835A (ja) 画像処理システム、及び画像処理プログラム
Li et al. Extended rule-based genetic network programming
Khoshdel et al. A new hybrid learning-based algorithm for data clustering
Alibabaie et al. A variant of genetic algorithm for non-homogeneous population
Ripon et al. Bi-level multi-objective image segmentation using texture-based color features
Oliva et al. Distance-constrained data clustering by combined k-means algorithms and opinion dynamics filters
Joo et al. Fuzzy systems modeling: an introduction
EP3852024A1 (en) Information processing apparatus, information processing method, and information processing program
Nguyen et al. GMKIT2-FCM: a genetic-based improved multiple kernel interval type-2 fuzzy C-means clustering
CN114821206B (zh) 基于对抗互补特征的多模态图像融合分类方法与系统

Legal Events

Date Code Title Description
A201 Request for examination
E902 Notification of reason for refusal
E701 Decision to grant or registration of patent right
GRNT Written decision to grant