mvdr appear The matrix is singular error about athena-signal HOT 3 CLOSED

int dios_ssp_mvdr_init_diffuse_rnn(objMVDR ptr_mvdr)
{
int i, j, k;
float f = 0, temp = 0;
for (k = 1; k < ptr_mvdr->m_sp_size; k++)
{
f = ptr_mvdr->m_deta_fsk;
for (i = 0; i < ptr_mvdr->m_channels; i++)
{
ptr_mvdr->m_sd_rnn_re[kptr_mvdr->m_rxx_size+iptr_mvdr->m_channels+i] = 1.0f*(1+ptr_mvdr->m_sd_factor)+ptr_mvdr->m_sd_eps;
for (j = i+1; j < ptr_mvdr->m_channels; j++)
{
temp = 2 * PI * f * ptr_mvdr->dist[i * ptr_mvdr->m_channels + j] / VELOCITY;
ptr_mvdr->m_sd_rnn_re[kptr_mvdr->m_rxx_size+iptr_mvdr->m_channels+j] = (float)(sin(temp)/temp);
ptr_mvdr->m_sd_rnn_re[kptr_mvdr->m_rxx_size+jptr_mvdr->m_channels+i] = ptr_mvdr->m_sd_rnn_re[kptr_mvdr->m_rxx_size+iptr_mvdr->m_channels+j];
}
}
for (i = 0; i < ptr_mvdr->m_channels; i++)
{
for (j = 0; j < ptr_mvdr->m_channels; j++)
{
ptr_mvdr->m_rxx_in[iptr_mvdr->m_channels2+2j] = ptr_mvdr->m_sd_rnn_re[kptr_mvdr->m_rxx_size+iptr_mvdr->m_channels+j];
ptr_mvdr->m_rxx_in[iptr_mvdr->m_channels2+2j+1] = 0;
}
}
dios_ssp_matrix_inv_process(ptr_mvdr->mvdrinv, ptr_mvdr->m_rxx_in, ptr_mvdr->m_irxx_out);
for (i = 0; i < ptr_mvdr->m_channels; ++i)
{
for (j = 0; j < ptr_mvdr->m_channels; ++j)
{
ptr_mvdr->m_sd_irnn_re[kptr_mvdr->m_rxx_size+iptr_mvdr->m_channels+j] = ptr_mvdr->m_irxx_out[iptr_mvdr->m_channels2+2j];
ptr_mvdr->m_sd_irnn_im[kptr_mvdr->m_rxx_size+iptr_mvdr->m_channels+j] = ptr_mvdr->m_irxx_out[iptr_mvdr->m_channels2+2j+1];
}
}
}

return 0;

}

我debug了 一下，原因是因为temp出现了0，导致ptr_mvdr->m_rxx_in中出现 nan 数据，最后导致求逆矩阵出现了问题。

from athena-signal.

ptr_mvdr->dist[i * ptr_mvdr->m_channels + j] 这个值出现了0值

from athena-signal.

float mic_coord[2][3]={ {-0.05, 0.0, 0.0}, {-0.05, 0.0, 0.0} }; // 麦克风阵列的位置
知道原因了，麦克风坐标设置错了

from athena-signal.