帮我写个MVDR的程序

MVDR（Minimum Variance Distortionless Response）是一种常见的自适应波束形成算法，用于在有噪声的环境中估计信号源的方向。下面是一个简单的MVDR算法的Python示例代码：

import numpy as np

def mvdr_beamforming(data, steering_vectors):
    """
    MVDR波束形成算法
    :param data: 麦克风阵列采集到的数据，shape为(N, T)，其中N为麦克风数量，T为采样点数
    :param steering_vectors: 麦克风阵列的阵列指向，shape为(N, D)，其中N为麦克风数量，D为阵列指向的数量
    :return: 估计的信号源方向
    """
    # 计算协方差矩阵
    cov = np.dot(data, data.conj().T) / data.shape[1]

    # 计算干扰噪声协方差矩阵的逆
    inv_cov = np.linalg.inv(cov)

    # 计算MVDR滤波器权重
    w = np.dot(inv_cov, steering_vectors) / np.dot(np.dot(steering_vectors.conj().T, inv_cov), steering_vectors)

    # 计算输出
    output = np.dot(w.conj().T, data)

    # 估计信号源方向
    source_direction = np.argmax(np.abs(output))

    return source_direction

其中，`data`为麦克风阵列采集到的数据，`steering_vectors`为麦克风阵列的阵列指向。函数返回估计的信号源方向。请注意，这只是一个简单的示例代码，实际应用中可能需要进行更多的优化和调整。