生物医学信号处理第二章.pptVIP

【例2-10】已知一连续信号为 试用DFT计算其幅度谱，并且与原信号进行比较。 解：由于 是周期信号，首先要确定该信号的周期，如果截断信号时没有完整的一个周期的信号则会发生泄漏现象。 第一部分的周期为0.8ms，第二部分的周期为1ms，第三部分周期为1.25ms，整个信号的周期为20ms，因此截取信号的持续时间最好为20ms的整数倍。 取采样频率和采样间隔分别为： 频率分辨率为 ,很容易满足该条件，因此DFT的点长N从持续时间中算得： 信号的最高频率为 ，因此采样频 率或采样间隔 令 ，代入原信号，利用FFT计算100点的DFT，幅度谱如图2.12所示。或者为了FFT计算，令N为2的幂次，但是要保证NT是20ms的整数倍,比如N取128，T取0.1563ms。只看前半图，从图2.12可以看出信号的频率成分包含k＝16、20、25，即包含三个数字频率成分 幅度谱的大小和频率和原信号完全一致。如果N取得不合适就会发生泄漏现象，整个频率上都会有一定幅度大小。因此要注意取法，当然，正确的采样频率和N的组合是多种的。 图2.12 信号的幅度谱 转换成模拟频率为 或者利用等式： 来求相应的Hz为单位的模拟频率。 第五节 幅值平方相干函数 （Magnitude-Squared Coherent Function） 相干函数表示两个函数之间的振幅、频率和相 位角的关联性，在脑电研究领域中该函数表示大脑不 同部位活动的相干性，这一节学习频域相干函数，也 称为幅值平方相干函数。 频域相干函数，也称为幅值平方相干函数，设有 两个信号 ，它们的幅值平方相干函数定义如 下： (2-40) 其中 表示两个信号的互功率谱,即这两个信 号的循环互相关函数的离散傅立叶变换. 为各 自的功率谱,这里的k代表含义即频率 。 的取值范围为0~1之间. =1, 说明两个信号是完全相干的，即一个信号 可以完全由另外一个信号决定； =0,这两个信号不相干,即这两个信号是完全 独立的； 在(0,1),说明这两个信号存在部分相干性,即 非线性关系或者有外界的干扰存在 可见，频域相干函数可以从频域上表示两个信号 各频率成分互相关联的程度。 例如：频域相干函数可以刻画大脑不同部位的活 动在节律上的一致性，但不能取得大脑活动的瞬态特 征。 【例2-11】设有两个信号 和 ,测量这两个信号 时，假设含有了不相关的噪声 和 ,即测量得到的两个信号为 和 ,比较理想信号和测量信号的幅值相干函数。 解：设 和 的功率谱分别为 ,互功率谱为 , 利用式(2-40)，理想的相干函数为 设噪声的功率谱为 ,由于信号和噪声不相 关，则噪声和信号的互相关函数为零,因此测量信号的 功率谱和互功率谱为： 因此测量信号的幅值相干函数： (2-41) 由于功率谱是非负的，因此上式满足 例2-11就是在有不相关噪声干扰下，频域相干函 数将小于1。一般我们能得到的信号都是测量信号而非 理想信号，因此要从测量信号的相干曲线来判断理想信 号各频率成分互相关联的程度。 下面用实际算例来看相干曲线。 【例2-12】设信号 假设观测时引入的噪声均为白噪声，它们的功率谱密度都为1；观测记录N＝10点，比较理想信号和观测值的幅值相干函数。 解：先利用FFT求出和的功率谱，再利用式（2-40）、（2-41）计算相干函数，MATLAB程序如下： N=10;n=0:N-1; x=N-1:-1:0;y=0:N-1; px=abs(fft(x)).^2/N; py=abs(fft(y)).^2/N; pxy=abs(conj(fft(x)).*fft(y)/N).^2; rxy=pxy./(px.*py); r=1./(1+1./px+1./py+1./(px.*py)); figure(1); subplot(1,2,1); stem(n,px);xlabel(k);title(x(n)功率谱) subplot(1,2,2);