1. 系统的表示方法
定义: 对于某个语音刺激(比如说音素), 将大脑假设成一个线性时不变系统, 大脑的整体响应是对语音刺激的响应加上其他的噪声(实验中应该尽量减小噪声)
$$
R(t)=TRF(t)S(t)+Noise(t)
$$
$𝑅(𝑡)$表示大脑总的神经响应Response , $𝑇𝑅𝐹$表示时间响应函数, $$ 表示卷积, $𝑆(𝑡)$表示语音刺激Stimulus, $Noise(𝑡)$表示测量过程中引入的噪声干扰
2. 已知条件
我们的最终目的是求出$TRF(t)$函数, 现在已知的有
- $S(t)$, 通过对播放的语音信号进行处理, 获取各个音素开始的时间, 我们可以知道语音刺激的函数
- $R(t)$, 通过脑电传感器我们获得了被试者的脑电信号
- 通过滤波函数对脑电信号进行去噪, 获得$R(t)-Noise(t)$
所以我们可以计算出$TRF(t)$, 计算的方法是最小二乘法, 也就是实际脑电信号和通过$TRF$计算出的脑电信号误差要最小.
因为我们获取的脑电信号和播放的语音刺激都是数字信号, 所以我们用离散形式来表示关系. 为了便于表示, 经过滤波之后的脑电信号记为$y(n)$, 语音刺激记为$x(n)$, 时间响应函数记为$h(n)$, 误差项记为$e(n)$, 可以得到
$$
y(n)=x(x)*h(n)+e(n)
$$
为了计算h, 可以将卷积写成矩阵乘积形式:
$$
\left[\begin{array}{c}
{y_{0}} \
{y_{1}} \
{\vdots} \
{y_{N}}
\end{array}\right]=\left[\begin{array}{cccc}
{x_{0}} & {x_{-1}} & {\cdots} & {x_{-k}} \
{x_{1}} & {x_{0}} & {\cdots} & {x_{-k+1}} \
{\vdots} & {\vdots} & {\ddots} & {\vdots} \
{x_{N}} & {x_{N-1}} & {\cdots} & {x_{0}}
\end{array}\right]\left[\begin{array}{c}
{h_{0}} \
{h_{1}} \
{\vdots} \
{h_{k}}
\end{array}\right]+\left[\begin{array}{c}
{e_{0}} \
{e_{1}} \
{\vdots} \
{e_{N}}
\end{array}\right]
$$
其中矩阵的行数$N$表示神经响应的长度
矩阵的列数$k$表示系统的阶数, 也就是神经响应函数的长度
用向量表示为
$$
\vec{y}=X h+\vec{e}
$$