3医学信号处理
本文关键词：医学信号处理
医学信号处理1第三章 线性时不变系 统对随机信号的响应要解决的问题：已知输入信号的统计特性， 其经过线性移不变系统后输出的统计特性类型有连续、 离散系统， 有二端、 多端系统， 要分别进行分析。 医学信号处理2确定信号通过LTI系统的响应：时域为卷积关系y(t)=x(t)*h(t)频域为乘积关系Fy(ω)=H(ω)Fx(ω)x(t) Fx(ω) y(t) Fy(ω) h(t) H(ω)h(t)H(ω)x(t)Fx(ω)y(t)Fy(ω) 医学信号处理3h(t)H(ω)Rx(τ)Sx(ω)Ry(τ)Sy(ω)Rxy(τ)Sxy(ω)当输入信号x(t)是随机信号时， 由于随机信号不存在傅立叶变换， 就不能对时域卷积做傅立叶变换。所以， 研究随机信号通过LTI系统的任务就是： 研究输入和输出间的相关函数和功率谱密度函数间的关系。具体为输出、 输入信号之间自 相关函数和功率谱的关系、 输入和输出间互相关函数和互谱， 并把它们和冲激响应与频率特性联系起来。h(t)H(ω)x(t)Fx(ω)y(t)Fy(ω) 医学信号处理4§ 3.1连续随机信号通过LTI 系统§ 3.1.1二端LTI 系统一、二端LTI 系统本章所有结论的适用条件为输入为平稳的， 而系统为稳定的LTI 系统。 则输出为平稳输出信号。 医学信号处理51、 输出均值式中H(0)是系统的频率特性在ω=0时的值。输出信号的均值等于输入随机信号均值乘以系统频率特性的直流分量。 医学信号处理62、 输入输出互相关函数 3、 输出随机信号的自相关函数输出随机信号的自相关函数等于输入输出互相关函数与系统冲激响应的负卷积， 还等于输入随机信号的自相关函数与系统冲激响应的正负两次卷积。 医学信号处理84、 功率谱 医学信号处理9h(τ)Rx(τ)Rxy(τ)h(-τ)Ry(τ)输入自相关函数、 输出自相关函数、输入输出互相关函数三者之间关系 医学信号处理105、 对协方差有类似公式 医学信号处理11二、两个随机过程的相干函数两个随机过程x(t)和y(t)的相干函数定义一为：如果x(t)和y(t)是某一线性系统的输入和输出如果x(t)和y(t)不相关 医学信号处理12可见， 相干函数是从频域来描述两个随机信号各频率成分之间互相关联的程度。如果输入输出随机信号的相干函数大于0而小于1，则存在两种可能：1、 联系x(t)和y(t)的系统是非线性的；2、 测量值中含有外界噪声， 即x(t)和y(t)是信号和噪声的叠加。 医学信号处理13两个随机序列x(n)和y(n)的相干函数定义二为：用rxy(0)、 rxy(m)代入上式得到： 医学信号处理14§ 3.1.2多端LTI 系统生理信号的测量和处理往往要同时记录多通道信号。 例如， 记录脑电图， 商业脑电图机有128导联， 记录心电图常用3导联和12导联。处理这些信号需要考虑各导信号间的相互作用。这就需要研究多输入、 多输出情况。先从二输入、 二输出开始。 医学信号处理15h11，H11x1(t)y11(t) ++∑y1(t)h22，H22+∑+x2(t)y2(t)y21(t)y12(t)y22(t) 医学信号处理16输入信号的功率谱和互谱所组成的2×2谱阵：输出信号的功率谱和互谱所组成的2×2谱阵：输入、 输出信号之间的互谱所组成的2×2谱阵： 医学信号处理17系统各端口间各频率特性所组成的矩阵：则：（2.1） 医学信号处理18§ 3.2离散随机信号通过LTI 系统∑−∞=k∞−=∗=knhkxnhnxny)()()()()(LTI系统的冲激响应为h(n)， 输入离散随机信号是一个平稳随机过程的一个取样序列， 系统产生的输出信号y(n)也是一个离散随机信号， 它是另一随机过程的一个取样序列。 对于系统来讲， 输入、 输出信号和系统冲激响应存在下列关系：下面研究输出随机序列的特征参数，， 以及输入输出随机序列这些参数的关系。 医学信号处理19∑−∞=k∑−∞=k∞∞−=−==yknxkhEknxkhEnyEm)]()([] )()([)]([xmknxEnxE=−=)]([)]([)()(0 jxxkyeHmmkhm==∑∞−∞=1、 输出随机序列均值H(ej0)为系统的频率特性在ω =0时的值由于输入随机过程为平稳随机过程 医学信号处理20( ,n n)[ ( ) (E y n y n)][( ) (h k x n)( ) (h r x n)]( )h k( ) [ (h r E x n) (k x n)]yykrkrRmmEkm r−m r−∞∞=−∞∑=−∞∑∞∞=−∞∑=−∞∑+=+=−+=−+2、 输出随机过程的自相关序列[ () (k x n)]()xxE x nmrRmkr−+−=R+ −m+由于：( ,n n)( )h k( )h r R()yyxxkrmkr∞∞=−∞∑=−∞∑=+−所以： 医学信号处理21从上式看出： 输出随机序列的自相关序列只与时间差m有关， 而与时间起点的选择（即n的选取） 无关，故有Ryy(n,n+m)=Ryy(m)， 则上式可化简为：( )m( )h k( )h r R()yyxxkrRmkr∞∞=−∞∑=−∞∑=+−( )m()( ) (h k h l)()( )=l( )m( )myyxxlkxxhhxxhhlRRmlkRm l R−RR∞∞=−∞∑=−∞∑∞=−∞∑=−+=∗( )l( ) (h k h l)( )()hhkRkh lhl∞=−∞∑=+=∗−其中，令l = r-k， 可表示为 （可考虑其他方法推导） 医学信号处理22结论： 确定信号的输出等于输入信号与系统的单位冲激响应的卷积；随机信号的输出过程的自相关函数， 等于输入信号的自相关函数与系统的单位冲激响应的自相关函数的线性卷积；确定性离散时间信号作用 于LTI系统与随机性离散时间信号作用 于LTI系统非常类似。( )m()( )=l( )m( )myyxxhhxxhhlRRml RRR∞=−∞∑=−∗ 医学信号处理233、 输出随机过程的功率谱( )m( )m( )myyxxhhRRR=∗上面式子等号两端进行Z变换， 得到：( )z( )z S( )zyyxxhhSS=Sxx(z)和Syy(z)分别等于Rxx(m)和Ryy(m)的Z变换( )z( )m z( )z( )m zmyyyymmxxxxmSRSR+∞∑−=−∞+∞∑−=−∞== 医学信号处理24Shh(z)是Rhh(z)的Z变换， 得到：1( )z( )m z( )()mhhhhmSRH z H z+∞∑−−=−∞==21z( )z( )z H z H( )()( )h n( )z( )z H z( )yyxxyyxxSSSS∗∗==于是， 有：如果为实序列， 则有：( )1( )h n( )z( )()hhH zSH z Hz∗∗=是系统的系统函数如果是复序列， 则( )z( )z S( )zyyxxhhSS=由可得 医学信号处理25物理意义： 一个随机信号通过系统， 从频域看其输出功率谱密度等于输入功率谱密度与系统频率响应的模平方的乘积。1( )z( )()hhSH z H z−=2( )z( )z H z( )yyxxSS=于是得到：|H(z)|为H(z)的模。如果系统稳定， 那么Syy(z)的收敛域包含单位圆， 有2()()()jωjωjωyyxxSeSeH e= 医学信号处理264、 输入随机过程和输出随机过程的互相关函数和互功率谱密度( )m[ ( ) (E x n y n)][ ( )E x n( ) (h k x n)]( ) [ ( ) (h k E x n x n)]( )h k R()( )m( )h mxykkxxxxkRmmkmkmkR∞=−∞∑∞=−∞∞∑=−∞=+=+−=+−=−=∗∑系 统的输入信号与输出信号之间的互相关序列， 等于输入信号自 相关序列与系统冲激响应的线性卷积。 医学信号处理27( )m( )m()yyxyRRhm=∗−( )l( ) (h k h l)( )h m()hhk===Rkhm∞=−∞∑=+=∗−餛( )m( )m( )m( )m( )xxm( )h m( )h m∗()yyxxhhxxRRRRRhm∗∗∗−( )mxyR∴h(m)Rxx(m)Rxy(m)h(-m)Ryy(m)输入自相关函数、 输出自相关函数、输入输出互相关函数三者之间关系 医学信号处理28设（自相关函数的z变换存在） ， 将式(2.8)与式(2.9)转换到z域有，( )( )xySzH z S=( )(yySzH z=用功率谱表示有( )z( )z( )z( )( )( )( )z( )z( )zxyxyxyxxxxxxSSSH z SH z SH z S===( )zxx1)( )zxyS−( )ω()( )ωjωxyxxSH eS=( )ω()( )ωjωyyxySH eS−= 医学信号处理292( )m( )mxxxxRσσ δ=2应用举例：如果输入为零均值平稳白噪声信号， 其方差为,自相关序列为冲激函数则12x( )m[( )]ω2xσ=1( )mxxxxRFSσ δ=( )h m−=( )mxyR2x()()jωjωxyH eSeσ=结论： 输入一个零均值平稳白噪声， 就可以通过计算系统输入输出之间的互相关序列或互功率谱， 来求系统的冲激响应或频率特性。求输入零均值平稳白噪声的LTI系统的单位冲激响应。 医学信号处理30ωaωb1如图为一理想带通滤波器的频率特性H(ejω)， 当输入端作用 一个功率谱为Sxx(ejω)、 均值为0的平稳随机过程，求输出 响应的平均功率。2()()()jωjωjωyyxxSeSeH e=1( )m()2πjωj mωyyyyRSeedωππ−=∫ 医学信号处理311(0)()2πjωyyyyRSedωππ−=∫2n1[](0)()2πjωyyyyE yRSedωππ−==∫2a2n20()()1[](0)()2π1=()()1=()bjωjωyyjωyyyyjωjωxxjωxxSeH eE yRSedωSeH edωSedωππ−πωω∫ωππ==∫∫由于和都是的偶函数， 医学信号处理32a2n1[]=()bjωxxE ySedωωωπ∫上面式子等号右边积分代表输入信号无衰减地通过带通滤波器的总功率， 该总功率除以π， 代表输入信号在该频带内的平均功率， 式子左边为输出信号的平均功率。 医学信号处理332n11[](0)=( )z z dz2π1()2πSyyyycjωyyE yRSjSedωππ−−==∫鹧∫z H z5、 输出随机过程均方值：2( )z( )( )yyxxS=22n11[]( )z H z( )2πxxcE ySz dzj−=∫ 鹧 医学信号处理34例题： 有一LTI系 统， 其系统函数为11( )=10.25H zz−−当 输入端作用 一个白噪声w(n)时， 在其输出端得到一个随机信号x(n)。 设白噪声 方差等于1， 即机信号x(n)的自 相关序列。=1， 求随2xσ 医学信号处理35解：2w111z1114 /15−11( )z( )()(10.25)(0.25)16 /150.2516 /150.25−10.2516 /1514−1xxSH z H zzzzzzzσ=−−−−−−−−=−−=+−=−将上式进行z反变换， 得：16 116( )m( )u m4(1)15 41516 115 4mmxxmRum =+×−− = 医学信号处理36作业：1、 令y(t)=x(t-D)+v(t)为接收信号， 式中D表示信号x(t)的传播延时， 已知x(t)和v(t)是相互独立的零均值的平稳随机过程， 且它们的自相关函数分别为：求互相关函数和互功率谱密度( )τ,( )τ( )xvRAeRBδ τβ τ−==( )τxyR( )xyP ω 医学信号处理37• 2、 将式（2.1） 展开推出的表达式。)(ω)(ω)(ω)(ω22122111yxyxyxyxssss、、、
本文编号：1328342
本文链接：https://www.wllwen.com/wenshubaike/dxkc/1328342.html
