肌组织实验报告
竭诚为您提供优质文档/双击可除
肌组织实验报告
篇一：表面肌实验报告
武汉理工大学
现代数字信号处理在前沿学科中的应用实验报告
基于semg时域特征的动作识别
学院：信息工程学院
学号：姓名：
班级：电子154
实验基于semg时域特征特的动作识别
一、实验目的
1.了解肌电信号常用的时域分析方法；
2.利用mATLAb对肌电信号进行去噪、特征提取及动作识别；
二、实验设备
1.wi-Fi表面肌电信号采集卡；
2.32位windowsxp台式机（matlab7.0软件）；
3.802.11b/g无线网卡；
三、实验内容
（1）学习信号的基本去噪方法，并用mATLAb实现；
（2）学习肌电信号常用的时域特征并利用matlab来进行波形长度（wL）符号改变数（ssc）、过零点（Zc）、威尔
逊赋值（wAmp）等特征的提取；
（3）学习神经网络信号处理方法，掌握bp神经网络的用法，将其用于肌电信号的动作识别。
学习以上三个部分，最终完成一整套肌电信号去噪、特征提取（选取一种特征）、基于特征的动作识别的mATLAb程序。
四、实验原理
（1）小波去噪
小波去噪方法是一种建立在小波变换基础上的新兴算法，基本思想是根据噪声在不同频带上的小波分解系数具有不同强度分布的特点，将各频带上的噪声对应的小系数去除，保留原始信号的小波分解系数，然后对处理后系数进行小波重构，得到纯净信号。
小波去噪的基本原理图如下
（2）特征提取
时域分析是将肌电信号看成均值为零，而方差随着信号强度的变化而变化的随机信号。时域特征的计算复杂度低，提取比较方便。
最常用的方法有：方差，过零点数（Zerocrossing,Zc），willison幅值（willisonAmplitude,wAmp），绝对值平均值(meanAbsoluteValue,mAV）和波形长度（wavelength，wL）等。在实际应用中，为了让特征可以包含更多的信息，往往选择用不同的时域特征组合形成联合特征向量。我们主要介绍一下几种方法：
过零率（Zc）：为波形通过零线的次数,从一定程度上反映了信号的频率特性。为了降低零点引入的噪声，往往会引入一个阈值δ。计算方式如下：
sgn(?xk?xk?1),(xk?xk?1??)（1）willison幅值：是由willison提出一种对表面肌电信号的幅值变化数量进行计
算的方法，经过后人的研究，对willison幅值的阈值有了明确的范围限定，目前认为50~100?V是最合适的阈值范围。其数学表示公式如公式（3-3）。
wAmp??fxi?xi?1
t?1n（2）
?1f(x)???0其中：ifx?阈值otherwise
波形长度（wL）：它是对某一分析窗中的波形长度的统计，波长可以体现该样本的持续时间、幅值、频率的特征。
1n?1
wL??x(i?1)?x(i)ni?1（3）符号改变斜率（ssc）：为信号的的频率性能提供了一些附加信息，对于3个连续的采样
点，给定阈值ω，通过下面的公式计算波峰波谷的个数。
?xi?xi?1???xi?xi?1???,?i?1,?,n?
（3）神经网络（4）
bp神经网络又称误差反向传播（backpropagation），它是一种多层的前向型神经网络。在bp网络中，信号是前向传播的，而误差是反向传播的。所谓的反向传播是指误差的调整过程是从最后的输出层依次向之前各层逐渐进行的。标准的bp网络采用梯度下降算法，与widrow-hoff学习规则相似，网络权值沿着性能函数的梯度反向调整。
前向型神经网络通常具有一个或多个由sigmoid神经元构成的隐层，以及一个由线性神经元构成的输出层。多个具有非线性传递函数的神经元层使得网络可以学习输入和输
出之间的非线性关系，而线性输出层使得网络可以产生[-1，+1]之外的输出值。
输入
隐层输出层
a1?tansig(Iw1,1p1?b1)a2?purelin(Iw2,1a1?b2)
由两层神经元构成的bp网络结构
（1）bp网络的训练算法
①bp算法
bp算法沿着误差函数减小最快的方向，也就是梯度的反方向改变权值和偏差，这一点与线性网络的学习算法是一致
的。bp算法的迭代计算公式可以表示为：
xk?1?xk?akgk（1）
其中，xk代表当前权值和偏差，xk?1代表迭代产生的下一次的权值与偏差，gk为当前误差函数的梯度，ak代表学习速率。
②有动量的梯度下降算法
标准的梯度下降法在调整权值时，仅仅按照当前时刻的负梯度方向进行调整，并没有考虑以前各次运算步骤中的梯度方向，因此新的样本对迭代过程影响太大，可能会导致训练过程中调整方向发生震荡，导致不稳定和收敛速度慢的问题，有动量的梯度下降算法则考虑了往前时刻的贡献，其权值迭代算法为：
wij(n?1)?wij(n)??[(1??)D(n)??D(n?1)（2）
其中，D(n),D(n?1)分别表示n时刻，n-1时刻的负梯度。由于加入了以前时刻梯度的贡献，相当于给迭代过程添加了一个低通滤波器，使得网络忽略误差曲面上细节特征，避免了陷入局部极小点的问题。
③共轭梯度算法
尽管标准的bp算法采用梯度下降算法，权值和偏差沿误差函数下降最快的方向调整，但却并不一定是收敛最快的算法。在改进的bp训练算法中，有一大类的算法称为共轭梯度算法。在这一类算法中，权值和偏差沿着共轭梯度方向
进行调整，通常能够获得比标准的梯度算法更快的收敛速度。
共轭梯度算法的第一次迭代都是从最陡下降的梯度方
向开始。梯度向量为：
p0??g0（3）
沿着此方向进行权值和偏差的调整，公式为：
xk?1?xk?akgk（4）
下一次搜索方向则由前两次搜索方向的共轭方向决定，表达式为：
pk??gk??kpk?1（5）
对于系数?k不同计算方法产生不同的共轭梯度算法。
a）F-R共轭梯度算法采取的系数确定方法为：
Tgkg?k?Tk（6）gk?1gk?1
即本次迭代梯度相对于上一次迭代梯度的归一化值。
b）p-R共轭梯度算法采取的系数确定方法为：
T?gkg?k?T?1k（7）gk?1gk?1
即上次迭代梯度与本次迭代梯度的内积对本次梯度的
归一化值。
c）scaled共轭梯度算法
到目前为止，讨论过的所有共轭梯度算法都需要在每一步迭代过程中对搜索方向进行计算，这样的计算量是比较大。对此moller提出了scaled梯度搜索算法[4]，在每一步迭
代过程中不计算搜索方向，以减少训练过程的计算量。其基
本原理是利用下面介绍的L-m算法与共轭梯度法相结合产生的。
④L-m算法
L-m算法其权值和阈值的更新过程为：
xk?1?xk?[JTJ??I]?1JTe（8）
其中，e为期望输出与实际输出的误差；J为误差对权值微分的Jacobi矩阵；?为标量因子。如果训练成功，误差性能函数减小，那么就减小?的值；反之就减小其值。
五、实验步骤
（1）数据格式转换：
n=4;%通道的个数
m=512;%每个通道的采样数
篇二：肌肉收缩实验报告
骨骼肌收缩实验
一．实验目的
１．肌肉标本收缩现象的描记及单收缩的分析，获得该肌肉收缩的阈值。
２．了解刺激强度对骨骼肌收缩的影响。
３．学习掌握刺激器和张力换能器的使用。
４．加强对神经和肌肉了解，熟练解剖。、
二．实验原理
１．肌肉标本收缩现象的描记
利用刺激器可诱发蛙的离体神经肌肉标本发生兴奋收
缩现象，可利用适当的参数和图形，客观、详细、准确地描述收缩的生理过程与现象。
骨骼肌受到一次短促的阈上刺激时，先是产生一次动作电位，紧接着出现一次机械收缩，称为单收缩。收缩的全过程可分为潜伏期、收缩期和舒张期。在一次单收缩中，肌峰电位的时程（相当于绝对不应期）仅1～2毫秒，而收缩过
程可达几十甚至上百毫秒（蛙的腓肠肌可达100毫秒以上）。
2.张力换能器
换能器是一种能将机械能、化学能、光能等非电量形式的能量转换为电能的器件或装置，并线性相关。利用物理性质和物理效应制成的物理换能器种类繁多，原理各异。张力换能器是一种能把非电量的
生理参数如力、位移等转换为电阻变化的间接型传感器，属于电阻应变式传感器。通常由弹性元件、电阻应变片和其他附件组成。弹性元件采用金属弹性悬梁，可根据机械力的大小选用不同厚度的弹性金属。弹性悬梁的厚度不同，张力换能器的量程亦不同。两组应变片R1、R4及R2、R3分别贴于梁的两面。两组应变片中间接一只调零电位器，并用5～
6V直流电源供电，组成差动式的惠斯登桥式电路（非平衡式电桥）输出电压值与应变片所受力的大小成正比，即力的变化转换成电桥输出电压的变化。此电信号经过记录仪器的放