信號與系統(tǒng)在生物醫(yī)學(xué)中應(yīng)用

信號與系統(tǒng)在生物醫(yī)學(xué)中的應(yīng)用

利用信號與系統(tǒng)研究生物神經(jīng)系統(tǒng)的方法及意義

生物神經(jīng)網(wǎng)絡(luò)模型及等效電路

神經(jīng)網(wǎng)絡(luò)的數(shù)學(xué)建模

數(shù)值計算方法1信號與系統(tǒng)在生物醫(yī)學(xué)中應(yīng)用

利用信號與系統(tǒng)研究生物神經(jīng)系統(tǒng)的方法

BNN數(shù)據(jù)采集BNN數(shù)據(jù)分析BNN等效電路BNN數(shù)學(xué)模型BNN數(shù)值計算BNN信號處理實驗數(shù)據(jù)與仿真結(jié)果相比較①

②③⑦⑥

BNN仿真輸出④⑤

根據(jù)實驗數(shù)據(jù)和概念模型建立相應(yīng)計算模型，從而逐步理解和推斷生物神經(jīng)系統(tǒng)的生理結(jié)構(gòu)和生化組成，以及內(nèi)部作用機理。2信號與系統(tǒng)在生物醫(yī)學(xué)中應(yīng)用

研究生物神經(jīng)系統(tǒng)的意義

驗證生物神經(jīng)網(wǎng)絡(luò)的模型和假設(shè)

系統(tǒng)化和理論化生物神經(jīng)系統(tǒng)的實驗數(shù)據(jù)

拓展實驗的范圍并延伸實驗的功能

探索生物神經(jīng)系統(tǒng)的未知領(lǐng)域

為其它學(xué)科模擬生物神經(jīng)網(wǎng)絡(luò)奠定必要基礎(chǔ)

為其它生物系統(tǒng)的仿真奠定了良好的基礎(chǔ)

3信號與系統(tǒng)在生物醫(yī)學(xué)中應(yīng)用

生物神經(jīng)系統(tǒng)的模型及等效電路神經(jīng)元生理結(jié)構(gòu)和生化組成靜息狀態(tài)下單個神經(jīng)元等效電路激勵狀態(tài)下單個神經(jīng)元等效電路神經(jīng)網(wǎng)絡(luò)中神經(jīng)元等效電路

4信號與系統(tǒng)在生物醫(yī)學(xué)中應(yīng)用神經(jīng)元生理結(jié)構(gòu)

非線性時變系統(tǒng)——神經(jīng)元5信號與系統(tǒng)在生物醫(yī)學(xué)中應(yīng)用1.神經(jīng)元生理結(jié)構(gòu)及生化組成

神經(jīng)元生理結(jié)構(gòu)一般由三個主要部分組成，即細胞體(soma)，軸梢(dendrite)，軸突(axon)。細胞體位于神經(jīng)元的中心部分，它包含細胞核，軸梢是從細胞核發(fā)射出的許多根狀物，軸突也是從細胞核發(fā)射出的一根管狀纖維。其中軸梢主要功能是從其它神經(jīng)元接受電信號；細胞體主要功能是積累來自許多軸梢的電位；軸突的主要功能是傳導(dǎo)電信號，并傳遞電信號至其它神經(jīng)元。

沿著神經(jīng)元軸突膜(membrane)分布的膜電位是描述神經(jīng)元內(nèi)信息傳遞的重要物理量。

6信號與系統(tǒng)在生物醫(yī)學(xué)中應(yīng)用1.神經(jīng)元生理結(jié)構(gòu)及生化組成

神經(jīng)膜內(nèi)外的離子濃度分布離子(Ion)

膜內(nèi)濃度(InsideCell)

膜外濃度(Outsidecell)

K+

397mM/l

20mM/l

Na+

49mM/l

440nM/l

Cl-

48mM/l

480mM/l

7信號與系統(tǒng)在生物醫(yī)學(xué)中應(yīng)用1.神經(jīng)元生理結(jié)構(gòu)及生化組成

這些離子主要為三種單元素離子，鉀離子(K+)、鈉離子(Na+)、和氯離子(Cl-)以及某些復(fù)合離子。其中，正極性的鉀離子(K+)主要分布在神經(jīng)膜內(nèi)，而正極性的鈉離子(Na+)和負極性的氯離子(Cl-)主要分布在神經(jīng)膜外。正是由于神經(jīng)膜內(nèi)外的這些離子的存在以及它們在膜內(nèi)外的濃度分布不同，形成了膜電位。膜電位的外在特性可分為明顯的兩個階段，即靜息膜電位(restingmembranepotential)階段和動作電位(actionpotential)階段。靜息膜電位為負極性，一般在-60與-70mV之間，神經(jīng)細胞在大多數(shù)情形下一直處于此平衡狀態(tài)。

8信號與系統(tǒng)在生物醫(yī)學(xué)中應(yīng)用2.靜息狀態(tài)下單個神經(jīng)元等效電路其中：gK,EK分別為鉀離子K+的等效電導(dǎo)和靜息電位；gNa,Ena分別為鈉離子Na+的等效電導(dǎo)和靜息電位；gCl,Ecl分別為氯離子Cl-的等效電導(dǎo)和靜息電位；

Vm為神經(jīng)細胞靜息膜電位。gKgNa

gcl

EKENa

ECl

IKINa

ICl

Vm

靜息膜電位為Vm=-60.2mv，鉀離子電流為IK

=5.1mA/cm2,鈉離子電流為gNa=-4.68mA/cm2,氯離子電流為gCl

=-0.25mA/cm2。9信號與系統(tǒng)在生物醫(yī)學(xué)中應(yīng)用3.激勵狀態(tài)下單個神經(jīng)元等效電路其中：gK,EK分別為鉀離子K+的等效電導(dǎo)和靜息電位；gNa,Ena分別為鈉離子Na+的等效電導(dǎo)和靜息電位；gCl,Ecl分別為氯離子Cl-的等效電導(dǎo)和靜息電位；

Vm為神經(jīng)細胞靜息膜電位。

CM為細胞膜電容；

IS為外部觸發(fā)gKgNa

gcl

EKENa

ECl

IKINa

ICl

Vm

CMMIC

IS

10信號與系統(tǒng)在生物醫(yī)學(xué)中應(yīng)用4.神經(jīng)系統(tǒng)中神經(jīng)元等效電路神經(jīng)元等效電路

Gion1Gion2GionmEion1Eion2EionmCM

IexGes1Ges2Gesn

V1V2Vn

Gcs1,1Gcs1,2Gcs1,pEcs1,1Ecs1,2Ecs1,pGcsn,1Gcsn,2Gcsn,pEcsn,1Ecsn,2Ecsn,pIonicconductances

Electricalsynapses(es)

Chemicalsynapses(cs)

+++++++++++++11信號與系統(tǒng)在生物醫(yī)學(xué)中應(yīng)用4.神經(jīng)系統(tǒng)中神經(jīng)元等效電路神經(jīng)元電化學(xué)模型

12信號與系統(tǒng)在生物醫(yī)學(xué)中應(yīng)用

神經(jīng)網(wǎng)絡(luò)的數(shù)學(xué)建模數(shù)學(xué)建模的基本原理神經(jīng)膜動作電位模型離子電導(dǎo)模型化學(xué)突觸電導(dǎo)模型電突觸電導(dǎo)模型離子電導(dǎo)調(diào)節(jié)因子函數(shù)模型電導(dǎo)隨機特性模型13信號與系統(tǒng)在生物醫(yī)學(xué)中應(yīng)用1.數(shù)學(xué)建模的基本原理

首先對在一定控制條件下得到的生物數(shù)據(jù)進行定性地分析，確定該實驗數(shù)據(jù)可能的概念模型。然后依照該概念模型確定相關(guān)的數(shù)學(xué)模型，按照最小均方誤差準則從這些數(shù)據(jù)中得到對應(yīng)的模型參數(shù)。為歸一化整個數(shù)學(xué)建模過程，通過簡單變換的方法，將其它模型都映射為直線模型，將數(shù)據(jù)也按相同的變換映射為新的數(shù)據(jù)集合，在得到直線模型參數(shù)后，再反變換成原來模型。14信號與系統(tǒng)在生物醫(yī)學(xué)中應(yīng)用1.數(shù)學(xué)建模的基本原理

其中：分別為直線參數(shù)a,b的估計； 為實驗數(shù)據(jù)序列對；

N

為數(shù)據(jù)序列的長度。15信號與系統(tǒng)在生物醫(yī)學(xué)中應(yīng)用2.神經(jīng)膜動作電位模型

其中：i為網(wǎng)絡(luò)中神經(jīng)元膜電位V的下標；

j為離子電流Iion的下標；

k為電突觸電流Ies的下標；

l為化學(xué)突觸電流Ics的下標；

r為外部激勵電流Iex的下標；

m

為各神經(jīng)元中電壓依賴電導(dǎo)的數(shù)目；

p

為每個電突觸伴隨的化學(xué)突觸的數(shù)目；16信號與系統(tǒng)在生物醫(yī)學(xué)中應(yīng)用3.離子電導(dǎo)模型

離子電流：17信號與系統(tǒng)在生物醫(yī)學(xué)中應(yīng)用3.離子電導(dǎo)模型

18信號與系統(tǒng)在生物醫(yī)學(xué)中應(yīng)用4.化學(xué)突觸電導(dǎo)模型

化學(xué)突觸電流：其中：f(At)是時間依賴性函數(shù)f(Av,t)是時間和電壓依賴性函數(shù)19信號與系統(tǒng)在生物醫(yī)學(xué)中應(yīng)用5.電突觸電導(dǎo)模型

電突觸電流：其中：是最大的電導(dǎo)是電突觸電導(dǎo)的隨機波動是突觸前神經(jīng)元(presynapticneuron)的膜電位是突觸后神經(jīng)元(postsynapticneuron)的膜電位

20信號與系統(tǒng)在生物醫(yī)學(xué)中應(yīng)用6.離子電導(dǎo)調(diào)節(jié)因子函數(shù)模型

其中：b為常數(shù)

gbr

反映離子濃度或第二信使?jié)舛?/p>

對膜電導(dǎo)G增強或衰減的程度離子池濃度第二信使?jié)舛?1信號與系統(tǒng)在生物醫(yī)學(xué)中應(yīng)用7.電導(dǎo)隨機特性模型

正態(tài)分布函數(shù)根據(jù)提供的四個參數(shù)來產(chǎn)生對應(yīng)的隨機序列，即均值，標準方差，種子值(seed)和更新速率(refreshrate)。相同的均值和標準方差，由于種子值的不同而產(chǎn)生不同的隨機序列。22信號與系統(tǒng)在生物醫(yī)學(xué)中應(yīng)用

數(shù)值計算方法問題的提出等間隔步長數(shù)值計算方法自適應(yīng)步長數(shù)值計算方法混合數(shù)值計算方法(HNCM)HNCM算法的應(yīng)用各數(shù)值計算方法比較

23信號與系統(tǒng)在生物醫(yī)學(xué)中應(yīng)用1.問題的提出

生物神經(jīng)網(wǎng)絡(luò)是一個異常復(fù)雜的非線性系統(tǒng)，難以通過數(shù)學(xué)的方法求解其非線性時變微分方程組，從而得到解析表達式，只能采用數(shù)值計算的方法得到微分方程的數(shù)值解。在數(shù)值計算過程中，為保證所需的計算精度，計算步長必須足夠小。在某特定的數(shù)值計算方法中，較小的計算步長可以得到較高的計算精度，但必然導(dǎo)致較低的計算效率。仿真系統(tǒng)的計算精度和計算效率是仿真系統(tǒng)的核心指標。為實現(xiàn)所需精度下的最佳效率，需要設(shè)計合適的計算方法以提高仿真系統(tǒng)的計算精度和效率。24信號與系統(tǒng)在生物醫(yī)學(xué)中應(yīng)用2.等間隔步長數(shù)值計算方法

①②

③

f(t)t1

t2

t3

t

…

Error=O(

2)Euler數(shù)值計算方法

25信號與系統(tǒng)在生物醫(yī)學(xué)中應(yīng)用2.等間隔步長數(shù)值計算方法

Midpoint數(shù)值計算方法

①④

⑤

f(t)t1

t2

t3

t

…

②

③

Error=O(

3)26信號與系統(tǒng)在生物醫(yī)學(xué)中應(yīng)用2.等間隔步長數(shù)值計算方法

Four-orderR-K數(shù)值計算方法Error=O(

5)

計算精度增高，計算復(fù)雜度不斷增大。27信號與系統(tǒng)在生物醫(yī)學(xué)中應(yīng)用3.自適應(yīng)步長數(shù)值計算方法

v1(t)

v2(t)

v3(t)

t

t

t

t2

t1

ti

t3

tn

…

…

…

R-K自適應(yīng)步長數(shù)值計算方法

各變量的計算時序不同由于變量之間存在依賴性，自適應(yīng)步長方法無法適用。28信號與系統(tǒng)在生物醫(yī)學(xué)中應(yīng)用4.混合數(shù)值計算方法

HybridNumericalComputationMethod(HNCM)

傳統(tǒng)數(shù)值計算方法的不足：

等間隔步長方法:

對所有的信號應(yīng)用相同的計算步長，這將導(dǎo)致較低的運算精度或較高的計算復(fù)雜度。因為若考慮到變化較快的信號，需要設(shè)置較小的步長，這對變化較慢的信號就會產(chǎn)生計算冗余；反之，若考慮到變化較慢的信號，需要設(shè)置較大的步長，這對變化較快的信號量就會產(chǎn)生計算誤差；

變步長方法:

由于仿真系統(tǒng)中的各信號之間存在依賴性，這種計算方法無法適用。29信號與系統(tǒng)在生物醫(yī)學(xué)中應(yīng)用4.