Nothing Special   »   [go: up one dir, main page]

TWI467184B - Spectrum analysis method and disease examination method - Google Patents

Spectrum analysis method and disease examination method Download PDF

Info

Publication number
TWI467184B
TWI467184B TW100147346A TW100147346A TWI467184B TW I467184 B TWI467184 B TW I467184B TW 100147346 A TW100147346 A TW 100147346A TW 100147346 A TW100147346 A TW 100147346A TW I467184 B TWI467184 B TW I467184B
Authority
TW
Taiwan
Prior art keywords
spectrum
energy
bandwidth
signals
frequency
Prior art date
Application number
TW100147346A
Other languages
Chinese (zh)
Other versions
TW201326834A (en
Inventor
Ling Sheng Jang
Ming Kun Chen
Fu Yu Chang
Hong Yi Chang
Original Assignee
Univ Nat Cheng Kung
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 Univ Nat Cheng Kung filed Critical Univ Nat Cheng Kung
Priority to TW100147346A priority Critical patent/TWI467184B/en
Publication of TW201326834A publication Critical patent/TW201326834A/en
Application granted granted Critical
Publication of TWI467184B publication Critical patent/TWI467184B/en

Links

Landscapes

  • Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)

Description

頻譜分析方法及疾病篩檢方法Spectrum analysis method and disease screening method

本發明係關於一種頻譜分析方法及疾病篩檢方法。The present invention relates to a spectrum analysis method and a disease screening method.

能量醫學(Energy Medicine,EM)是利用科學儀器,以能量來對人體進行非侵入性的疾病檢測,進而可提供病人一種藥品,甚至非藥品的治療方式。Energy Medicine (EM) is a scientific instrument that uses energy to detect non-invasive diseases in humans, thereby providing patients with a drug or even non-pharmaceutical treatment.

生物體本身都會發出一種特有的微能量訊號,而該微能量訊號內可包含該生物體的內在資訊,例如人體是否具有某種疾病的訊息。不過,於直觀上,此微能量訊號是一種不易分析的雜亂訊號,而且目前並沒有任何一種能量訊號的頻譜分析方法,可將該微能量的頻譜進行分析,進而可將此分析應用在病理資訊上的快速篩檢。另外,也沒有一種疾病的篩檢方法,可根據人體本身的微能量訊號來進行各種疾病之可能性與潛在性之探討與分析。The organism itself emits a unique micro-energy signal, and the micro-energy signal can contain the internal information of the organism, such as whether the human body has a certain disease message. However, intuitively, this micro-energy signal is a kind of messy signal that is difficult to analyze, and there is no spectrum analysis method of energy signal at present, which can analyze the spectrum of the micro-energy, and then apply the analysis to pathological information. Quick screening on the top. In addition, there is no screening method for diseases, which can be based on the microenergy signal of the human body to explore and analyze the possibility and potential of various diseases.

因此,如何提供一種頻譜分析方法及一種疾病篩檢方法,可針對(微)能量訊號的頻譜進行分析,並可應用此頻譜分析方法於疾病的快速篩檢工作,已成為重要課題之一。Therefore, how to provide a spectrum analysis method and a disease screening method can analyze the spectrum of the (micro) energy signal, and can apply this spectrum analysis method to the rapid screening of diseases, which has become one of the important topics.

有鑑於上述課題,本發明之目的為提供一種可針對(微)能量訊號的頻譜進行分析,並可應用此頻譜分析方法於疾病的快速篩檢工作之頻譜分析方法及疾病篩檢方法。In view of the above problems, an object of the present invention is to provide a spectrum analysis method and a disease screening method capable of analyzing a spectrum of a (micro) energy signal and applying the spectrum analysis method to a rapid screening operation of a disease.

為達上述目的,依據本發明之一種頻譜分析方法係用以分析一能量訊號,能量訊號具有複數第一頻譜訊號及複數第二頻譜訊號,頻譜分析方法包括:一第一步驟:分別找出各第一頻譜訊號及各第二頻譜訊號之較大能量的頻譜,並依據該等較大能量的頻譜分別對應建立其第一頻寬矩陣及其第二頻寬矩陣,其中各第一頻寬矩陣及各第二頻寬矩陣分別具有各第一頻譜訊號及各第二頻譜訊號之複數主能量頻率;一第二步驟:依據各第一頻寬矩陣及各第二頻寬矩陣,分別計算各第一頻譜訊號及各第二頻譜訊號之該等主能量頻率的能量;一第三步驟:分別整合該等第一頻寬矩陣及該等第二頻寬矩陣,以分別建立該等第一頻譜訊號之一第一整合頻寬矩陣及該等第二頻譜訊號之一第二整合頻寬矩陣;一第四步驟:依據第一整合頻寬矩陣及第二整合頻寬矩陣,以分別將各第一頻譜訊號及各第二頻譜訊號之相同的主能量頻率的能量累加,進而分別建立一第一能量頻譜及一第二能量頻譜;以及一第五步驟:依據第一能量頻譜及第二能量頻譜,建立能量訊號的主能量頻率之一主能量頻譜。To achieve the above objective, a spectrum analysis method according to the present invention is for analyzing an energy signal, the energy signal having a plurality of first spectrum signals and a plurality of second spectrum signals, and the spectrum analysis method includes: a first step: separately finding each a spectrum of a larger energy of the first spectrum signal and each of the second spectrum signals, and correspondingly establishing a first bandwidth matrix and a second bandwidth matrix thereof according to the spectra of the larger energy, wherein each of the first bandwidth matrix And each of the second bandwidth matrices respectively has a plurality of primary energy frequencies of the first spectral signal and each of the second spectral signals; and a second step: calculating each of the first bandwidth matrix and each of the second bandwidth matrix The energy of the main energy frequencies of a spectrum signal and each of the second spectrum signals; a third step: integrating the first bandwidth matrix and the second bandwidth matrix to respectively establish the first spectrum signals a first integrated bandwidth matrix and a second integrated bandwidth matrix of the second spectral signals; a fourth step: according to the first integrated bandwidth matrix and the second integrated bandwidth matrix, respectively The energy of the same main energy frequency of each of the first spectrum signal and each of the second spectrum signals is accumulated, thereby respectively establishing a first energy spectrum and a second energy spectrum; and a fifth step: according to the first energy spectrum and the second The energy spectrum establishes the main energy spectrum of one of the main energy frequencies of the energy signal.

在一實施例中,於第一步驟中,各第一頻寬矩陣及各第二頻寬矩陣分別具有各第一頻譜訊號及各第二頻譜訊號之該等主能量頻率及其頻寬。In an embodiment, in the first step, each of the first bandwidth matrix and each of the second bandwidth matrices respectively have the main energy frequencies and the bandwidths of the first spectrum signals and the second spectrum signals.

在一實施例中,於第三步驟中,第一整合頻寬矩陣及第二整合頻寬矩陣係分別整合該等第一頻譜訊號及該等第二頻譜訊號之各主要能量頻率,並記錄該等第一頻譜訊號及該等第二頻譜訊號之各主要能量頻率的最大左頻寬及最大右頻寬。In an embodiment, in a third step, the first integrated bandwidth matrix and the second integrated bandwidth matrix respectively integrate the primary energy signals of the first spectral signals and the second spectral signals, and record the The maximum left bandwidth and the maximum right bandwidth of each of the primary energy signals of the first spectral signal and the second spectral signals.

在一實施例中,於第四步驟中,係依據第一整合頻寬矩陣及第二整合頻寬矩陣分別整合該等第一頻寬矩陣及該等第二頻寬矩陣,以分別將該等第一頻寬矩陣及該等第二頻寬矩陣中同一主能量頻率的能量累加。In an embodiment, in the fourth step, the first bandwidth matrix and the second bandwidth matrix are separately integrated according to the first integrated bandwidth matrix and the second integrated bandwidth matrix, respectively, to respectively The first bandwidth matrix and the energy of the same primary energy frequency in the second bandwidth matrix are accumulated.

在一實施例中,於第五步驟中,係將第一能量頻譜中,與第二能量頻譜具有相同頻率之能量頻譜去除,以建立能量訊號的主能量頻率之主能量頻譜。In an embodiment, in the fifth step, the energy spectrum of the first energy spectrum having the same frequency as the second energy spectrum is removed to establish a main energy spectrum of the main energy frequency of the energy signal.

在一實施例中,於第五步驟中,當建立能量訊號的主能量頻率之主能量頻譜時,需以一固定資料庫對主能量頻譜進行修正。In an embodiment, in the fifth step, when the main energy spectrum of the main energy frequency of the energy signal is established, the main energy spectrum needs to be corrected by a fixed database.

在一實施例中,當建立能量訊號的主能量頻率之主能量頻譜時,若第一能量頻譜上的某個主能量頻率在第二能量頻譜上找不到對應的頻率時,則將固定資料庫之主能量頻率的頻寬導入,以取代並修正主能量頻率之主能量頻譜的頻寬。In an embodiment, when the main energy spectrum of the main energy frequency of the energy signal is established, if a certain main energy frequency in the first energy spectrum cannot find a corresponding frequency on the second energy spectrum, the fixed data is fixed. The bandwidth of the main energy frequency of the library is introduced to replace and correct the bandwidth of the main energy spectrum of the main energy frequency.

在一實施例中,頻譜分析方法更包括一第六步驟:依據第一整合頻寬矩陣及第二整合頻寬矩陣,分別建立一第一累加頻寬頻譜及一第二累加頻寬頻譜。In an embodiment, the spectrum analysis method further includes a sixth step of: establishing a first accumulated bandwidth spectrum and a second accumulated bandwidth spectrum according to the first integrated bandwidth matrix and the second integrated bandwidth matrix, respectively.

在一實施例中,於第六步驟中,係分別將各第一頻譜訊號及各第二頻譜訊號之各主能量頻率的左頻寬平方加右頻寬平方,以形成一新頻寬指標,並分別將該等第一頻譜訊號及該等第二頻譜訊號之頻譜中,各主能量頻率的新頻寬指標相加,以分別建立該等第一頻譜訊號及該等第二頻譜訊號的第一累加頻寬頻譜及第二累加頻寬頻譜,其中,第一累加頻寬頻譜及第二累加頻寬頻譜分別具有新頻寬指標相加的頻寬值。In an embodiment, in the sixth step, the square of the left bandwidth and the square of the right bandwidth of each main energy frequency of each of the first spectrum signal and each of the second spectrum signals are respectively squared to form a new bandwidth index. And adding the new bandwidth indices of the primary energy signals and the spectrums of the second spectral signals to the first spectrum signals and the second spectrum signals respectively An accumulated bandwidth spectrum and a second accumulated bandwidth spectrum, wherein the first accumulated bandwidth spectrum and the second accumulated bandwidth spectrum respectively have a bandwidth value added by the new bandwidth indicator.

在一實施例中,頻譜分析方法更包括一第七步驟:找出第一累加頻寬頻譜之新頻寬指標相加的頻寬值大於或等於2倍的第二累加頻寬頻譜之主能量頻率之新頻寬指標相加的頻寬值,以建立一主能量頻寬比例頻譜。In an embodiment, the spectrum analysis method further comprises a seventh step of: finding a main energy of the second accumulated bandwidth spectrum whose bandwidth value added by the new bandwidth index of the first accumulated bandwidth spectrum is greater than or equal to 2 times The new bandwidth of the frequency is added to the bandwidth value to establish a main energy bandwidth proportional spectrum.

為達上述目的,依據本發明之一種疾病篩檢方法係用以篩檢一疾病,疾病篩檢方法包括如上所述之頻譜分析方法,疾病篩檢方法包括:對具有疾病的複數病患及複數正常人分別進行檢測,以分別得到該等病患及該等正常人之該等能量訊號;分別以上述之頻譜分析方法分析該等病患及該等正常人之該等能量訊號,以分別得到該等病患及該等正常人之該等主能量頻譜及該等主能量頻寬比例頻譜;分別對應比較各病患與各正常人之各主能量頻譜及各主能量頻寬比例頻譜,以分別找出該等主能量頻譜及該等主能量頻寬比例頻譜中,只存在該等病患之疾病的複數頻率;以及依據該等主能量頻譜及該等主能量頻寬比例頻譜中只存在該等病患之疾病的該等頻率,分別計算疾病之各頻率的重覆比例。To achieve the above object, a disease screening method according to the present invention is for screening a disease, and the disease screening method includes the spectrum analysis method as described above, and the disease screening method includes: a plurality of diseases and plurals having diseases The normal persons are separately tested to obtain the energy signals of the patients and the normal persons respectively; the energy signals of the patients and the normal persons are respectively analyzed by the above-mentioned spectrum analysis methods to obtain respectively The main energy spectrum of the patients and the normal persons and the spectrum of the main energy bandwidth ratios; respectively, comparing the main energy spectrum of each patient and each normal person with the spectrum of the main energy bandwidth ratio, Finding the main energy spectrum and the main energy bandwidth ratio spectrum separately, only the complex frequencies of the diseases of the patients exist; and only the spectrum exists according to the main energy spectrum and the ratio of the main energy bandwidths The frequencies of the diseases of the patients are calculated as the repetition ratios of the frequencies of the diseases, respectively.

在一實施例中,於檢測步驟中,係檢測一天或複數天,且每天檢測至少一次。In one embodiment, in the detecting step, one or more days are detected and detected at least once a day.

在一實施例中,疾病篩檢方法更包括:當一被檢測者進行檢測,並以上述之頻譜分析方法分析後,被檢測者具有疾病之該等頻率的至少其中之一,且頻率的重覆比例越高者,則被檢測者得到疾病之機率越高。In an embodiment, the disease screening method further comprises: when a test subject performs the test and analyzes by the above-mentioned spectrum analysis method, the test subject has at least one of the frequencies of the disease, and the frequency is heavy. The higher the proportion, the higher the probability that the subject will get the disease.

在一實施例中,疾病篩檢方法更包括:當一被檢測者進行檢測,並以上述之頻譜分析方法分析後,被檢測者具有疾病之該等頻率的至少其中之一,且頻率的重覆比例越低者,則被檢測者得到疾病之機率越低。In an embodiment, the disease screening method further comprises: when a test subject performs the test and analyzes by the above-mentioned spectrum analysis method, the test subject has at least one of the frequencies of the disease, and the frequency is heavy. The lower the ratio, the lower the probability that the subject will get the disease.

在一實施例中,疾病篩檢方法更包括:當一被檢測者進行檢測,並以上述之頻譜分析方法分析後,被檢測者沒有疾病之該等頻率的其中之一,則被檢測者得到疾病之機率越低。In one embodiment, the disease screening method further comprises: when a subject is detected and analyzed by the above-described spectrum analysis method, and the subject does not have one of the frequencies of the disease, the subject obtains The lower the chance of disease.

承上所述,本發明提供一種頻譜分析方法可針對不易分析、且雜亂的(微)能量訊號的頻譜進行分析,並可將背景雜訊盡可能地排除(即將背景雜訊之能量頻譜去除),以得到真正可代表該能量訊號的主能量頻率及其能量。藉由本分析方法,可得到該能量訊號的主能量頻率之主能量頻譜及主能量頻寬比例頻譜。另外,可將本發明之頻譜分析方法得到的主能量頻譜及主能量頻寬比例頻譜應用於疾病篩檢,以進行人體病理資訊上的快速篩檢,並可進行各種疾病之可能性與潛在性的分析。In view of the above, the present invention provides a spectrum analysis method for analyzing the spectrum of a (micro) energy signal that is difficult to analyze and disorder, and can eliminate background noise as much as possible (ie, remove the energy spectrum of the background noise). To get the main energy frequency and its energy that can truly represent the energy signal. By the analysis method, the main energy spectrum and the main energy bandwidth ratio spectrum of the main energy frequency of the energy signal can be obtained. In addition, the main energy spectrum and the main energy bandwidth ratio spectrum obtained by the spectrum analysis method of the present invention can be applied to disease screening for rapid screening of human pathological information, and the possibility and potential of various diseases can be performed. Analysis.

以下將參照相關圖式,說明依本發明較佳實施例之一種頻譜分析方法及疾病篩檢方法,其中相同的元件將以相同的參照符號加以說明。DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS Hereinafter, a spectrum analysis method and a disease screening method according to a preferred embodiment of the present invention will be described with reference to the accompanying drawings, wherein the same elements will be described with the same reference numerals.

請參照圖1所示,其為本發明較佳實施例之一種頻譜分析方法的流程示意圖。先說明的是,本發明之頻譜分析方法係用以分析一能量訊號,該能量訊號可為一微弱之能量訊號,並例如可為一生物體或一非生物體的微能量場所散發出的微能量。其中,本實施例之能量訊號係以一人體所發出之微弱能量訊號為例(可為正常人或患有某種疾病的人所發出之微能量訊號)。由於以儀器量測人體微弱的能量訊號時,因環境雜訊的緣故,量測時會得到一能量訊號及一背景雜訊。於此,係將頭髮(也可為人體的其它組織)利用一光子波動能量分析儀(photon code subtle energy analyzers,PCSEA)來進行量測。其中,當將頭髮放置於光子波動能量分析儀而未調整其共振旋鈕時,得到的訊號即為背景雜訊。另外,當將頭髮放置並調整光子波動能量分析儀之共振旋鈕以達到共振的情況時,所量測到的訊號即為該能量訊號。Please refer to FIG. 1, which is a schematic flowchart of a spectrum analysis method according to a preferred embodiment of the present invention. First, the spectrum analysis method of the present invention is used to analyze an energy signal, which can be a weak energy signal, and can be, for example, a micro energy emitted by a living body or a non-living micro energy field. . The energy signal of the embodiment is exemplified by a weak energy signal emitted by a human body (a micro energy signal emitted by a normal person or a person suffering from a certain disease). When the human body's weak energy signal is measured by the instrument, due to environmental noise, an energy signal and a background noise are obtained during the measurement. Here, the hair (which may also be other tissues of the human body) is measured using a photon code subtle energy analyzer (PCSEA). Among them, when the hair is placed on the photon fluctuation energy analyzer without adjusting its resonance knob, the obtained signal is background noise. In addition, when the hair is placed and the resonance knob of the photon fluctuation energy analyzer is adjusted to reach resonance, the measured signal is the energy signal.

另外,得到該能量訊號與背景雜訊後,再將量測到的能量訊號與背景雜訊先分別進行切割,例如1秒鐘的時間即切割出一筆,共切割100次,即可獲得100筆主要能量訊號,圖2A顯示為其中一筆,其中,圖2A之縱座標為振幅(即能量),而橫座標為時間(秒)。同樣地,也可得到100筆的背景雜訊。當然,也可1秒取2次,共分別取200次或為其它切割方式,使用者可依其需求取不同的切割時間及數量。再者,以本發明之頻譜分析方法進行分析前需先將該等主要能量訊號分別以快速傅立葉轉換(Fourier Transform),以得到複數第一頻譜訊號(共100筆),圖2B即顯示其中一筆第一頻譜訊號。其中,圖2B之縱座標為振幅(即能量),橫座標為頻率(Hz),而頻率範圍為1Hz至1025Hz。當然,於其它的實施態樣中,使用不同的儀器量測,其頻率範圍也可為不同,於此並不加以限定。另外,再將該等背景雜訊分別以快速傅立葉轉換後,得到複數第二頻譜訊號(共100筆,圖未顯示)。之後,再將傅立葉轉換後的複數第一頻譜訊號及複數第二頻譜訊號進行本發明的分析工作。In addition, after obtaining the energy signal and the background noise, the measured energy signal and the background noise are separately cut first, for example, one stroke is cut in one second, and 100 times are cut, and 100 pens are obtained. The main energy signal, shown in Figure 2A, is one of them, where the ordinate of Figure 2A is the amplitude (i.e., energy) and the abscissa is time (seconds). Similarly, 100 background noises are also available. Of course, you can take 2 times in 1 second, take 200 times or other cutting methods. Users can take different cutting times and quantities according to their needs. Furthermore, before performing the analysis by the spectrum analysis method of the present invention, the main energy signals are respectively subjected to Fourier Transform to obtain a plurality of first spectrum signals (100 words in total), and FIG. 2B shows one of them. First spectrum signal. Wherein, the ordinate of FIG. 2B is the amplitude (ie, energy), the abscissa is the frequency (Hz), and the frequency range is 1 Hz to 1025 Hz. Of course, in other implementations, the frequency range may be different using different instrument measurements, which is not limited herein. In addition, the background noises are respectively converted by fast Fourier transform to obtain a plurality of second spectrum signals (100 in total, not shown). Then, the Fourier-converted complex first spectrum signal and the complex second spectrum signal are subjected to the analysis work of the present invention.

本發明之頻譜分析方法包括以下第一步驟S01至第五步驟S05。The spectrum analysis method of the present invention includes the following first step S01 to fifth step S05.

第一步驟S01係為:分別找出各第一頻譜訊號及各第二頻譜訊號之較大能量的頻譜,並依據該等較大能量的頻譜分別對應建立其第一頻寬矩陣(bandwidth matrix,BM)及其第二頻寬矩陣。在本實施例中,係分別找出各第一頻譜訊號(共100筆)及各第二頻譜訊號(共100筆)之前100大(Top 100)能量的頻譜,且一個第一頻譜訊號建立一個第一頻寬矩陣(共100個第一頻寬矩陣),而一個第二頻譜訊號建立一個第二頻寬矩陣(共100個第二頻寬矩陣)。當然,也可取其它數量的較大能量的頻譜,例如前50大或前200大或其它數量。其中,各第一頻寬矩陣及各第二頻寬矩陣分別具有各第一頻譜訊號及各第二頻譜訊號之複數主能量頻率(harmonic),也就是說,這些主能量頻率所具有的能量較高(各主能量頻率有其對應的能量)。The first step S01 is: separately finding a spectrum of a larger energy of each of the first spectrum signal and each of the second spectrum signals, and respectively establishing a first bandwidth matrix according to the spectrum of the larger energy. BM) and its second bandwidth matrix. In this embodiment, the spectrum of the top 100 energy of each of the first spectrum signals (100 lines in total) and the second spectrum signals (100 lines in total) is respectively found, and a first spectrum signal is used to establish a spectrum. A first bandwidth matrix (a total of 100 first bandwidth matrices) and a second spectral signal establishes a second bandwidth matrix (a total of 100 second bandwidth matrices). Of course, other quantities of larger energy spectrum may be taken, such as the top 50 or the first 200 or other quantities. Each of the first bandwidth matrix and each of the second bandwidth matrix respectively has a complex primary energy frequency of each of the first spectral signal and each of the second spectral signals, that is, the energy of the primary energy frequencies is greater than High (each main energy frequency has its corresponding energy).

請參照圖3A所示,其為一頻譜示意圖。為了計算某一頻率的頻寬,於此係取該頻率的最大振幅(能量)中,小於3db時之振幅的另一頻率與該頻率之差為頻寬。例如圖3A中,頻率f之最大能量為X2,則比能量X2小3db之能量X1的頻率為f1,其頻率差為f-f1即為左頻寬,而比能量X2小3db之能量X3的頻率為f2,其頻率差為f2-f即為右頻寬。於此,係利用相位雜訊的概念分析前100大能量的頻譜,並找出各第一頻譜訊號及各第二頻譜訊號之主能量頻率為何。換言之,即區別出各第一頻譜訊號及各第二頻譜訊號中,何為主能量頻率(即能量較大者)、何為雜訊(即能量較小者)後,並建立各自的第一頻寬矩陣及各第二頻寬矩陣,而各第一頻寬矩陣及各第二頻寬矩陣係分別具有各第一頻譜訊號及各第二頻譜訊號之該等主能量頻率及其左、右頻寬。Please refer to FIG. 3A, which is a schematic diagram of a spectrum. In order to calculate the bandwidth of a certain frequency, in the maximum amplitude (energy) of the frequency, the difference between the other frequency of the amplitude less than 3 db and the frequency is the bandwidth. For example, in FIG. 3A, the maximum energy of the frequency f is X2, and the frequency of the energy X1 which is 3 db less than the energy X2 is f1, and the frequency difference is f-f1 is the left bandwidth, and the energy X3 is 3 db less than the energy X2. The frequency is f2, and the frequency difference is f2-f, which is the right bandwidth. Here, the spectrum of the first 100 energy is analyzed by the concept of phase noise, and the main energy frequencies of the first spectrum signals and the second spectrum signals are found out. In other words, distinguishing between the first spectrum signal and each of the second spectrum signals, what is the main energy frequency (ie, the greater energy), and what is the noise (ie, the less energy), and establishes the first a bandwidth matrix and each of the second bandwidth matrix, wherein each of the first bandwidth matrix and each of the second bandwidth matrix respectively have the primary energy frequencies of the first spectrum signal and each of the second spectrum signals and their left and right bandwidth.

如表一所示,其為某一個第一頻譜訊號之第一頻寬矩陣。例如表一之頻率為803Hz的主能量頻率,其左頻寬為2(Hz),而其右頻寬為3(Hz)。另外,頻寬值為0代表該頻率沒有頻寬,也代表該第一頻譜訊號的能量中,並不存在有該頻率及其能量。As shown in Table 1, it is the first bandwidth matrix of a certain first spectrum signal. For example, the frequency of the primary energy of Table 803 is 803 Hz, the left bandwidth is 2 (Hz), and the right bandwidth is 3 (Hz). In addition, a bandwidth value of 0 means that the frequency has no bandwidth, and also represents the energy of the first spectrum signal, and the frequency and its energy do not exist.

再者,第二步驟S02係為:依據各第一頻寬矩陣及各第二頻寬矩陣,分別計算各第一頻譜訊號及各第二頻譜訊號之該等主能量頻率的能量。於此,能量的計算係可正比於能量X1、X2及X3的平方和,即該等主能量頻率的能量可正比於(X12 +X22 +X32 )。故依照此公式,分別計算各100筆之第一頻譜訊號及各100筆之第二頻譜訊號的各主能量頻率的能量(並非能量的絕對值,而是正比的相對值),因此,第一頻譜訊號可得到100個能量頻譜圖,如圖3B所示,其為某一筆第一頻譜訊號的能量頻譜示意圖,而第二頻譜訊號也會得到100個能量頻譜圖(圖未顯示),其中,橫座標為頻率,而縱座標為能量的強度(magnitude)。Furthermore, the second step S02 is: calculating the energy of the main energy frequencies of the first spectrum signal and each of the second spectrum signals according to each of the first bandwidth matrix and each of the second bandwidth matrix. Here, the energy calculation can be proportional to the sum of the squares of the energies X1, X2, and X3, that is, the energy of the main energy frequencies can be proportional to (X1 2 + X2 2 + X3 2 ). Therefore, according to the formula, the energy of each main energy frequency of each of the first 100 spectral signals and the second spectral signals of each of the 100 pens (not the absolute value of the energy, but a proportional relative value) is calculated, and therefore, the first The spectrum signal can obtain 100 energy spectrum diagrams, as shown in FIG. 3B, which is a schematic diagram of the energy spectrum of a certain first spectrum signal, and the second spectrum signal also obtains 100 energy spectrum diagrams (not shown), wherein The abscissa is the frequency and the ordinate is the intensity of the energy.

另外,第三步驟S03係為:分別整合該等第一頻寬矩陣及該等第二頻寬矩陣,以分別建立該等第一頻譜訊號之一第一整合頻寬矩陣及該等第二頻譜訊號之一第二整合頻寬矩陣。其中,第三步驟S03係分別將該等第一頻寬矩陣及該等第二頻寬矩陣進行整合,於此,係以較大能量為判斷標準,判斷不同筆頻譜中,哪些主能量頻率係屬於同一組頻率。因此,第一整合頻寬矩陣及第二整合頻寬矩陣可分別整合該等第一頻譜訊號及該等第二頻譜訊號之各主要能量頻率,並記錄該等第一頻譜訊號及該等第二頻譜訊號之各主要能量頻率的最大左頻寬及最大右頻寬。In addition, the third step S03 is: separately integrating the first bandwidth matrix and the second bandwidth matrix to respectively establish a first integrated bandwidth matrix of the first spectrum signals and the second spectrum. One of the signals is a second integrated bandwidth matrix. The third step S03 is to integrate the first bandwidth matrix and the second bandwidth matrix respectively, and determine which main energy frequency systems in different pen spectra are determined by using a larger energy as a criterion. Belong to the same set of frequencies. Therefore, the first integrated bandwidth matrix and the second integrated bandwidth matrix can respectively integrate the first spectrum signals and the main energy frequencies of the second spectrum signals, and record the first spectrum signals and the second The maximum left bandwidth and the maximum right bandwidth of each of the main energy frequencies of the spectrum signal.

第四步驟S04係為:依據步驟S03之第一整合頻寬矩陣及第二整合頻寬矩陣,以分別將各第一頻譜訊號及各第二頻譜訊號之相同的主能量頻率的能量累加,進而分別建立一第一能量頻譜及一第二能量頻譜。於此,係依據第一整合頻寬矩陣及第二整合頻寬矩陣分別整合該等第一頻寬矩陣及該等第二頻寬矩陣,以分別將該等第一頻寬矩陣及該等第二頻寬矩陣中同一主能量頻率的能量累加。因為第一整合頻寬矩陣及第二整合頻寬矩陣係分別將該等第一頻譜訊號及該等第二頻譜訊號的所有頻譜進行整合,並判斷不同筆頻譜中哪些頻率係屬於同一組。故於步驟S04中,係分別將該等第一頻譜訊號及該等第二頻譜訊號的所有頻譜中,同一組頻率的能量累加,以分別建立該等第一頻譜訊號及該等第二頻譜訊號的之第一能量頻譜(如圖4A所示)及第二能量頻譜(圖4B所示)。於此,第一能量頻譜及第二能量頻譜係分別稱為第一頻譜訊號及第二頻譜訊號之前一百大能量頻譜(Top100 energy spectrum,TES)。The fourth step S04 is: accumulating the energy of the same main energy frequency of each of the first spectrum signal and each of the second spectrum signals according to the first integrated bandwidth matrix and the second integrated bandwidth matrix of step S03, respectively A first energy spectrum and a second energy spectrum are respectively established. In this case, the first bandwidth matrix and the second bandwidth matrix are respectively integrated according to the first integrated bandwidth matrix and the second integrated bandwidth matrix, respectively, to respectively use the first bandwidth matrix and the first The energy of the same main energy frequency in the two-bandwidth matrix is accumulated. The first integrated bandwidth matrix and the second integrated bandwidth matrix respectively integrate the first spectrum signals and all the spectrums of the second spectrum signals, and determine which frequencies in the different pen spectrum belong to the same group. Therefore, in step S04, the energy of the same group of frequencies in the first spectrum signals and the second spectrum signals are respectively accumulated to establish the first spectrum signals and the second spectrum signals respectively. The first energy spectrum (shown in Figure 4A) and the second energy spectrum (shown in Figure 4B). Here, the first energy spectrum and the second energy spectrum are respectively referred to as a first spectrum signal and a top 100 energy spectrum (TES) before the second spectrum signal.

另外,第五步驟S05係為:依據第一能量頻譜(即第一頻譜訊號之TES)及第二能量頻譜(即第二頻譜訊號之TES),建立該能量訊號的主能量頻率之一主能量頻譜。於此,係將第一能量頻譜中,與第二能量頻譜具有相同頻率之能量頻譜去除(即將背景雜訊之能量頻譜去除),以建立該能量訊號的主能量頻率之主能量頻譜。如圖5所示,其為該能量訊號的主能量頻率之主能量頻譜示意圖。其中,圖5已將背景雜訊(排除了第二頻譜訊號),也就是只顯示該能量訊號的主能量頻率之主能量頻譜。換言之,係從第一頻譜訊號之Top100能量頻譜中,挑出同時出現於第二頻譜訊號的主能量頻率之頻率並去除之(即將第一頻譜訊號之TES中有第二頻譜訊號的主能量頻率之頻率及其能量去除),以建立該能量訊號的主能量頻率之主能量頻譜。In addition, the fifth step S05 is: establishing one main energy of the main energy frequency of the energy signal according to the first energy spectrum (ie, the TES of the first spectrum signal) and the second energy spectrum (ie, the TES of the second spectrum signal) Spectrum. In this case, the energy spectrum of the first energy spectrum having the same frequency as the second energy spectrum is removed (ie, the energy spectrum of the background noise is removed) to establish a main energy spectrum of the main energy frequency of the energy signal. As shown in FIG. 5, it is a schematic diagram of the main energy spectrum of the main energy frequency of the energy signal. Among them, Figure 5 has background noise (excluding the second spectrum signal), that is, only the main energy spectrum of the main energy frequency of the energy signal. In other words, from the Top 100 energy spectrum of the first spectrum signal, the frequency of the main energy frequency simultaneously appearing at the second spectrum signal is picked and removed (ie, the main energy frequency of the second spectrum signal in the TES of the first spectrum signal) The frequency and its energy removal) to establish the main energy spectrum of the main energy frequency of the energy signal.

值得注意的是,為了使本發明更可準確分析,需要有另一修正步驟,以修正第五步驟S05之該能量訊號的主能量頻率之主能量頻譜。It is worth noting that in order to make the invention more accurate, another correction step is needed to correct the main energy spectrum of the main energy frequency of the energy signal of the fifth step S05.

因此,於第五步驟S05中,當建立該能量訊號的主能量頻率之主能量頻譜時,需以一固定資料庫(fixing database)對主能量頻譜進行修正。而修正原因為:由於量測的誤差可能使主能量的頻率之變動大於主能量頻率本身的頻寬(變動超過左、右頻寬),以至於無法以相位雜訊的概念找出正確的主能量頻率,故建構一固定資料庫來解決由量測所造成的頻率變動(即誤差)。而此誤差是因人為操作中,原有雜訊的頻率飄移而變成主要能量的頻率。Therefore, in the fifth step S05, when the main energy spectrum of the main energy frequency of the energy signal is established, the main energy spectrum needs to be corrected by a fixing database. The reason for the correction is that the error of the measurement may cause the frequency of the main energy to be greater than the bandwidth of the main energy frequency itself (variation exceeds the left and right bandwidths), so that the correct master cannot be found by the concept of phase noise. The energy frequency, so a fixed database is constructed to solve the frequency variation (ie, error) caused by the measurement. This error is due to the frequency at which the frequency of the original noise drifts and becomes the main energy during human operation.

在本實施例中,固定資料庫的建構方法是:以上述的方法取得複數組能量訊號之複數第一頻譜訊號(於此係以5組為例,當然也可為其它數量,而每組能量訊號可分別具有100筆第一頻譜訊號及100筆第二頻譜訊號),於此稱為群組R1~R5,並依照上述的分析步驟S01~S05分別得到群組R1~R5各自的能量頻譜(於此稱為TES1~TES5)。從群組R1~R5中任取兩個不同的能量頻譜(TES)來做比較,如此,可組成10對的能量頻譜(TES)比較組。例如TESR1與TESR2比較(或TESR2與TESR3比較等等),找出在TESR1與TESR2中沒有彼此相關連的頻率,即只單獨存在TESR1或TESR2中,而非兩者共有的頻率,這些頻率的頻寬即是由於測量誤差所造成的,導致無法在別的能量頻譜(TES)上找到相關連的頻率,觀察並紀錄這些頻率的頻寬變異。於此,係以每100 Hz為間隔將10組比較後所得之頻寬變異作成固定資料庫(即找出誤差最大之左、右頻寬的最大者),也就是固定資料庫具有各100 Hz之頻率間隔的最大之左、右頻寬。因此,固定資料庫可得到下表二所示。In this embodiment, the method for constructing the fixed database is to obtain a plurality of first spectrum signals of the complex array energy signals by the above method (for example, five groups are used as examples, and of course, other quantities may be used, and each group of energy is used. The signals may have 100 first spectrum signals and 100 second spectrum signals respectively, which are referred to herein as groups R1 to R5, and the respective energy spectra of the groups R1 to R5 are obtained according to the above analysis steps S01 to S05 ( This is referred to herein as TES1 to TES5). Two different energy spectra (TES) are taken from the groups R1 to R5 for comparison, so that 10 pairs of energy spectrum (TES) comparison groups can be formed. For example, TESR1 is compared with TESR2 (or TESR2 is compared with TESR3, etc.), and the frequencies that are not associated with each other in TESR1 and TESR2 are found, that is, only TESR1 or TESR2 exist alone, not the frequencies shared by the two, and the frequencies of these frequencies. The width is due to measurement errors, making it impossible to find the associated frequencies on other energy spectra (TES), observing and recording the bandwidth variation of these frequencies. Here, the bandwidth variation obtained by comparing 10 groups is made into a fixed database at intervals of 100 Hz (that is, the largest one of the left and right bandwidths with the largest error is found), that is, the fixed database has 100 Hz each. The maximum left and right bandwidth of the frequency interval. Therefore, the fixed database can be obtained as shown in Table 2 below.

接著,當建立該能量訊號的主能量頻率之主能量頻譜時,若第一能量頻譜上的某個主能量頻率在第二能量頻譜上找不到對應的頻率時,則將固定資料庫之主能量頻率的左、右頻寬導入,以取代並修正該能量訊號的主能量頻率之主能量頻譜的頻寬值。如此,可修正因量測誤差而使第二頻譜訊號的頻率移至第一頻譜訊號中,可使本發明之分析方法更準確。Then, when the main energy spectrum of the main energy frequency of the energy signal is established, if a certain main energy frequency in the first energy spectrum cannot find a corresponding frequency on the second energy spectrum, the owner of the fixed database will be fixed. The left and right bandwidths of the energy frequency are introduced to replace and correct the bandwidth value of the main energy spectrum of the main energy frequency of the energy signal. In this way, the frequency of the second spectrum signal can be corrected to be shifted into the first spectrum signal due to the measurement error, so that the analysis method of the present invention can be more accurate.

另外,請參照圖6所示,其為本發明之頻譜分析方法的另一流程示意圖。In addition, please refer to FIG. 6, which is another schematic flowchart of the spectrum analysis method of the present invention.

除了上述之第一步驟S01至第五步驟S05之外,本發明之頻譜分析方法更可包第六步驟S06及第七步驟S07。In addition to the first step S01 to the fifth step S05 described above, the spectrum analysis method of the present invention may further include the sixth step S06 and the seventh step S07.

第六步驟S06係為:依據第三步驟S03之第一整合頻寬矩陣及第二整合頻寬矩陣,分別建立一第一累加頻寬頻譜及一第二累加頻寬頻譜(於此,累加頻寬頻稱為Top100 cumulative bandwidth spectrum,TCBS)。在本實施例中,係分別將各第一頻譜訊號及各第二頻譜訊號之各主能量頻率的左頻寬平方加上右頻寬平方,以形成一新頻寬指標,並分別將該等第一頻譜訊號及該等第二頻譜訊號之頻譜中,各主能量頻率的新頻寬指標再相加,以分別建立該等第一頻譜訊號及該等第二頻譜訊號的第一累加頻寬頻譜(如圖7A所示)及第二累加頻寬頻譜(如圖7B所示)。其中,第一累加頻寬頻譜及第二累加頻寬頻譜分別具有該新頻寬指標相加的頻寬值。The sixth step S06 is: according to the first integrated bandwidth matrix and the second integrated bandwidth matrix of the third step S03, respectively establishing a first accumulated bandwidth spectrum and a second accumulated bandwidth spectrum (here, accumulating frequency Broadband is called Top100 cumulative bandwidth spectrum, TCBS). In this embodiment, the square of the left bandwidth of each main energy frequency of each of the first spectrum signal and each of the second spectrum signals is respectively added to the square of the right bandwidth to form a new bandwidth index, and respectively In the spectrum of the first spectrum signal and the second spectrum signals, the new bandwidth indices of the main energy frequencies are further added to respectively establish the first spectrum signal and the first accumulated frequency broadband of the second spectrum signals. The spectrum (shown in Figure 7A) and the second accumulated bandwidth spectrum (as shown in Figure 7B). The first accumulated bandwidth spectrum and the second accumulated bandwidth spectrum respectively have bandwidth values added by the new bandwidth indicator.

此外,第七步驟S07係為:找出步驟S06之第一累加頻寬頻譜之該新頻寬指標相加的頻寬值大於或等於2倍的第二累加頻寬頻譜之主能量頻率之該新頻寬指標相加的頻寬值,以建立一主能量頻寬比例頻譜(bandwidth ratio spectrum)。換言之,即找出該等第一頻譜訊號中,新頻寬指標相加的頻寬值,比第二累加頻寬頻譜之該新頻寬指標相加的頻寬值2倍(含)以上的頻寬值後,以建立該能量訊號之主能量頻寬比例頻譜,其結果如圖8所示。In addition, the seventh step S07 is: finding the main energy frequency of the second accumulated bandwidth spectrum whose bandwidth value of the first added bandwidth spectrum of the first accumulated bandwidth spectrum of step S06 is greater than or equal to 2 times. The bandwidth value of the new bandwidth index is added to establish a bandwidth ratio spectrum. In other words, the bandwidth value of the new bandwidth index added in the first spectrum signal is found to be more than 2 times (inclusive) the bandwidth value of the new bandwidth index of the second accumulated bandwidth spectrum. After the bandwidth value, the main energy bandwidth ratio spectrum of the energy signal is established, and the result is shown in FIG. 8.

承上,本發明之頻譜分析方法係用以分析一能量訊號(尤其是微能量訊號),並將其背景雜訊盡可能地排除(即將背景雜訊之能量頻譜去除),以得到真正可代表該能量訊號的主能量頻率及其能量。經過本分析方法後,會得到第五步驟S05之主能量頻譜(如圖5所示,於此,若經過上述之固定資料庫修正後可較為準確)及第七步驟S07之主能量頻寬比例頻譜(如圖8所示),並將上述之頻譜分析方法及得到的主能量頻譜及主能量頻寬比例頻譜應用於以下的疾病篩檢方法,以應用在病理資訊上的快速篩檢。The spectrum analysis method of the present invention is for analyzing an energy signal (especially a micro energy signal) and excluding its background noise as much as possible (that is, removing the energy spectrum of the background noise) to obtain a truly representative The main energy frequency of the energy signal and its energy. After the analysis method, the main energy spectrum of the fifth step S05 is obtained (as shown in FIG. 5, which may be more accurate after being corrected by the fixed database described above) and the main energy bandwidth ratio of the seventh step S07. The spectrum (as shown in Fig. 8), and the above spectrum analysis method and the obtained main energy spectrum and main energy bandwidth ratio spectrum are applied to the following disease screening methods for rapid screening of pathological information.

請參照圖9所示,其為本發明之一種疾病篩檢方法的流程示意圖。Please refer to FIG. 9, which is a schematic flow chart of a disease screening method of the present invention.

本發明之疾病篩檢方法係以一種非侵入式的方式對某一種疾病進行快速篩檢,可提供一被檢測者是否得到該種疾病可能比例的高低。先說明的是,本發明之疾病篩檢方法係可應用上述的頻譜分析方法分析生物體所具有的能量訊號(微能量訊號)。The disease screening method of the present invention provides a rapid screening of a disease in a non-invasive manner, providing a possible proportion of whether or not the subject is getting the disease. First, the disease screening method of the present invention can analyze the energy signal (micro energy signal) possessed by the living body by using the above spectrum analysis method.

本發明對之疾病篩檢方法包括以下步驟P01~步驟P04。The disease screening method of the present invention includes the following steps P01 to P04.

步驟P01係為:對具有某一疾病的複數病患及複數正常人分別進行檢測,以分別得到該等病患及該等正常人之該等能量訊號。其中,係可分別檢測一天或複數天,且每天檢測至少一次,以得到一天或複數天數的複數病患及複數正常人的微能量訊號。在本實施例中,篩檢的疾病係以糖尿病為例(於此只是舉例,當然也可以本方法篩檢其它的疾病,例如肝病、高血壓等等)。另外,再分別檢測10個糖尿病患者及10個正常人的微能量訊號(共20人),且一天檢測一次,共檢測三天,故總共得到60筆微能量訊號(30筆正常人的微能量訊號,30筆糖尿病患者的微能量訊號),且每筆微能量訊號又可分別具有複數第一頻譜訊號及複數第二頻譜訊號(背景雜訊)。其中,上述的人數、天數及每天的檢測次數均為舉例,當然,也可使用不同的人數、不同的天數及每天不同的檢測次數。其中,當人數、天數及檢測次數越多時(即樣本數越多時),其後續的疾病篩檢結果將越準確。Step P01 is: detecting a plurality of patients having a certain disease and a plurality of normal persons to obtain the energy signals of the patients and the normal persons respectively. Among them, one or more days can be detected separately, and at least once a day to obtain a micro energy signal of a plurality of patients and a plurality of normal persons in one day or a plurality of days. In the present embodiment, the disease to be screened is exemplified by diabetes (here is merely an example, and of course, other diseases such as liver disease, hypertension, etc.) can be screened by the method. In addition, 10 micro-energy signals (20 people) were detected in 10 diabetic patients and 10 normal people, respectively, and tested once a day for a total of three days, so a total of 60 micro-energy signals (30 normal human micro-energy) were obtained. Signal, 30 micro-energy signals of diabetic patients, and each micro-energy signal can have a plurality of first spectrum signals and a plurality of second spectrum signals (background noise). The number of people, the number of days, and the number of tests per day are all examples. Of course, different numbers of people, different days, and different number of tests per day can be used. Among them, the more the number of people, the number of days and the number of tests (ie, the more the number of samples), the more accurate the subsequent disease screening results will be.

再者,步驟P02係為:分別以上述之頻譜分析方法分析該等病患及該等正常人之該等能量訊號,以分別得到該等病患及該等正常人之該等主能量頻譜及該等主能量頻寬比例頻譜。於此,係以上述的頻譜分析方法分析,並可得到60個的主能量頻譜及60個主能量頻寬比例頻譜。In addition, the step P02 is to analyze the energy signals of the patients and the normal persons by using the above-mentioned spectrum analysis method to obtain the main energy spectrums of the patients and the normal persons, respectively. The main energy bandwidth ratio spectrum. Here, the spectrum analysis method described above is used, and 60 main energy spectra and 60 main energy bandwidth ratio spectra are obtained.

另外,步驟P03係為:分別對應比較各病患與各正常人之各主能量頻譜及各主能量頻寬比例頻譜,以分別找出該等主能量頻譜及該等主能量頻寬比例頻譜中,只存在該等病患之疾病的複數頻率。在本實施例中,係將一天的檢測資料中(即主能量頻譜及主能量頻寬比例頻譜),找出一個糖尿病患者與一個正常人之間只存在於此糖尿病患者的頻率,且重覆三天的此步驟,再從這三天中找出出現於糖尿病患者兩天以上的頻率,並製作成糖尿病患者之主能量頻譜出現頻率的表格。於此,可得到下表三及下表四所示的頻率表格。因此,20人共可分為10組,共可得到10個表三的表格及10個表四的表格。其中,表三及表四中出現的數字代表這三天的檢測資料中,具有糖尿病的病患之(微)能量訊號中出現兩天以上之主能量頻譜的頻率(表示可能是糖尿病患者身上特有的頻率)。In addition, the step P03 is: respectively comparing the main energy spectrum and the main energy bandwidth ratio spectrum of each patient and each normal person to respectively find the main energy spectrum and the main energy bandwidth ratio spectrum. There are only multiple frequencies of the disease in these patients. In this embodiment, the frequency of one day of detection (ie, the main energy spectrum and the main energy bandwidth ratio spectrum) is used to find out the frequency of a diabetic patient and a normal person only present in the diabetic patient, and repeat In this three-day step, the frequency of occurrence of diabetes patients for more than two days is found out from these three days, and a table showing the frequency of occurrence of the main energy spectrum of diabetic patients is made. Here, the frequency tables shown in Table 3 below and Table 4 below can be obtained. Therefore, 20 people can be divided into 10 groups, and a total of 10 tables in Table 3 and 10 tables in Table 4 can be obtained. Among them, the figures appearing in Tables 3 and 4 represent the frequency of the main energy spectrum of more than two days in the (micro) energy signal of patients with diabetes in the three-day test data (indicating that it may be specific to diabetic patients) Frequency of).

此外,步驟P04係為:依據該等主能量頻譜及該等主能量頻寬比例頻譜中只存在該等病患之該疾病的該等頻率,分別計算該疾病之各頻率的重覆比例。在本實施例中,係從上述之10個表三的表格及10個表四的表格10中,計算其主能量頻率(於此係為糖尿病的主能量頻率)出現的重覆比例百分比,並製作成下表五所示。其中,表格五的重覆比例百分比代表這一個頻率在上述所有表三及表四的資料中,出現的比例。例如8Hz的重覆出現的比例為20%,代表所有的表三及表四的頻率中,100次中出現的比例為20次,以此類推。In addition, step P04 is: calculating the repetition ratio of each frequency of the disease according to the main energy spectrum and the frequencies of the diseases in which only the patients are present in the main energy bandwidth spectrum. In this embodiment, the percentage of the repetition rate of the main energy frequency (the main energy frequency of the diabetes) is calculated from the table of the above 10 Table 3 and the table 10 of the 10 Table 4, and Make the following table five. Among them, the percentage of repeated proportions in Table 5 represents the proportion of this frequency appearing in all the data in Tables 3 and 4 above. For example, the repetition rate of 8Hz appears to be 20%, which represents the frequency of all the three tables and four, the ratio of occurrence in 100 times is 20, and so on.

另外,疾病篩檢方法更可包括:當一被檢測者進行檢測,並以上述之頻譜分析方法分析後,被檢測者具有糖尿病之該等頻率(即表五中所列之頻率)的至少其中之一,且該頻率的重覆比例越高者,則被檢測者得到糖尿病的機率越高。In addition, the disease screening method may further include: when a test subject performs the test and analyzes by the above-mentioned spectrum analysis method, the test subject has at least one of the frequencies of diabetes (ie, the frequencies listed in Table 5). One, and the higher the repetition ratio of the frequency, the higher the probability that the subject gets diabetes.

舉例而言,例如一被檢測者經過上述的頻譜分析方法得到的主能量頻譜及主能量頻寬比例頻譜之頻率中有1073Hz(其重覆比例為60%)的話,則該被檢測者得到糖尿病的機率將相當高(可能有60%)。因此,被檢測者需再進一步進行其它的傳統檢測方式,以進一步釐清是否真的患有糖尿病。For example, if a subject is subjected to the above-described spectrum analysis method and the frequency of the main energy spectrum and the main energy bandwidth ratio spectrum is 1073 Hz (the repetition ratio is 60%), the subject receives diabetes. The chances are quite high (possibly 60%). Therefore, the testee needs to further carry out other traditional detection methods to further clarify whether or not there is diabetes.

反之,疾病篩檢方法更可包括:當一被檢測者進行檢測,並以上述之頻譜分析方法分析後,被檢測者具有糖尿病之該等頻率的至少其中之一,且該頻率的重覆比例越低者,則被檢測者得到疾病之機率將越低。例如表五之125Hz的頻率,其重覆比例只有10%,表示被檢測者得到糖尿病的機率較低(較低不表示一定沒有糖尿病)。Conversely, the disease screening method may further include: when a test subject performs the test and analyzes by the above-mentioned spectrum analysis method, the test subject has at least one of the frequencies of diabetes, and the repetition ratio of the frequency The lower the rate, the lower the chance that the subject will get the disease. For example, the frequency of 125 Hz in Table 5 is only 10%, indicating that the probability of getting diabetes in the testee is low (lower means no diabetes).

此外,疾病篩檢方法更可包括:當一被檢測者進行檢測,並以上述之頻譜分析方法分析後,被檢測者沒有糖尿病之該等頻率的其中之一,則被檢測者得到糖尿病的機率也相當低(初步可判定沒有糖尿病)。In addition, the disease screening method may further include: when a test subject performs the test and analyzes by the above-mentioned spectrum analysis method, and the test subject does not have one of the frequencies of diabetes, the probability of the testee getting diabetes It is also quite low (previously it can be judged that there is no diabetes).

綜上所述,本發明提供一種頻譜分析方法可針對不易分析、且雜亂的(微)能量訊號的頻譜進行分析,並可將背景雜訊盡可能地排除(即將背景雜訊之能量頻譜去除),以得到真正可代表該能量訊號的主能量頻率及其能量。藉由本分析方法,可得到該能量訊號的主能量頻率之主能量頻譜及主能量頻寬比例頻譜。另外,可將本發明之頻譜分析方法得到的主能量頻譜及主能量頻寬比例頻譜應用於疾病篩檢,以進行人體病理資訊上的快速篩檢,並可進行各種疾病之可能性與潛在性的分析。In summary, the present invention provides a spectrum analysis method for analyzing the spectrum of a (micro) energy signal that is difficult to analyze and disorder, and can eliminate background noise as much as possible (ie, removing the energy spectrum of the background noise). To get the main energy frequency and its energy that can truly represent the energy signal. By the analysis method, the main energy spectrum and the main energy bandwidth ratio spectrum of the main energy frequency of the energy signal can be obtained. In addition, the main energy spectrum and the main energy bandwidth ratio spectrum obtained by the spectrum analysis method of the present invention can be applied to disease screening for rapid screening of human pathological information, and the possibility and potential of various diseases can be performed. Analysis.

以上所述僅為舉例性,而非為限制性者。任何未脫離本發明之精神與範疇,而對其進行之等效修改或變更,均應包含於後附之申請專利範圍中。The above is intended to be illustrative only and not limiting. Any equivalent modifications or alterations to the spirit and scope of the invention are intended to be included in the scope of the appended claims.

S01~S07、P01~P04‧‧‧步驟S01~S07, P01~P04‧‧‧ steps

圖1為本發明較佳實施之一種頻譜分析方法的流程示意圖;圖2A為一主要能量訊號示意圖;圖2B為一第一頻譜訊號的示意圖;圖3A為一頻譜的示意圖;圖3B為第一頻譜訊號的能量頻譜示意圖;圖4A及圖4B分別為第一能量頻譜示意圖及第二能量頻譜示意圖;圖5為能量訊號的主能量頻率之主能量頻譜示意圖;圖6為本發明之頻譜分析方法的另一流程示意圖;圖7A及圖7B分別為第一累加頻寬頻譜示意圖及第二累加頻寬頻譜示意圖;圖8為能量訊號之主能量頻寬比例頻譜示意圖;以及圖9為本發明之一種疾病篩檢方法的流程示意圖。1 is a schematic flowchart of a spectrum analysis method according to a preferred embodiment of the present invention; FIG. 2A is a schematic diagram of a main energy signal; FIG. 2B is a schematic diagram of a first spectrum signal; FIG. 3A is a schematic diagram of a spectrum; Schematic diagram of the energy spectrum of the spectrum signal; FIG. 4A and FIG. 4B are respectively a schematic diagram of the first energy spectrum and a second energy spectrum; FIG. 5 is a schematic diagram of the main energy spectrum of the main energy frequency of the energy signal; FIG. 6 is a spectrum analysis method of the present invention; FIG. 7A and FIG. 7B are respectively a schematic diagram of a first accumulated bandwidth spectrum and a second accumulated bandwidth spectrum; FIG. 8 is a schematic diagram of a main energy bandwidth proportional spectrum of an energy signal; and FIG. 9 is a schematic diagram of the present invention A schematic diagram of a process for a disease screening method.

S01~S05...步驟S01~S05. . . step

Claims (15)

一種頻譜分析方法,係用以分析一能量訊號,該能量訊號具有複數第一頻譜訊號及複數第二頻譜訊號,該頻譜分析方法包括:一第一步驟:分別找出各該第一頻譜訊號及各該第二頻譜訊號之較大能量的頻譜,並依據該等較大能量的頻譜分別對應建立其第一頻寬矩陣及其第二頻寬矩陣,其中各該第一頻寬矩陣及各該第二頻寬矩陣分別具有各該第一頻譜訊號及各該第二頻譜訊號之複數主能量頻率;一第二步驟:依據各該第一頻寬矩陣及各該第二頻寬矩陣,分別計算各該第一頻譜訊號及各該第二頻譜訊號之該等主能量頻率的能量,其中該等主能量頻率的能量分別正比於能量X1、X2及X3的平方和,X2為各該第一頻譜訊號或各該第二頻譜訊號的最大能量值,且X1及X3分別為比X2小一特定值的能量值;一第三步驟:分別整合該等第一頻寬矩陣及該等第二頻寬矩陣,以分別建立該等第一頻譜訊號之一第一整合頻寬矩陣及該等第二頻譜訊號之一第二整合頻寬矩陣;一第四步驟:依據該第一整合頻寬矩陣及該第二整合頻寬矩陣,以分別將各該第一頻譜訊號及各該第二頻譜訊號之相同的主能量頻率的能量累加,進而分 別建立一第一能量頻譜及一第二能量頻譜;以及一第五步驟:依據該第一能量頻譜及該第二能量頻譜,建立該能量訊號的主能量頻率之一主能量頻譜。 A spectrum analysis method for analyzing an energy signal having a plurality of first spectrum signals and a plurality of second spectrum signals, the spectrum analysis method comprising: a first step: separately finding each of the first spectrum signals and a spectrum of a larger energy of each of the second spectral signals, and correspondingly establishing a first bandwidth matrix and a second bandwidth matrix thereof according to the spectra of the larger energy, wherein each of the first bandwidth matrix and each of the first bandwidth matrix The second bandwidth matrix has a complex primary energy frequency of each of the first spectral signal and each of the second spectral signals; and a second step: calculating respectively according to each of the first bandwidth matrix and each of the second bandwidth matrix The energy of the primary energy signals of the first spectral signal and each of the second spectral signals, wherein the energy of the primary energy frequencies is proportional to the sum of the squares of the energy X1, X2, and X3, respectively, and X2 is the first frequency spectrum. The maximum energy value of the signal or each of the second spectral signals, and X1 and X3 are respectively smaller than X2 by a specific value; a third step: respectively integrating the first bandwidth matrix and the second bandwidth matrix, Establishing a first integrated bandwidth matrix of the first spectrum signals and a second integrated bandwidth matrix of the second spectrum signals respectively; a fourth step: according to the first integrated bandwidth matrix and the second integration a bandwidth matrix for accumulating the energy of the same main energy frequency of each of the first spectrum signal and each of the second spectrum signals, and further dividing Do not establish a first energy spectrum and a second energy spectrum; and a fifth step: establishing a primary energy spectrum of the primary energy frequency of the energy signal according to the first energy spectrum and the second energy spectrum. 如申請專利範圍第1項所述之頻譜分析方法,其中於該第一步驟中,各該第一頻寬矩陣及各該第二頻寬矩陣分別具有各該第一頻譜訊號及各該第二頻譜訊號之該等主能量頻率及其頻寬。 The spectrum analysis method of claim 1, wherein in the first step, each of the first bandwidth matrix and each of the second bandwidth matrix has each of the first spectrum signals and each of the second The main energy frequencies of the spectrum signal and their bandwidth. 如申請專利範圍第1項所述之頻譜分析方法,其中於該第三步驟中,該第一整合頻寬矩陣及該第二整合頻寬矩陣係分別整合該等第一頻譜訊號及該等第二頻譜訊號之各主要能量頻率,並記錄該等第一頻譜訊號及該等第二頻譜訊號之各該主要能量頻率的最大左頻寬及最大右頻寬。 The spectrum analysis method of claim 1, wherein in the third step, the first integrated bandwidth matrix and the second integrated bandwidth matrix respectively integrate the first spectrum signals and the first The main energy frequencies of the two spectral signals, and recording the maximum left bandwidth and the maximum right bandwidth of the primary energy signals and the primary energy frequencies of the second spectral signals. 如申請專利範圍第1項所述之頻譜分析方法,其中於該第四步驟中,係依據該第一整合頻寬矩陣及該第二整合頻寬矩陣分別整合該等第一頻寬矩陣及該等第二頻寬矩陣,以分別將該等第一頻寬矩陣及該等第二頻寬矩陣中同一主能量頻率的能量累加。 The spectrum analysis method of claim 1, wherein in the fourth step, the first bandwidth matrix is integrated according to the first integrated bandwidth matrix and the second integrated bandwidth matrix, respectively. And a second bandwidth matrix to accumulate the energy of the same primary energy frequency in the first bandwidth matrix and the second bandwidth matrix, respectively. 如申請專利範圍第1項所述之頻譜分析方法,其中於該第五步驟中,係將該第一能量頻譜中,與該第二能量頻譜具有相同頻率之能量頻譜去除,以建立該能量訊號的主能量頻率之該主能量頻譜。 The spectrum analysis method of claim 1, wherein in the fifth step, the energy spectrum of the first energy spectrum having the same frequency as the second energy spectrum is removed to establish the energy signal. The main energy spectrum of the main energy frequency. 如申請專利範圍第1項所述之頻譜分析方法,其中於該第五步驟中,當建立該能量訊號的主能量頻率之該 主能量頻譜時,需以一固定資料庫對該主能量頻譜進行修正。 The spectrum analysis method according to claim 1, wherein in the fifth step, when the main energy frequency of the energy signal is established In the main energy spectrum, the main energy spectrum needs to be corrected with a fixed database. 如申請專利範圍第6項所述之頻譜分析方法,其中當建立該能量訊號的主能量頻率之該主能量頻譜時,若該第一能量頻譜上的某個主能量頻率在該第二能量頻譜上找不到對應的頻率時,則將該固定資料庫之該主能量頻率的頻寬導入,以取代並修正該主能量頻率之該主能量頻譜的頻寬。 The spectrum analysis method according to claim 6, wherein when the main energy spectrum of the main energy frequency of the energy signal is established, if a main energy frequency in the first energy spectrum is in the second energy spectrum When the corresponding frequency is not found, the bandwidth of the main energy frequency of the fixed database is introduced to replace and correct the bandwidth of the main energy spectrum of the main energy frequency. 如申請專利範圍第1項所述之頻譜分析方法,更包括:一第六步驟:依據該第一整合頻寬矩陣及該第二整合頻寬矩陣,分別建立一第一累加頻寬頻譜及一第二累加頻寬頻譜。 The method for spectrum analysis according to claim 1, further comprising: a sixth step: respectively establishing a first accumulated bandwidth spectrum and a first integrated bandwidth matrix and the second integrated bandwidth matrix The second accumulated bandwidth spectrum. 如申請專利範圍第8項所述之頻譜分析方法,其中於該第六步驟中,係分別將各該第一頻譜訊號及各該第二頻譜訊號之各主能量頻率的左頻寬平方加右頻寬平方,以形成一新頻寬指標,並分別將該等第一頻譜訊號及該等第二頻譜訊號之頻譜中,各主能量頻率的新頻寬指標相加,以分別建立該等第一頻譜訊號及該等第二頻譜訊號的該第一累加頻寬頻譜及該第二累加頻寬頻譜,其中,該第一累加頻寬頻譜及該第二累加頻寬頻譜分別具有該新頻寬指標相加的頻寬值。 The method of spectrum analysis according to claim 8, wherein in the sixth step, the square of the left bandwidth of each of the first spectrum signal and each of the second spectrum signals is respectively added to the right The bandwidth is squared to form a new bandwidth index, and the new bandwidth indices of the primary energy signals and the spectrums of the second spectral signals are respectively added to each other to establish the first a first accumulated bandwidth spectrum of the spectral signal and the second spectral signals and the second accumulated bandwidth spectrum, wherein the first accumulated bandwidth spectrum and the second accumulated bandwidth spectrum respectively have the new bandwidth The bandwidth value to which the metrics are added. 如申請專利範圍第8項所述之頻譜分析方法,更包括:一第七步驟:找出該第一累加頻寬頻譜之該新頻寬指標相加的頻寬值大於或等於2倍的該第二累加頻寬 頻譜之主能量頻率之該新頻寬指標相加的頻寬值,以建立一主能量頻寬比例頻譜。 The method for spectrum analysis according to claim 8 of the patent application, further comprising: a seventh step: finding that the added bandwidth value of the new bandwidth indicator of the first accumulated bandwidth spectrum is greater than or equal to 2 times Second accumulated bandwidth The bandwidth value of the new bandwidth index of the main energy frequency of the spectrum is added to establish a main energy bandwidth proportional spectrum. 一種疾病篩檢方法,係用以篩檢一疾病,該疾病篩檢方法包括如申請專利範圍第1項至第10項任一項所述之頻譜分析方法,該疾病篩檢方法包括:對具有該疾病的複數病患及複數正常人分別進行檢測,以分別得到該等病患及該等正常人之該等能量訊號;分別以該頻譜分析方法分析該等病患及該等正常人之該等能量訊號,以分別得到該等病患及該等正常人之該等主能量頻譜及該等主能量頻寬比例頻譜;分別對應比較各該病患與各該正常人之各該主能量頻譜及各該主能量頻寬比例頻譜,以分別找出該等主能量頻譜及該等主能量頻寬比例頻譜中,只存在該等病患之該疾病的複數頻率;以及依據該等主能量頻譜及該等主能量頻寬比例頻譜中只存在該等病患之該疾病的該等頻率,分別計算該疾病之各該頻率的重覆比例。 A disease screening method for screening a disease, the disease screening method comprising the method of spectrum analysis according to any one of claims 1 to 10, wherein the disease screening method comprises: The multiple patients and the normal persons of the disease are separately tested to obtain the energy signals of the patients and the normal persons respectively; the spectrum analysis method is used to analyze the patients and the normal persons respectively. Equal energy signals to obtain the main energy spectrum of the patients and the normal persons and the spectrum of the main energy bandwidth ratios; respectively, respectively comparing the main energy spectrum of each of the patients and each of the normal persons And each of the main energy bandwidth ratio spectrums to separately find the main energy spectrum and the main energy bandwidth ratio spectrum, only the complex frequencies of the diseases of the patients are present; and according to the main energy spectrum And the frequencies of the diseases in the main energy bandwidth ratio spectrum are only present in the patients, and the repetition ratios of the frequencies of the diseases are respectively calculated. 如申請專利範圍第11項所述之疾病篩檢方法,其中於該檢測步驟中,係檢測一天或複數天,且每天檢測至少一次。 The disease screening method according to claim 11, wherein in the detecting step, detecting one day or multiple days, and detecting at least once a day. 如申請專利範圍第11項所述之疾病篩檢方法,更包括:當一被檢測者進行檢測,並以該頻譜分析方法分析 後,該被檢測者具有該疾病之該等頻率的至少其中之一,且該頻率的重覆比例越高者,則該被檢測者得到該疾病之機率越高。 The method for screening diseases according to claim 11 of the patent application further includes: when a test subject performs the test, and analyzes by the spectrum analysis method Thereafter, the subject has at least one of the frequencies of the disease, and the higher the repetition ratio of the frequency, the higher the probability that the subject gets the disease. 如申請專利範圍第11項所述之疾病篩檢方法,更包括:當一被檢測者進行檢測,並以該頻譜分析方法分析後,該被檢測者具有該疾病之該等頻率的至少其中之一,且該頻率的重覆比例越低者,則該被檢測者得到該疾病之機率越低。 The method for screening a disease according to claim 11, further comprising: when a subject is tested and analyzed by the spectrum analysis method, the subject has at least one of the frequencies of the disease. First, the lower the repetition ratio of the frequency, the lower the probability that the subject will get the disease. 如申請專利範圍第11項所述之疾病篩檢方法,更包括:當一被檢測者進行檢測,並以該頻譜分析方法分析後,該被檢測者沒有該疾病之該等頻率的其中之一,則該被檢測者得到該疾病之機率越低。The method for screening a disease according to claim 11 further includes: when a test subject performs the test and analyzes by the spectrum analysis method, the test subject does not have one of the frequencies of the disease. The lower the chance that the subject will get the disease.
TW100147346A 2011-12-20 2011-12-20 Spectrum analysis method and disease examination method TWI467184B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
TW100147346A TWI467184B (en) 2011-12-20 2011-12-20 Spectrum analysis method and disease examination method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
TW100147346A TWI467184B (en) 2011-12-20 2011-12-20 Spectrum analysis method and disease examination method

Publications (2)

Publication Number Publication Date
TW201326834A TW201326834A (en) 2013-07-01
TWI467184B true TWI467184B (en) 2015-01-01

Family

ID=49224968

Family Applications (1)

Application Number Title Priority Date Filing Date
TW100147346A TWI467184B (en) 2011-12-20 2011-12-20 Spectrum analysis method and disease examination method

Country Status (1)

Country Link
TW (1) TWI467184B (en)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1151684A (en) * 1994-05-26 1997-06-11 卡罗琳娜丝心脏研究所 Microwave tomographic spectroscopy system and method
US20020007122A1 (en) * 1999-12-15 2002-01-17 Howard Kaufman Methods of diagnosing disease
US20100228145A1 (en) * 2009-03-06 2010-09-09 Bo-Jau Kuo Device and method for evaluating sympathetic function using electrooculography
TW201132336A (en) * 2010-03-17 2011-10-01 Univ Nat Chiao Tung Electrode simulating method and system thereof

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1151684A (en) * 1994-05-26 1997-06-11 卡罗琳娜丝心脏研究所 Microwave tomographic spectroscopy system and method
US20020007122A1 (en) * 1999-12-15 2002-01-17 Howard Kaufman Methods of diagnosing disease
US20100228145A1 (en) * 2009-03-06 2010-09-09 Bo-Jau Kuo Device and method for evaluating sympathetic function using electrooculography
TW201132336A (en) * 2010-03-17 2011-10-01 Univ Nat Chiao Tung Electrode simulating method and system thereof

Also Published As

Publication number Publication date
TW201326834A (en) 2013-07-01

Similar Documents

Publication Publication Date Title
Sharanreddy et al. EEG signal classification for epilepsy seizure detection using improved approximate entropy
Neto et al. Regularized linear discriminant analysis of EEG features in dementia patients
França et al. Fractal and multifractal properties of electrographic recordings of human brain activity: toward its use as a signal feature for machine learning in clinical applications
Robinson et al. Age effect in hearing-a comparative analysis of published threshold data
Rezaeezadeh et al. Attention Deficit Hyperactivity Disorder Diagnosis using non-linear univariate and multivariate EEG measurements: a preliminary study
Hartings et al. Direct current electrocorticography for clinical neuromonitoring of spreading depolarizations
Bladh et al. Psychometric performance of a generic walking scale (Walk-12G) in multiple sclerosis and Parkinson’s disease
Cleworth et al. Differences in head impulse test results due to analysis techniques
Vianna et al. How to detect progression in glaucoma
Niu et al. Abnormal dynamic functional connectivity and brain states in Alzheimer’s diseases: functional near-infrared spectroscopy study
KR20090024808A (en) Assessing dementia and dementia-type disorders
US20110230781A1 (en) Apparatus and method for dementia diagnosis through eeg (electroencephalogram) analysis
CN111739647B (en) Magnetic resonance image-based bipolar disorder suicide risk prediction method and related device
Parks et al. Bootstrap signal-to-noise confidence intervals: an objective method for subject exclusion and quality control in ERP studies
Luke et al. Characterization of Mayer-wave oscillations in functional near-infrared spectroscopy using a physiologically informed model of the neural power spectra
Grossman et al. The noisy brain: power of resting-state fluctuations predicts individual recognition performance
Marín-Franch et al. Choice of statistical method influences apparent association between structure and function in glaucoma
Sadus et al. An explorative multiverse study for extracting differences in P3 latencies between young and old adults
CN116509417B (en) Consumer preference consistency prediction method based on nerve similarity
US9618597B2 (en) Method and magnetic resonance apparatus for automated analysis of the raw data of a spectrum
TWI467184B (en) Spectrum analysis method and disease examination method
Baboukani et al. Classifying single-trial eeg during motor imagery using a multivariate mutual information based phase synchrony measure
Davoudi et al. Inter-individual variability during neurodevelopment: an investigation of linear and nonlinear resting-state EEG features in an age-homogenous group of infants
Reichenheim et al. Re-assessing the dimensional structure of the Informant Questionnaire on Cognitive Decline in the Elderly (IQCODE): empirical evidence for a shortened Brazilian version
Bhattacharyya et al. Integration of electroencephalogram (EEG) and motion tracking sensors for objective measure of attention-deficit hyperactivity disorder (MAHD) in pre-schoolers

Legal Events

Date Code Title Description
MM4A Annulment or lapse of patent due to non-payment of fees