CN102657526A - Method for evaluating R value and power spectrums of electroencephalogram signals causing discomfort when people watch three-dimensional (3D) images - Google Patents

Abstract

Provided is a method for evaluating an R value and power spectrums of electroencephalogram signals causing discomfort when people watch three-dimensional (3D) images. The method comprises arranging an experimental environment and experimental conditions, arranging experimental sequence of a subject, explaining the experiment process, arranging an electroencephalogram electrical electrode, collecting electroencephalogram signals of the subject before watching images for 2 minutes, and performing a first questionnaire survey; collecting electroencephalogram signals when the subject watches the whole movie and performing a second questionnaire survey after the movie is over; collecting electroencephalogram signals for 5 minutes when the subject rests for 10-20 minutes after watching the movie and performing a third questionnaire survey; preprocessing the collected electroencephalogram signals; analyzing the power spectrums in a plurality of wave bands according to the frequency range obtained by preprocessing the electroencephalogram signals; and analyzing the R value and performing statistics t-inspection under 2D experimental conditions and 3D experimental conditions. The method can perform comparison of 2D images and the 3D images, analyzes centralized brain regions of the electroencephalogram signals from the energy angle, provides an experimental basis for follow-up researches, and then provides thoughts for technical standards for evaluating the 3D images.

Description

Watch the 3D image to cause the EEG signals power spectrum and the evaluation methodology of R value of sense of discomfort

Technical field

The present invention relates to a kind of collection and feature extracting method of EEG signals.Particularly relate to EEG signals power spectrum and the evaluation methodology of R value that a kind of 3D of watching image causes sense of discomfort

Background technology

The 3D TV is the abbreviation of 3 D stereoscopic image TV.It is slightly variant that the 3 D stereoscopic image TV utilizes people's eyes to observe the angle of object; Therefore can distinguish that object is far and near; Produce this principle of spatial vision, the image of seeing right and left eyes separates, thereby makes the user need not can bore hole experience stereoperception by anaglyph spectacles.The improvement to picture quality is not only in the generation of 3D TV, especially to change vivid in the image.Though stereo display technique has third dimension and the very strong advantage of feeling of immersion; If but watch stereoscopic image (more than 30 minutes) for a long time; Just produce visual fatigue, uncomfortable symptom such as dizzy probably, this has influenced popularizing and development of stereo display technique to a certain extent.Health and safety also is the difficult problem that 3D stereo display urgent need is captured.Result of study shows, spectators are when watching stereoscopic image, because eyes can promptly move around, thereby causes sense of discomfort such as eyestrain easily.

Human background brain wave such as Li and event related potential are estimated visual fatigue, and the result shows that the frequency spectrum of background brain wave and the preclinical delay of P700 all depend on visual disparity and viewing time, and this point also obtains checking in the measuring method of subjectivity.Therefore, the delay of electrical wave measurement's visual information in transmission course of requiring mental skill is a kind of effective ways of estimating visual fatigue.To the senses of discomfort such as visual fatigue of watching 3-dimensional image to cause, corresponding research work is carried out also fewerly.

For a new technology product, freely watch and need be guaranteed, therefore, picture quality and visual comfort all have advantage than conventional television.If can find a kind of three-dimensional stereo display technique; When can reduce existing stereoscopic display device cost greatly; Can be simple and easy to again with, take up space limited; Of paramount importance is can alleviate or solve present stereo display technique greatly to give numerous physiology and the psychological problem that people brought, and the development that will inevitably promote stereo display technique lets 3 D stereo show that science service is in popular life.

Summary of the invention

Technical problem to be solved by this invention is; Cause the similarities and differences of human body electroencephalogram's signal when providing a kind of basis to watch 2D, 3D image; Thereby a kind of method of assessing visual fatigue is provided, promotes the 3D image of watching of 3D image technology development to cause the EEG signals power spectrum and the evaluation methodology of R value of sense of discomfort.

The technical scheme that the present invention adopted is: a kind of 3D of watching image causes the EEG signals power spectrum and the evaluation methodology of R value of sense of discomfort, comprises the steps:

(1) experimental situation and condition are set, arrange experimenter's experimental sequence;

(2) explain orally experiment flow, settle electrode for encephalograms, gather the preceding EEG signals 2min of experimenter's viewing, and carry out the questionnaire survey first time, described questionnaire comprises visual discomfort sense, the uncomfortable problem of body part;

(3) watch whole film and gather EEG signals, carry out the questionnaire survey second time behind the wrapped film, described questionnaire comprises visual discomfort sense, the uncomfortable problem of body part;

(4) experimenter behind the 10～20min that has a rest, gathers EEG signals 5min after finishing the viewing stage, carries out questionnaire survey for the third time, and described questionnaire comprises visual discomfort sense, the uncomfortable problem of body part;

(5) EEG signals of being gathered are carried out pretreatment, comprising:

(a) EEG signals of gathering being reseted reference, filtering de-noising and normalization handles;

(b) utilize the ICA method to remove the eye electricity, soon signal decomposition becomes each independent components and component is analyzed;

(6), carry out the power spectrumanalysis of a plurality of wave bands to the frequency range of pretreatment gained EEG signals:

To the experimental data that collects, on MATLAB 7.11 platforms, handle; At first; Eeg data is averaged power spectrum find the solution, choose corresponding data point and carry out Fourier transformation, get frequency spectrum and the conjugated product of this frequency spectrum then; Obtain the estimation of power spectrum; Calculate the percent value that each wave band accounts for general power spectrum, and to whole experiment duration totally 85～100min carry out fitting a straight line, seek overall variation trend; Secondly, on the SPSS18.0 platform, respectively power spectrum numerical value and percent value are carried out one factor analysis of variance, in order to judge the difference under 2D and two kinds of experiment conditions of 3D;

(7) the R value is analyzed, and carries out the statistics t-check of two kinds of experiment condition 2D and 3D;

Step 1) comprises: the bias light illumination that is provided with in the room is 2.5～3.5Lux; Select the experimenter to watch the film that to show for the first time; The 3D television set is positioned at the dead ahead of experimenter 2.5～3.5m; When watching the 3D image, wear the fast gate-type 3D glasses of active of 120/second of field frequencies, whole experimental period is 85～100min.

Step 1) also comprises: the experimenter is divided into two groups, watches the 3D image mode earlier for one group, after watch the 2D image mode, another group is watched the 2D pattern earlier, after watch the 3D pattern.

Step 2) the preceding EEG signals 2min of described collection viewing is to gather the experimenter to test the EEG signals under the preceding quiescent condition, the 1min that closes one's eyes, and the 1min that opens eyes, the sampling of EEG signals guarantees that electrode impedance is less than 5000 ohm.

The described interval of watching each experimenter of film requirement to watch 2D and two images of 3D to test of step 3) will be more than 10 days.

The described ICA method of utilizing of step 5) is removed the eye electricity, specifically:

Under the situation of unknown source signal S and this unknown source signal S aliasing characteristic A, only from the mixed signal X=AS that receives, extract source signal, that is: one group of source signal s (t)=[s independently ₁(t) ..., s _n(t)] ^TMix through hybrid system A, obtain observation signal x (t)=[x ₁(t) ..., x _n(t)] ^T, represent as follows:

$\left(\begin{array}{c}{x}_1\left(t\right)\\ x_2\left(t\right)\\ \·\\ \·\\ \·\\ x_n\left(t\right)\end{array}\right)=A\left(\begin{array}{c}{s}_1\left(t\right)\\ s_2\left(t\right)\\ \·\\ \·\\ \·\\ s_n\left(t\right)\end{array}\right)$

Wherein, source signal s (t) and hybrid system A are unknown, have only mixed x (t) to observe, and the task of ICA is under the prerequisite of A and S the unknown, finds and separates mixed matrix W, makes output matrix

U＝W·X＝W·A·S

Wherein, x (t)=[x ₁(t) ..., x _n(t)] ^TBe observation signal, s (t)=[s ₁(t) ..., s _n(t)] ^TBe source signal, A is a hybrid system, and S is a source signal, and U is an output matrix, and W is for separating mixed matrix, and X is a mixed signal

Step 7) uses alpha to cause an index of sense of discomfort as the average power spectra ratio of two frequency bands of 18-22Hz as the 3D image as 10-13Hz and beta, establishes signal and in the mean power spectrum density of frequency band h is:

$G\left(h\right)=\frac{{\∫}_{f_d}^{f_u}p\left(f\right)\mathrm{df}}{f_u-f_d}$

$R=\frac{G\left(h_1\right)}{G\left(h_2\right)}$

Wherein, G is the mean power spectrum density, and h is the frequency band of signal, f _uBe the upper limit of frequency band h, f _dBe the lower limit of frequency band h, p (f) is the power spectral density of signal, h ₁And h ₂Be the different frequency bands of EEG signals, R is the ratio of result of calculation definition;

All experimenters are carried out the calculating of R value, and the R value of 0～10min and 35～45min is carried out the t-check.

The 3D of watching image of the present invention causes the EEG signals power spectrum and the evaluation methodology of R value of sense of discomfort; Can carry out the contrast of 2D and 3D image; Analyze the concentrated brain district of EEG signals from energy point of view; For follow-up research provides experiment basis, thereby, promote the development of 3D image technology for the technical standard of estimating the 3D image provides thinking.

Description of drawings

Fig. 1 is the experiment figure that leads that adopts among the present invention;

Fig. 2 is the flow chart of the inventive method;

Fig. 3 is the lead average power spectra ratio figure of delta wave band of FP2;

Fig. 4 is under 2D and two kinds of experiment conditions of 3D, and Oz leads and locates the average power spectra of alpha wave band;

Fig. 5 (a) is the power spectrum ratio figure of the gamma wave band of 2D experimental group;

Fig. 5 (b) is the power spectrum ratio figure of the gamma wave band of 3D experimental group;

Fig. 6 is that the R value behind 0-10min and the 35-45min is t-assay figure under the 3D condition.

The specific embodiment

Below in conjunction with embodiment and accompanying drawing detailed description is made in EEG signals power spectrum and the evaluation methodology of R value that the 3D of watching image of the present invention causes sense of discomfort.

The 3D of watching image of the present invention causes the EEG signals power spectrum and the evaluation methodology of R value of sense of discomfort, and is as shown in Figure 2, comprises the steps:

(1) experimental situation and condition are set, arrange experimenter's experimental sequence;

Comprise: the bias light illumination that is provided with in the room is 2.5～3.5Lux; Select the experimenter to watch the image that to show for the first time; The 3D television set is positioned at the dead ahead of experimenter 2.5～3.5m; When watching the 3D image, wear the fast gate-type 3D glasses of active of 120/second of field frequencies, whole experimental period is 85～100min.

Also comprise: the experimenter is divided into two groups, watches the 3D pattern earlier for one group, after watch the 2D pattern, another group is watched the 2D pattern earlier, after watch the 3D pattern.

The 64 conducts digital eeg recording appearance that present embodiment adopts U.S. Neuroscan company to produce are gathered eeg data.The experimenter is arranged in the room that electromagnetic shielding is good, soundproof effect is good, temperature humidity is suitable and experimentizes.Block external light source through lighttight curtain, the bias light illumination in the room is 3Lux.The experimenter to feel comfortable but the posture that does not influence data acquisition sit in an armchair.Apart from the dead ahead about experimental subject 3m is 55 inches 3D TVs of a Samsung 55C8000XF, and the image of watching is the 3D animation: " Monster Aliens " and " Megamind ".When watching the 3D image, the experimenter wears the fast gate-type 3D glasses of active of 120/second of field frequencies.Whole experimental period is about 95min.The experimenter is 10 people in this experiment, and 5 people watch the 3D pattern earlier, after watch the 2D pattern, 5 people watch the 2D pattern earlier, after watch the 3D pattern.

(2) explain orally experiment flow, settle electrode for encephalograms, gather the preceding EEG signals 2min of experimenter's viewing, and carry out the questionnaire survey first time, described questionnaire comprises the problem of aspects such as visual discomfort sense, body part discomfort;

EEG signals 2min before the described collection viewing is to gather the experimenter to test the EEG signals under the preceding quiescent condition, the 1min that closes one's eyes, and the 1min that opens eyes, the sampling of EEG signals guarantees that electrode impedance is less than 5000 ohm.

This period is the preceding preparatory stage of experiment, and this period is introduced experiment flow by the experimenter to the experimenter.

The placement location of electrode adopts international 10/20 system standard, and is as shown in Figure 1.

Electrodes use Ag/AgCl electrode, and hang down (M2) as bipolar reference with left ear-lobe (M1) and auris dextra, the sample frequency of EEG signals is 1000Hz, and adopts the 50Hz wave trap to remove power frequency and disturb.In the experiment, guarantee that electrode impedance is less than 5000 ohm.

(3) watch whole image and gather EEG signals, image carries out questionnaire survey second time after finishing, and described questionnaire comprises the problem of visual discomfort sense, body part discomfort;

The described interval of watching each experimenter of image requirement to watch 2D and two images of 3D to test will be more than 10 days.

This period is the viewing period, and this window duration of present embodiment 80min requires experimenter's peace and quiet to be sitting on the armchair, watches image.According to the experiment progress, watch the image of 2D or 3D pattern respectively, when image finishes, carry out the questionnaire survey second time.At whole viewing experimental session, the experimenter watches image for the first time, does not have priori, requires the experimenter to keep relaxation state in the experimentation, does not allow any actual act.Whole experimental program requires each experimenter to accomplish 2D and two kinds of experiments of 3D, between two kinds of experiments, guarantees that each experimenter is at interval more than 10 days, so that carry out fatigue recovery.

(4) experimenter behind the 10～20min that has a rest, gathers EEG signals 5min after finishing the viewing stage, carries out questionnaire survey for the third time, and described questionnaire comprises visual discomfort sense, the uncomfortable problem of body part;

This period is convalescent period, and in the present embodiment, the experimenter is after finishing the viewing stage, and rest 15min carries out questionnaire survey once more, and gathers the eeg data after having a rest, duration 5min.

(5) EEG signals of being gathered are carried out pretreatment, comprising:

(a) EEG signals of gathering being reseted reference, filtering de-noising and normalization handles;

After data acquisition finishes, be the accuracy that guarantees that later data is analyzed, adopted and reseted reference, filtering de-noising and normalization processing.

(b) utilize the ICA method to remove the eye electricity, soon signal decomposition becomes each independent components and component is analyzed;

Independent component analysis (ICA) decomposes signal, is broken down into each independent components and component is analyzed.Under the situation of unknown source signal S and this unknown source signal S aliasing characteristic A, only from the mixed signal X=AS that receives, extract source signal.That is: one group of source signal s (t)=[s independently ₁(t) ..., s _n(t)] ^TMix through linear system A, obtain observation signal x (t)=[x ₁(t) ..., x _n(t)] ^TRepresent as follows:

$\left(\begin{array}{c}{x}_1\left(t\right)\\ x_2\left(t\right)\\ \·\\ \·\\ \·\\ x_n\left(t\right)\end{array}\right)=A\left(\begin{array}{c}{s}_1\left(t\right)\\ s_2\left(t\right)\\ \·\\ \·\\ \·\\ s_n\left(t\right)\end{array}\right)---\left(1\right)$

Wherein, source signal s (t) and hybrid system A are unknown, have only mixed x (t) to observe.The task of ICA is under the prerequisite of A and S the unknown, finds and separates mixed matrix W, makes output matrix

U＝W·X＝W·A·S (2)

In the formula, x (t)=[x ₁(t) ..., x _n(t)] ^TBe observation signal, s (t)=[s ₁(t) ..., s _n(t)] ^TBe source signal, A is a hybrid system, and S is a source signal, and U is an output matrix, and W is for separating mixed matrix, and X is a mixed signal

The EEG signals of being gathered in the experiment because of repeatedly blinking, obviously influence the calculating of average power spectra value.We use the ICA method to separate pretreated EEG signals; Find component and the noise component(s) the strongest with the signal of blinking dependency; Remove this two components; The data of separating after mixing matrix W and removing component to having preserved are inverted, thereby reconstruct EEG signals of removing the eye electricity, for follow-up analysis to energy lays the first stone.

(6), carry out the power spectrumanalysis of a plurality of wave bands to the frequency range of pretreatment gained EEG signals:

To the experimental data that collects, on MATLAB 7.11 platforms, handle; At first; Eeg data is averaged power spectrum find the solution, choose corresponding data point and carry out Fourier transformation, get frequency spectrum and the conjugated product of this frequency spectrum then; Obtain the estimation of power spectrum; Calculate the percent value that each wave band accounts for general power spectrum, and to whole experiment duration totally 85～100min carry out fitting a straight line, seek overall variation trend; Secondly, on the SPSS18.0 platform, respectively power spectrum numerical value and percent value are carried out one factor analysis of variance, in order to judge the difference under 3D and two kinds of experiment conditions of 2D;

Because duration of experiment is longer; The method that adopts in the present embodiment is that signal is carried out segment processing; Every separated 20min gets the 10s data, and alpha (8-13Hz), beta (13-30Hz), delta (1-4Hz), four wave bands of theta (4-8Hz) are carried out power spectrumanalysis, and the result sees from power spectrumanalysis; Subjects's sense of discomfort is apparent in view behind about 50-60min, shown in concrete outcome such as Fig. 3 to Fig. 5.

Shown in the viewing process result that record data analyze that leads among Fig. 3 to FP2; Can see at the average power spectra ratio that this leads from each wave band; FP2 is led after the data characteristics amount carries out fitting a straight line, and the power spectrum ratio of delta wave band shows that 3D image experimental group presents ascendant trend; 2D image experimental group presents slow ascendant trend.Because the delta wave band shows the adult and extremely tired and lethargy occurs, this presentation of results is compared to watches the 2D image, and the brain ammeter reveals more tangible sense of discomfort in the frontal region to watch the experimenter of 3D image.

For power spectrum situation over time, compare with the average profile diagram.Shown in Figure 4 for to the EEG signals that write down in the experimentation segmentation in per 10 minutes calculate average power spectral density value, carry out the variance analysis of repeated measure, the coordinate points of abscissa time represents 0,10,20,30,40,50,60,70 respectively, 80min.In the calculating that averages power spectrum, the alpha wave band that Oz leads and locates, after viewing began, increase trend all appearred in the situation of the image of 2D and two kinds of patterns of 3D, and the growth rate of 3D image is apparently higher than the 2D image.

Gamma wave band to two experimental grouies compares, and can see from Fig. 5, and 2 FPz that lead of frontal lobe, FP2 wave form varies trend are consistent, and 3 POz that lead of occipital lobe, Oz, O1 wave form varies trend are consistent; 3D image experimental group changes more obvious than 2D image experimental group, behind viewing 50min, reach peak value.

The people is under the situation that feels the stress; High band EEG signals (like the gamma ripple) are easy to observe; Among Fig. 5 (a) the power spectrum ratio figure of the gamma wave band of 2D image experimental group; Fig. 5 (b) is the power spectrum ratio figure of the gamma wave band of 3D image experimental group, gamma wave band analysis result is illustrated in the long period viewing after, by try psychentonia to occur, feel the stress the increase situation.

(7) the R value is analyzed, and carries out the statistics t-check of two kinds of experiment condition 2D and 3D;

The average power spectra ratio that uses alpha (10-13Hz) and two frequency bands of beta (18-22Hz) causes an index of sense of discomfort as the 3D image, establishes signal and in the mean power spectrum density of frequency band h is:

$G\left(h\right)=\frac{{\∫}_{f_d}^{f_u}p\left(f\right)\mathrm{df}}{f_u-f_d}---\left(3\right)$

$R=\frac{G\left(h_1\right)}{G\left(h_2\right)}---\left(4\right)$

Wherein, G is the mean power spectrum density, and h is the frequency band of signal, f _uBe the upper limit of frequency band h, f _dBe the lower limit of frequency band h, p (f) is the power spectral density of signal, h ₁And h ₂Be the different frequency bands of EEG signals, R is the ratio of result of calculation definition.

10 experimenters are carried out the calculating of R value, and the R value of 0～10min and 35～45min is carried out t-check, the result is as shown in Figure 6, and wherein, in order to simplify analysis, present embodiment has only adopted 21 in the standard lead to lead.Behind viewing 30min, the notable difference of R value has appearred in frontal region (1-3 number lead).

Final result shows, under the identical situation of other experiment conditions, with watching 3D image time lengthening, causes that easily the subjects produces fatigue state.Watch the difference of 2D and two kinds of EEG signals that image causes of 3D, comparatively obvious with the reaction in occipital lobe and two brain districts of frontal lobe.EEG signals power spectrum and the evaluation methodology of R value that the 3D image causes sense of discomfort can be used as a research direction and deeply inquire into, for the evaluation index that the 3D image further is discussed has been established strong basis.

Claims (7)

1. EEG signals power spectrum and the evaluation methodology of R value of watching the 3D image to cause sense of discomfort is characterized in that, comprises the steps:

(1) experimental situation and condition are set, arrange experimenter's experimental sequence;

(2) explain orally experiment flow, settle electrode for encephalograms, gather the preceding EEG signals 2min of experimenter's viewing, and carry out the questionnaire survey first time, described questionnaire comprises visual discomfort sense, the uncomfortable problem of body part;

(3) watch whole film and gather EEG signals, carry out the questionnaire survey second time behind the wrapped film, described questionnaire comprises visual discomfort sense, the uncomfortable problem of body part;

(4) experimenter behind the 10～20min that has a rest, gathers EEG signals 5min after finishing the viewing stage, carries out questionnaire survey for the third time, and described questionnaire comprises visual discomfort sense, the uncomfortable problem of body part;

(5) EEG signals of being gathered are carried out pretreatment, comprising:

(a) EEG signals of gathering being reseted reference, filtering de-noising and normalization handles;

(b) utilize the ICA method to remove the eye electricity, soon signal decomposition becomes each independent components and component is analyzed;

(6), carry out the power spectrumanalysis of a plurality of wave bands to the frequency range of pretreatment gained EEG signals:

To the experimental data that collects, on MATLAB 7.11 platforms, handle; At first; Eeg data is averaged power spectrum find the solution, choose corresponding data point and carry out Fourier transformation, get frequency spectrum and the conjugated product of this frequency spectrum then; Obtain the estimation of power spectrum; Calculate the percent value that each wave band accounts for general power spectrum, and to whole experiment duration totally 85～100min carry out fitting a straight line, seek overall variation trend; Secondly, on the SPSS18.0 platform, respectively power spectrum numerical value and percent value are carried out one factor analysis of variance, in order to judge the difference under 2D and two kinds of experiment conditions of 3D;

(7) the R value is analyzed, and carries out the statistics t-check of two kinds of experiment condition 2D and 3D.

2. the 3D of watching image according to claim 1 causes the EEG signals power spectrum and the evaluation methodology of R value of sense of discomfort; It is characterized in that step 1) comprises: the bias light illumination that is provided with in the room is 2.5～3.5Lux, selects the experimenter to watch the film that will show for the first time; The 3D television set is positioned at the dead ahead of experimenter 2.5～3.5m; When watching the 3D image, wear the fast gate-type 3D glasses of active of 120/second of field frequencies, whole experimental period is 85～100min.

3. the 3D of watching image according to claim 1 causes the EEG signals power spectrum and the evaluation methodology of R value of sense of discomfort; It is characterized in that; Step 1) also comprises: the experimenter is divided into two groups, watches the 3D image mode earlier for one group, after watch the 2D image mode; Another group is watched the 2D pattern earlier, after watch the 3D pattern.

4. the 3D of watching image according to claim 1 causes the EEG signals power spectrum and the evaluation methodology of R value of sense of discomfort; It is characterized in that; Step 2) the preceding EEG signals 2min of described collection viewing is to gather the experimenter to test the EEG signals under the preceding quiescent condition, and 1min closes one's eyes; The 1min that opens eyes, the sampling of EEG signals guarantees that electrode impedance is less than 5000 ohm.

5. the 3D of watching image according to claim 1 causes the EEG signals power spectrum and the evaluation methodology of R value of sense of discomfort, it is characterized in that, the described interval of watching each experimenter of film requirement to watch 2D and two images of 3D to test of step 3) will be more than 10 days.

6. the 3D of watching image according to claim 1 causes the EEG signals power spectrum and the evaluation methodology of R value of sense of discomfort, it is characterized in that, the described ICA method of utilizing of step 5) is removed the eye electricity, specifically:

Under the situation of unknown source signal S and this unknown source signal S aliasing characteristic A, only from the mixed signal X=AS that receives, extract source signal, that is: one group of source signal s (t)=[s independently ₁(t) ..., s _n(t)] ^TMix through hybrid system A, obtain observation signal x (t)=[x ₁(t) ..., x _n(t)] ^T, represent as follows:

Wherein, source signal s (t) and hybrid system A are unknown, have only mixed x (t) to observe, and the task of ICA is under the prerequisite of A and S the unknown, finds and separates mixed matrix W, makes output matrix

U＝W·X＝W·A·S

Wherein .x (t)=[x ₁(t) ..., x _n(t)] ^TBe observation signal, s (t)=[s ₁(t) ..., s _n(t)] ^TBe source signal, A is a hybrid system, and S is a source signal, and U is an output matrix, and W is for separating mixed matrix, and X is a mixed signal.

7. the 3D of watching image according to claim 1 causes the EEG signals power spectrum and the evaluation methodology of R value of sense of discomfort; It is characterized in that; Step 7) uses alpha to cause an index of sense of discomfort as the average power spectra ratio of two frequency bands of 18-22Hz as the 3D image as 10-13Hz and beta, establishes signal and in the mean power spectrum density of frequency band h is:

Wherein, G is the mean power spectrum density, and h is the frequency band of signal, f _uBe the upper limit of frequency band h, f _dBe the lower limit of frequency band h, p (f) is the power spectral density of signal, h ₁And h ₂Be the different frequency bands of EEG signals, R is the ratio of result of calculation definition;

All experimenters are carried out the calculating of R value, and the R value of 0～10min and 35～45min is carried out the t-check.

Images