KR101948417B1 - Multiple neuro-physiological signal based animal behavior control system and method thereof - Google Patents

Multiple neuro-physiological signal based animal behavior control system and method thereof Download PDF

Info

Publication number
KR101948417B1
KR101948417B1 KR1020170039216A KR20170039216A KR101948417B1 KR 101948417 B1 KR101948417 B1 KR 101948417B1 KR 1020170039216 A KR1020170039216 A KR 1020170039216A KR 20170039216 A KR20170039216 A KR 20170039216A KR 101948417 B1 KR101948417 B1 KR 101948417B1
Authority
KR
South Korea
Prior art keywords
signal
animal
equation
quot
vector
Prior art date
Application number
KR1020170039216A
Other languages
Korean (ko)
Other versions
KR20180109451A (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 KR1020170039216A priority Critical patent/KR101948417B1/en
Publication of KR20180109451A publication Critical patent/KR20180109451A/en
Application granted granted Critical
Publication of KR101948417B1 publication Critical patent/KR101948417B1/en

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A01AGRICULTURE; FORESTRY; ANIMAL HUSBANDRY; HUNTING; TRAPPING; FISHING
    • A01KANIMAL HUSBANDRY; AVICULTURE; APICULTURE; PISCICULTURE; FISHING; REARING OR BREEDING ANIMALS, NOT OTHERWISE PROVIDED FOR; NEW BREEDS OF ANIMALS
    • A01K15/00Devices for taming animals, e.g. nose-rings or hobbles; Devices for overturning animals in general; Training or exercising equipment; Covering boxes
    • A01K15/02Training or exercising equipment, e.g. mazes or labyrinths for animals ; Electric shock devices ; Toys specially adapted for animals
    • A01K15/029Electric or similar shock devices, e.g. prods
    • AHUMAN NECESSITIES
    • A01AGRICULTURE; FORESTRY; ANIMAL HUSBANDRY; HUNTING; TRAPPING; FISHING
    • A01KANIMAL HUSBANDRY; AVICULTURE; APICULTURE; PISCICULTURE; FISHING; REARING OR BREEDING ANIMALS, NOT OTHERWISE PROVIDED FOR; NEW BREEDS OF ANIMALS
    • A01K29/00Other apparatus for animal husbandry
    • A01K29/005Monitoring or measuring activity, e.g. detecting heat or mating
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/02Agriculture; Fishing; Forestry; Mining

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Environmental Sciences (AREA)
  • Animal Husbandry (AREA)
  • Engineering & Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Biodiversity & Conservation Biology (AREA)
  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Marine Sciences & Fisheries (AREA)
  • Marketing (AREA)
  • Animal Behavior & Ethology (AREA)
  • Agronomy & Crop Science (AREA)
  • Physical Education & Sports Medicine (AREA)
  • Mining & Mineral Resources (AREA)
  • Biophysics (AREA)
  • Economics (AREA)
  • Human Resources & Organizations (AREA)
  • Zoology (AREA)
  • Primary Health Care (AREA)
  • Strategic Management (AREA)
  • Tourism & Hospitality (AREA)
  • Physics & Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)

Abstract

본 발명은 동물 행동 제어를 위한 시스템 및 방법에 관한 것으로, 본 발명의 따른 동물 행동 제어 시스템은 동물의 뇌신경신호 및 생리신호를 감지 및 기록하는 신호 감지 및 기록부; 상기 신호 감지 및 기록부로부터 전달받은 뇌신경신호 및 생리신호를 통해 각각의 뇌신경신호 특징벡터 및 생리신호 특징벡터를 추출하고, 추출된 특징벡터를 해독하여 동물의 운동 의도 및 생체 상태를 파악하는 신호 분석부; 파악된 동물의 운동 의도 및 생체 상태에 따라 전기 자극 파라미터를 조절하여 그에 상응하는 보상 혹은 처벌 자극을 주는 제어부; 및 조절된 전기 자극 파라미터를 이용하여 특정 대뇌 영역에 전기 자극을 출력하는 행동제어 명령 지시부;를 포함한다. 이를 이용하여 동물의 다중 신경생리신호를 기록 및 해독하여 동물의 운동 의도 및 생체 상태를 파악할 수 있고, 파악된 동물의 운동 의도 및 생체 상태를 통해 동물에게 적절한 전기 자극을 가해줌으로써 동물의 움직임을 효율적으로 제어하고 생체 상태를 조절할 수 있는 효과가 있다.The present invention relates to a system and a method for controlling animal behavior, and an animal behavior control system according to the present invention comprises a signal sensing and recording unit for sensing and recording a neural signal and a physiological signal of an animal; A signal analyzing unit for extracting each of the cranial nerve signal feature vector and physiological signal feature vector from the signal sensing and recording unit and receiving the cranial nerve signal and the physiological signal, ; A control unit for controlling electric stimulation parameters according to a motion intention and a living state of the identified animal to give corresponding compensation or punishment stimulation; And a behavior control instruction instructing unit to output electric stimulation to a specific cerebral region using the adjusted electric stimulation parameter. By using this, it is possible to record and decode animal's multi-nervous physiological signal to understand the motion intention and living state of the animal, and to apply the appropriate electric stimulation to the animal through the motion intention and the living state of the animal, And the living body state can be controlled.

Figure R1020170039216
Figure R1020170039216

Description

다중 신경생리신호 기반의 동물 행동 제어 시스템 및 그 방법{MULTIPLE NEURO-PHYSIOLOGICAL SIGNAL BASED ANIMAL BEHAVIOR CONTROL SYSTEM AND METHOD THEREOF}BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to an animal behavior control system based on multiple neural physiological signals,

본 발명은 동물의 행동 제어를 위한 시스템에 관한 것으로, 보다 구체적으로는 다중 신경생리신호의 기록 및 해독을 통해 동물의 생체의 상태와 행동 상태를 확인할 수 있게 하여 보다 효과적으로 동물의 행동을 제어할 수 있는 동물 행동 제어 시스템 및 이를 이용한 동물 행동의 제어 방법에 관한 것이다.The present invention relates to a system for controlling the behavior of an animal, and more particularly, by recording and deciphering multiple neuronal physiological signals, it is possible to check the state and behavior of an animal's living body, And an animal behavior control system using the same.

생체는 주변의 환경 변화에 적응하고 몸의 항상성을 유지하기 위하여 몸에 적절한 변화를 일으킨다. 또한, 인지 활동을 포함한 운동을 실행하기 위해 생체는 주변의 변화를 감지하고 움직임에 대한 계획을 설립하고 계획을 실천하는 복잡한 생체 활동을 일으킨다. 그리고 그 과정에서 다양한 종류의 생체 신호를 획득하고 해석할 수 있으며 이를 통해 생체가 현재 어떤 일을 진행하고 있는지에 대해 간접적으로 확인할 수도 있다.The living body adapts itself to changes in the surrounding environment and causes appropriate changes in the body to maintain the body's homeostasis. In addition, in order to carry out the exercise including the cognitive activity, the living body senses the change around, establishes a plan of movement, and creates a complex bio-activity that implements the plan. In the process, various kinds of bio-signals can be acquired and analyzed, and indirectly checking what the living body is currently doing.

한편, 전 세계적으로 동물을 훈련하여 특수 목적에 사용하는 많은 케이스가 있고 국내의 경우에도 인명구조견, 마약탐지견 등 동물을 이용 및 제어하여 특수 목적에 사용하고 있다. 하지만, 동물을 훈련하는 데 있어도 동물 역시 생물인 한계로 강제적인 명령에 대한 복종은 불가능하고, 또 동물의 생체 상태를 수치적으로 확인하고 적절한 자극을 이용해 생체 상태를 완화시켜 임무 완수에 도움을 주기가 힘들다는 한계점이 분명하다.On the other hand, there are many cases in which animals are trained and used for special purposes all over the world, and even in Korea, animals are used and controlled for special purposes such as rescue dogs and drug detection dogs. However, even in the training of animals, animals are also limited to living creatures, and it is impossible to obey mandatory orders, and numerical confirmation of the animal's state of life and mitigation of vital conditions by using appropriate stimuli It is clear that the limit is difficult.

한국에 출원한 특허출원번호 제10-2015-0022655호에서는 마이크, 동물의 활동량, 체온 등을 제어 단말로 전송해 동물의 상태를 알 수 있게 하고 사용자의 목소리 등을 전달할 수 있으나 실질적인 자극을 통한 강제적 명령 복종 또는 생체 상태의 완화를 시킬 수 있게 하는 것은 아니다.Patent Application No. 10-2015-0022655 filed in Korea can transmit the microphone, the activity amount of the animal, and the body temperature to the control terminal so as to know the state of the animal and to transmit the voice of the user. However, It does not allow instruction obedience or relaxation of the living state.

또한, 한국에 출원한 특허출원번호 제10-2011-0033068호에서는 보상 및 처벌 자극을 이용해서 동물의 움직임을 유도하는 시스템이긴 하나 피드백의 방식이 수동 형식이고 가상 환경 장치 내에서 훈련하고 적용할 수 있는 한계점이 있는 문제점이 있다.Patent Application No. 10-2011-0033068 filed in Korea also discloses a system for inducing animal movement using compensation and punishment stimulation, but the feedback method is a manual type and can be trained and applied in a virtual environment apparatus There is a problem with the limit.

한국 특허출원번호 제10-2015-0022655호Korean Patent Application No. 10-2015-0022655 한국 특허출원번호 제10-2011-0033068호Korean Patent Application No. 10-2011-0033068

이에 상기와 같은 점을 감안한 본 발명은 동물의 뇌신경신호 및 다중신경생리신호의 기록(recording) 및 해독(decoding) 기술을 이용함으로써 동물의 운동 의도, 생체 상태를 파악하고, 동물에게 적절한 전기 전극을 가해줌으로써 동물 움직임을 효과적 제어하고 생체 상태를 조절할 수 있는 동물 행동 제어 시스템 및 방법의 제공에 목적이 있다.In view of the above, the present invention provides a method of recording an animal's movement intention and a living state by using a recording and decoding technique of an animal's brain signal and multiple nerve physiological signals, And an animal behavior control system and method capable of effectively controlling an animal movement and controlling a living body by applying the animal behavior control system.

상기와 같은 목적을 달성하기 위해 본 발명의 일 양태에 따른 다중 신경생리신호 기록 및 해석 기반의 동물행동 제어 시스템은, 동물의 뇌신경신호 및 생리신호를 감지 및 기록하는 신호 감지 및 기록부(100), 상기 신호 감지 및 기록부로부터 전달받은 뇌신경신호 및 생리신호를 통해 각각의 뇌신경신호 특징벡터 및 생리신호 특징벡터를 추출하고, 추출된 뇌신경신호 특징벡터 및 생리신호 특징벡터를 해독하여 동물의 운동 의도 및 생체 상태를 파악하는 신호 분석부(110), 파악된 동물의 운동 의도 및 생체 상태에 따라 전기 자극 파라미터를 조절하여 그에 상응하는 보상 혹은 처벌 자극을 주는 제어부(120) 및 조절된 전기 자극 파라미터를 이용하여 특정 대뇌 영역에 전기 자극을 출력하는 행동제어 명령 지시부(130)를 포함하여 구성될 수 있다.According to an aspect of the present invention, there is provided an animal behavior control system for recording and analyzing multiple neural physiological signals, comprising: a signal sensing and recording unit (100) for sensing and recording a neural signal and a physiological signal of an animal; And extracts the cranial nerve signal feature vector and the physiological signal feature vector through the cranial nerve signal and the physiological signal received from the signal sensing and recording unit, decodes the extracted cranial nerve signal feature vector and physiological signal feature vector, A control unit 120 for adjusting electrical stimulation parameters according to the intention of the animal and physical conditions of the animal, and providing a corresponding compensation or punishment stimulus, And a behavior control command instruction unit 130 for outputting electric stimulation to a specific cerebral region.

신호 감지 및 기록부(100)에서 측정되는 뇌신경신호는 신경세포로부터 발생하는 활동전위인 세포 외 기록 신호(extracellular recording signal), 국소뇌전위(local field potential), 뇌피질 뇌파(subdural electrocortigram, ECoG) 및 뇌전도(electroencephalography, EEG) 중에서 선택된 어느 하나 이상이다.The cranial nerve signal measured in the signal sensing and recording unit 100 is an extracellular recording signal, a local field potential, a subdural electrocortigram (ECoG), and a cerebral cortex And electroencephalography (EEG).

상기 생리신호는 피부표면온도(skin temperature, SKT), 광혈류량(photo-plethysmography, PPG), 심전도(electrocardiography, ECG), 피부 전도도 반응(skin conductance response, SCR) 및 호흡 신호 중에서 선택되는 어느 하나 이상이다.The physiological signal may be at least one selected from skin temperature (SKT), photo-plethysmography (PPG), electrocardiography (ECG), skin conductance response to be.

상기 신호 감지 및 기록부(100)는 동물의 뇌신경신호를 감지 및 기록하는 뇌신경신호 기록 전극(100a), 동물의 생리신호를 감지 및 기록하는 생리신호 기록 전극(100b), 및 감지된 뇌신경신호 및 생리신호를 증폭하는 증폭부(100c)를 포함한다.The signal sensing and recording unit 100 includes a cranial nerve signal recording electrode 100a for sensing and recording a cranial nerve signal of an animal, a physiological signal recording electrode 100b for sensing and recording a physiological signal of an animal, And an amplification section 100c for amplifying the signal.

상기 신호 분석부(110)는 뇌신경신호로부터 뇌신경신호 특징벡터를 추출하는 뇌신경신호 특징벡터 추출부(110a-1)와, 생리신호로부터 생리신호 특징벡터를 추출하는 생리신호 특징벡터 추출부(110a-2)로 이루어진 특징벡터 추출부(110a), 및 상기 뇌신경신호 특징벡터와 상기 생리신호 특징벡터를 해독하여 동물의 운동 의도와 생체 상태를 파악하는 해독부(110b)를 포함한다.The signal analyzer 110 includes a cranial nerve signal feature vector extracting unit 110a-1 for extracting a cranial nerve signal feature vector from a cranial nerve signal, a physiological signal feature vector extracting unit 110a-c for extracting a physiological signal feature vector from the physiological signal, 2), and a decryption unit 110b for deciphering the cranial nerve signal feature vector and the physiological signal feature vector to determine the motion intention and the biological state of the animal.

상기 해독부(110b)는 서포트 벡터 머신(Support Vector Machine,이하 'SVM'라고도 함)에 의해 동물의 운동 의도를 이진 값으로 출력하는 이진 결정 의도 인식기(110b-1), 및 서포트 벡터 회귀(Support Vector Regression, 이하 'SVR'라고도함)에 의해 동물의 생체 상태를 연속된 실수 값으로 출력하는 회귀부(110b-2)를 포함한다.The decryption unit 110b includes a binary decision intention recognizer 110b-1 for outputting the motion intention of the animal as a binary value by a support vector machine (SVM), and a support vector regression And a regression unit 110b-2 that outputs the biological state of the animal as a continuous real value by Vector Regression (SVR).

상기 제어부(120)는 파악된 동물의 운동 의도 및 생체 상태를 토대로 보상 혹은 처벌 여부를 결정하는 보상 혹은 처벌 결정부(120a), 및 상기 보상 혹은 처벌 결정부에서 수신받은 결정으로부터 동물의 행동을 제어하기 위한 전기 자극 파라미터를 생성하는 자극 생성 및 제시부(120b)를 포함한다.The control unit 120 includes a compensation or punishment decision unit 120a for determining whether compensation or punishment is based on the detected motion intention and the living state of the animal, and a control unit for controlling the behavior of the animal from the decision received from the compensation or punishment decision unit And a stimulus generation and presentation unit 120b for generating an electric stimulation parameter for performing the stimulation.

상기 행동제어 명령 지시부(130)는 동물의 행동 및 생체 상태의 제어에 대한 사용자 명령을 입력받는 입력부(130a), 상기 신호 분석부(110)로부터 수신된 운동 의도와 생체 상태를 디스플레이부로 출력하는 출력부(130b), 상기 출력부에서 출력된 운동 의도와 생체 상태를 표시하는 디스플레이부(130c), 및 상기 입력부로부터 전달받은 사용자 명령을 토대로 움직임 방향에 관한 큐(cue) 자극 명령을 받아 행동 제어 명령과 생체 상태를 조절할 수 있도록 동물의 특정 대뇌 영역에 전기 자극을 출력하는 명령 자극 제시부(130d)를 포함한다.The behavior control instruction instructing unit 130 includes an input unit 130a for receiving a user command for controlling behavior and anatomy of an animal, an output unit 130a for outputting motion intention and biological condition received from the signal analyzing unit 110 to the display unit, A display unit 130c for displaying the motion intention and the living body state outputted from the output unit, and a cue stimulation command for the direction of movement based on the user command received from the input unit, And a command stimulus presentation unit 130d for outputting an electric stimulation to a specific cerebral region of the animal so as to control the living body state.

또 다른 본 발명의 목적을 달성하기 위한 본 발명의 동물 행동 제어 시스템을 이용한 동물 행동의 제어 방법은 도 8을 참조하면, (a) 대상 동물로부터 뇌신경신호 및 생리신호를 측정 및 저장하는 단계(S110), (b) 뇌신경신호 및 생리신호를 분석하여 각각의 뇌신경신호 특징벡터 및 생리신호 특징벡터를 추출하는 단계(S120), (c) 추출된 뇌신경신호 특징벡터 및 생리신호 특징벡터를 해독하여 동물의 운동 의도 및 생체 상태를 파악하는 단계(S130), (d) 파악된 동물의 운동 의도 및 생체 상태를 기반으로 보상 혹은 처벌여부를 결정하거나 동물의 임무지속 가능여부를 결정하여 동물의 행동을 제어하기 위한 전기 자극 파라미터를 조절하는 단계(S140), (e) 결정된 전기 자극 파라미터에 따른 전기 자극 파형을 특정 대뇌 영역에 전달하여 동물의 행동을 제어하는 단계(S150)를 포함하여 이루어진다.According to another aspect of the present invention, there is provided a method of controlling animal behavior using an animal behavior control system, comprising the steps of: (a) measuring and storing a neural signal and a physiological signal from a subject animal; (b) analyzing a cranial nerve signal and a physiological signal to extract a cranial nerve signal characteristic vector and a physiological signal characteristic vector, respectively, (c) extracting the extracted cranial nerve signal characteristic vector and physiological signal characteristic vector, (D) determining whether the animal is to be compensated or punished based on the intention and the biological state of the animal, or determining whether the animal can sustain the animal's behavior, (S140) of controlling the electric stimulation parameters for the animal, (e) controlling the behavior of the animal by transmitting the electric stimulation waveform according to the determined electric stimulation parameter to a specific cerebral region It comprises a (S150).

상기 (a) 단계(S110)는 신호 감지 및 기록부(100)를 통해 동물의 뇌신경신호 및 생리신호를 측정하고 저장할 수 있다.The step (a) S110 may measure and store the cerebral nerve signal and the physiological signal of the animal through the signal sensing and recording unit 100.

여기서 측정되는 동물의 뇌신경신호 및 생리신호는 앞서 설명한 동물행동 제어 시스템에서 측정된 뇌신경신호 및 생리신호와 동일하므로, 중복되는 설명은 생략하기로 한다.The cranial nerve signal and the physiological signal of the animal measured here are the same as the cranial nerve signal and the physiological signal measured in the animal behavior control system described above, and therefore, a duplicate description will be omitted.

상기 (b) 단계는 상기 (a) 단계 측정된 뇌신경신호 및 생리신호를 신호 분석부(110)의 특징벡터 추출부(110a)에서 각각의 뇌신경신호 및 생리신호 특징벡터를 추출하는 단계이다.The step (b) extracts the cranial nerve signals and the physiological signals from the feature vector extraction unit 110a of the signal analysis unit 110 in step (a).

구체적으로 (b) 단계는 상기 신호 감지 및 기록부(100)에서 측정된 뇌신경신호를 이용하여 특징벡터 추출부(110a)의 뇌신경신호 특징벡터 추출부(110a-1)에서 스파이크 트레인(spike train), 파워 스펙트럼 밀도(power spectrum density), 대뇌 영역 간 기능적 연결도(functional connectivity), 및 대뇌 영역 간 효과적 연결도(effective connectivity) 중의 어느 하나 이상을 통해 뇌신경신호 특징벡터를 추출할 수 있다.More specifically, in step (b), a cranial nerve signal characteristic vector extraction unit 110a-1 of the feature vector extraction unit 110a extracts a spike train, It is possible to extract the cranial nerve signal characteristic vector through at least one of power spectrum density, functional connectivity between cerebral regions, and effective connectivity between cerebral regions.

상기 파워 스펙트럼 밀도(power spectrum density)는 측정된 뇌신경신호 중 어느 하나의 신호를 M차원(M은 1 이상의 정수)의 벡터

Figure 112017030378610-pat00001
로 하여, 하기 수학식 1을 통해 분석하는 이산 푸리에 변환(discrete Fourier transform), 하기 수학식 2를 통해 분석하는 단기간 푸리에 변환(short-time Fourier transform) 및 하기 수학식 3을 통해 분석하는 웨이브릿 변환(wavelet transform) 중의 어느 하나를 이용하여 특정 시간-주파수 대역의 파워를 얻어 뇌신경신호 특징벡터를 추출할 수 있다.The power spectral density is obtained by dividing one of the measured nerve signals into M dimensions (M is an integer of 1 or more)
Figure 112017030378610-pat00001
A discrete Fourier transform for analysis through the following equation (1), a short-time Fourier transform for analysis through the following equation (2) and a wavelet transform for analysis through the following equation (3) and a wavelet transform to obtain a power of a specific time-frequency band and to extract a cranial nerve signal feature vector.

[수학식 1][Equation 1]

Figure 112017030378610-pat00002
Figure 112017030378610-pat00002

수학식 1에서

Figure 112017030378610-pat00003
는 주파수를 지칭하는 변수
Figure 112017030378610-pat00004
에 대하여
Figure 112017030378610-pat00005
를 의미하는 각주파수를 나타낸다.In Equation (1)
Figure 112017030378610-pat00003
Is a variable that refers to frequency
Figure 112017030378610-pat00004
about
Figure 112017030378610-pat00005
≪ / RTI >

[수학식 2]&Quot; (2) "

Figure 112017030378610-pat00006
Figure 112017030378610-pat00006

수학식 2에서

Figure 112017030378610-pat00007
은 특정한 길이의 특정한 모양의 시간 윈도우 벡터를 나타낸다.In Equation 2,
Figure 112017030378610-pat00007
Represents a time window vector of a particular shape of a specific length.

[수학식 3]&Quot; (3) "

Figure 112017030378610-pat00008
Figure 112017030378610-pat00008

수학식 3에서

Figure 112017030378610-pat00009
는 웨이블릿 함수의 공액복소수(complex conjugate),
Figure 112017030378610-pat00010
는 정규화 파라미터,
Figure 112017030378610-pat00011
는 특정 스케일,
Figure 112017030378610-pat00012
는 목표 레이턴시를 나타낸다.In Equation 3,
Figure 112017030378610-pat00009
Is a complex conjugate of the wavelet function,
Figure 112017030378610-pat00010
Is a normalization parameter,
Figure 112017030378610-pat00011
Lt; / RTI >
Figure 112017030378610-pat00012
Represents the target latency.

상기 대뇌 영역 간 기능적 연결도(functional connectivity)는 뇌신경신호 중 어느 둘의 뇌신경신호를 한 쌍의 M차원(M은 1 이상의 정수)의 벡터

Figure 112017030378610-pat00013
,
Figure 112017030378610-pat00014
로 하여, 하기 수학식 4를 통해 분석하는 간섭성(coherence), 하기 수학식 5를 통해 분석하는 위상 동기 값(phase locking value, PLV), 및 하기 수학식 6을 통해 분석하는 weighted phase lag index(wPLI) 중의 어느 하나를 이용하여 한 쌍의 뇌신경 영역 간의 기능적 연결성을 확인하여 뇌신경신호 특징벡터를 추출할 수 있다.The functional connectivity between the cerebral regions is determined by connecting two neurons of the cerebral nerve signal to a pair of M dimensions (M is an integer of 1 or more)
Figure 112017030378610-pat00013
,
Figure 112017030378610-pat00014
, A phase locking value (PLV) analyzed through Equation (5), and a weighted phase lag index (PLV) analyzed through Equation (6) wPLI) to identify the functional connectivity between a pair of neuronal regions to extract the neural signal characteristic vector.

[수학식 4]&Quot; (4) "

Figure 112017030378610-pat00015
Figure 112017030378610-pat00015

수학식 4에서

Figure 112017030378610-pat00016
Figure 112017030378610-pat00017
Figure 112017030378610-pat00018
Figure 112017030378610-pat00019
의 진폭,
Figure 112017030378610-pat00020
는 두 시계열의 상대적 위상차를 나타낸다.In Equation 4,
Figure 112017030378610-pat00016
Wow
Figure 112017030378610-pat00017
The
Figure 112017030378610-pat00018
Wow
Figure 112017030378610-pat00019
The amplitude,
Figure 112017030378610-pat00020
Represents the relative phase difference between two time series.

[수학식 5]&Quot; (5) "

Figure 112017030378610-pat00021
Figure 112017030378610-pat00021

수학식 5에서,

Figure 112017030378610-pat00022
Figure 112017030378610-pat00023
는 각각 시계열 벡터
Figure 112017030378610-pat00024
Figure 112017030378610-pat00025
에 힐버트 변환을 적용하여 얻은 순간 위상(instantaneous phase)을 나타낸다.In Equation (5)
Figure 112017030378610-pat00022
and
Figure 112017030378610-pat00023
Respectively,
Figure 112017030378610-pat00024
and
Figure 112017030378610-pat00025
And shows the instantaneous phase obtained by applying the Hilbert transform to the signal.

[수학식 6]&Quot; (6) "

Figure 112017030378610-pat00026
Figure 112017030378610-pat00026

수학식 6에서,

Figure 112017030378610-pat00027
은 기대값 연산,
Figure 112017030378610-pat00028
은 복소수의 허수 항만을 취하는 연산자,
Figure 112017030378610-pat00029
는 두 시계열 벡터
Figure 112017030378610-pat00030
Figure 112017030378610-pat00031
의 교차 스펙트럼(cross spectrum)을 나타낸다.In Equation (6)
Figure 112017030378610-pat00027
The expected value calculation,
Figure 112017030378610-pat00028
Is an operator that takes an imaginary number of complex numbers,
Figure 112017030378610-pat00029
The two time series vectors
Figure 112017030378610-pat00030
and
Figure 112017030378610-pat00031
Of the cross spectrum.

상기 대뇌 영역 간 효과적 연결도(effective connectivity)는 뇌신경신호 중 어느 둘의 뇌신경신호를 한 쌍의 M차원(M은 1 이상의 정수)의 벡터

Figure 112017030378610-pat00032
,
Figure 112017030378610-pat00033
로 하여, 하기 수학식 7을 통해 분석하는 그랜저 인과관계(Granger's causality), 하기 수학식 8을 통해 분석하는 PDC(Partial Directed Coherence), 및 하기 수학식 9를 통해 분석하는 DTF(Directed Transfer Function) 중 어느 하나 이상을 이용하여 한 쌍의 뇌신경 영역 간의 방향성이 있는 연결성을 확인하여 뇌신경신호 특징벡터로 대뇌 영역 간 효과적 연결도(effective connectivity)를 추출할 수 있다.The effective connectivity between the cerebral regions is obtained by connecting two neurons of the cerebral nerve signal to a pair of M dimensions (M is an integer of 1 or more)
Figure 112017030378610-pat00032
,
Figure 112017030378610-pat00033
(Granger's causality) analyzed by the following equation (7), PDC (Partial Directed Coherence) analyzed through the following equation (8), and DTF (Directed Transfer Function) analyzed through the following equation Using one or more of the above, effective connectivity between the cerebral regions can be extracted with the neural signal characteristic vector by confirming the directional connectivity between the pair of neuronal regions.

[수학식 7]&Quot; (7) "

Figure 112017030378610-pat00034
Figure 112017030378610-pat00034

수학식 7에서,

Figure 112017030378610-pat00035
Figure 112017030378610-pat00036
Figure 112017030378610-pat00037
의 과거 값만을 이용하여 자동 회귀 모델을 이용해 추정할 때 생기는 에러,
Figure 112017030378610-pat00038
Figure 112017030378610-pat00039
Figure 112017030378610-pat00040
Figure 112017030378610-pat00041
의 과거값을 모두 사용하여 자동 회귀 모델을 이용해
Figure 112017030378610-pat00042
Figure 112017030378610-pat00043
의 현재 값을 추정할 때 생기는 에러이다.In Equation (7)
Figure 112017030378610-pat00035
silver
Figure 112017030378610-pat00036
of
Figure 112017030378610-pat00037
The error of the estimation using the automatic regression model using only the past value of the error,
Figure 112017030378610-pat00038
silver
Figure 112017030378610-pat00039
of
Figure 112017030378610-pat00040
and
Figure 112017030378610-pat00041
Using all of the past values of the regression model
Figure 112017030378610-pat00042
and
Figure 112017030378610-pat00043
Is an error that occurs when estimating the current value of < RTI ID = 0.0 >

[수학식 8]&Quot; (8) "

Figure 112017030378610-pat00044
Figure 112017030378610-pat00044

상기 수학식 8에서,

Figure 112017030378610-pat00045
Figure 112017030378610-pat00046
Figure 112017030378610-pat00047
을 자동 회귀 모델의 계수 행렬을 푸리에 변환한 것이며,
Figure 112017030378610-pat00048
Figure 112017030378610-pat00049
행렬의
Figure 112017030378610-pat00050
번째 열벡터이다.In Equation (8)
Figure 112017030378610-pat00045
The
Figure 112017030378610-pat00046
and
Figure 112017030378610-pat00047
Is the Fourier transform of the coefficient matrix of the autoregressive model,
Figure 112017030378610-pat00048
The
Figure 112017030378610-pat00049
Matrix
Figure 112017030378610-pat00050
Th column vector.

[수학식 9]&Quot; (9) "

Figure 112017030378610-pat00051
Figure 112017030378610-pat00051

상기 수학식 9에서,

Figure 112017030378610-pat00052
는 상기 수학식 8에서 나온
Figure 112017030378610-pat00053
의 역행렬이다.In Equation (9)
Figure 112017030378610-pat00052
Is expressed by Equation (8)
Figure 112017030378610-pat00053
/ RTI >

그리고 (b) 단계에서는 신호 감지 및 기록부(100)에서 측정된 생리신호를 통해 특징벡터 추출부(110a)의 생리신호 특징벡터 추출부(110a-2)에서 피부표면온도(skin temperature, SKT)를 이용한 특징벡터, 광혈류량(photo-plethysmography, PPG)을 이용한 특징벡터, 심전도(electrocardiography, ECG)을 이용한 특징벡터, 피부 전도도 반응(skin conductance response, SCR)을 이용한 특징벡터 및 호흡 신호를 이용한 특징벡터 중에서 선택되는 어느 하나 이상으로 생리신호 특징벡터를 추출할 수 있다.In step (b), the skin temperature (SKT) is calculated in the physiological signal feature vector extraction unit 110a-2 of the feature vector extraction unit 110a through the physiological signal measured by the signal detection and recording unit 100 A feature vector using the feature vector using the photo-plethysmography (PPG), a feature vector using the electrocardiography (ECG), a feature vector using the skin conductance response (SCR) The physiological signal feature vector may be extracted from the physiological signal feature vector.

일 구체예로 생리신호 특징벡터 추출부(110a-2)에서 추출하는 생리신호 특징 벡터는 아래에 설명하는 방법으로 추출할 수 있다.In one embodiment, the physiological signal feature vector extracted by the physiological signal feature vector extraction unit 110a-2 can be extracted by the method described below.

예를 들어 획득된 생리신호 중 피부 전도도(SCR)를 이용한 특징벡터의 경우, 피부 전도도 반응의 진폭 평균값(mean of SCR amplitude), 피부 전도도 반응의 음수 값의 평균값, 피부 전도도 반응의 표준 편차, 피부 전도도 반응의 1차 미분 신호의 평균값(mean of 1st derivative of SCR), 피부 전도도 반응의 2차 미분 신호의 평균값(mean of 2nd derivative of SCR), 피부 전도도 반응의 전체 샘플 개수 대비 음의 값을 가지는 샘플의 비율을 포함하는 특징 벡터를 추출하거나, 또는 피부 전도도가 느린 반응 값(skin conductance slow response, SCSR)을 이용하여 SCSR의 평균값(Mean of SCSR), SCSR 진폭의 평균값(mean of SCSR ampliutde), SCSR의 표준편차(standard deviation of SCSR), 1차 미분 신호의 평균값(mean of 1st derivative of SCSR), 2차 미분 신호의 평균값(mean of 2nd derivative of SCSR), SCSR occurrence 등을 포함하는 특징벡터를 추출할 수 있다.For example, in the case of feature vectors using skin conductivity (SCR) among acquired physiological signals, mean values of SCR amplitude of skin conductance response, mean value of negative values of skin conductivity response, standard deviation of skin conductivity response, The mean value of the first derivative signal of the conductivity response, the mean derivative of the second derivative signal of the skin conductivity response, and the negative value of the total number of samples of the skin conductivity response. The mean of the SCSR amplitude, the mean of the SCSR amplitude, and the mean value of the SCSR using the skin conductance slow response (SCSR) A feature vector including a standard deviation of the SCSR, a mean of the first derivative of the SCSR, a mean of the second derivative of the SCSR, Can be extracted.

또 다른 하나의 생리신호 특징벡터로 광혈류량을 이용한 특징벡터의 경우는, 광혈류량을 이용하여 얻을 수 있는 심박수의 평균, 심박 수 변이(Heart rate variability)의 평균, 광혈류량의 표준편차, 광혈류량의 1차 미분의 평균, 광혈류량의 1차 미분의 표준편차, 및 광혈류량의 저주파 성분과 고주파 성분의 비율중 어느 하나 이상을 포함하는 특징 벡터를 추출할 수 있다.In the case of the feature vector using the optical blood flow as the other physiological signal feature vector, the average of the heart rate obtained by using the optical blood flow, the average of the heart rate variability, the standard deviation of the optical blood flow, A standard deviation of the first derivative of the optical blood flow, and a ratio of the low-frequency component and the high-frequency component of the optical blood flow can be extracted.

또 다른 하나의 생리신호 특징벡터로 피부표면온도를 이용한 특징벡터 경우는, 피부 표면 온도 변화를 나타내는 1차 미분 값의 평균, 피부표면온도의 최대값, 및 0.1 내지 0.2 헤르츠(Hz)의 주파수 대역 스펙트럼 파워 중 어느 하나 이상을 포함하는 특징벡터를 추출할 수 있다.The feature vector using the skin surface temperature as the other physiological signal feature vector includes the average of the first derivative showing the skin surface temperature change, the maximum value of the skin surface temperature, and the frequency band of 0.1 to 0.2 hertz (Hz) The spectral power, and the spectral power.

(c) 단계(S130)는 상기 (b) 단계에서 추출된 뇌신경신호 특징벡터 및 생리신호 특징벡터를 신호 분석부(110)의 해독부(110b)를 통해 상기 뇌신경신호 특징벡터 및 생리신호 특징벡터를 이용하는 기계학습 알고리즘을 사용하여 동물의 운동 의도 및 생체 상태를 파악하는 단계이다.(c) In operation S130, the cranial nerve signal feature vector and the physiological signal feature vector extracted in the step (b) are input to the deciphering unit 110b of the signal analyzer 110 to convert the cranial nerve signal feature vector and the physiological signal feature vector Is a step of grasping the motion intention and the biological state of the animal using a machine learning algorithm using the algorithm.

여기서 상기 운동 의도는 동물의 움직임, 정지 및 움직임의 좌우 방향의 정보 중 어느 하나 이상이며, 상기 생체 상태는 두려움, 경계심, 정신적 워크로드, 피로도, 수면상태, 감정상태, 각성 상태 및 해마의 공간 지도 형성의 정보 중 어느 하나 이상이다.Here, the motion intention is at least one of information on the movement, stop, and movement of the animal, and the biological condition includes at least one of fear, vigilance, mental workload, fatigue, sleep state, emotional state, Forming information.

상기 기계학습 알고리즘은, 구체적으로 획득한 스파이크 파형을 주성분 분석법(principal component analysis, PCA)과 독립성분 분석법(independent component analysis, ICA)을 수행하여 뇌신경신호 및 생리신호에 관한 기저벡터를 추출하는 단계; 상기 기저 벡터로부터 특징벡터를 추출하는 단계; 상기 특징벡터를 서브클래스(subclass)로 뇌신경신호 특징벡터 및 생리신호 특징벡터로 분류하는 클러스터링(clustering) 단계; 및 상기 클러스터링 방법에 의하여 분류된 특징벡터를 서포터 벡터 머신(Support Vector Machine, SVM) 및 서포트 벡터 회귀(Support Vector Regression, SVR)을 기반으로 동물의 운동 의도 및 생체 상태를 파악하는 단계;를 포함할 수 있다.The machine learning algorithm includes extracting a basis vector related to a cranial nerve signal and a physiological signal by performing principal component analysis (PCA) and independent component analysis (ICA) on a spike waveform obtained in detail; Extracting a feature vector from the basis vector; A clustering step of classifying the feature vector into a cranial nerve signal feature vector and a physiological signal feature vector in a subclass; And a step of grasping the feature vector classified by the clustering method based on a support vector machine (SVM) and a support vector regression (SVR) .

상기 동물의 운동 의도 및 생체 상태를 파악하는 단계는 추출된 뇌신경신호 특징벡터 및 생리신호 특징벡터가 각각 해독부(110b)의 동물의 움직임 및 정지 여부와 왼쪽 및 오른쪽에 대한 방향성에 관한 운동 의도를 이진(binary) 값으로 출력하는 이진 결정 의도 인식기(110b-1)와 공포 수준, 경계심 수준, 피로도, 워크로드(스트레스), 졸림 수준 및 통증 수준 등을 포함한 생체 상태를 연속 실수 값으로 출력하는 회귀부(110b-2)에 입력되어 동물의 운동 의도 및 생체 상태를 파악할 수 있다.In the step of grasping the motion intention and the biological state of the animal, it is preferable that the extracted cranial nerve signal characteristic vector and physiological signal characteristic vector respectively indicate motion intention and stop motion of the animal of the decryption unit 110b and motion intention about the left and right direction A binary decision intention recognizer 110b-1 for outputting a binary value and a regression line for outputting a continuous biometric state including a fear level, an alertness level, a fatigue level, a workload (stress), a sleepiness level and a pain level, (110b-2) so that the intention of the animal and the biological state of the animal can be grasped.

구체적으로 상기 (c) 단계는 상기 (b) 단계에서 분류된 특징벡터를 통해 이진 결정 의도 인식기(110b-1)에서 서포터 벡터 머신(SVM)을 기반으로 서포터 벡터 머신 모델을 생성하여 동물의 운동 의도를 출력하는 단계; 및 상기 (b) 단계에서 분류된 특징벡터를 통해 회귀부(110b-2)에서 서포트 벡터 회귀(SVR)을 기반으로 하여 동물의 생체 상태를 출력하는 단계;를 포함할 수 있다.Specifically, in the step (c), a supporter vector machine model is generated based on a supporter vector machine (SVM) in the binary decision intention recognizer 110b-1 through the feature vectors classified in the step (b) ; And outputting the biological state of the animal based on the support vector regression (SVR) in the regression unit 110b-2 through the feature vectors classified in the step (b).

상기 서포터 벡터 머신(SVM)은 분류된 특징벡터를

Figure 112017030378610-pat00054
이라고 취급할 때, 아래의 수학식 10 및 수학식 11에 의해 최적의 가중벡터
Figure 112017030378610-pat00055
와 최적의 바이어스
Figure 112017030378610-pat00056
를 구하여 서포트 벡터 머신 모델을 생성하는 것이다.The supporter vector machine (SVM)
Figure 112017030378610-pat00054
, The optimal weight vector < RTI ID = 0.0 >
Figure 112017030378610-pat00055
And optimum bias
Figure 112017030378610-pat00056
To generate a support vector machine model.

[수학식 10]&Quot; (10) "

Figure 112017030378610-pat00057
Figure 112017030378610-pat00057

상기 수학식 10에서,

Figure 112017030378610-pat00058
는 벡터의 전치(transpose) 연산이고,
Figure 112017030378610-pat00059
는 커널 함수이다.In Equation (10)
Figure 112017030378610-pat00058
Is a transpose operation of the vector,
Figure 112017030378610-pat00059
Is a kernel function.

상기 수학식 10에서 상기 커널 함수는 다음의 수학식 11와 같이 결정되는 RBF(radial-basis function)이다.In Equation (10), the kernel function is a radial-basis function (RBF) determined according to Equation (11).

[수학식 11]&Quot; (11) "

Figure 112017030378610-pat00060
Figure 112017030378610-pat00060

상기 수학식 11에서,

Figure 112017030378610-pat00061
는 특징벡터
Figure 112017030378610-pat00062
의 표준편차이다.In Equation (11)
Figure 112017030378610-pat00061
Is a feature vector
Figure 112017030378610-pat00062
Lt; / RTI >

상기 서포트 벡터 회귀(SVR)는 분류된 특징벡터를

Figure 112017030378610-pat00063
이라고 하고, 상기 분류된 특징벡터에 해당하는 뇌신경신호 및 신경신호의 관계로부터 얻게 되는 동물의 생체 상태의 지표 값을
Figure 112017030378610-pat00064
이라 할 때 종속변수
Figure 112017030378610-pat00065
로부터 최대 적절한 양수
Figure 112017030378610-pat00066
만큼의 편차(deviation)를 갖는 함수
Figure 112017030378610-pat00067
를 구하는 알고리즘이다.The support vector regression (SVR)
Figure 112017030378610-pat00063
, And an indicator value of an animal's biological condition obtained from a relationship between a cranial nerve signal and a nerve signal corresponding to the classified characteristic vector
Figure 112017030378610-pat00064
In this case,
Figure 112017030378610-pat00065
The most appropriate positive number from
Figure 112017030378610-pat00066
A function having a deviation of
Figure 112017030378610-pat00067
.

(d) 단계(S140)는 상기 (c) 단계에서 파악된 동물의 운동 의도 및 생체 상태를 보상 혹은 처벌여부를 결정하거나 생체 상태를 통해 동물의 체력 소모 정도와 임무 지속 가능 여부를 파악하여 임무지속 가능 여부를 결정하며, 결정된 사항에 맞춰 동물의 행동을 제어하기 위한 전기 자극 파라미터를 조절하는 단계이다.(d) Step (S140) determines whether the animal's exercise intention and / or the biological state determined in the step (c) is compensated or punished, , And adjusting the electrical stimulation parameters to control the behavior of the animal in accordance with the decision.

상기 전기 자극 파라미터는 전기 펄스의 폭, 진폭, 주파수, 지속시간 중 선택되는 어느 하나 이상이며, 이에 한정되지 않는다.The electric stimulation parameter may be any one or more selected from among width, amplitude, frequency, and duration of the electric pulse, but is not limited thereto.

상기 (e) 단계(S150)는 (d) 단계(S140)에서 조절하여 결정된 전기 자극 파라미터를 통해 동물의 특정 대뇌 영역에 전기 신호 자극을 가하여 동물의 행동을 제어하는 단계로서 동물의 자발적인 방향성 있는 움직임 자극, 강제적인 방향성 있는 움직임 자극, 보상 자극 및 처벌 자극 중 어느 하나의 자극을 가하여 동물의 행동을 제어할 수 있다.The step (e) (S150) comprises the steps of (d) controlling an animal's behavior by applying electrical signal stimulation to a specific cerebral region of the animal through the electrical stimulation parameter determined in step (S140), wherein the spontaneous directional movement The behavior of the animal can be controlled by applying any one of stimulation, forced directional motion stimulation, compensation stimulation, and punitive stimulation.

상기 동물의 특정 대뇌 영역은 자발적인 방향성 있는 움직임을 제시하기 위한 대뇌 체성 감각 영역, 강제적인 방향성 있는 움직임을 시키기 위한 바닥핵(basal ganglia)의 흑색선조체 경로(nigrostriatal pathway) 또는 시상(thalamus)의 체성감각 통로인 복측후내측핵(ventral posteromedial nucleus), 보상 자극을 가하기 위한 뇌측전다발 영역, 및 공포심을 이용한 처벌 자극을 가하기 위한 편도체(amygdala) 중 어느 하나 이상을 포함할 수 있으며, 이에 한정되지 않고 동물의 행동을 제어할 수 있는 대뇌 영역이라면 본 발명의 실시 범위에 포함될 수 있다.The specific cerebral region of the animal may include a cerebral somatosensory region to present spontaneous directional movement, a nigrostriatal pathway of basal ganglia for forced directional movement, or a somatosensory sense of thalamus A ventral posteromedial nucleus which is a pathway, an area of a brain-conduction bundle for applying a compensatory stimulus, and an amygdala for applying a punitive stimulus by using a fear. The amygdala may include, but is not limited to, The cerebral region capable of controlling the behavior of the brain can be included in the scope of the present invention.

본 발명의 다중 신경생리신호 기반의 동물 행동 제어 시스템 및 방법은 실시간으로 동물의 운동 의도 및 생체 상태를 신속하게 파악하고, 파악된 동물의 운동 의도 및 생체상태를 토대로 동물에게 적절한 전기 전극을 가해줌으로써 동물 움직임을 효과적 제어하고 생체 상태를 조절할 수 있는 효과가 있다.The animal behavior control system and method based on the multiple neuronal physiological signals of the present invention can quickly identify an animal's exercise intention and a living state in real time and apply an appropriate electric electrode to the animal based on the animal's exercise intention and living state There is an effect that can control the animal movement effectively and control the biological state.

도 1은 본 발명에 따른 다중 신경생리신호 기록 및 해독 기반의 동물 행동 제어 시스템의 모식도이다.
도 2는 본 발명의 동물 행동 제어 시스템에서 신호 기록부(100)가 뇌신경 활동과 자율신경 활동을 검출하는 과정을 간략하게 보여주는 도면이다.
도 3은 본 발명의 동물 행동 제어 시스템에서 신호 분석부(110)가 다중 신경생리신호로부터 특징 벡터를 추출하고 운동 의도와 생체 상태를 디코딩하는 과정을 보여주는 도면이다.
도 4는 본 발명의 동물 행동 제어 시스템에서 신호 분석부(110)의 특징 벡터 추출부(110a)가 다중 뇌신경신호와 다중 생리신호로부터 뇌신경신호의 특징 벡터와 생리신호의 특징 벡터를 추출하는 과정을 보여주는 도면이다.
도 5는 본 발명의 동물 행동 제어 시스템에서 신호 분석부(110)의 해독부(110b)가 뇌신경신호와 생리신호로부터 추출된 특징 벡터를 통해 운동 의도와 생체 상태를 해독하는 과정을 보여주는 도면이다.
도 6은 본 발명의 동물 행동 제어 시스템에서 제어부(120)가 운동 의도와 생체 상태에 따른 보상 또는 처벌 여부를 결정하고 자극을 생성하는 과정을 보여주는 도면이다.
도 7은 본 발명의 동물 행동 제어 시스템에서 행동 제어 명령 지시부(130)가 운동 의도와 상태 상태를 출력하고 사용자로부터 수동으로 행동 제어 명령과 생체 상태 조절에 관한 자극을 입력받고 명령 자극을 제시하는 과정을 보여주는 도면이다.
도 8은 본 발명의 동물 행동 제어 시스템을 이용한 동물 행동 제어 방법의 순서도이다.
FIG. 1 is a schematic diagram of an animal behavior control system based on recording and decoding of multiple neuronal physiological signals according to the present invention.
FIG. 2 is a schematic view showing a process in which the signal recording unit 100 detects cranial nerve activity and autonomic nervous activity in the animal behavior control system of the present invention.
FIG. 3 is a diagram illustrating a process in which the signal analysis unit 110 in the animal behavior control system of the present invention extracts a feature vector from multiple neural physiological signals and decodes a motion intention and a biological condition.
4 is a flowchart illustrating a process of extracting a feature vector of a cranial nerve signal and a feature vector of a physiological signal from a multiple cranial nerve signal and a multiple physiological signal in a feature vector extraction unit 110a of the signal analysis unit 110 in the animal behavior control system of the present invention Fig.
FIG. 5 is a diagram illustrating a process in which the decoding unit 110b of the signal analysis unit 110 in the animal behavior control system of the present invention decodes the motion intention and the biological state through a feature vector extracted from a cranial nerve signal and a physiological signal.
FIG. 6 is a diagram illustrating a process in which the controller 120 in the animal behavior control system of the present invention determines whether to compensate or punish according to a motion intention and a living body state, and generates a stimulus.
FIG. 7 is a flowchart illustrating an operation of the behavior control command instructing unit 130 in the animal behavior control system of the present invention. FIG. 7 is a flowchart illustrating a process of outputting a motion intention and a state state, manually inputting a behavior control command, Fig.
8 is a flowchart of an animal behavior control method using the animal behavior control system of the present invention.

이하, 본 발명을 도면을 참조하여 상세히 설명하기로 한다. 본 명세서 및 청구범위에 사용된 용어나 단어는 통상적이거나 사전적인 의미로 한정해서 해석되어서는 아니 되며, 발명자는 그 자신의 발명을 가장 최선의 방법으로 설명하기 위해 용어의 개념을 적절하게 정의할 수 있다는 원칙에 입각하여 본 발명의 기술적 사상에 부합하는 의미와 개념으로 해석되어야만 한다.Hereinafter, the present invention will be described in detail with reference to the drawings. The terms and words used in the present specification and claims should not be construed as limited to ordinary or dictionary terms and the inventor may appropriately define the concept of the term in order to best describe its invention It should be construed as meaning and concept consistent with the technical idea of the present invention.

도 1은 본 발명의 동물행동 제어 시스템의 모식도이다.1 is a schematic diagram of an animal behavior control system of the present invention.

도 1을 참조하면, 본 발명의 동물행동 제어 시스템은 동물의 다중 뇌신경신호 및 다중 생리신호를 감지 및 기록하는 신호 감지 및 기록부(100), 상기 획득된 다중 뇌신경신호 및 다중 생리신호를 분석하고, 다중 뇌신경신호 및 다중 생리신호의 시간, 시간-주파수, 위상 영역에서 동물의 상태를 반영하는 특징을 수집하는 신호 분석부(110), 상기 수집된 특징을 분류하여 동물의 상태를 파악하고 방향성 있는 움직임을 제어하기 위하여 전기 자극 파라미터를 조절하는 제어부(120), 및 상기 조절된 자극 파라미터를 이용하여 특정 대뇌 영역에 전기 자극을 출력하는 행동제어 명령 지시부(130)를 포함하여 구성된다.Referring to FIG. 1, the animal behavior control system of the present invention includes a signal sensing and recording unit 100 for sensing and recording multiple neural signals and multiple physiological signals of an animal, analyzing the obtained multiple neuronal signals and multiple physiological signals, A signal analysis unit 110 for collecting features reflecting the state of the animal in the time, time-frequency, and phase regions of the multiple cranial nerve signals and the multiple physiological signals, And a behavior control instruction instructing unit 130 for outputting electric stimulation to a specific cerebral region using the adjusted stimulation parameters.

상기 신호 감지 및 기록부(100)는 동물의 생리신호를 신호의 발생되는 해부학적 위치에 따라 중추 신경계(central nervous system) 신호와 말초 신경계(peripheral nervous system) 신호로 두 종류의 기록 전극을 이용하여 측정하며, 도 2에 도시된 바와 같이, 신호 감지 및 기록부(100)는 뇌신경신호 기록 전극(100a), 생리신호 기록 전극(100b) 및 증폭부(100c)를 포함한다.The signal sensing and recording unit 100 measures a physiological signal of an animal using two kinds of recording electrodes as a central nervous system signal and a peripheral nervous system signal according to an anatomical position where a signal is generated As shown in FIG. 2, the signal sensing and recording unit 100 includes a cranial nerve signal recording electrode 100a, a physiological signal recording electrode 100b, and an amplification unit 100c.

좀 더 상세하게 상기 뇌신경 신호 기록 전극(100a)은 생체 신호 중에서 중추 신경계(central nervous system) 반응에 의한 신경세포로부터 발생하는 활동전위, 국소뇌전위(local field potential), 뇌피질 뇌파(subdural electrocortigram, ECoG) 및 뇌전도(electroencephalography, EEG) 등을 포함한 다중 뇌신경신호를 획득하기 위해 복수의 미세전극을 이용한 세포외 기록(extracellular recording) 방법으로 측정되며, 여기서 사용되는 뇌신경 신호 기록 전극으로 미세전극 어레이(multi-electrode array, MEA), 경막하 전극, 심부 전극(deep electrode) 또는 뇌파 전극 중에서 선택되는 어느 하나 이상을 사용할 수 있다.In more detail, the nerve signal recording electrode 100a may be used to measure an action potential, a local field potential, a subdural electrocortigram, or the like generated from nerve cells in response to a central nervous system reaction, (E. G., ECoG), electroencephalography (EEG), etc., is measured by extracellular recording using a plurality of microelectrodes to obtain a multi-nerve signal, an electrodial array (MEA), a subdural electrode, a deep electrode, or an EEG electrode.

상기 생리신호 기록 전극(100b)은 생체 신호 중에서 말초 신경계 반응에 의한 동물의 피부표면온도(skin temperature, SKT), 광혈류량(photo-plethysmography, PPG), 심전도(electrocardiography, ECG), 피부 전도도 반응(skin conductance response, SCR) 및 호흡 신호 등을 포함하는 다중 생리신호를 측정하기 위한 것으로, 동물의 체표면에 부착되는 다양한 생체 상대 감지 수단을 이용하여 상기 제시된 다중 생리신호를 얻을 수 있다.The physiological signal recording electrode 100b is an electrode for physiologic signal recording in which the skin temperature (SKT), photo-plethysmography (PPG), electrocardiography (ECG), skin conductivity skin conductance response (SCR), respiration signal, etc., and the multiple physiological signals can be obtained using various biological relative sensing means attached to the body surface of an animal.

이처럼 신호 감지 및 기록부(100)에서는 뇌신경신호 기록 전극(100a)과 생리신호 기록 전극(100b)을 통해 동물의 생체 신호로서 다중 뇌신경신호 및 다중 생리신호를 감지하고, 감지된 신호들을 기초로 동물의 생체신호를 기록한다.As described above, the signal sensing and recording unit 100 senses multiple neural signals and multiple physiological signals as biological signals of an animal through the nerve signal recording electrode 100a and the physiological signal recording electrode 100b, Record the biomedical signal.

본 발명에서 뇌신경신호이나 생리신호를 측정할 때에는 머리를 포함한 신체의 여러 부분에서 신호를 측정하므로 본 발명의 뇌신경신호 및 생리신호는 단일의 신호 뿐만 아니라 둘 이상의 다중 신호를 측정한다. In the present invention, when measuring a cranial nerve signal or a physiological signal, since the signal is measured in various parts of the body including the head, the cranial and physiological signals of the present invention measure not only a single signal but also two or more multiple signals.

그러므로 신호 감지 및 기록부(100)에서 측정된 뇌신경신호 및 생리신호를 또 다른 용어로 다중 뇌신경신호 및 다중 생리신호라 명칭 할 수 있으며, 이하에서는 이러한 다중 뇌신경신호 및 다중 생리신호를 통틀어 "다중 신경생리신호"이라 명칭하기도 한다.Therefore, the cranial nerve signals and physiological signals measured by the signal sensing and recording unit 100 may be referred to as multipath cranial signals and multiple physiological signals in another term. Hereinafter, such multiple cranial nerve signals and multiple physiological signals are collectively referred to as " Signal "

그리고 상기 신호 감지 및 기록부(100)에서 구성된 증폭부(100c)는 뇌신경신호 기록 전극(100a)과 생리신호 기록 전극(100b)을 통해 감지된 각각의 다중 신경생리신호를 증폭하고, 증폭된 다중 신경 생리신호를 신호 분석부(110)로 전송하는 역할을 한다.The amplification unit 100c configured in the signal sensing and recording unit 100 amplifies each of the multiple neuro physiological signals sensed through the brain signal recording electrode 100a and the physiological signal recording electrode 100b, And transmits a physiological signal to the signal analyzer 110.

도 3 내지 도 5는 신호 분석부(110)에서 다중 신경생리신호로부터 추출된 특징벡터를 통해 운동 의도와 생체 상태를 해독하는 과정을 나타낸 것이다.FIGS. 3 to 5 illustrate a process of deciphering a motion intention and a living body through a feature vector extracted from a multi-nervous physiological signal in the signal analyzer 110. FIG.

도 3를 참조하면, 신호 분석부(110)는 상기 신호 감지 및 기록부(100)에서 전송받은 다중 신경생리신호로부터 특징벡터를 추출하고, 추출된 특징벡터를 통해 대상 동물의 운동 의도와 생체 상태를 해독하는 것으로, 특징벡터 추출부(110a)와 해독부(110b)를 포함하여 구성될 수 있다.3, the signal analyzer 110 extracts a feature vector from the multi-nervous physiological signals transmitted from the signal sensing and recording unit 100, and outputs the motion intention and the biological state of the target animal through the extracted feature vector And may include a feature vector extraction unit 110a and a decryption unit 110b.

도 4에 도시된 바와 같이 특징벡터 추출부(110a)는 뇌신경신호 특징 벡터 추출부(110a-1)와 생리신호 특징벡터 추출부(110a-2)를 포함하여 구성되며, 상기 신호 감지 및 기록부(100)에서 측정된 다중 뇌신경신호 및 다중 생리신호를 입력신호로 하여 각 해당 입력 신호로부터 특징벡터를 추출한다. 4, the feature vector extraction unit 110a includes a cranial nerve signal feature vector extraction unit 110a-1 and a physiological signal feature vector extraction unit 110a-2. The signal detection and recording unit 100), and extracts a feature vector from each corresponding input signal using the multiple physiological signals as input signals.

뇌신경신호 특징벡터 추출부(110a-1)는 상기 뇌신경신호 기록 전극(100a)에서 얻은 뇌신경신호로부터 스파이크 트레인(spike train), 파워 스펙트럼 밀도(power spectrum density), 뇌 영역 간의 방향성을 포함하지 않는 연결정도를 표현하는 기능적 연결도(functional connectivity) 및 뇌 영역 간 방향성을 가지는 연결정도를 표현하는 효과적 연결도(effective connectivity)를 추출하여 뇌신경신호 특징벡터를 추출할 수 있다.The cranial nerve signal feature vector extraction unit 110a-1 extracts a spike train, a power spectrum density, and a connection including no directionality between brain regions from the cranial nerve signal obtained from the cranial nerve signal recording electrode 100a And the effective connectivity expressing the degree of connection having the directionality between the brain regions can be extracted to extract the cranial nerve signal feature vector.

구체적으로 뇌신경신호 특징 벡터 추출부(110a-1)는 아래와 같은 방법으로 계산되어 뇌신경신호 특징벡터를 추출할 수 있다.Specifically, the cranial nerve signal feature vector extraction unit 110a-1 can calculate the cranial nerve signal feature vector by the following method.

스파이크 트레인(spike train)은 측정된 세포외 기록 신호(extracellular recording signal)를 300 헤르츠(Hz)의 컷오프 주파수를 가지는 고대역 통과 필터 (highpass filter)를 통과시켜 주변의 저주파 잡음과 국소뇌전위(local field potential)를 제거하는 단계; 상기 고대역 통과 필터를 통과한 출력 신호에 음의 값을 가지는 문턱치(threshold)를 적용하여 주변 잡음보다 값이 큰 피크를 잡는 피크설정 단계; 상기 피크설정 단계로 부터 얻어진 피크를 중심으로 좌우로 15개의 샘플을 획득하여 신경신호의 스파이크 파형을 획득하는 단계; 상기 획득된 스파이크 파형을 주성분 분석법(principal component analysis, PCA) 및 K-평균 알고리즘(K-means algorithm)을 이용해 스파이크들의 형태에 따라 두 그룹 내지는 세 그룹으로 분류하는 단계; 상기 분류된 스파이크로부터 단위시간에 들어있는 스파이크의 갯수를 요소(element)로 하는 뇌신경신호 특징벡터를 추출할 수 있다.The spike train passes the measured extracellular recording signal through a highpass filter with a cutoff frequency of 300 Hertz to generate low frequency noise and localized local removing a field potential; A peak setting step of applying a threshold value having a negative value to an output signal passed through the high pass filter to obtain a peak having a value larger than the ambient noise; Acquiring fifteen samples to the left and right about a peak obtained from the peak setting step to obtain a spike waveform of a neural signal; Classifying the obtained spike waveform into two or three groups according to the types of spikes using principal component analysis (PCA) and K-means algorithm; A cranial nerve signal feature vector having the number of spikes included in the unit time as an element can be extracted from the classified spikes.

뇌신경신호 특징벡터중 하나로 파워 스펙트럼 밀도(power spectrum density)는 측정된 뇌신경신호 중 어느 하나의 신호를 M차원의 벡터

Figure 112017030378610-pat00068
로 하여, 하기 수학식 1을 통해 분석하는 이산 푸리에 변환(discrete Fourier transform), 하기 수학식 2를 통해 분석하는 단기간 푸리에 변환(short-time Fourier transform) 및 하기 수학식 3을 통해 분석하는 웨이브릿 변환(wavelet transform) 중의 어느 하나를 이용하여 특정 시간-주파수 대역의 파워를 얻어 뇌신경신호 특징벡터를 추출할 수 있다.The power spectrum density, which is one of the cranial nerve signal characteristic vectors, is a function of a signal of one of the measured nerve signals as an M-dimensional vector
Figure 112017030378610-pat00068
A discrete Fourier transform for analysis through the following equation (1), a short-time Fourier transform for analysis through the following equation (2) and a wavelet transform for analysis through the following equation (3) and a wavelet transform to obtain a power of a specific time-frequency band and to extract a cranial nerve signal feature vector.

이하, 본 발명에서 설명하는 모든 수학식의 M차원의 벡터에서 변수 M은 이산 시계열의 길이를 의미하고, 바람직하게 M은 1 이상의 정수이며, 1 내지 1000의 정수일 수 있다.Hereinafter, the variable M in the M-dimensional vector of all the equations described in the present invention means the length of discrete time series, and preferably M is an integer of 1 or more, and may be an integer of 1 to 1000.

[수학식 1][Equation 1]

Figure 112017030378610-pat00069
Figure 112017030378610-pat00069

수학식 1에서

Figure 112017030378610-pat00070
는 주파수를 지칭하는 변수
Figure 112017030378610-pat00071
에 대하여
Figure 112017030378610-pat00072
를 의미하는 각주파수를 나타낸다.In Equation (1)
Figure 112017030378610-pat00070
Is a variable that refers to frequency
Figure 112017030378610-pat00071
about
Figure 112017030378610-pat00072
≪ / RTI >

[수학식 2]&Quot; (2) "

Figure 112017030378610-pat00073
Figure 112017030378610-pat00073

수학식 2에서

Figure 112017030378610-pat00074
은 특정한 길이의 특정한 모양의 시간 윈도우 벡터를 나타낸다.In Equation 2,
Figure 112017030378610-pat00074
Represents a time window vector of a particular shape of a specific length.

[수학식 3]&Quot; (3) "

Figure 112017030378610-pat00075
Figure 112017030378610-pat00075

수학식 3에서

Figure 112017030378610-pat00076
는 웨이블릿 함수의 공액복소수(complex conjugate),
Figure 112017030378610-pat00077
는 정규화 파라미터,
Figure 112017030378610-pat00078
는 특정 스케일,
Figure 112017030378610-pat00079
는 목표 레이턴시를 나타낸다.In Equation 3,
Figure 112017030378610-pat00076
Is a complex conjugate of the wavelet function,
Figure 112017030378610-pat00077
Is a normalization parameter,
Figure 112017030378610-pat00078
Lt; / RTI >
Figure 112017030378610-pat00079
Represents the target latency.

본 발명의 일 구체예에 따르면 파워스펙트럼 밀도(power spectrum density)를 얻는 것은, 상기 측정된 뇌신경신호로 세포외 기록 신호(extracellular recording) 또는 뇌전도 신호를 저대역 통과 필터(lowpass filter)를 통과시켜 주변의 고주파 잡음을 제거하는 단계와, 상기 수학식 1로 이산 푸리에 변환(discrete Fourier transform)과 자동회귀 모델링(autoregressive modeling)을 이용하여 시간정보가 없는 파워스펙트럼 밀도를 얻거나, 상기 수학식 2로 단시간 푸리에 변환(short-time Fourier transform)를 계산하거나 또는 상기 수학식 3으로 웨이블릿 변환(wavelet transform)를 계산하여 시간 정보를 갖는 파워스펙트럼 밀도를 얻는 단계를 포함할 수 있다.According to one embodiment of the present invention, obtaining the power spectrum density is performed by passing the extracellular recording signal or the electroencephalogram signal through the low pass filter with the measured cranial nerve signal Frequency noise without time information is obtained by using discrete Fourier transform and autoregressive modeling in Equation (1), or the power spectral density without time information is obtained by Equation (1) Calculating a short-time Fourier transform or calculating a wavelet transform in Equation (3) to obtain a power spectral density having time information.

상기 대뇌 영역 간 기능적 연결도(functional connectivity)는 뇌신경신호 중 어느 둘의 뇌신경신호를 한 쌍의 M차원의 벡터

Figure 112017030378610-pat00080
,
Figure 112017030378610-pat00081
로 하여, 하기 수학식 4를 통해 분석하는 간섭성(coherence), 하기 수학식 5를 통해 분석하는 위상 동기 값(phase locking value, PLV), 및 하기 수학식 6을 통해 분석하는 weighted phase lag index(wPLI) 중의 어느 하나를 이용하여 한 쌍의 뇌신경 영역 간의 기능적 연결성을 확인하여 뇌신경신호 특징벡터를 추출할 수 있다.The functional connectivity between the cerebral regions is determined by connecting two neurons of the cerebral nerve signal to a pair of M-dimensional vectors
Figure 112017030378610-pat00080
,
Figure 112017030378610-pat00081
, A phase locking value (PLV) analyzed through Equation (5), and a weighted phase lag index (PLV) analyzed through Equation (6) wPLI) to identify the functional connectivity between a pair of neuronal regions to extract the neural signal characteristic vector.

[수학식 4]&Quot; (4) "

Figure 112017030378610-pat00082
Figure 112017030378610-pat00082

수학식 4에서

Figure 112017030378610-pat00083
Figure 112017030378610-pat00084
Figure 112017030378610-pat00085
Figure 112017030378610-pat00086
의 진폭,
Figure 112017030378610-pat00087
는 두 시계열의 상대적 위상차를 나타낸다.In Equation 4,
Figure 112017030378610-pat00083
Wow
Figure 112017030378610-pat00084
The
Figure 112017030378610-pat00085
Wow
Figure 112017030378610-pat00086
The amplitude,
Figure 112017030378610-pat00087
Represents the relative phase difference between two time series.

[수학식 5]&Quot; (5) "

Figure 112017030378610-pat00088
Figure 112017030378610-pat00088

수학식 5에서

Figure 112017030378610-pat00089
Figure 112017030378610-pat00090
는 각각 시계열 벡터
Figure 112017030378610-pat00091
Figure 112017030378610-pat00092
에 힐버트 변환을 적용하여 얻은 순간 위상(instantaneous phase)을 나타낸다.In Equation (5)
Figure 112017030378610-pat00089
and
Figure 112017030378610-pat00090
Respectively,
Figure 112017030378610-pat00091
and
Figure 112017030378610-pat00092
And shows the instantaneous phase obtained by applying the Hilbert transform to the signal.

[수학식 6]&Quot; (6) "

Figure 112017030378610-pat00093
Figure 112017030378610-pat00093

수학식 6에서

Figure 112017030378610-pat00094
은 기대값 연산,
Figure 112017030378610-pat00095
은 복소수의 허수 항만을 취하는 연산자,
Figure 112017030378610-pat00096
는 두 시계열 벡터
Figure 112017030378610-pat00097
Figure 112017030378610-pat00098
의 교차 스펙트럼(cross spectrum)을 나타낸다.In Equation (6)
Figure 112017030378610-pat00094
The expected value calculation,
Figure 112017030378610-pat00095
Is an operator that takes an imaginary number of complex numbers,
Figure 112017030378610-pat00096
The two time series vectors
Figure 112017030378610-pat00097
and
Figure 112017030378610-pat00098
Of the cross spectrum.

보다 구체적으로 대뇌 영역 간 기능적 연결도는 뇌신경신호 중 어느 둘의 신호로 이루어진 한 쌍의 뇌신경신호에서 DC 성분을 제거하는 단계와, DC 성분이 제거된 한 쌍의 뇌신경신호를 관찰하고자 하는 주파수 대역 통과 필터로 통과시켜주는 단계와, 힐버트 변환(Hilbert transform)하여 순간 위상(instantaneous phase) 신호를 생성하는 과정과, 상기 수학식 8로 한 쌍의 뇌신경신호의 순간 위상(instantaneous phase)의 차이를 계산하는 단계와, 상기 수학식 9로 한 쌍의 뇌신경신호의 순간 위상의 차이,

Figure 112017030378610-pat00099
를 위상으로 가지는
Figure 112017030378610-pat00100
함수의 시간 평균 또는 트라이얼(trial) 평균을 계산하여 얻은 위상 동기 값(phase locking value, PLV)을 통해 기능적 연결도를 구하는 단계를 포함할 수 있다.More specifically, functional connectivity between the cerebral regions includes a step of removing a DC component from a pair of neuronal signals consisting of a signal of either one of the signals of the brain, a step of passing through a frequency band for observing a pair of neuronal signals from which the DC component is removed Passing the signal through a filter, generating an instantaneous phase signal by Hilbert transform, and calculating a difference between instantaneous phases of a pair of cranial nerve signals in Equation (8) The difference of the instantaneous phase of a pair of cranial nerve signals in Equation (9)
Figure 112017030378610-pat00099
In phase
Figure 112017030378610-pat00100
And obtaining a functional connectivity through a phase locking value (PLV) obtained by calculating a time average or a trial average of the function.

또 다른 대뇌 영역 간 기능적 연결도를 구하는 방법은 뇌신경신호 중 어느 둘의 신호로 이루어진 한 쌍의 뇌신경신호에서 DC 성분을 제거하는 과정, DC 성분이 제거된 한 쌍의 뇌신경신호의 상호스펙트럼(cross spectrum)을 구해주는 과정 및 상기 상호스펙트럼(cross spectrum)을 상기 수학식 9를 통해 얻을 수 있는 weighted phase lag index(wPLI)를 통해 기능적 연결도를 구해주는 과정을 포함할 수 있다.Another method for obtaining functional connectivity between the cerebral regions is to remove DC components from a pair of neuronal signals consisting of either of two signals of the nervous signal, to measure the cross spectrum of a pair of neuronal signals from which the DC component has been removed And a process of obtaining the cross spectrum of the functional connectivity through a weighted phase lag index (wPLI) obtained by Equation (9).

상기 대뇌 영역 간 효과적 연결도(effective connectivity)의 경우는 뇌신경신호 중 어느 둘의 뇌신경신호를 한 쌍의 M차원의 벡터

Figure 112017030378610-pat00101
,
Figure 112017030378610-pat00102
로 하여, 하기 수학식 7을 통해 분석하는 그랜저 인과관계(Granger's causality), 하기 수학식 8을 통해 분석하는 PDC(Partial Directed Coherence), 및 하기 수학식 9를 통해 분석하는 DTF(Directed Transfer Function) 중 어느 하나 이상을 이용하여 한 쌍의 뇌신경 영역 간의 방향성이 있는 연결성을 확인하여 뇌신경신호 특징벡터로 대뇌 영역 간 효과적 연결도(effective connectivity)를 추출할 수 있다.In the case of the effective connectivity between the cerebral regions, two of the cerebral nerve signals among the cerebral nerve signals are divided into a pair of M-dimensional vectors
Figure 112017030378610-pat00101
,
Figure 112017030378610-pat00102
(Granger's causality) analyzed by the following equation (7), PDC (Partial Directed Coherence) analyzed through the following equation (8), and DTF (Directed Transfer Function) analyzed through the following equation Using one or more of the above, effective connectivity between the cerebral regions can be extracted with the neural signal characteristic vector by confirming the directional connectivity between the pair of neuronal regions.

[수학식 7]&Quot; (7) "

Figure 112017030378610-pat00103
Figure 112017030378610-pat00103

수학식 7에서

Figure 112017030378610-pat00104
Figure 112017030378610-pat00105
Figure 112017030378610-pat00106
의 과거 값만을 이용하여 자동 회귀 모델을 이용해 추정할 때 생기는 에러,
Figure 112017030378610-pat00107
Figure 112017030378610-pat00108
Figure 112017030378610-pat00109
Figure 112017030378610-pat00110
의 과거값을 모두 사용하여 자동 회귀 모델을 이용해
Figure 112017030378610-pat00111
Figure 112017030378610-pat00112
의 현재 값을 추정할 때 생기는 에러이다.In Equation (7)
Figure 112017030378610-pat00104
silver
Figure 112017030378610-pat00105
of
Figure 112017030378610-pat00106
The error of the estimation using the automatic regression model using only the past value of the error,
Figure 112017030378610-pat00107
silver
Figure 112017030378610-pat00108
of
Figure 112017030378610-pat00109
and
Figure 112017030378610-pat00110
Using all of the past values of the regression model
Figure 112017030378610-pat00111
and
Figure 112017030378610-pat00112
Is an error that occurs when estimating the current value of < RTI ID = 0.0 >

[수학식 8]&Quot; (8) "

Figure 112017030378610-pat00113
Figure 112017030378610-pat00113

상기 수학식 8에서

Figure 112017030378610-pat00114
Figure 112017030378610-pat00115
Figure 112017030378610-pat00116
을 자동 회귀 모델의 계수 행렬을 푸리에 변환한 것이며,
Figure 112017030378610-pat00117
Figure 112017030378610-pat00118
행렬의
Figure 112017030378610-pat00119
번째 열벡터이다.In Equation (8)
Figure 112017030378610-pat00114
The
Figure 112017030378610-pat00115
and
Figure 112017030378610-pat00116
Is the Fourier transform of the coefficient matrix of the autoregressive model,
Figure 112017030378610-pat00117
The
Figure 112017030378610-pat00118
Matrix
Figure 112017030378610-pat00119
Th column vector.

[수학식 9]&Quot; (9) "

Figure 112017030378610-pat00120
Figure 112017030378610-pat00120

상기 수학식 9에서

Figure 112017030378610-pat00121
는 상기 수학식 8에서 나온
Figure 112017030378610-pat00122
의 역행렬이다.In Equation (9)
Figure 112017030378610-pat00121
Is expressed by Equation (8)
Figure 112017030378610-pat00122
/ RTI >

생리신호 특징 벡터 추출부(110a-2)는 상기 생리신호 기록 전극(100b)로부터 얻은 다중 생리신호로부터 피부 전도도 반응 값(SCR)을 이용한 특징벡터, 광혈류량(PPG) 이용한 특징 벡터 및 피부표면온도(SKT)를 이용한 특징벡터를 계산하여 추출할 수 있다.The physiological signal feature vector extraction unit 110a-2 extracts, from the physiological signal recording electrode 100b, a feature vector using the skin conductivity response value (SCR), a feature vector using the optical blood flow rate (PPG) (SKT) can be calculated and extracted.

좀 더 상세하게 생리신호 특징벡터 추출부(110a-2)에서는 생리신호 특징 벡터를 추출하는 방법로서, 아래와 같이 설명하는 방법을 통해 생리신호 특징벡터를 추출할 수 있다.In more detail, the physiological signal feature vector extracting unit 110a-2 extracts the physiological signal feature vector by the following method.

일예로 피부 전도도 반응(SCR)의 경우, 피부 전도도 반응 값(SCR)을 이용하여 피부 전도도 진폭의 평균값(mean of SCR amplitude), 피부 전도도 반응의 음수 값 평균값, 피부 전도도 반응의 표준 편차, 피부 전도도 반응의 1차 미분 신호의 평균값(mean of 1st derivative of SCR), 2차 미분 신호의 평균값(mean of 2nd derivative of SCR), 및 전체 샘플 개수 대비 음의 값을 가지는 샘플의 비율을 포함하는 특징 벡터를 구하거나, 또는 피부 전도도가 느린 반응 값(skin conductance slow response, SCSR)을 이용하여 SCSR의 평균값(Mean of SCSR), SCSR 진폭의 평균값(mean of SCSR ampliutde), SCSR의 표준편차(standard deviation of SCSR), 1차 미분 신호의 평균값(mean of 1st derivative of SCSR), 2차 미분 신호의 평균값(mean of 2nd derivative of SCSR), SCSR occurrence 등을 포함하는 특징 벡터를 검출할 수 있다.For example, in the case of skin conductivity response (SCR), mean value of SCR amplitude of skin conductance, mean value of skin conductivity response, standard deviation of skin conductivity response, skin conductivity A characteristic vector including a mean of a first derivative signal of the response (SCR), a mean of a second derivative signal (SCR), and a ratio of a sample having a negative value relative to the total number of samples (SCSR), the mean of the SCSR amplitudes, and the standard deviation of the SCSR using the skin conductance slow response (SCSR) (SCSR), a mean value of the first derivative signal (SCSR), a mean value of the second derivative signal (SCSR), and a SCSR occurrence.

예를 들어, 상기 피부 전도도 반응(skin conductance response, SCR)을 이용한 특징벡터는 측정한 피부전도도 반응 신호에서 분석하고자 하는 주파수 대역을 제한하기 위하여 샘플링 주파수를 변환하는 과정과, 피부전도도 변화를 관찰하기 위하여 미분 및 bartlett 윈도우(삼각형 윈도우)를 통한 컨볼루션(convolution) 과정과, 피부전도도의 상 활성(phasic activity)을 관찰하기 위하여 한계점을 초과하는 신호를 분리하는 과정과, 피부전도도의 긴장성 활성(tonic activity)을 관찰하기 위하여 피부전도도 반응 신호를 저대역 통과 필터(lowpass filter)에 통과시켜 주변의 고주파 잡음을 제거하는 과정과, 상대적인 변화 정도를 정규화하기 위하여 표준화(normalization)하는 과정과, 기저 수준에서 변화하는 생체 활동을 관찰하기 위하여 선형 트렌드를 제거해주는 디트렌드(detrend) 과정을 포함하여 추출할 수 있다.For example, the characteristic vector using the skin conductance response (SCR) may include a process of converting a sampling frequency to limit the frequency band to be analyzed in the measured skin conductivity response signal, A convolution process through differential and bartlett windows (a triangular window), a process of separating signals exceeding the threshold to observe the phasic activity of the skin conductivity, the process of tonic activity is filtered through a lowpass filter to remove the high frequency noise of the skin, a process of normalization to normalize the degree of relative change, Detrend, which removes linear trends to observe changing biological activity, A can be extracted, including.

또 다른 생리신호로 광혈류량의 경우는, 광혈류량을 이용하여 얻을 수 있는 심박수의 평균, 표준편차, 심박변이율의 평균, 표준편차, 1차 미분의 평균, 1차 미분의 표준편차, 저주파 성분과 고주파 성분의 비율을 포함하는 특징 벡터를 검출할 수 있다.In the case of photoperiod, the mean, standard deviation, mean, standard deviation, mean of the first derivative, standard deviation of the first derivative, low frequency component And a ratio of high-frequency components.

예를 들어, 상기 광혈류량(photo-plethysmography, PPG)를 이용한 특징벡터를 추출하는 방법은 측정한 광혈류량 측정 신호에서 지터(jitter)에 비하여 주기가 긴 변동을 제거해 주기 위한 기준선 변동 제거(baseline wander removal)하는 과정과, 심전도 신호 R파의 진폭 값을 추출하기 위하여 광혈류량 측정 신호를 두 번 미분한 후 부호(sign)의 변화에 따른 R파 진폭 부분을 분리해주는 R-peak 검출(R-peak detection) 과정과, 분리된 R파 진폭 간의 구간길이(R-R interval)를 구한 후 평활(smoothing)하여 심박 수 변이(heart rate variability) 신호를 얻는 보간(interpolation)하는 과정과, 심박 수 변이(heart rate variability) 신호를 자동 회귀 모델(autoregressive model)을 통하여 파워 스펙트럼(power spectrum)을 추정하는 과정을 포함하여 추출할 수 있다.For example, a method of extracting a feature vector using the photo-plethysmography (PPG) is a method of extracting a feature vector from a baseline wander to remove a long-period variation in a measured optical blood flow measurement signal from jitter (R-peak), which separates the R wave amplitude part according to the sign change after differentiating the optical blood flow measurement signal twice in order to extract the amplitude value of the R wave of the electrocardiogram signal, a step of interpolating a heart rate variability signal by smoothing the RR interval between the separated R wave amplitudes and a heart rate variability signal, and estimating the power spectrum of the variability signal through an autoregressive model.

또 다른 생리신호로 피부표면온도 경우는, 피부표면온도를 통해 얻을 수 있는 피부표면온도의 평균, 1차 미분 값의 평균, 최대값, 0.1 내지 0.2 헤르츠(Hz)의 주파수 성분의 스펙트럼 파워를 포함하는 특징벡터를 검출할 수 있다.As the other physiological signals, the case of the skin surface temperature includes the average of the skin surface temperature obtained through the skin surface temperature, the average of the first derivative, the maximum value, and the spectral power of the frequency component of 0.1 to 0.2 hertz (Hz) Can be detected.

예를 들어, 상기 피부표면온도(skin temperature, SKT)를 이용한 특징벡터를 추출하는 방법은 측정된 피부표면온도 신호에서 생체 상태를 반영하는 주파수 대역 분석을 위해 샘플링 주파수를 변환하는 과정, 생체 활동 변화를 관찰하기 위해 전체 DC를 제거해주는 기준선 제거(baseline removal) 과정, 긴장성(tonic) 변화를 관찰하기 위해 자동 회귀 모델(autoregressive model)을 이용하여 추정한 저 주파수 대역 스펙트럼 파워를 얻는 과정, 및 통상적인 피부 표면 온도 변화를 관찰하기 위해 미분을 하는 과정을 포함하여 특징벡터를 추출할 수 있다.For example, a method of extracting a feature vector using the skin temperature (SKT) includes a process of converting a sampling frequency to analyze a frequency band that reflects a living body in a measured skin surface temperature signal, A baseline removal process that removes the entire DC to observe the tonic change, a process of obtaining the low-frequency band spectral power estimated using an autoregressive model to observe the tonic change, A feature vector can be extracted including a process of differentiating to observe a skin surface temperature change.

도 5는 신호 분석부(110)의 해독부(110b)가 다중 신경생리신호로부터 추출된 특징 벡터를 통해 운동 의도와 생체 상태를 해독하는 과정을 나타낸 것으로, 도 5에 도시된 바와 같이 해독부(110b)는 이진 결정 의도 인식기(110b-1)와 회귀부(110b-2)를 포함하며, 해독부(110b)는 상기 특징벡터 추출부(110a)에서 추출된 상기 뇌신경신호의 특징 벡터와 생리 신호의 특징 벡터를 입력받아 이를 해독하여 동물의 운동 의도와 생체 상태를 파악한다.5 shows a process in which the decoding unit 110b of the signal analysis unit 110 decodes the motion intention and the living state through the feature vector extracted from the multiple neurophysiological signals. As shown in FIG. 5, 110b includes a binary decision intent recognizer 110b-1 and a regression unit 110b-2, and the decoding unit 110b includes a feature vector extracting unit 110a extracting a feature vector of the cranial nerve signal, And the animal's motion intention and the biological state are determined.

이때, 해독부(110b)에서 이진 결정 의도 인식기(110b-1)와 회귀부(110b-2)에 입력되는 특징 벡터는 뇌신경신호 특징벡터와 생리신호 특징벡터 모두를 입력하여 사용할 수 있다.At this time, the feature vector input to the binary decision intention recognizer 110b-1 and the regression unit 110b-2 in the decryption unit 110b may be used by inputting both the nerve signal feature vector and the physiological signal feature vector.

이진 결정 의도 인식기(110b-1)는 동물의 움직임 및 정지 여부와 왼쪽 및 오른쪽에 대한 방향성에 관한 운동 의도를 이진(binary) 값으로 출력하게 되고, 상기 회귀부(110b-2)는 공포 수준, 경계심 수준, 피로도, 워크로드(스트레스), 졸림 수준 및 통증 수준 등을 포함한 생체 상태를 연속 실수 값으로 출력하게 된다.The binary decision intent recognizer 110b-1 outputs a binary value of the motion intention regarding the movement and stop of the animal and the directionality to the left and right, and the regression unit 110b-2 outputs the motion intention, The body condition including the level of beating, fatigue, workload (stress), sleepiness level and pain level are outputted as continuous real values.

좀 더 상세하게 이진 결정 의도인식기(110b-1)는 서포트 벡터 머신(SVM)을 기반으로 하여 이진 결정을 실시하게 된다. 여기서 서포트 벡터 머신(SVM)은 상기 분류된 특징 벡터를

Figure 112017030378610-pat00123
이라고 취급할 때, 아래 제시된 수학식 5에 의해 최적의 가중벡터
Figure 112017030378610-pat00124
와 최적의 바이어스
Figure 112017030378610-pat00125
를 구하여 서포트 벡터 머신 모델을 생성하는 기계 학습 알고리즘이다.In more detail, the binary decision intent recognizer 110b-1 performs a binary decision based on a support vector machine (SVM). Wherein the support vector machine (SVM)
Figure 112017030378610-pat00123
, The optimal weight vector < RTI ID = 0.0 >
Figure 112017030378610-pat00124
And optimum bias
Figure 112017030378610-pat00125
To generate a support vector machine model.

그러므로 이진 결정 의도인식기(110b-1)는 상기 다중 신경생리신호와 동물의 운동 의도의 관계에 대해 기록한 소정의 운동정보 데이터베이스를 통해 학습 될 수 있으며, 이것을 바탕으로 동물의 운동 의도를 인식할 수 있게 된다.Therefore, the binary decision intention recognizer 110b-1 can be learned through a predetermined exercise information database recorded about the relationship between the multiple nerve physiological signals and the animal's motion intention, and based on this, do.

여기서, 운동정보 데이터베이스는 서포트 벡터 머신을 훈련하기 위하여 실험을 통해 제작된 뇌신경신호, 생리신호의 특징 벡터와 해당 특징 벡터에 대응하는 동물의 움직임과 정지 여부와, 좌우 움직임 방향 등을 스칼라 값으로 기록한 데이터 벡터이다.Herein, the exercise information database is a database for recording support information such as a cranial nerve signal, a physiological signal feature vector, an animal movement and a stop motion, and a left / right movement direction corresponding to the feature vector, Data vector.

[수학식 10]&Quot; (10) "

Figure 112017030378610-pat00126
Figure 112017030378610-pat00126

상기 수학식 10에서,

Figure 112017030378610-pat00127
는 벡터의 전치(transpose) 연산이고,
Figure 112017030378610-pat00128
는 커널 함수이다.In Equation (10)
Figure 112017030378610-pat00127
Is a transpose operation of the vector,
Figure 112017030378610-pat00128
Is a kernel function.

상기 수학식 10에서 상기 커널 함수는 다음의 수학식 11과 같이 결정되는 RBF(radial-basis function)이다.In Equation (10), the kernel function is a radial-basis function (RBF) determined by Equation (11).

[수학식 11]&Quot; (11) "

Figure 112017030378610-pat00129
Figure 112017030378610-pat00129

상기 수학식 11에서,

Figure 112017030378610-pat00130
는 특징벡터
Figure 112017030378610-pat00131
의 표준편차이다.In Equation (11)
Figure 112017030378610-pat00130
Is a feature vector
Figure 112017030378610-pat00131
Lt; / RTI >

해독부(110b)에서 회귀부(110b-2)는 서포트 벡터 회귀(SVR)를 기반으로 하여 공포 수준, 경계심 수준, 피로도, 워크로드(workload), 졸림 수준, 통증 수준을 연속 실수 값으로 출력할 수 있게 된다.In the decryption unit 110b, the regression unit 110b-2 outputs the fear level, the alertness level, the fatigue, the workload, the sleepiness level, and the pain level as continuous real values based on the support vector regression .

여기서 서포트 벡터 회귀(SVR)는 상기 분류된 특징 벡터를

Figure 112017030378610-pat00132
이라고 하고, 그에 해당하는 다중 신경생리신호의 관계에 대한 소정의 생체정보 데이터베이스로부터 얻게 되는 공포 수준, 경계심 수준, 피로도, 워크로드, 졸림 수준, 통증 수준의 지표 값을
Figure 112017030378610-pat00133
이라 할 때 종속변수
Figure 112017030378610-pat00134
로 부터 최대 적절한 양수
Figure 112017030378610-pat00135
만큼의 편차(deviation)를 갖는 함수
Figure 112017030378610-pat00136
를 구하는 알고리즘이다.Wherein the support vector regression (SVR)
Figure 112017030378610-pat00132
And the index values of fear level, beating level, fatigue, workload, sleepiness level, and pain level obtained from a predetermined biometric information database on the relationship of the corresponding multiple neurophysiological signals
Figure 112017030378610-pat00133
In this case,
Figure 112017030378610-pat00134
The most appropriate positive number from
Figure 112017030378610-pat00135
A function having a deviation of
Figure 112017030378610-pat00136
.

여기서, 생체정보 데이터베이스는 서포트 벡터 회귀 알고리즘을 훈련하기 위하여 실험을 통해 제작된 뇌신경신호, 생리신호의 특징 벡터와 해당 특징벡터에 대응하는 공포 수준, 경계심 수준, 피로도, 워크로드, 졸림 수준, 통증 수준의 지표 값 등을 기록한 데이터 벡터이다.Here, the biometric information database includes a brain nerve signal, a feature vector of a physiological signal, a fear level corresponding to the feature vector, an alertness level, fatigue, a workload, a sleepiness level, a pain level And the like.

따라서 회귀부(110b-2)는 상기 다중 신경생리신호와 동물의 운동 의도의 관계에 대해 기록한 소정의 생체정보 데이터베이스를 통해 학습될 수 있으며, 이것을 바탕으로 동물의 생체 상태를 점수화 할 수 있게 된다.Therefore, the regression unit 110b-2 can be learned through a predetermined biometric information database recorded about the relationship between the multi-nervous physiological signals and the motion intention of the animal, and based on this, the biological state of the animal can be scored.

도 6은 상기 신호 분석부(110)로부터 해독된 운동 의도와 생체 상태를 입력받아 보상 혹은 처벌 자극의 종류를 결정하게 되는 제어부(120)를 도시된 것으로, 제어부(120)는 보상 혹은 처벌 결정부(120a)와 자극생성 및 제시부(120b)를 포함한다.6 is a block diagram of a controller 120 that receives a decoded motion intention and a living body state from the signal analyzer 110 to determine a type of compensation or punishment stimulus. (120a) and a stimulus generation and presentation unit (120b).

상기 보상/처벌 결정부(120a)에서는 사용자가 행동제어 명령 지시부(130)를 통해 입력한 명령 자극과 일치하는 운동 의도를 보였을 경우 보상을 준다고 결정할 수 있으며, 또한 상기 명령자극과 불일치하는 운동 의도를 보였을 경우 처벌을 준다고 결정할 수 있다.The compensation / punishment determining unit 120a can determine that compensation is provided when the user shows the motion intention that matches the command stimulus inputted through the behavior control instruction instructing unit 130, and also determines the motion intention that does not match the command stimulus If you see it, you can decide to punish.

또한 보상 혹은 처벌 결정부(120a)에서는 생체 상태에 따라 보상을 결정할 수 있으며, 예를 들어 피로도 또는 워크로드 또는 통증 수준이 문턱치(threshold)보다 높았을 경우 보상을 준다고 결정할 수 있다.In addition, the compensation or punishment decision unit 120a can determine compensation according to the living body state, for example, it can be determined that the compensation is provided when the fatigue or workload or the pain level is higher than the threshold.

자극생성 및 제시부(120b)는 상기 보상 혹은 처벌 결정부(120a)에서 결정된 보상 혹은 처벌 사항에 따라 자극 신호 파라미터로 전기 펄스의 폭, 진폭, 주파수, 지속시간 중 선택되는 어느 하나 이상을 조절하여 동물에게 보상 자극 및 처벌 자극을 준다. The stimulus generation and presentation unit 120b adjusts one or more selected among the width, amplitude, frequency, and duration of the electric pulse as the stimulus signal parameter according to the compensation or punishment determined by the compensation or punishment decision unit 120a, To stimulate reward and punishment.

도 7은 동물의 운동의도와 생체 상태를 사용자가 파악할 수 있도록 하는 출력부와 사용자가 동물에게 행동 제어 명령을 지시할 수 있도록 하는 행동제어 명령 지시부(130)를 도시한 것이며, 이 행동제어 명령 지시부(130)는 입력부(130a), 출력부(130b), 디스플레이부(130c), 명령 자극 제시부(130d)로 구성되어 있다.7 shows an output unit for allowing the user to grasp the movement state and the living body state of the animal and a behavior control instruction instruction unit 130 for allowing the user to designate a behavior control instruction to the animal, The display unit 130 includes an input unit 130a, an output unit 130b, a display unit 130c, and a command stimulus suggestion unit 130d.

상기 출력부(130b)에서는 신호 분석부(110)로부터 운동 의도와 생체 상태에 대한 정보를 입력받아 디스플레이부(130c)를 통해 사용자에게 출력할 수 있다. 예를 들어 사용자는 출력부(130b)와 디스플레이부(130c)를 통해 동물의 운동 의도에 대한이진 값 결과와 생체 상태에 대한 점수를 확인할 수 있다.The output unit 130b receives information on the motion intention and the living body state from the signal analysis unit 110 and outputs the information to the user through the display unit 130c. For example, the user can check the result of the binary value of the animal's motion intention and the score of the living body through the output unit 130b and the display unit 130c.

이렇게 디스플레이부(130c)를 통해 동물의 운동 의도 및 생체 상태에 정보를 확인한 사용자는 동물의 행동 및 생체 상태의 제어에 대한 명령을 입력부(130a)를 통해 입력한다. 구체적으로 상기 입력부(130a)에서는 사용자로부터 행동 명령에 관한 입력을 수동으로 행동 및 생체 상태 제어에 대한 지시 받으며, 예를 들어 조이스틱, 리모컨 등의 형태로 동물 행동에 대한 사용자 명령을 수동으로 입력 받을 수 있다.The user who has confirmed the information on the motion intention and the living state of the animal through the display unit 130c inputs a command for controlling the behavior and the living state of the animal through the input unit 130a. Specifically, the input unit 130a receives an instruction for manually inputting a behavior command from a user, and receives a user command for an animal behavior in the form of, for example, a joystick or a remote controller. have.

상기 입력부(130a)에서 입력받은 사용자 명령은 명령 자극 제시부(130d)로 전달되어 사용자로부터 받은 행동 및 상체상태 제어 명령을 전기 자극을 통해 동물에게 제시할 수 있다. 예를 들어 왼쪽, 오른쪽의 방향에 관한 명령을 사용자로부터 받았다고 하면 명령 자극 제시부(130d)는 각각의 명령을 동물의 대뇌의 감각 영역의 오른쪽, 왼쪽에 전기 자극을 줌으로써 명령 자극을 제시할 수 있다.The user command received from the input unit 130a may be transmitted to the command stimulus presentation unit 130d to present the behavior and the body condition control command received from the user to the animal through electric stimulation. For example, if a command regarding the left and right directions is received from the user, the command stimulus suggestion unit 130d can present command stimuli by giving electrical stimulation to the right and left sides of the sensory region of the animal's cerebrum.

또한 명령 자극 제시부 (130d)는 사용자로부터 받은 생체 상태 조절 자극 명령을 전기 자극을 통해 제시하는 기능을 가질 수 있다. 예를 들어, 동물의 공포 수준이 높아져 있다고 판단되면 편도체(amygdala)의 억제를 조절하는 변연계 아래 피질(infralimbic cortex)에 전기 자극을 가함으로써 공포 수준을 낮춰주는 조절을 할 수 있다.Also, the command stimulus suggestion unit 130d may have a function of presenting a biological condition adjustment stimulation command received from the user through electrical stimulation. For example, if the level of fear of an animal is determined to be elevated, electrical stimulation may be applied to the infralimbic cortex to regulate the inhibition of the amygdala, thereby lowering the level of fear.

또한, 척수(spinal cord)의 등쪽뿔(dorsal horn)에서 통각 신호를 실시간으로 인지하여 심한 통증이 유발된다고 판단되면 수관주변 회색질(periaqueductal gray, PAG)에 전기 자극을 가함으로써 통각을 억제할 수 있다.In addition, when the sensory signal is recognized in real time in the dorsal horn of the spinal cord and it is judged that severe pain is caused, the pain can be suppressed by applying electrical stimulation to the periaqueductal gray (PAG) .

본 발명에 따른 동물 행동 제어 시스템은 컴퓨터로 읽을 수 있는 기록매체에 컴퓨터가 읽을 수 있는 코드로서 구현하는 것이 가능하다. 상기 컴퓨터가 읽을 수 있는 기록매체는 컴퓨터 시스템에 의하여 읽힐 수 있는 데이터가 저장되는 모든 종류의 기록 장치를 포함한다. 기록매체의 예로는 ROM(Read Only Memory), RAM(Random Access Memory), 광학 디스크, 자기 테이프, 플로피 디스크, 하드 디스크, 비휘발성 메모리 등을 포함한다.The animal behavior control system according to the present invention can be implemented as computer-readable code on a computer-readable recording medium. The computer-readable recording medium includes all kinds of recording apparatuses in which data that can be read by a computer system is stored. Examples of the recording medium include ROM (Read Only Memory), RAM (Random Access Memory), optical disk, magnetic tape, floppy disk, hard disk, nonvolatile memory and the like.

또한 컴퓨터가 읽을 수 있는 기록매체는 네트워크로 연결된 컴퓨터 시스템에 분산되어 분산방식으로 컴퓨터가 읽을 수 있는 코드가 저장되고 실행될 수 있다.The computer-readable recording medium may also be distributed over a networked computer system so that computer readable code can be stored and executed in a distributed manner.

이상 본 발명의 동물 행동 제어 시스템의 구성 및 이를 이용한 동물 행동의 제어 방법에 대해 구체적으로 설명하였으나, 이는 예시에 불과한 것으로 본 발명의 기술적 사상이 벗어나지 않는 범주 내에서 여러 가지 변형이 가능함을 물론이다.While the present invention has been particularly shown and described with respect to the preferred embodiments thereof, it will be understood by those skilled in the art that various changes and modifications may be made without departing from the spirit and scope of the invention.

100 : 신호 감지 및 기록부
100a : 뇌신경신호 기록 전극
100b : 생리신호 기록 전극
100c : 증폭부
110 : 신호 분석부
110a : 특징벡터 추출부
110a-1 : 뇌신경신호 특징벡터 추출부
110a-2 : 신경신호 특징벡터 추출부
110b : 해독부
110b-1 : 이진 결정 의도 인식기
110b-2 : 회귀부
120 : 제어부
120a : 보상 혹은 처벌 결정부
120b : 자극 생성 및 제시부
130 : 행동제어 명령 지시부
130a : 입력부
130b : 출력부
100: Signal detection and recording unit
100a: cranial nerve signal recording electrode
100b: physiological signal recording electrode
100c:
110: Signal Analysis Unit
110a: a feature vector extracting unit
110a-1: brain nerve signal feature vector extracting unit
110a-2: a neural signal characteristic vector extracting unit
110b:
110b-1: binary decision intention recognizer
110b-2:
120:
120a: compensation or punishment decision section
120b: stimulus generation and presentation unit
130: Behavior control command indicator
130a:
130b:

Claims (24)

삭제delete 삭제delete 삭제delete 삭제delete 삭제delete 삭제delete 삭제delete 삭제delete (a) 대상 동물로부터 동물의 뇌신경신호 및 생리신호를 획득하는 단계;
(b) 획득된 뇌신경신호 및 생리신호를 분석하여 뇌신경신호 특징벡터 및 생리신호 특징벡터를 추출하는 단계;
(c) 추출된 뇌신경신호 특징벡터 및 생리신호 특징벡터를 해독하여 동물의 운동 의도 및 생체 상태를 파악하는 단계;
(d) 파악된 동물의 운동 의도 및 생체 상태를 기반으로 보상 혹은 처벌여부를 결정하거나 임무 지속 가능 여부를 결정하여 동물의 행동을 제어하기 위한 전기 자극 파라미터를 조절하는 단계; 및
(e) 결정된 전기 자극 파라미터에 따른 전기자극 파형을 특정 대뇌 영역에 전달하여, 동물의 행동을 제어하는 단계;를 포함하되,
상기 (c) 단계는,
상기 (b) 단계에서 분류된 특징벡터를 통해 서포터 벡터 머신(Support Vector Machine, SVM)을 기반으로 서포터 벡터 머신 모델을 생성하여 동물의 운동 의도를 출력하는 단계; 및
서포트 벡터 회귀(Support Vector Regression, SVR)을 기반으로 하여 동물의 생체 상태를 출력하는 단계;를 포함하며,
상기 서포터 벡터 머신은 분류된 특징벡터를
Figure 112018094740120-pat00212
이라고 할 때, 아래의 수학식 10과 수학식 11에 의해 최적의 가중벡터
Figure 112018094740120-pat00213
와 최적의 바이어스
Figure 112018094740120-pat00214
를 구하여 서포트 벡터 머신 모델을 생성하는 것을 특징으로 하는 동물 행동 제어 방법.
[수학식 10]
Figure 112018094740120-pat00215

(수학식 10에서,
Figure 112018094740120-pat00216
는 벡터의 전치(transpose) 연산이고,
Figure 112018094740120-pat00217
는 커널 함수임)
상기 수학식 10에서 상기 커널 함수는 다음의 수학식 11와 같이 결정되는 RBF(radial-basis function)임
[수학식 11]
Figure 112018094740120-pat00218

(수학식 11에서,
Figure 112018094740120-pat00219
는 특징벡터
Figure 112018094740120-pat00220
의 표준편차임)
(a) obtaining a neural signal and a physiological signal of an animal from a subject animal;
(b) extracting a cranial nerve signal characteristic vector and a physiological signal characteristic vector by analyzing the obtained cranial nerve signal and physiological signal;
(c) decoding the extracted cranial nerve signal characteristic vector and physiological signal characteristic vector to determine the motion intention and the biological state of the animal;
(d) adjusting electrical stimulation parameters to control the behavior of the animal by determining whether compensation or punishment based on the animal's exercise intention and the biological state of the identified animal, or determining the sustainability of the mission; And
(e) delivering an electric stimulation waveform according to the determined electric stimulation parameter to a specific cerebral region to control the behavior of the animal,
The step (c)
Generating a supporter vector machine model based on a support vector machine (SVM) through the feature vectors classified in the step (b) and outputting the intention of the animal; And
Outputting the biological state of the animal on the basis of a support vector regression (SVR)
The supporter vector machine uses the sorted feature vector
Figure 112018094740120-pat00212
, The optimal weight vector < RTI ID = 0.0 >
Figure 112018094740120-pat00213
And optimum bias
Figure 112018094740120-pat00214
To generate a support vector machine model.
&Quot; (10) "
Figure 112018094740120-pat00215

(In equation (10)
Figure 112018094740120-pat00216
Is a transpose operation of the vector,
Figure 112018094740120-pat00217
Is a kernel function)
In Equation (10), the kernel function is a radial-basis function (RBF) determined by Equation (11)
&Quot; (11) "
Figure 112018094740120-pat00218

(In the equation (11)
Figure 112018094740120-pat00219
Is a feature vector
Figure 112018094740120-pat00220
Lt; / RTI >
제9항에 있어서,
상기 뇌신경신호는 세포 외 기록 신호(extracellular recording signal), 국소뇌전위(local field potential), 뇌피질 뇌파(subdural electrocortigram, ECoG) 및 뇌전도(electroencephalography, EEG) 중에서 선택된 어느 하나 이상인 것을 특징으로 하는 동물 행동 제어 방법.
10. The method of claim 9,
Wherein the brain signal is at least one selected from an extracellular recording signal, a local field potential, a subdural electrocortigram (ECoG), and an electroencephalography (EEG) Control method.
제9항에 있어서,
상기 생리신호는 피부표면온도(skin temperature, SKT), 광혈류량(photo-plethysmography, PPG), 심전도(electrocardiography, ECG), 피부 전도도 반응(skin conductance response, SCR) 및 호흡 신호 중에서 선택되는 어느 하나 이상인 것을 특징으로 하는 동물 행동 제어 방법.
10. The method of claim 9,
The physiological signal may be any one or more selected from skin temperature (SKT), photo-plethysmography (PPG), electrocardiography (ECG), skin conductance response ≪ / RTI >
제9항 또는 제10항에 있어서,
상기 (b) 단계는,
상기 획득된 뇌신경신호를 통해 스파이크 트레인(spike train), 파워 스펙트럼 밀도(Power Spectrum Density), 대뇌 영역 간 기능적 연결도(functional connectivity), 및 대뇌 영역 간 효과적 연결도(effective connectivity) 중 선택되는 어느 하나의 뇌신경신호 특징벡터를 추출하는 것을 특징으로 하는 동물 행동 제어 방법.
11. The method according to claim 9 or 10,
The step (b)
The obtained neural signal is used to select either one of spike train, power spectral density, functional connectivity between cerebral regions, and effective connectivity between cerebral regions. And extracting a cranial nerve signal feature vector of the animal.
제12항에 있어서,
상기 파워 스펙트럼 밀도는 뇌신경신호 중 어느 하나의 신호를 M차원(M은 1 이상의 정수)의 벡터
Figure 112017030378610-pat00137
로 하여,
하기 수학식 1을 통해 분석하는 푸리에 변환(Fourier transform);
하기 수학식 2를 통해 분석하는 단기간 푸리에 변환(short-time Fourier transform); 및
하기 수학식 3을 통해 분석하는 웨이브릿 변환(wavelet transform); 중의 어느 하나를 이용하여 특정 시간-주파수 대역의 파워를 얻어 뇌신경신호 특징벡터를 추출하는 것을 특징으로 하는 동물 행동 제어 방법.
[수학식 1]
Figure 112017030378610-pat00138

(수학식 1에서,
Figure 112017030378610-pat00139
는 주파수를 지칭하는 변수
Figure 112017030378610-pat00140
에 대하여
Figure 112017030378610-pat00141
를 의미하는 각주파수임)
[수학식 2]
Figure 112017030378610-pat00142

(수학식 2에서,
Figure 112017030378610-pat00143
은 특정한 길이의 특정한 모양의 시간 윈도우 벡터임)
[수학식 3]
Figure 112017030378610-pat00144

(수학식 3에서
Figure 112017030378610-pat00145
는 웨이블릿 함수의 공액복소수(complex conjugate),
Figure 112017030378610-pat00146
는 정규화 파라미터,
Figure 112017030378610-pat00147
는 특정 스케일,
Figure 112017030378610-pat00148
는 목표 레이턴시임)
13. The method of claim 12,
The power spectral density is obtained by multiplying one of the signals of the cerebral nerve signal by an M-dimensional (M is an integer of 1 or more)
Figure 112017030378610-pat00137
As a result,
A Fourier transform analyzing by the following equation (1);
A short-time Fourier transform analyzing by the following equation (2); And
A wavelet transform to be analyzed through Equation (3); To obtain a power of a specific time-frequency band and to extract a brain nerve signal feature vector.
[Equation 1]
Figure 112017030378610-pat00138

(In the equation (1)
Figure 112017030378610-pat00139
Is a variable that refers to frequency
Figure 112017030378610-pat00140
about
Figure 112017030378610-pat00141
≪ / RTI >
&Quot; (2) "
Figure 112017030378610-pat00142

(In the equation (2)
Figure 112017030378610-pat00143
Is a time window vector of a particular shape of a particular length)
&Quot; (3) "
Figure 112017030378610-pat00144

(3)
Figure 112017030378610-pat00145
Is a complex conjugate of the wavelet function,
Figure 112017030378610-pat00146
Is a normalization parameter,
Figure 112017030378610-pat00147
Lt; / RTI >
Figure 112017030378610-pat00148
Is the target latency)
제12항에 있어서,
상기 대뇌 영역 간 기능적 연결성은 뇌신경신호 중 어느 둘의 뇌신경신호를 한 쌍의 M차원(M은 1 이상의 정수)의 벡터
Figure 112017030378610-pat00149
,
Figure 112017030378610-pat00150
로 하여,
하기 수학식 4을 통해 분석하는 간섭성(coherence);
하기 수학식 5를 통해 분석하는 위상 동기 값(phase locking value, PLV); 및
하기 수학식 6을 통해 분석하는 weighted phase lag index(wPLI); 중의 어느 하나 이상을 이용하여 한 쌍의 뇌신경 영역 간의 기능적 연결성을 확인하여 뇌신경신호 특징벡터를 추출하는 것을 특징으로 하는 동물 행동 제어 방법.
[수학식 4]
Figure 112017030378610-pat00151

(수학식 4에서,
Figure 112017030378610-pat00152
Figure 112017030378610-pat00153
Figure 112017030378610-pat00154
Figure 112017030378610-pat00155
의 진폭,
Figure 112017030378610-pat00156
는 두 시계열의 상대적 위상차임)
[수학식 5]
Figure 112017030378610-pat00157

(수학식 5에서,
Figure 112017030378610-pat00158
Figure 112017030378610-pat00159
는 각각 시계열 벡터
Figure 112017030378610-pat00160
Figure 112017030378610-pat00161
에 힐버트 변환을 적용하여 얻은 순간 위상(instantaneous phase)임)
[수학식 6]
Figure 112017030378610-pat00162

(수학식 6에서,
Figure 112017030378610-pat00163
은 기대값 연산,
Figure 112017030378610-pat00164
은 복소수의 허수 항만을 취하는 연산자,
Figure 112017030378610-pat00165
는 두 시계열 벡터
Figure 112017030378610-pat00166
Figure 112017030378610-pat00167
의 교차 스펙트럼(cross spectrum)임)
13. The method of claim 12,
The functional connectivity between the cerebral regions may be determined by connecting either of two neural signals of the cerebral nerve signal to a pair of M dimensions (M is an integer of 1 or more)
Figure 112017030378610-pat00149
,
Figure 112017030378610-pat00150
As a result,
Coherence to be analyzed through Equation (4);
A phase locking value (PLV) analyzed through Equation (5); And
A weighted phase lag index (wPLI) analyzed through Equation (6); And extracting a cranial nerve signal feature vector by confirming the functional connectivity between the pair of cranial nerve regions.
&Quot; (4) "
Figure 112017030378610-pat00151

(In Equation 4,
Figure 112017030378610-pat00152
Wow
Figure 112017030378610-pat00153
The
Figure 112017030378610-pat00154
Wow
Figure 112017030378610-pat00155
The amplitude,
Figure 112017030378610-pat00156
Is the relative phase of the two time series)
&Quot; (5) "
Figure 112017030378610-pat00157

(In equation (5)
Figure 112017030378610-pat00158
and
Figure 112017030378610-pat00159
Respectively,
Figure 112017030378610-pat00160
and
Figure 112017030378610-pat00161
Which is the instantaneous phase obtained by applying the Hilbert transform to the < RTI ID = 0.0 >
&Quot; (6) "
Figure 112017030378610-pat00162

(In Equation 6,
Figure 112017030378610-pat00163
The expected value calculation,
Figure 112017030378610-pat00164
Is an operator that takes an imaginary number of complex numbers,
Figure 112017030378610-pat00165
The two time series vectors
Figure 112017030378610-pat00166
and
Figure 112017030378610-pat00167
≪ / RTI >
제12항에 있어서,
상기 대뇌 영역 간 방향성 있는 연결성 분석은 뇌신경신호 중 어느 둘의 뇌신경신호를 한 쌍의 M차원(M은 1 이상의 정수)의 벡터
Figure 112017030378610-pat00168
,
Figure 112017030378610-pat00169
로 하여,
수학식 7을 통해 분석하는 그랜저 인과관계(Granger's causality);
수학식 8을 통해 분석하는 PDC(Partial Directed Coherence); 및
수학식 9를 통해 분석하는 DTF(Directed Transfer Function); 중 어느 하나를 이용하여 한 쌍의 뇌신경 영역 간의 방향성이 있는 연결성을 확인하여 뇌신경신호 특징벡터를 추출하는 것을 특징으로 하는 동물 행동 제어 방법.
[수학식 7]
Figure 112017030378610-pat00170

(수학식 7에서,
Figure 112017030378610-pat00171
Figure 112017030378610-pat00172
Figure 112017030378610-pat00173
의 과거 값만을 이용하여 자동 회귀 모델을 이용해 추정할 때 생기는 에러,
Figure 112017030378610-pat00174
Figure 112017030378610-pat00175
Figure 112017030378610-pat00176
Figure 112017030378610-pat00177
의 과거값을 모두 사용하여 자동 회귀 모델을 이용해
Figure 112017030378610-pat00178
Figure 112017030378610-pat00179
의 현재 값을 추정할 때 생기는 에러임)
[수학식 8]
Figure 112017030378610-pat00180

(수학식 8에서,
Figure 112017030378610-pat00181
Figure 112017030378610-pat00182
Figure 112017030378610-pat00183
을 자동 회귀 모델의 계수 행렬을 푸리에 변환한 것이며,
Figure 112017030378610-pat00184
Figure 112017030378610-pat00185
행렬의
Figure 112017030378610-pat00186
번째 열벡터임)
[수학식 9]
Figure 112017030378610-pat00187

(수학식 9에서,
Figure 112017030378610-pat00188
는 상기 수학식 8에서 나온
Figure 112017030378610-pat00189
의 역행렬임)
13. The method of claim 12,
The directional connectivity analysis between the cerebral regions can be performed by connecting any two of the cerebral nerve signals of the cerebral nerve signal to a pair of M-dimensional (M is an integer of 1 or more)
Figure 112017030378610-pat00168
,
Figure 112017030378610-pat00169
As a result,
Granger's causality analyzing through equation (7);
PDC (Partial Directed Coherence) to be analyzed through Equation (8); And
DTF (Directed Transfer Function) analyzed through Equation (9); And extracting a cranial nerve signal characteristic vector by confirming directional connectivity between the pair of cranial nerve regions.
&Quot; (7) "
Figure 112017030378610-pat00170

(In equation (7)
Figure 112017030378610-pat00171
silver
Figure 112017030378610-pat00172
of
Figure 112017030378610-pat00173
The error of the estimation using the automatic regression model using only the past value of the error,
Figure 112017030378610-pat00174
silver
Figure 112017030378610-pat00175
of
Figure 112017030378610-pat00176
and
Figure 112017030378610-pat00177
Using all of the past values of the regression model
Figure 112017030378610-pat00178
and
Figure 112017030378610-pat00179
≪ / RTI >< RTI ID = 0.0 >
&Quot; (8) "
Figure 112017030378610-pat00180

(In equation (8)
Figure 112017030378610-pat00181
The
Figure 112017030378610-pat00182
and
Figure 112017030378610-pat00183
Is the Fourier transform of the coefficient matrix of the autoregressive model,
Figure 112017030378610-pat00184
The
Figure 112017030378610-pat00185
Matrix
Figure 112017030378610-pat00186
Th column vector)
&Quot; (9) "
Figure 112017030378610-pat00187

(9)
Figure 112017030378610-pat00188
Is expressed by Equation (8)
Figure 112017030378610-pat00189
Lt; / RTI >
제9항 또는 제10항에 있어서,
상기 (b) 단계는,
상기 측정된 생리신호를 통해 피부표면온도(skin temperature, SKT)를 이용한 특징벡터, 광혈류량(photo-plethysmography, PPG)을 이용한 특징벡터, 심전도(electrocardiography, ECG)를 이용한 특징벡터, 피부 전도도 반응(skin conductance response, SCR)을 이용한 특징벡터 및 호흡 신호를 이용한 특징벡터 중에서 선택되는 어느 하나의 생리신호 특징벡터를 추출하는 것을 특징으로 하는 동물 행동 제어 방법.
11. The method according to claim 9 or 10,
The step (b)
A feature vector using skin temperature (SKT), a feature vector using photo-plethysmography (PPG), a feature vector using electrocardiography (ECG), skin conductivity response skin behavior response (SCR), and a feature vector using a respiration signal are extracted from the physiological signal feature vector.
제16항에 있어서,
상기 피부표면온도를 이용한 특징벡터는 피부표면온도의 평균, 피부 표면 온도 변화를 나타내는 1차 미분 값의 평균, 피부표면온도의 최대값, 및 0.1 내지 0.2 헤르츠(Hz)의 주파수 대역 스펙트럼 파워 중 어느 하나 이상을 포함하는 생리신호 특징벡터를 추출하는 것을 특징으로 하는 동물 행동 제어 방법.
17. The method of claim 16,
The characteristic vector using the surface temperature of the skin may be any one of an average of the skin surface temperature, an average of the first derivative indicating the skin surface temperature change, a maximum value of the skin surface temperature, and a spectral power of 0.1 to 0.2 hertz (Hz) And extracting physiological signal feature vectors including at least one of the physiological signal feature vectors.
제16항에 있어서,
상기 피부 전도도 반응을 이용한 특징벡터는 피부 전도도 반응 진폭의 평균값(mean of SCR amplitude), 피부 전도도 반응의 음수 값 평균값, 피부 전도도 반응의 표준 편차, 피부 전도도 반응의 1차 미분 신호 평균값(mean of 1st derivative of SCR), 피부 전도도 반응의 2차 미분 신호 평균값(mean of 2nd derivative of SCR), 및 피부 전도도 반응의 전체 샘플 개수 대비 음의 값을 가지는 샘플의 비율 중 어느 하나 이상을 포함하는 생리신호 특징벡터를 추출하는 것을 특징으로 하는 동물 행동 제어 방법.
17. The method of claim 16,
The feature vector using the skin conductivity response includes a mean of SCR amplitude, a mean value of the skin conductivity response, a standard deviation of the skin conductance response, a mean of the first derivative signal of the skin conductance response a derivative of SCR, a mean of a second derivative of the skin conductivity response, and a ratio of a sample having a negative value relative to the total number of samples of the skin conductivity response. And the vector is extracted.
제16항에 있어서,
상기 광혈류량을 이용한 특징벡터는 심박수의 평균, 표준편차, 심박 수 변이(Heart rate variability)의 평균, 광혈류량의 표준편차, 광혈류량의 1차 미분의 평균, 광혈류량의 1차 미분의 표준편차, 및 광혈류량의 저주파 성분과 고주파 성분의 비율중 어느 하나 이상을 포함하는 생리신호 특징벡터를 추출하는 것을 특징으로 하는 동물 행동 제어 방법.
17. The method of claim 16,
The feature vector using the optical blood flow can be expressed as the mean of the heart rate, the standard deviation, the average of the heart rate variability, the standard deviation of the optical blood flow, the average of the first derivative of the optical blood flow, , And a ratio of a low-frequency component and a high-frequency component of the optical blood flow to the physiological signal feature vector.
제9항에 있어서,
상기 (b) 단계는,
획득된 동물의 뇌신경신호 및 생리신호를 주성분 분석법(principal component analysis, PCA)과 독립성분 분석법(independent component analysis, ICA)을 수행하여 뇌신경신호 및 생리신호에 관한 기저벡터를 추출하는 단계;
상기 기저벡터로부터 특징벡터를 추출하는 단계; 및
상기 특징벡터를 서브클래스(subclass)로 분류하는 클러스터링(clustering) 단계;를 더 포함하여 특징벡터를 추출하는 것을 특징으로 하는 동물 행동 제어 방법.
10. The method of claim 9,
The step (b)
Extracting a base vector related to a cranial nerve signal and a physiological signal by performing principal component analysis (PCA) and independent component analysis (ICA) on the nerve signal and physiological signal of the obtained animal;
Extracting a feature vector from the basis vector; And
And a clustering step of classifying the feature vectors into subclasses, wherein the feature vectors are extracted.
삭제delete 삭제delete 제9항에 있어서,
상기 서포트 벡터 회귀는 분류된 특징벡터를
Figure 112018094740120-pat00199
이라고 하고, 상기 분류된 특징벡터에 해당하는 뇌신경신호 및 신경신호의 관계로부터 얻게 되는 동물의 생체 상태의 지표 값을
Figure 112018094740120-pat00200
이라 할 때, 종속변수
Figure 112018094740120-pat00201
로부터 최대 양수
Figure 112018094740120-pat00202
만큼의 편차(deviation)를 갖는 함수
Figure 112018094740120-pat00203
를 구하는 것을 특징으로 하는 동물 행동 제어 방법.
10. The method of claim 9,
The support vector regression may be performed by
Figure 112018094740120-pat00199
, And an indicator value of an animal's biological condition obtained from a relationship between a cranial nerve signal and a nerve signal corresponding to the classified characteristic vector
Figure 112018094740120-pat00200
, The dependent variable
Figure 112018094740120-pat00201
The maximum positive
Figure 112018094740120-pat00202
A function having a deviation of
Figure 112018094740120-pat00203
Of the animal behavior.
제9항에 있어서,
상기 (e) 단계에서 상기 동물의 특정 대뇌 영역은,
자발적인 방향성 있는 움직임 자극을 가하기 위한 대뇌 체성 감각 영역;
강제적인 방향성 있는 움직임 자극을 가하기 위한 바닥핵(basal ganglia)의 흑색선조체 경로(nigrostriatal pathway) 또는 시상(thalamus)의 체성감각 통로인 복측후내측핵 (ventral posteromedial nucleus);
보상 자극을 가하기 위한 뇌측전다발 영역; 및
처벌 자극을 가하기 위한 편도체(amygdala);를 포함하는 것을 특징으로 하는 동물 행동 제어 방법.
10. The method of claim 9,
The specific cerebral region of the animal in the step (e)
A cerebral somatosensory area for applying spontaneous directional motion stimulation;
Ventral posteromedial nucleus, a somatic sensory pathway of the nigrostriatal pathway or thalamus of the basal ganglia for applying forced directional motor stimulation;
A brainstem bundle region for applying compensatory stimulation; And
And an amygdala for applying a punitive stimulus.
KR1020170039216A 2017-03-28 2017-03-28 Multiple neuro-physiological signal based animal behavior control system and method thereof KR101948417B1 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
KR1020170039216A KR101948417B1 (en) 2017-03-28 2017-03-28 Multiple neuro-physiological signal based animal behavior control system and method thereof

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
KR1020170039216A KR101948417B1 (en) 2017-03-28 2017-03-28 Multiple neuro-physiological signal based animal behavior control system and method thereof

Publications (2)

Publication Number Publication Date
KR20180109451A KR20180109451A (en) 2018-10-08
KR101948417B1 true KR101948417B1 (en) 2019-02-14

Family

ID=63864778

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1020170039216A KR101948417B1 (en) 2017-03-28 2017-03-28 Multiple neuro-physiological signal based animal behavior control system and method thereof

Country Status (1)

Country Link
KR (1) KR101948417B1 (en)

Families Citing this family (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2015145422A1 (en) 2014-03-26 2015-10-01 Scr Engineers Ltd Livestock location system
US11071279B2 (en) 2014-09-05 2021-07-27 Intervet Inc. Method and system for tracking health in animal populations
US10986817B2 (en) 2014-09-05 2021-04-27 Intervet Inc. Method and system for tracking health in animal populations
AU2019261293A1 (en) 2018-04-22 2020-12-10 Vence, Corp. Livestock management system and method
CN112911927B (en) 2018-10-10 2023-06-27 世亚工程设备有限公司 Method and device for drying milk of livestock
CN109566446A (en) * 2018-12-07 2019-04-05 中国人民解放军军事科学院军事医学研究院 The training device and method grasped for primate moving target object
KR102245768B1 (en) * 2019-06-28 2021-04-28 국방과학연구소 Animal Behavior Control Apparatus and System including the same
KR102296838B1 (en) 2019-08-12 2021-09-02 단국대학교 천안캠퍼스 산학협력단 Motion Intention Discrimination Method and Prosthesis Control Method by Linking EEG Signal and Peripheral Nerve Signal
CN110743101B (en) * 2019-10-29 2023-08-25 中国电子科技集团公司第五十四研究所 Animal robot movement behavior recording and analyzing method
KR102439639B1 (en) * 2020-03-24 2022-09-02 비웨이브 주식회사 Method for providing information of major depressive disorders and device for providing information of major depressive disorders using the same
USD990063S1 (en) 2020-06-18 2023-06-20 S.C.R. (Engineers) Limited Animal ear tag
USD990062S1 (en) 2020-06-18 2023-06-20 S.C.R. (Engineers) Limited Animal ear tag
IL275518B (en) 2020-06-18 2021-10-31 Scr Eng Ltd An animal tag
CN113421141A (en) * 2020-07-29 2021-09-21 阿里巴巴集团控股有限公司 Shopping processing method and device based on brain-computer and brain nerve signals and electronic equipment
WO2022113062A1 (en) 2020-11-25 2022-06-02 Scr Engineers Ltd. A system and method for tracing members of an animal population

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002345862A (en) * 2001-05-22 2002-12-03 Hokkaido Technology Licence Office Co Ltd Individually adaptive biological signal passive-movement appliance control system and its control method
JP2006079321A (en) * 2004-09-09 2006-03-23 Nippon Telegr & Teleph Corp <Ntt> Parameter estimation device, method, and program, and recording medium of parameter estimation program
JP2012055679A (en) * 2010-09-10 2012-03-22 Jiaotong Univ Interactive biosensing device and method thereof

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002345862A (en) * 2001-05-22 2002-12-03 Hokkaido Technology Licence Office Co Ltd Individually adaptive biological signal passive-movement appliance control system and its control method
JP2006079321A (en) * 2004-09-09 2006-03-23 Nippon Telegr & Teleph Corp <Ntt> Parameter estimation device, method, and program, and recording medium of parameter estimation program
JP2012055679A (en) * 2010-09-10 2012-03-22 Jiaotong Univ Interactive biosensing device and method thereof

Also Published As

Publication number Publication date
KR20180109451A (en) 2018-10-08

Similar Documents

Publication Publication Date Title
KR101948417B1 (en) Multiple neuro-physiological signal based animal behavior control system and method thereof
CA2524617C (en) System and method of assessment of the efficacy of treatment of neurological disorders using the electroencephalogram
US7706871B2 (en) System and method of prediction of response to neurological treatment using the electroencephalogram
Blinowska et al. Electroencephalography (eeg)
US9302111B2 (en) Modulation and analysis of cerebral perfusion in epilepsy and other neurological disorders
EP2575608B1 (en) Detector for identifying physiological artifacts from physiological signals and method
JP3769023B2 (en) A system for predicting, early detection, warning, prevention or control of changes in brain activity
EP3217869B1 (en) Scoring method based on improved signals analysis
JP2019528104A (en) In-ear sensing system and method for monitoring biological signals
JP2007515200A5 (en)
WO2020132315A1 (en) Systems and methods to detect and treat obstructive sleep apnea and upper airway obstruction
US11154234B2 (en) Method for the management of nociception
CN115486818B (en) Method, system and device for motor nerve pathway detection and quantification
JP6304626B2 (en) Rehabilitation equipment
Paulraj et al. Fractal feature based detection of muscular and ocular artifacts in EEG signals
Saidatul et al. The assessment of developed mental stress elicitation protocol based on heart rate and EEG signals
Dubey et al. Digital analysis of EEG brain signal
Prakash et al. A study on artifacts removal from physiological signals
KR20230136991A (en) Method and device for neuro/biofeedback to stabilize motor cortex of tic-disorders based on biomedcal signals
Hasan et al. Artifact Reduction
Zhang et al. Estimation of mental fatigue based on wavelet packet parameters and kernel learning algorithms
Hapuarachchi Feature selection and artifact removal in sleep stage classification
Qaraqe Patient-Specific Epileptic Seizure Onset Detection via Fused Eeg and Ecg Signals

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