3.算法 Prlt-1=AP:-tA +Q (7) Kg,=P+-1H'/(HP:lt.1H'+R) (9) P::=(I-KgH)XP1:.1 (10)
3. 算法
3.算法 ·初始值设定 在卡尔曼滤波器开始工作之前，应该设置合理 的初始值。假设脑电源强度矩阵X0=0,随着卡尔 曼滤波器的运行，X逐渐收敛接近真实值。P是估 计值X的误差协方差矩阵：P=E{(X·E(X)》(X· E(X)},以协方差矩阵P表示对估计值的置信度， 若P越接近零，则表示估计值越接近真实值。因此 可根据先验知识估计脑电源实际强度X然后由P= E{(X-E(X》(X·E(X)》T确定Po1o
3. 算 法 • 初 始 值 设 定
4.实验结果与讨论 4.1仿真实验结果 (1)模拟实验 本文采用3层同心球的头模型进行仿真实验 该模型的半径分别为79.983m.84.427m, 91.181mm.其对应的电导率分别为0.33Sm, 0.0042S/1m.0.33S/m 源模型：选择位于大脑白质的单个电流偶极子 作为源模型。在头模型的直角坐标系中，源模型的 位置坐标为(21.4,27.3,101.9)
4. 实 验 结 果 与 讨 论 4.1 仿真实验结果 （ 1）模拟实验
4.实验结果与讨论 4.1仿真实验结果 (1)模拟实验 测量电极：在头皮表面安放256个电极.这些电 极按照国际标准10/20导联系统均一地分布在头皮 表面。先通过正向解析求解就可以得到每个电极的 电位值[8】再对测量数据添加一定的高斯噪声（从 5%到50%不等)。实验的任务就是利用含有噪声 的测量数据来反演脑电源分布.实验结果如图1~ 图3所示。并用以下两个指标对实验结果进行评价 和分析：
4. 实 验 结 果 与 讨 论 4.1 仿真实验结果 （ 1）模拟实验
4.实验结果与讨论 4.1仿真实验结果 (2)实验结果评估指标 (1)定位误差（工E):计算得到的最大强度活动 源的位置与实际脑电源位置的偏差。 (2)残差(Residua):反演求得的脑电源通过正 向计算得到的电位分布与测量
4. 实验结果与讨论 4.1 仿真实验结果 （2）实验结果评估指标