def getanmval(pvec, B, Y, W):
    """
    Return the SHD for a single time-frequency bin

    :param pvec: Vector containing M-channel STFTs of a single time-frequency bin
    :param B: Response equalisation matrix
    :param Y: SHD matrix
    :param W: Cubature matrix
    :return: SHD for a single time-frequency bin; (N+1)^2 by 1
    """
    anm = B * Y * W * pvec
    return anm