JPS62148642A - System for analysis of living body signal - Google Patents

Abstract

(57) [Summary] This bulletin contains application data before electronic filing, so abstract data is not recorded.

Description

DETAILED DESCRIPTION OF THE INVENTION [Field of Industrial Application] The present invention relates to a biological signal analysis method suitable for the publication of biological phenomena.

[Summary of the Invention] The present invention provides a method for quantifying the temporal correlation of a plurality of living organisms (, No. This makes it easy to determine whether the device is active or active.

(Prior art) Conventionally, action potentials of a living body are directly recorded, such as electrocardiograms and electroencephalograms, or various medical quantities, such as +111 intraluminal pressure, are recorded as electrical signals through appropriate transusers. It is widely practiced to convert and record the recorded biological signals for diagnostics, examinations, etc. In addition, such biological signals are processed by a computer and the old parameters are converted to f4. ? ζ, biological viscosity υ) and diagnosis 1tl[ have been conventionally performed.

For example, in an electroencephalograph, C is the amount of brain waves.
(1Jti, depending on the person's personality and direction゛(、英１・″
The line 11 is marked with a blank color, and is published in JP-A No. 6 o-+.
63t; proposed as No. a8. In addition, as another example, "C" can be used to display various parameters such as the amplitude of one brain wave, frequency, hiss I-gram and amplitude distribution of each frequency, and the appearance rate of average amplitude 11Vj. City 1
There is a slope.

Historically, attempts have been made to detect brain waves at two 'd1 poles and perform brain analysis JJ+ from the correlation between the two signals.

[Problem that the invention seeks to solve]

However, since the source of brain waves is thought to be a multi-dipole, it has been difficult to treat the brain as a volumetric conductor and analyze it using Heckler's method. In addition, in order to judge whether a brain is abnormal or not based on a wide variety of parameters, a certain level of specialized knowledge and experience is required, and it is extremely difficult for non-experts to analyze brain waves. Met. Historically, since the frequency and waveform of brain waves were not constant, it was difficult to set the correlation constant to 1 and obtain parameters that were effective for analysis.

In view of these problems, the purpose of the present invention is to
We are about to provide a biological fK number analysis method that can easily analyze each issue.

[Means for solving problems]

The present invention separates combinations having a predetermined polarity at each instant (11° 1) at the same point in time from a plurality of biological signals of a predetermined length of time, and combines the separated polar set combinations. This is a biological signal analysis method that expresses the temporal correlation of a plurality of biological signals by the appearance rate of at least one of the following.

[Effect]

According to the present invention, the temporal correlation of a plurality of biological signals is quantified using the appearance rate of a predetermined polarity combination.

〔Example〕

Hereinafter, an embodiment in which the biological fg analysis method according to the present invention is applied to brain waves will be described with reference to FIGS. 1 to 5.

The configuration of one embodiment of the present invention is shown in FIG.

In Fig. 4, ζ, a plurality of electrodes (11), (12),
... (1n) is supplied to the multichannel amplifier (3) via a buffer in the electrode box (2). A plurality of output signals of the amplifier (3) are temporally serialized by a multiplexer (4) and are sequentially supplied to an A/D converter (5).

The output data of the A/D converter (5) is supplied to the computer's central processing unit (processor) (6). 10 The memo (January 7) is connected to the setter (6), and the display section (10), recording section (11), and the figure H are connected via the display recording control section (8) and the keyboard control section (9). A magnetic recording section and keyboard 1'' (12), which are omitted, are connected.

The electrode arrangement part (37 (detection part)) on the LJU skin for electroencephalogram detection is internationally standardized, but in this example, the FPI of the left hemisphere, as shown in Fig. 5, was used.

T], 11 sites of C3 and 01, and 4 sites of FP2゜'T''4.C4 and 02 of the right hemisphere are selected. Then ζ
, as shown in Figure A, the names of both the left and right hemispheres, and the names of the posterior three
By selecting the detection site, 11 detection patterns P1 to F) are formed by the fIl frame guidance method. Also, between 113
As shown in Figure 4, one channel of each detection signal is obtained from both the '1゛3 and C3 detection sites in the left hemisphere using the jJi pole induction method (indicated by the broken line), and one channel of detection signals are obtained from both the FPI and O] detection sites. One channel of detection signals is obtained by the bipolar induction method (indicated by solid line C), and a detection pattern P5 of three channels of detection signals is formed.

Similarly, a detection pattern P6 is formed in the right hemisphere.

Next, the functions of the computer in this embodiment will be explained with reference to FIGS. 1 to 3.

The functional block diagram of this embodiment is shown in Figure 1. In FIG. 1, 8 channels of detection signals are inputted to the preprocessing calculation means (22) from the 4-input terminals (211) to (21s) corresponding to the above-mentioned eight detection parts (step 2),
Unnecessary components such as noise (including biological signals from non-target organs, etc.) and fluctuation of the line are removed by filtering, etc., and each detection 1,
7. Preprocessing calculations are performed such as setting a reference DC level or an appropriate reference DC level such that the sum of I+ becomes zero (Step 2). l1ii processing calculation means (2
The output of 8,000 yen in 2) is the pattern setting hand 1 (23
), six detection patterns P1 to P6 as shown in FIG.
It is supplied to the correlation calculation means (24) (suasop(0)).

In this embodiment, the calculation of the temporal correlation of the 3,000-yen contradictory signals of each detection pattern is performed as follows.

For example, three detection sites T3 . . . forming the first detection pattern P1. Each detection signal (i?ij processed) of C3 and FPt (i?ij processed) SA, SR, and Sc has their instantaneous width and polarity changing over time, as shown in Figure 3 A, 13, and C. shall be. These 4g SA, SR
and S. For example, in the case of screening, over a measurement period of several seconds to several minutes.

When sampling is performed simultaneously at a period of ζ, 4 mS,
Outside the frame of each signal at each sampling time point are eight combinations (, 11 to Q8) as shown in Table 1.

Table 1 Correlation '6? Data sorting means (25) of JW' means (24)
) The odor ζ is separated into each set of sample values at arbitrary sampling points of the 13 detection signals SA, SB, and Sc, or into each polarity combination Q[~Q8 at the sampling point in question (step ■). The Kaji set of these separated sample values is (a*, 7. bij, clj) (i=
l −N, j = 1-8. However, N may differ depending on the classification group), and each wave 1qti calculation is performed (
In 261) to (26s), the interpupillary distance (root power) Kij of the detection signal of the detection pattern P1 is first allowed to be 6Q according to the following equation (1).

...(11) From the rule 1-power Kij, the polarity area ratio Mj for each polarity combination Qj determined by the following equation (2) is calculated (step 2).

Depending on the target of analysis, the root power at each sampling point may be calculated as 1'' (ζ with K i j as a constant),
The appearance rate for each polarity combination may be determined.

The output of the correlation calculation means (24) as described above, that is, the polar plane, 4Mj, is supplied to the convergence rate/liJ calculation means (27), and the pattern P is determined by the following equation (3).
A convergence rate U[ of 1 is calculated.

[J1-Ml+Mt...(3
) As is clear from equation (3), the convergence rate is determined when all the detection signals of the detection pattern are 1 and the polarity combination b (Y 1
and] h tendency 'I of polarity combination 07 where all 1 are negative
If the measurement time is p, I, for example, several seconds, the polar area ratios M1 and M1 in the 1111 polarity combinations Ql and Q? are generally not equal.However, if the measurement time is, for example, several seconds When the average efficiency of the pre-processing operation (step ■) described above appears, the bipolar conversion rates M1 and Ml become approximately equal, and the convergence rate U1g in this case is as follows (4 ) is expressed as the formula.

[JiJ! "-2M1'= 2M7...
・(4) Similarly, other five detection patterns P2 to P
The convergence rates U2 to LJ6 of 6 are determined, and the total convergence rate U is calculated as the average of the convergence rates [Jt to U6 of each detection pattern (step 2).

The output of the convergence rate operator 1 (27), that is, the convergence rate U1 of each detection pattern P1 to P ) For example, bar graphs, histodarams, 4Ir line graphs, etc. are displayed on the screen, and the list I format and plot (in the case of changes over time) format are recorded in the record (11) etc. (Step ■) .

In addition, depending on the subject of analysis in 11, the left side of the brain, E I half I
, each 1111 part of C, the ratio/intersection between each rear part, or the ratio/intersection between each front and rear are performed, but in such a case, the same polarity is used as in the case of convergence with Ql and Q7, which are combinations of the same polarity. ,
For example, it is possible to effectively use a specific combination of polar area ratios, etc., such as Q3 and Q5, in which the poles do not necessarily match.

- Analyze the state of brain activity according to the level of mental activity that changes during sleep, wakefulness, etc., based on the convergence rate of each detection pattern obtained as described above, the total convergence rate, etc. can do.

It is generally believed that a normal brain is active under the control of its central part, and in this case, both the convergence rate and the total convergence rate will be large values, and the distribution (frequency) will be in the quotient convergence rate region. Concentrate on. On the other hand, in the case of abnormal electroencephalogram propagation, ectopic dipole, cerebral cortical damage, etc., the convergence rate decreases (the distribution concentrates in the Iζ convergence rate region, and a specific pattern of convergence occurs depending on the abnormal location). The low rate becomes noticeable.

Therefore, according to this embodiment, by using the convergence rate as a parameter, even people who are not necessarily specialists can improve the brain activity. J) normality/
Abnormal rJ+ can be determined easily in a short time.

Historically, abnormal areas can be detected by comparing the convergence rate of each detection pattern, and it is also possible to quantitatively understand changes in υ, discernment health, changes during various activations, medication effects, etc., and make judgments easily. be able to.

In addition, in the above explanation, the case where six detection patterns are obtained from eight electrode placement sites was explained, but if necessary, combinations of patterns of 4:F: of the shortest number of electrodes may be obtained according to the desired arrangement. can be used and the desired 1'
i. It is also possible to set a detection area for the body.

Next, another embodiment of the present invention will be described with reference to FIG.

The configuration of the main parts of another embodiment of the present invention is shown in FIG. This FIG. 6 corresponds to the correlation calculation means (24) and convergence fX calculation means (27) of FIG. 1, and has three input terminals (30a), (30b) and (30c).
The preprocessed 3-channel input signal from j°
amplifiers (31a), (31b) and (31
c) and the corresponding inverting amplifier (41
a).

(41,b) and (41c), respectively.

These inverting amplifiers are provided corresponding to the polarity combinations shown in Table 1 above.

Amplifiers (31a), (31b) and (31c)
), each output of which has a square characteristic (32a),
(32b) and (32c) as well as a comparator (33a).

(331+') and (33c), respectively. The outputs of the comparators (33a) to (33c) are both supplied to an AND circuit (34). ;) When the signal poles 111 of the two input terminals (30a) to (30c) are all stopped, the output of the tent circuit (34) becomes "1" and the switch (35++) + (35b)
and (35c) are closed, and the stretchers (,'(2a)~(
Each output of 32 (:) is passed through an adder (36) to
It is supplied to the condenser (37) of the next generation H. Compressor (37
) is the root barry expressed by the above equation (1), which is integrated by the integrator (38) and becomes equivalent to the polarity surface 4 ratio M1 of the above polarity combination Ql. .

In addition, input terminal (30, +) ~ (30 (: ) 0 shoulder i
If the polarity of -X is negative, the inversion increase 111i
The outputs of the li devices (41a) to (4Lc) are supplied to the integral z): (48) in exactly the same manner as described above, and a polarity combination (corresponding to the near polarity area ratio M7) is obtained.

The outputs of the inner integrators (38) and (48) are added to the W adder (39)
The convergence rate of equation (3) above is obtained.

Although the diagram is omitted, by appropriately providing an amplifier and an inverting amplifier on the input side for other polarity combinations, components corresponding to the respective polarity area ratios can be obtained with the same configuration as above. . Therefore, in this example, −(also, +
' + ij statement) AI! Convergence; ↑・t can be found as in the example, and the same effect is achieved.

The embodiments in which the present invention was applied to electroencephalograms have been described above, but it is also possible to apply the invention to induced electroencephalograms, etc.
When the present invention is applied to other biological signal waveforms, such as electrocardiograms, the analysis results are quantitatively highlighted, making it easy to determine whether the IS is normal or abnormal. It is effective in detecting abnormal areas such as infarction.

(Effects of the Invention) As described below, according to the present invention, a plurality of biological signals are temporally determined by using the appearance rate of a combination in which each instantaneous value at the same time has a predetermined polarity. By quantifying the correlation, a 4F body signal analysis method can be obtained that can easily determine the state of the body.

[Brief explanation of drawings]

FIG. 1 is a functional block diagram showing the configuration of an embodiment of the biological signal 1+'l folding method according to the present invention, and FIGS. 2 and 3 are diagrams for explaining the operation of an embodiment of the present invention. Flowcharts and waveform diagrams, Figure 4 is a neo-block diagram of the configuration of an embodiment of the present invention, Figure 5 is a plan view of the route for agricultural use according to the present invention,
FIG. 6 is a block diagram showing the configuration of another embodiment of the present invention. (24) is a correlation calculation means, (25) is a data classification stage, (26) is an amplitude calculation means, and (27) is a convergence rate calculation means.

Claims (1)

[Claims] 1. From a plurality of biological signals of a predetermined length of time, a combination in which each instantaneous value at the same point in time of the plurality of biological signals has a predetermined polarity is separated, and at least one of the separated polarity combinations. A biological signal analysis method characterized in that the temporal correlation of the plurality of biological signals is expressed by the appearance rate of the biological signals. 2. The biological signal analysis method according to claim 1, wherein the polarity area ratio of at least one of the classified polarity combinations is obtained.