信号处理方法分析脑血流自动调节功能文档资料

1、内内 容容 研究意义研究意义1 相关知识介绍相关知识介绍2方法介绍方法介绍3工作内容工作内容4研究意义研究意义大脑大脑： 脑组织几乎没有能量储备，大脑依靠连续不断血液循环来提供氧和脑组织几乎没有能量储备，大脑依靠连续不断血液循环来提供氧和葡萄糖等营养物质来维持正常的功能和结构的完整性。葡萄糖等营养物质来维持正常的功能和结构的完整性。 大脑的重量约占人体体重的大脑的重量约占人体体重的2%，但是却要消耗人体，但是却要消耗人体25%的氧和的氧和1520%的供血（安静平躺的情况下）。的供血（安静平躺的情况下）。 大脑对缺血十分敏感，当大脑的供血不足或中断时，由于没有氧、大脑对缺血十分敏感，当大脑的供血

2、不足或中断时，由于没有氧、葡萄糖和其他营养物质，在葡萄糖和其他营养物质，在35分钟后脑细胞就开始软化甚至死亡。为分钟后脑细胞就开始软化甚至死亡。为保护脑组织不受缺血的影响，大脑就需要一种能够保证相对稳定的供血保护脑组织不受缺血的影响，大脑就需要一种能够保证相对稳定的供血机能，这就是脑血流调节机能。机能，这就是脑血流调节机能。具有正常的自动调节能力对于脑机体是至关重要的。但是，当发生脑卒中、高血压、脑外伤、脑肿瘤等严重的脑损伤时会使脑血流自动调节功能下降。目前在临床上虽然可以控制动脉血压，但由于不了解每一个体的脑血流自动调节状态，无法实时的调整脑灌注压，也就无法保证脑血流量的恒定。因此，如果能及

3、时、准确的评价患者的脑血流自动调节能力，对脑血管疾病的诊断及治疗有重要的指导意义。研究意义研究意义研究意义研究意义大脑自动调节功能可以简单的描述成在当灌注压（cpp）变化时大脑维持脑血流在恒定的范围的机制 。cbf = cpp / cvrcbf = cpp/(8*l)/r4相关知识介绍相关知识介绍1、肌源性学说2、代谢学说 3、神经源性学说 4、内皮细胞源学说 脑血流自动调节功能的生理机制还没有完全的被解释.目前主要有4种学说以解释脑血流自动调节功能的生理机制。数据采集连续血压测量仪 多谱勒超声测量脑血流速度 (tcd)分析dca的数据1、 外界模拟abp改变 束带法通过使受试者周期的呼吸引起

4、abp周期波动等、自发abp改变 休息状态的abp改变 方法介绍方法介绍静态脑血流自动调节分析：研究血压变化引起的血流变化脑血流自动调节能力，在没有考虑时间因素条件下。脑血流动脉血压静态脑血流自动调节自动调节的上限：当动脉灌注压升高到维持脑血流量恒定的最高值时，脑血管阻力增高到最大限度，如果动脉灌注压超过该点时，脑血流则呈线性增高，脑血管阻力反而降低。维持脑血流恒定的最高灌注压值，即为脑血流自动调节的上限。脑血流自动调节上限，是机体保护脑组织免受因灌注压过高而使脑血流增高的最大耐受能力。若超过自动调节的上限，则可立即导致脑的过渡灌注，从而引起脑水肿和颅内压增高等症状。自动调节的下限：若动脉灌注

5、压降低到维持脑血流恒定的最低值时，脑血流量固定不变，脑血管阻力达到最低，此值即为脑血流自动调节的下限。若动脉灌注压下降至低于自动调节的下限时，脑血流则呈线性减少。动脉灌注压降低到出现脑缺血症状时，此灌注压值即为机体的最低耐受压。脑血流动态自动调节（dca）动脉血压脑血流速度传统动态分析方法1、tiecks model tiecks1995年提出一个简单的二阶微分方程公式去描述脑血压和脑年提出一个简单的二阶微分方程公式去描述脑血压和脑血流的关系，它是第一个试图去量化大脑自动调节的动态特性。血流的关系，它是第一个试图去量化大脑自动调节的动态特性。 dp = (mabp - iabp)/(iabp

6、- ccp) (4.1) x2 = x2 + (x1 - 2d * x2)/(f * t) (4.2) x1 = x1 + (dp- x2)/(f * t) (4.3)mcav = imcav * (1 + dp - k * x2) (4.4) matlab仿真结果:在这个数学模型中,定义脑自动调节指数ari,从上图底到上依次定ari=0-9,9个等级,9代表脑血流调节的最佳状态,0代表自动调节机能完全消失。 通过下肢束带法实际测得的曲线与之比较,利用最小二乘方法选择最符合的曲线，即为受试者实际的自动调节指数.在算法中t、d、k分别表示时间常量、衰减因子、动态自动调节因数。 由于算法是最先被引

7、入的方法，较为经典，被各国研究者广泛应用。但是tiecks模型完全根据数学公式，引入t、d、k分别表示时间、衰减因子、动态调节因数三个参数，并没有考虑脑血流自动调节的生理机制，同时由于该算法需要外界干扰血压和可重复性较差，所以不是“金”标准。2、transfer function analysis 利用脑血流自动调节的高频滤波特性，同步记录平均动脉压和脑血流速度的自然波动曲线，利用谱分析的方法来评价自动调节功能。把脑血流速度和动脉血压的时域自发曲线经离散傅立叶变换获得频域的动脉血压、脑血流速度功率谱和交叉谱，分析得出脑血流速度和动脉血压相位移。 （f）、v（f）是abp和cbfv的fft频域转

8、换。通过频域转换表明脑血流动态自动调节与频率相关，三个频率范围可能有不同的特性：1、低频部分(0.2hz)有高相关系数、较大幅值和较小相位3、中间频率（0.07-0.2hz）增加的相关系数0.5hz、增加的幅值、减小的相位。( )( )* ( )ppgfe p fp f( )( )* ( )ppgfe p fp f( )( )* ( )ppgfe p fp f( )( )* ( )ppgfe p fp f()()*()ppgfe p fp f()()*()pvgfe p fvf22()()()()p vv vp pgffgfgf()()()p vp pgfhfgf1/222( )( )( )r

9、ih fhfh f1()()tan()irhffhfabp自相关谱abp、cbfv互相关谱平方相关系数传递函数传递函数幅值是由传递函数实部和虚部求得相位3、physiological models生理模型是一种使用生理参数（abp、cbf、icp和cpp等）生理系统的数学表示法，这种模型被用于研究系统的生理和病理机制。 下面是ursinos的描述血流控制系统特征的一种非线性电子电路模型。通过电子原器件模拟脑血流循环动力学系统。使用ursinos提供的参数变量，模型能够仿真出相应的icp和cbfv在生理和病理条件下任意的abp变化。比起上面描述的模型，这个模型考虑了较全面的一系列脑血流调节系统参

10、数。但是这种模型也是在部分理想假设的生理条件下 。图 ursinos 模型4、linear adaptive filter model 在频率0.04hz-0.2hz范围abp和cbfv有较高相关系数表现更强线性关系,使用线性自适应滤波器有：算法简单参数相对较少的特点可以快速跟踪cbfv变换,量化脑血流较强抗干扰能力等优点并且易应用到临床中( ) ( ), (1), (2),., (1)x nx n x nx nx nm12()(),(), .,()mwnwnwnwn()()()tynxnwn()()()()()()tendnyndnxn wn(1 )()2()()wnwnenxn传统传统lm

11、s算法自适应滤波器算法自适应滤波器n时刻输入信号矢量：m是自适应滤波器的阶数，加权向量：滤波器输出：d （n）为系统的期望响应信号，以e（n）值控制权向量使得y （n）逼近d（n），权系数的修改算法。步 骤系统辩识5、nonlinear adaptive filter model由于多方面的原因,volterra级数模型是非线性系统最广泛采用的模型.非线性系统的volterra级数展开是非递归级数组成的, (针对二阶级数和n阶滤波器, volterra lms算法)输入信号和权系数.22( )(1).()( )( )( ) (1).( ) ().() (1)()x nx nx nnx nx n

12、x n x nx n x nnx nn x nnx nn010,00,10,1,( )( ).( )( )( )( ).( ).( )( )nnn nn nw nw nwnwnwnw nwnwnwn， ()()()tynwnxn22( )( ) 2 ( )( ) ( )( ) ( )( ) ( )ttte nd nd n w n x nw n x n x n w n1111(1)( )2( ) ()llwnwne n x nl1, 22212(1)( ) 2( ) () ()l llwnw ne n x n l x n l输出信号：均方误差（mse）目标函数的估计：采用lms算法使目标函数最小化：一阶和二阶项采用不同的收敛因子线性自适应滤波器： mse =6.3391非线性自适应滤波器： mse=2.2875异常信号检测与处理步骤：1、选窗 约1s 采样点 50 2、求pressure and velocity cross-correla