Fine-grained brain load determination method based on clustering algorithm and related device

Abstract

The invention provides a fine-grained brain load determination method based on a clustering algorithm and a related device, and belongs to the technical field of ergonomics. The method comprises the following steps: combining power spectrum densities of different frequency bands, dimensionality-reduced electroencephalogram signals and starting time of the electroencephalogram signals to obtain combined electroencephalogram characteristics; clustering the combined electroencephalogram features to obtain a plurality of clustering centers; performing cosine similarity scoring on the plurality of clustering centers to obtain a plurality of cosine similarity scoring results; and interpolating a plurality of cosine similarity scoring results to obtain a brain load value in the whole time period. According to the invention, the problem of low accuracy of determining the fine-grained brain load in the prior art is solved.

Description

Technical Field

[0001] The present invention belongs to the technical field of ergonomics, and in particular relates to a fine-grained brain load determination method based on a clustering algorithm and a related device. Background Art

[0002] Most of the existing fine-grained brain load determination methods are based on experimental paradigms, which are easy to classify brain load levels. First, the EEG data is labeled, and then deep learning or machine learning is used to complete the classification. The existing fine-grained brain load determination methods also directly calculate the brain load by refining EEG features (power spectral density value ratio).

[0003] The method based on the experimental paradigm requires the labeling of brain load, but the data in real scenarios is difficult to label. Moreover, the method based on the experimental paradigm is a coarse-grained brain load calculation, which cannot reflect the load change trend in detail. The method based on the experimental paradigm is not a real task scenario and is not practical. The method based on the experimental paradigm scores the load level and is not targeted at the experimental individual (the same task has different load levels for different individuals), making it difficult to implement applications such as load monitoring and reminders for individuals.

[0004] The brain load is directly calculated by refining EEG features (ratio of power spectral density values). This method only uses the power spectral density of EEG features, which does not fully utilize the information, resulting in insufficient load calculation, which leads to low accuracy in determining fine-grained brain load. Summary of the invention

[0005] The purpose of the present invention is to provide a fine-grained brain load determination method and related devices based on a clustering algorithm, so as to solve the problem of low accuracy in determining fine-grained brain load in the prior art.

[0006] In order to achieve the above object, the present invention adopts the following technical solution: In a first aspect, the present invention provides a fine-grained brain load determination method based on a clustering algorithm, comprising the following steps: Obtaining EEG signals, filtering, whole-brain averaging, independent component analysis, and artifact removal preprocessing of the obtained EEG signals; Perform wavelet packet decomposition on the preprocessed EEG signal to obtain continuous time-frequency information after wavelet packet decomposition; Calculate the power spectrum density of different frequency bands based on the continuous time-frequency information after wavelet packet decomposition; Performing dimension reduction on the obtained EEG signal to obtain a dimension-reduced EEG signal; The power spectrum density of different frequency bands, the EEG signal after dimensionality reduction and the start time of the EEG signal are combined to obtain the combined EEG feature; Cluster the merged EEG features to obtain several cluster centers; Scoring the cosine similarity of several cluster centers to obtain several cosine similarity scoring results; Interpolate several cosine similarity scoring results to obtain the brain load value for the entire time period.

[0007] A further improvement of the present invention is that, in the steps of acquiring EEG signals and filtering, whole-brain averaging, independent component analysis and artifact removal preprocessing the obtained EEG signals, the EEG signals are specifically acquired through an EEG device or a 32-channel EEG cap, and the obtained EEG signals are filtered, whole-brain averaged, independent component analysis and artifact removal preprocessing are performed.

[0008] A further improvement of the present invention is that, in the step of performing wavelet packet decomposition on the preprocessed EEG signal to obtain continuous time-frequency information after wavelet packet decomposition, the MATLAB wavelet packet decomposition toolkit is specifically used to perform wavelet packet decomposition on the preprocessed EEG signal to obtain continuous time-frequency information after wavelet packet decomposition.

[0009] A further improvement of the present invention is that, in the step of reducing the dimension of the obtained EEG signal to obtain the EEG signal after dimension reduction, a principal component analysis method is specifically used to reduce the dimension of the obtained EEG signal to obtain the EEG signal after dimension reduction.

[0010] A further improvement of the present invention is that, in the step of clustering the combined EEG features to obtain a plurality of cluster centers, a k-means clustering algorithm is specifically used to cluster the combined EEG features to obtain a plurality of cluster centers.

[0011] A further improvement of the present invention is that, in the step of interpolating a plurality of cosine similarity scoring results to obtain the brain load value for the entire time period, a weighted summation method is specifically used to interpolate a plurality of cosine similarity scoring results to obtain the brain load value for the entire time period.

[0012] In a second aspect, the present invention provides a fine-grained brain load determination system based on a clustering algorithm, comprising a data acquisition module, a data preprocessing module, a wavelet packet decomposition module, a power spectrum density acquisition module, a data dimension reduction module, a data merging module, a data clustering module, a cosine similarity scoring module and a data interpolation module; The data acquisition module is used to acquire EEG signals; The data preprocessing module is used to perform filtering, whole-brain averaging, independent component analysis and artifact removal preprocessing on the obtained EEG signal; The wavelet packet decomposition module is used to perform wavelet packet decomposition on the preprocessed EEG signal to obtain continuous time-frequency information after wavelet packet decomposition; The power spectrum density acquisition module is used to calculate the power spectrum density of different frequency bands based on the continuous time-frequency information after wavelet packet decomposition; The data dimension reduction module is used to reduce the dimension of the obtained EEG signal to obtain the EEG signal after dimension reduction; The data merging module is used to merge the power spectrum density of different frequency bands, the EEG signal after dimensionality reduction and the start time of the EEG signal to obtain a merged EEG feature; The data clustering module is used to cluster the combined EEG features to obtain a number of cluster centers; The cosine similarity scoring module is used to perform cosine similarity scoring on a plurality of cluster centers to obtain a plurality of cosine similarity scoring results; The data interpolation module is used to interpolate a number of cosine similarity scoring results to obtain the brain load value of the entire time period.

[0013] A further improvement of the present invention is that the data dimension reduction module uses principal component analysis to reduce the dimension of the obtained EEG signal to obtain a signal after dimension reduction.

[0014] In a third aspect, the present invention provides an electronic device including a memory and a processor, wherein the memory stores a computer program, and the processor implements the steps of the above-mentioned method for determining fine-grained brain load based on a clustering algorithm when executing the computer program.

[0015] In a fourth aspect, the present invention provides a storage medium having a computer program stored thereon, which, when executed by a processor, implements the steps of the above-mentioned method for determining fine-grained brain load based on a clustering algorithm.

[0016] Compared with the prior art, the present invention has the following beneficial effects: The present invention is an improved invention. Compared with the existing fine-grained brain load determination method based on clustering algorithm, on the one hand, the present invention combines the power spectrum density of different frequency bands, the EEG signal after dimensionality reduction and the start time of the EEG signal to obtain a combined EEG feature. It can be seen that the factors considered by the present invention are relatively comprehensive, which improves the accuracy of the fine-grained brain load determined later. On the other hand, the present invention clusters the combined EEG features to obtain several clustering centers. The clustering algorithm can fully and effectively utilize the multi-dimensional information of the EEG signal, making the subsequent analysis and processing (interpolation to obtain brain load) more efficient and accurate. In addition, interpolation of several cosine similarity scoring results to obtain the brain load value of the entire time period can ensure the integrity of the data. BRIEF DESCRIPTION OF THE DRAWINGS

[0017] Figure 1 It is a flow chart of the fine-grained brain load determination method based on clustering algorithm of the present invention; Figure 2 is a schematic diagram of a fine-grained brain load determination system based on a clustering algorithm of the present invention; Figure 3 The brain load curves of different human bodies of the present invention are Figure 1 ; Figure 4 The brain load curves of different human bodies of the present invention are Figure 2 ; Figure 5 It is a schematic diagram of the structure of the electronic device of the present invention. DETAILED DESCRIPTION

[0018] In order to further understand the content of the present invention, the present invention is described in detail below in conjunction with the accompanying drawings and specific embodiments. It should be understood that the embodiments are only for explaining the present invention and are not intended to limit it.

[0019] The fine-grained brain load determination method based on clustering algorithm proposed by the present invention combines the power spectrum density of different frequency bands, the EEG signal after dimensionality reduction and the start time of the EEG signal to obtain a combined EEG feature; clusters the combined EEG feature to obtain several cluster centers; performs cosine similarity scoring on several cluster centers to obtain several cosine similarity scoring results; interpolates several cosine similarity scoring results to obtain the brain load value of the entire time period. Compared with the prior art, the present invention effectively solves the problem of low accuracy in determining fine-grained brain load in the prior art.

[0020] Embodiment 1: The flowchart of the fine-grained brain load determination method based on clustering algorithm of the present invention is as follows: Figure 1 As shown, the fine-grained brain load determination method based on clustering algorithm of the present invention comprises the following steps: S1. Obtain EEG signals, and perform filtering, whole-brain averaging, independent component analysis, and artifact removal preprocessing on the obtained EEG signals.

[0021] S2. Perform wavelet packet decomposition on the preprocessed EEG signal to obtain continuous time-frequency information after wavelet packet decomposition.

[0022] S3. Calculate the power spectral density of different frequency bands based on the continuous time-frequency information after wavelet packet decomposition.

[0023] S4. Perform dimension reduction on the obtained EEG signal to obtain an EEG signal after dimension reduction.

[0024] S5. Combine the power spectral density of different frequency bands, the EEG signal after dimensionality reduction and the start time of the EEG signal to obtain a combined EEG feature.

[0025] S6. Cluster the merged EEG features to obtain several cluster centers.

[0026] S7. Perform cosine similarity scoring on several cluster centers to obtain several cosine similarity scoring results.

[0027] S8. Interpolate several cosine similarity scoring results to obtain the brain load value of the entire time period.

[0028] Embodiment 2: The schematic diagram of the fine-grained brain load determination system based on clustering algorithm of the present invention is as follows: Figure 2 As shown, the fine-grained brain load determination system based on clustering algorithm of the present invention includes a data acquisition module, a data preprocessing module, a wavelet packet decomposition module, a power spectrum density acquisition module, a data dimension reduction module, a data merging module, a data clustering module, a cosine similarity scoring module and a data interpolation module.

[0029] The data acquisition module is used to acquire EEG signals.

[0030] The data preprocessing module is used to filter, average the whole brain, perform independent component analysis and remove artifacts preprocessing on the obtained EEG signals.

[0031] The wavelet packet decomposition module is used to perform wavelet packet decomposition on the preprocessed EEG signal to obtain continuous time-frequency information after wavelet packet decomposition.

[0032] The power spectrum density acquisition module is used to calculate the power spectrum density of different frequency bands based on the continuous time-frequency information after wavelet packet decomposition.

[0033] The data dimension reduction module is used to reduce the dimension of the obtained EEG signal to obtain the EEG signal after dimension reduction.

[0034] The data merging module is used to merge the power spectral density of different frequency bands, the EEG signal after dimensionality reduction and the start time of the EEG signal to obtain the merged EEG features.

[0035] The data clustering module is used to cluster the merged EEG features to obtain several cluster centers.

[0036] The cosine similarity scoring module is used to perform cosine similarity scoring on a number of cluster centers to obtain a number of cosine similarity scoring results.

[0037] The data interpolation module is used to interpolate several cosine similarity scoring results to obtain the brain load value of the entire time period.

[0038] Embodiment 3: S1. Obtain EEG signals, and perform filtering, whole-brain averaging, independent component analysis, and artifact removal preprocessing on the obtained EEG signals.

[0039] First, the EEG signal is obtained through an EEG device or a 32-channel EEG cap, and then the obtained EEG signal is filtered, averaged across the whole brain, subjected to independent component analysis, and preprocessed to remove artifacts.

[0040] Specifically, the EEG signal is first filtered with a notch filter to eliminate 50Hz power line interference, then the EEG signal in the frequency range of 1 to 30Hz is bandpass filtered with a Butterworth filter, and then the EEG signal is averaged across the whole brain, followed by independent component analysis (also called independent component analysis) of the EEG signal to extract 20 independent components. Finally, artifacts such as electrocardiogram and electrooculogram interference are automatically identified and deleted.

[0041] S2. Perform wavelet packet decomposition on the preprocessed EEG signal to obtain continuous time-frequency information after wavelet packet decomposition.

[0042] In this step, the MATLAB wavelet packet decomposition toolkit (its function is sym6) is used to perform wavelet packet decomposition on the preprocessed EEG signal to obtain continuous time-frequency information after wavelet packet decomposition.

[0043] S3. Calculate the power spectral density of different frequency bands based on the continuous time-frequency information after wavelet packet decomposition.

[0044] In this step, the power spectral density of different frequency bands is calculated based on the continuous time-frequency information after wavelet packet decomposition (specifically, the power spectral density of different frequency bands is obtained by calling the wpspectrum function using MATLAB).

[0045] The calculation formula of wavelet packet decomposition is:

[0046] in n =0,1,2, indicating the wavelet packet determined by the orthogonal scaling function. and are two sets of coefficients, usually used to define filters or scaling functions and wavelet functions in wavelet packet decomposition. For the n The wavelet packet function of level , k Represents an index, and its domain is the integer domain. t is the independent variable of the wavelet packet function, representing the time domain. Specifically, this embodiment adopts 4-level wavelet packet decomposition, which can easily decompose different frequency bands.

[0047] S4. Perform dimension reduction on the obtained EEG signal to obtain an EEG signal after dimension reduction.

[0048] In this step, the principal component analysis method is specifically used to reduce the dimension of the obtained EEG signal to obtain the EEG signal after dimension reduction.

[0049] S5. Combine the power spectral density of different frequency bands, the EEG signal after dimensionality reduction and the start time of the EEG signal to obtain a combined EEG feature.

[0050] In this step, FEATURE is specifically used to represent the process of merging the power spectral density of different frequency bands, the EEG signal after dimensionality reduction, and the start time of the EEG signal to obtain the merged EEG feature.

[0051] The process of merging the power spectral density of different frequency bands, the EEG signal after dimensionality reduction and the start time of the EEG signal using FEATURE to obtain the merged EEG feature is specifically expressed as follows: FEATURE=[start time of EEG signal, EEG signal after dimensionality reduction, power spectrum density ratio of different frequency bands, power corresponding to different frequency bands].

[0052] The power spectral density ratios of different frequency bands include three cases: theta / alpha, theta / beta and beta / (alpha+theta).

[0053] The powers corresponding to different frequency bands are alpha, beta, theta, delta, and gamma.

[0054] S6. Cluster the merged EEG features to obtain several cluster centers.

[0055] In this step, the k-means clustering algorithm is used to cluster the merged EEG features to obtain several cluster centers.

[0056] S7. Perform cosine similarity scoring on several cluster centers to obtain several cosine similarity scoring results.

[0057] The calculation formula for cosine similarity is:

[0058] in, similarity is the cosine similarity, A is the cluster center feature vector, B are vectors of different scoring criteria.

[0059] S8. Interpolate several cosine similarity scoring results to obtain the brain load value of the entire time period.

[0060] In this step, a weighted sum method is used to interpolate several cosine similarity scoring results to obtain the brain load value of the entire time period. The calculation formula is as follows:

[0061] in, is the influence of different cluster centers on the load value in this time period ( Proportional to the weights of different cluster centers ), The calculation of is based on the exponential decrease of e, with the coefficient of decrease proportional to the distance of the point from the cluster center (the farther the distance, the smaller the impact), represents the weight, The influence of all cluster centers is normalized (the weights are converted to between 0 and 1). is the distance between the point (the feature vector of the time segment, one time segment every three seconds) and the cluster center (the Euclidean distance in the feature vector space), For all the cluster centers listed , is the MWL discreteness of all points in a certain class in the cluster, k The value of is [1-7], k For an index, k Indicates that the weight values ​​of all cluster centers are calculated.

[0062] In order to verify the effectiveness of the fine-grained brain load determination method based on clustering algorithm proposed in the present invention, this embodiment observes the trend of brain load curves of different human bodies to determine whether it can match the difficulty of the actual task (flight task). The results of the two brain load curves of different human bodies are as follows: Figure 3 and Figure 4 shown. Figure 3 and Figure 4 The horizontal axis represents time, and the vertical axis represents brain load. Figure 3 The curve data shows that the brain load of different human bodies decreases first and then increases, which is consistent with the difficulty of the actual mission (flight mission). Figure 4 and Figure 3 It can be seen that the changes in brain load of different groups of individuals (also called a collection of many individuals) are similar, which shows that the method of the present invention can reflect the actual situation of brain load, thereby verifying the effectiveness of the method of the present invention.

[0063] Embodiment 4: See also Figure 5 As shown, the present invention also provides an electronic device 100 for a fine-grained brain load determination method based on a clustering algorithm; the electronic device 100 includes a memory 101, at least one processor 102, a computer program 103 stored in the memory 101 and executable on the at least one processor 102, and at least one communication bus 104.

[0064] The memory 101 can be used to store the computer program 103. The processor 102 implements the steps of the fine-grained brain load determination method based on the clustering algorithm described in Example 1 by running or executing the computer program stored in the memory 101 and calling the data stored in the memory 101. The memory 101 can mainly include a program storage area and a data storage area, wherein the program storage area can store an operating system, an application required for at least one function (such as a sound playback function, an image playback function, etc.), etc.; the data storage area can store data (such as audio data) created according to the use of the electronic device 100, etc. In addition, the memory 101 can include a non-volatile memory, such as a hard disk, a memory, a plug-in hard disk, a smart memory card (Smart Media Card, SMC), a secure digital (Secure Digital, SD) card, a flash card (Flash Card), at least one disk storage device, a flash memory device, or other non-volatile solid-state storage devices.

[0065] The at least one processor 102 may be a central processing unit (CPU), or other general-purpose processors, digital signal processors (DSP), application-specific integrated circuits (ASIC), field-programmable gate arrays (FPGA) or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. The processor 102 may be a microprocessor or any conventional processor, etc. The processor 102 is the control center of the electronic device 100, and uses various interfaces and lines to connect various parts of the entire electronic device 100.

[0066] The memory 101 in the electronic device 100 stores a plurality of instructions to implement a fine-grained brain load determination method based on a clustering algorithm, and the processor 102 may execute the plurality of instructions to implement: Obtaining EEG signals, filtering, whole-brain averaging, independent component analysis, and artifact removal preprocessing of the obtained EEG signals; Perform wavelet packet decomposition on the preprocessed EEG signal to obtain continuous time-frequency information after wavelet packet decomposition; Calculate the power spectrum density of different frequency bands based on the continuous time-frequency information after wavelet packet decomposition; Performing dimension reduction on the obtained EEG signal to obtain a dimension-reduced EEG signal; The power spectrum density of different frequency bands, the EEG signal after dimensionality reduction and the start time of the EEG signal are combined to obtain the combined EEG feature; Cluster the merged EEG features to obtain several cluster centers; Scoring the cosine similarity of several cluster centers to obtain several cosine similarity scoring results; Interpolate several cosine similarity scoring results to obtain the brain load value for the entire time period.

[0067] Embodiment 5: If the module / unit integrated in the electronic device 100 is implemented in the form of a software functional unit and sold or used as an independent product, it can be stored in a computer-readable storage medium. Based on this understanding, the present invention implements all or part of the process in the above-mentioned embodiment method, and can also be completed by instructing the relevant hardware through a computer program. The computer program can be stored in a computer-readable storage medium, and the computer program can implement the steps of the above-mentioned method embodiments when executed by the processor. Among them, the computer program includes computer program code, and the computer program code can be in source code form, object code form, executable file or some intermediate form. The computer-readable medium may include: any entity or device that can carry the computer program code, recording medium, U disk, mobile hard disk, disk, optical disk, computer memory and read-only memory (ROM, Read-Only Memory).

[0068] It will be appreciated by those skilled in the art that embodiments of the present invention may be provided as methods, systems, or computer program products. Therefore, the present invention may take the form of a complete hardware embodiment, a complete software embodiment, or an embodiment combining software and hardware. Furthermore, the present invention may take the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.

[0069] The present invention is described with reference to flowcharts and / or block diagrams of methods, devices (systems), and computer program products according to embodiments of the present invention. It should be understood that each process and / or block in the flowchart and / or block diagram, as well as the combination of processes and / or blocks in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing device to produce a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the processes in the flowchart and / or block diagram. Figure 1 A process or multiple processes and / or boxes Figure 1A device that provides the functions specified in a block or multiple blocks.

[0070] These computer program instructions may also be stored in a computer-readable memory capable of directing a computer or other programmable data processing device to operate in a specific manner, so that the instructions stored in the computer-readable memory produce an article of manufacture comprising an instruction device, which implements the process Figure 1 A process or multiple processes and / or boxes Figure 1 A function specified in one or more boxes.

[0071] These computer program instructions can also be loaded onto a computer or other programmable data processing device so that a series of operating steps are executed on the computer or other programmable device to produce a computer-implemented process, thereby providing instructions for implementing the process. Figure 1 A process or multiple processes and / or boxes Figure 1 The steps for the functions specified in one or more boxes.

[0072] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention rather than to limit it. Although the present invention has been described in detail with reference to the above embodiments, ordinary technicians in the relevant field should understand that the specific implementation methods of the present invention can still be modified or replaced by equivalents. Any modification or equivalent replacement that does not depart from the spirit and scope of the present invention should be covered within the scope of protection of the claims of the present invention.

Claims

1. A fine-grained brain load determination method based on clustering algorithm, characterized in that: The following steps are involved: Obtaining EEG signals, filtering, whole-brain averaging, independent component analysis, and artifact removal preprocessing of the obtained EEG signals; Perform wavelet packet decomposition on the preprocessed EEG signal to obtain continuous time-frequency information after wavelet packet decomposition; Calculate the power spectrum density of different frequency bands based on the continuous time-frequency information after wavelet packet decomposition; Performing dimension reduction on the obtained EEG signal to obtain a dimension-reduced EEG signal; The power spectrum density of different frequency bands, the EEG signal after dimensionality reduction and the start time of the EEG signal are combined to obtain the combined EEG feature; Cluster the merged EEG features to obtain several cluster centers; Scoring the cosine similarity of several cluster centers to obtain several cosine similarity scoring results; Interpolate several cosine similarity scoring results to obtain the brain load value for the entire time period.

2. The fine-grained brain load determination method based on clustering algorithm according to claim 1 is characterized in that: In the steps of obtaining EEG signals, filtering, whole-brain averaging, independent component analysis and artifact removal preprocessing of the obtained EEG signals, the EEG signals are obtained specifically through an EEG device or a 32-channel EEG cap, and the obtained EEG signals are filtered, whole-brain averaged, independent component analysis and artifact removal preprocessing.

3. The fine-grained brain load determination method based on clustering algorithm according to claim 1 is characterized in that: In the step of performing wavelet packet decomposition on the preprocessed EEG signal to obtain continuous time-frequency information after wavelet packet decomposition, specifically, MATLAB wavelet packet decomposition toolkit is used to perform wavelet packet decomposition on the preprocessed EEG signal to obtain continuous time-frequency information after wavelet packet decomposition.

4. The fine-grained brain load determination method based on clustering algorithm according to claim 1 is characterized in that: In the step of reducing the dimension of the obtained EEG signal to obtain the EEG signal after dimension reduction, a principal component analysis method is specifically used to reduce the dimension of the obtained EEG signal to obtain the EEG signal after dimension reduction.

5. The fine-grained brain load determination method based on clustering algorithm according to claim 1 is characterized in that: In the step of clustering the combined EEG features to obtain a plurality of cluster centers, a k-means clustering algorithm is specifically used to cluster the combined EEG features to obtain a plurality of cluster centers.

6. The fine-grained brain load determination method based on clustering algorithm according to claim 1, characterized in that: In the step of interpolating a plurality of cosine similarity scoring results to obtain the brain load value for the entire time period, a weighted sum method is specifically used to interpolate a plurality of cosine similarity scoring results to obtain the brain load value for the entire time period.

7. A fine-grained brain load determination system based on clustering algorithm, characterized in that: It includes data acquisition module, data preprocessing module, wavelet packet decomposition module, power spectrum density acquisition module, data dimension reduction module, data merging module, data clustering module, cosine similarity scoring module and data interpolation module; The data acquisition module is used to acquire EEG signals; The data preprocessing module is used to perform filtering, whole-brain averaging, independent component analysis and artifact removal preprocessing on the obtained EEG signal; The wavelet packet decomposition module is used to perform wavelet packet decomposition on the preprocessed EEG signal to obtain continuous time-frequency information after wavelet packet decomposition; The power spectrum density acquisition module is used to calculate the power spectrum density of different frequency bands based on the continuous time-frequency information after wavelet packet decomposition; The data dimension reduction module is used to reduce the dimension of the obtained EEG signal to obtain the EEG signal after dimension reduction; The data merging module is used to merge the power spectrum density of different frequency bands, the EEG signal after dimensionality reduction and the start time of the EEG signal to obtain a merged EEG feature; The data clustering module is used to cluster the combined EEG features to obtain a number of cluster centers; The cosine similarity scoring module is used to perform cosine similarity scoring on a plurality of cluster centers to obtain a plurality of cosine similarity scoring results; The data interpolation module is used to interpolate a number of cosine similarity scoring results to obtain the brain load value of the entire time period.

8. The fine-grained brain load determination system based on clustering algorithm according to claim 7, characterized in that: The data dimension reduction module uses principal component analysis to reduce the dimension of the obtained EEG signal to obtain the EEG signal after dimension reduction.

9. An electronic device comprising a memory and a processor, wherein the memory stores a computer program, wherein: When the processor executes the computer program, the processor implements the steps of the fine-grained brain load determination method based on a clustering algorithm according to any one of claims 1 to 6.

10. A storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the steps of the fine-grained brain load determination method based on a clustering algorithm according to any one of claims 1 to 6 are implemented.

Images