ncepuwanglin
2009-03-07, 22:33
各位大侠，本人的课题关于声源空间定位，其中需要用matlab进行一些信号处理求估计时延值τ，其中需要求互功率谱（即公式（5）），具体过程如下：
假定声源信号为s(t),在时间段T中，两个传声器接收到的信号分别为：
x1(t)=h1(t)*s(t)+n1(t) 0≤t≤T （1）
x2(t)=h2(t)*s(t-τ)+n2(t) 0≤t≤T （2）
其中，h1(t)和h2(t)为房间脉冲响应，τ为待估计时延值，n1(t)和n2(t)为环境噪声。
在频域上，可得：
X1(ω)=H1(ω)S(ω)+N1(ω) （3）
X2(ω)=H2(ω)S(ω)exp(-jωτ)+N2(ω) （4）
可得信号互功率谱的相位为：
θ(ω)=arg[X1(ω)X2'(ω)]=ωτ+υ(ω) （5）
其中，υ(ω)是由于噪声引起的相位。
课题中声源信号s(t）、传声器接收到的信号x1(t)、x12(t)和噪声信号n1(t)、n2(t），均为已知量；公式（1）中的*为卷积符号；公式（5）中的X2'(ω)表示X2(ω)的共轭。
本人matlab不太懂，时间又比较紧张，在这里求助各位高手，请多多指教，不胜感激！
假定声源信号为s(t),在时间段T中，两个传声器接收到的信号分别为：
x1(t)=h1(t)*s(t)+n1(t) 0≤t≤T （1）
x2(t)=h2(t)*s(t-τ)+n2(t) 0≤t≤T （2）
其中，h1(t)和h2(t)为房间脉冲响应，τ为待估计时延值，n1(t)和n2(t)为环境噪声。
在频域上，可得：
X1(ω)=H1(ω)S(ω)+N1(ω) （3）
X2(ω)=H2(ω)S(ω)exp(-jωτ)+N2(ω) （4）
可得信号互功率谱的相位为：
θ(ω)=arg[X1(ω)X2'(ω)]=ωτ+υ(ω) （5）
其中，υ(ω)是由于噪声引起的相位。
课题中声源信号s(t）、传声器接收到的信号x1(t)、x12(t)和噪声信号n1(t)、n2(t），均为已知量；公式（1）中的*为卷积符号；公式（5）中的X2'(ω)表示X2(ω)的共轭。
本人matlab不太懂，时间又比较紧张，在这里求助各位高手，请多多指教，不胜感激！