【正文】
ith the change of loading signal frequency. It is a frequency domain measurement method to measure a relatively wide frequency domain impedance spectrum in researching the physiological feature of biological tissues. Bioelectrical impedance measurement has the requirement that of tiny wound or wound absence, rapid measuring and high accuracy. This thesis employs a fourelectrode method adopting current source excitation to measure bioelectrical impedance. It designs a highspeed acquisition system and afterwards processes data to acquire impedance value with a set of algorithm. The main research content is as follows: 1. The thesis firstly research and analyze properties of m sequence. Due to autocorrelation function of m sequence possesses the merits, such as close to impulse function, discrete power spectrum, strong antiinterference ability, convenient bandwidth adjustment and binary function is easy to realize, and the m sequence pseudo random signal is chosen as the current source excitation model and is implemented on FPGA platform. When using m sequence as bioelectricity impedance model excitation source, it can obtain system impact response through solving correlation function. Therefore, impedance spectrum is convenient to be acquired. 2. It designs the highspeed acquisition system that based on FPGA and ARM in order to achieve rapid measurement. The system use with FPGA logic programmable chip and STM32 microprocessor as a core, it designs several modules, such as power supply module, excitation signal source module, analogtodigital converter module, data buffer module, FPGA front control module and STM32 backend control module. It realizes synchronous sampling about m sequence excitation and response signal, and achieves correct sampling about experimental data. 西安理工大學(xué)碩士學(xué)位論文 IV 3. In the thesis, a FFT detection algorithm for calculating correlation function fast is designed for obtaining the correlation between the incentive current and response voltage signal which is the measured impedance of the time domain impulse response by means of the method which is extending sequence through adding zero, and it is based on rapid correlation algorithm and the AllPhase FFT calculating impedance value algorithm. The result obtained is the time domain impulse response of the measured impedance. In order to get the impedance spectrum, the measured impedance frequency response is get with adopting the AllPhase spectrum analysis method. So multifrequency measuring the electrical impedance model simultaneously and rapidly is realized 4. The thesis employs impedance to do measurement experiment for measuring system and studying algorithm, it demarcate measurement system and analysis the result of the measurement error. This thesis realizes simultaneous measurement of bioelectrical impedance many frequency, and provides a feasible method for bioelectricity impedance rapid measurement research. Key words: m sequence。它是一種 在頻域 測(cè)量 生物電阻抗的 方法,以測(cè)得頻域較寬的阻抗譜來(lái)研究生物組織的生理、病理等狀態(tài)。其一般的實(shí)現(xiàn)方法是在生物 體 的 表面加上激勵(lì)電極,然后對(duì)激勵(lì)電極施加微小的交變電壓信號(hào)或交變電流信號(hào),通過(guò)測(cè)量電極檢測(cè)對(duì)應(yīng)的電流信號(hào)或電壓信號(hào),進(jìn)而計(jì)算出相應(yīng)的電阻抗及其變化,最后根據(jù)不同的應(yīng)用目的獲取 所需要的相關(guān) 信息 【 2】 。如果能夠把掃頻時(shí)間減少,則有可能避免被測(cè) 生物組織 出現(xiàn)極化效應(yīng),在較短時(shí)間內(nèi)所測(cè)得的各個(gè)頻率點(diǎn)的阻抗信息更能準(zhǔn)確反映被測(cè) 生物組織 的真實(shí)阻抗信息。據(jù)相關(guān)資料表明, BIS 在 組織缺血 【 5】 、 生物 水合 作用估計(jì) 【 6】 、 肺水腫 【 7】 、 臨床診斷中 干重 的估計(jì) 【 8】 、 惡性皮膚癌的診斷 【 9】 、 病人 在 營(yíng)養(yǎng)不良方面 的 鑒別 【 10】 、腫瘤 的早期診斷和 監(jiān)測(cè) 【 1113】 、血液透析 【 1415】 、乳腺癌 【 1617】 、 病人手術(shù)后的 細(xì)胞變化監(jiān)測(cè) 【 18】 、 淋巴水腫探測(cè) 【 1920】 以及 輻射損傷檢測(cè) 【 21】 等 方面 有 一些成果。 西安理工大學(xué)碩士學(xué)位論文 2 國(guó)內(nèi)外研究現(xiàn)狀 縱觀國(guó)內(nèi)外 BIS 測(cè)量技術(shù)一百多年的發(fā)展歷程, 生物電阻抗研究的 激勵(lì)信號(hào)源的選擇 經(jīng)歷了從直流到交流 、被測(cè)對(duì)象模型經(jīng)歷了 從純電阻到復(fù)阻 抗 、測(cè)量原理經(jīng)歷了從單頻測(cè)量到多頻測(cè)量的發(fā)展過(guò)程,多頻測(cè)量包含了 分時(shí)單頻 率 掃描 測(cè)量 模式、多頻率混頻 測(cè)量模式、多 頻率 同步 測(cè)量 模式。 但 在 單一 的 頻率 成分 下, 獲得的生物信息不多,并不能滿足測(cè)量的 要求 【 22】 。但是在測(cè)量過(guò)程中 測(cè)量頻率 要 改變 多次,用時(shí)比較長(zhǎng)。 近 些 年 ,在 國(guó)內(nèi) 有 部分 學(xué)者 嘗試了 多路信號(hào)混頻 方式 。 混頻 模式主要是很難實(shí)現(xiàn)各頻率分量相位的精確控制, 很難同時(shí)將更多路 信號(hào)的 正弦波進(jìn)行混頻。國(guó)外學(xué)者 Min等 【 28】 對(duì) sinc 脈沖 、 高斯脈沖 等寬頻帶 的 脈沖做了一個(gè) 比較全面的 回顧,發(fā)現(xiàn) 與 高斯脈沖 、矩形脈沖 相比較, 啁啾脈沖 ,也叫 線性調(diào)頻脈沖 , 在 BIS 的 測(cè)量 中有更 好的性能,但其復(fù)雜的數(shù)學(xué)表達(dá)式使得其波形實(shí)現(xiàn)需要復(fù)雜的硬件設(shè)備(如 高端的 任意函數(shù)發(fā)生器),數(shù)據(jù)采集和處理系統(tǒng)必須有高速性能(在 us 內(nèi)獲取激勵(lì)脈沖及其對(duì)應(yīng)響應(yīng)),實(shí)際實(shí)驗(yàn)需要很高端 的 儀器。 測(cè)量平臺(tái)的構(gòu)建 有了激勵(lì)源模型 ,還要有 一個(gè)能對(duì)激勵(lì)源和生物電阻抗響應(yīng)信號(hào)同步采 樣 的采集系統(tǒng),要求 速度 快,并能夠?qū)?shù)據(jù)進(jìn)行正確的存儲(chǔ)與讀取 。 設(shè)計(jì)基于 FPGA 的數(shù)據(jù)采集系統(tǒng) 為了保證 能得到 可靠的數(shù)據(jù) , 要求 建立一套準(zhǔn)確的采集系統(tǒng) ,這步 非常關(guān)鍵,構(gòu)建基于 FPGA 和模數(shù)轉(zhuǎn)換器的采集部分,能對(duì)激勵(lì)信號(hào)和響應(yīng)信號(hào)實(shí)現(xiàn) 2 路同步采 樣 ,為后西安理工大學(xué)碩士學(xué)位論文 4 續(xù)的相關(guān)檢測(cè)算法提供數(shù)據(jù)。 本章小結(jié) 本文首先介紹了生物電阻抗的背景和研究意義, 生物電阻抗的研究對(duì)醫(yī)學(xué)方面有很大的貢獻(xiàn)。 m序列的研究和激勵(lì)信號(hào)源的設(shè)計(jì)思路 5 2 m 序列的研究和激勵(lì)信號(hào)源的設(shè)計(jì)思路 生物電阻抗的測(cè)量和研究中,最基本的是選擇一個(gè)合適的信號(hào)源模型,所以要研究一種信號(hào)源模型,本文選擇的是偽隨機(jī)信號(hào) m 序列,偽隨機(jī)信號(hào)有自己的利用價(jià)值, m 序列作為其中一種,它是一種偽隨機(jī)二電平信號(hào)。 就偽隨機(jī)信號(hào)和隨機(jī)信號(hào)相比較,有兩個(gè)比較明顯的特點(diǎn) 【 30】 : ( 1) 隨 機(jī)信號(hào)是沒(méi)有周期性的,而偽隨機(jī)信號(hào)在經(jīng)過(guò)一定的時(shí)間后是會(huì)重現(xiàn)的。 所以,足夠長(zhǎng)的偽隨機(jī)信號(hào)可以作為隨機(jī)信號(hào)用。 m 序列 m 序列的定義 及介 紹 m 序列是最長(zhǎng)線性反饋移位寄存器的簡(jiǎn)稱,它是由線性反饋的移位寄存器產(chǎn)生的周期最長(zhǎng)的一種序列。 隨著 m序列階數(shù)的增大, m序列的周期也增大。一般而言,在 m 序列中,游程長(zhǎng)度為 1 的數(shù)目占總游程數(shù)目的 1/2,游程長(zhǎng)度為 2 的數(shù)目占總游程數(shù)目的 1/4,游程長(zhǎng)度為 3 的數(shù)目占總游程的 1/8, 依次 類推,游程長(zhǎng)度為 k 的數(shù)目占總游程的 2k。任何一個(gè) m序列 信號(hào)與 對(duì)應(yīng)的 橫坐標(biāo)都有 [(N+1)/2]個(gè)交點(diǎn),因此可以得到上述 由于 錯(cuò)位 而 引起 減小的平均正面積 S+為 21()2NaS Nt?? ?? () 增加的平均負(fù)面積 S為 2_ 1( )( )2NaS Nt???? () 即可以得到 2221 1 1( 0) ( ) ( ) ( ) [ 1 ( ) ]220x x x xN a N a NR R aN t N t N tt?? ??? ? ?? ? ? ? ??? () 如 圖 21 所示。相應(yīng)的雙極性波形 可表示 為 ( ) ( )nc t a t nT??????? () 其中 , ( ) ( 1,1)at?? () 取 m序列的一個(gè)周期 , c(t)的歸一化自相關(guān)函數(shù)為 01( ) ( ) ( )Tr c c t d tT? ? ?? ? ?? () 如果令 1 1()0cT cN Tr NT? ???? ????? ??? ????? ??? ????, 其 他 () 則 ? ? ? ? ? ?111Tnr r n T rNN? ? ??? ? ?? ? ? ? ?? () 其中 , 1 ( ) ( )Tnr r n T????????? () 由于 c(t)的功率譜密度 G(ω)和 r(τ)互為傅里葉變換,所以有 西安理工大學(xué)碩士學(xué)位論文 10 ? ? ? ?1 1( ) ( )G F r F r F N? ? ? ??? ? ??? ???? ?? () 周期性函數(shù) r1(τ)可以展為傅里葉級(jí)數(shù) ? ?010 2,j n tnnr F e T? ????? ??? ? ?? () 其中, ? ??? ??? 22 101TT tjnn dterTF ?? ? ? ? ?0011j n tTT nr e d t F rTT? ????? ? ???? ? ? ????? () 可以計(jì)算 ? ?? ? ? ? ? ?? ? ? ?? ????? ???? ???? n n Tn nrFTnFrF 001 22 ?????????? ? ? ? ?20