MATLAB室内多径环境下直达波TDOA估计方法研究(4)
最后将通过DATEMM算法求得的直达波TDOA值与通过传统的几何TDOA算法求得的直达波TDOA值进行比较，得出误差。
分析传感器接收到的白噪声信号与单位脉冲响应卷积后信号的自相关函数图像，本文取其中所有的峰值（根据本文下文设定的参数一般不会超过10个）来代表白噪声通过直接路径和若干个反射路径被传感器接收到的信号。再通过分析两个传感器接收到的白噪声信号与单位脉冲响应卷积后信号的互相关函数图像，本文取其中所有的峰值（根据本文下文设定的参数一般不会超过20个）来代表白噪声通过各自的直接路径和若干个各自的反射路径被两个传感器接收到的信号。记录下取得这些峰值的采样点，再通过采样点之间的差值比较将采集到峰值的采样点与声源通过各自的直接路径传输到两个传感器的TDOA值和声源通过一个直接路径和一个反射路径传输到两个传感器的TDOA值相匹配，从而得到白噪声在混响环境中的直达波的TDOA值。
2.3 语音声源时直达波TDOA估计方法
对于语音声源直达波TDOA值的分析方法和白噪声直达波TDOA值的分析方法基本相同，唯一的区别在于，由于语音声源信号相比白噪声信号功率谱波动要大的多，在得到的传感器接收到的信号自相关函数和互相关函数图像中呈现的峰值也更多且重叠，所以本文在每个图像中取得峰值时只取其中最大的几个峰值。本文在每个自相关函数图像中只取最大的6个峰值来代表语音声源通过直接路径和5个反射路径传输到每个传感器的信号，在每个互相关函数图像中只取最大的10个峰值来代表语音声源通过直接路径或反射路径传输到两个传感器的信号。
3 混响条件下直达波TDOA估计方法的实现与仿真
论文 使用MATLAB软件模拟了一个真实的室内环境并在其中进行声源直达波的TDOA值计算仿真，其模拟出的室内环境要求有：（1）包含一个白声源或语音声源；（2）包含一个由多个传感器组成的传感器阵列组；（3）包含多个反射路径的混响环境。
本文将使用MATLAB软件模拟一个室内多路径的混响环境。本文使用的MATLAB模拟程序是Roomsimove程序。该程序使用MATLAB可以模拟室内移动声源或静止声源（白声源或语音信号声源皆可）通过直接路径和多个反射路径传输到传感器阵列组的情况，并可以自行设置房间尺寸，传感器阵列组尺寸及坐标，传感器阵列组角度，声源类型和声源位置，混响时间和采样频率等多个参数。
在本章中将先进行对白噪声在混响环境中的仿真模拟并将得到的结果进行分析估算其直达波的TDOA值，之后进行对语音声源在混响环境中的仿真模拟并将得到的结果进行分析估算其直达波的TDOA值。
3.1 白噪声声源时直达波TDOA估计方法的实现与仿真
混响环境中白噪声声源的模拟过程：
首先构建一个包含一个白噪声声源的混响环境模型。
（1）设置房间尺寸参数
为了使得到的信号数据集中且不过于分散，房间尺寸设置得比较小，设置为一个长为4米，宽为4米，高为3米的长方体。
（2）设置混响时间参数
为了使得到信号的峰值清晰且相对独立，设置的混响时间为100ms，即可以认为在模拟房间内声源通过较少次数的反射（一般为一次或两次）将信号通过反射路径传输到传感器阵列组。如果混响时间过长，那么模拟房间内声源就很有可能通过多次反射后被传感器阵列组接收到，得到信号的峰值就比较重叠，不利于通过分析识别出TDOA值。通过Roomsimove的M文件RT60toA将混响时间RT60换算为A，A即是每毫秒反射次数。