Pwelch

Welch方法是一种修正周期图功率谱密度估计方法，它通过选取的窗口对数据进行加窗处理，先分段求功率谱之后再进行平均。

函数接口声明如下：

HmppResult HMPPS_Pwelch_32f(const float* src, const float* window, float* dst, int32_t len, double overLap, int32_t nfft, int32_t windowLen);

HmppResult HMPPS_Pwelch_64f(const double* src, const double* window, double* dst, int32_t len, double overLap, int32_t nfft, int32_t windowLen);

HmppResult HMPPS_Pwelch_32fc(const Hmpp32fc* src, const float* window, float* dst, int32_t len, double overLap, int32_t nfft, int32_t windowLen);

HmppResult HMPPS_Pwelch_64fc(const Hmpp64fc* src, const double* window, double* dst, int32_t len, double overLap, int32_t nfft, int32_t windowLen);

参数

参数名	描述	取值范围	输入/输出
src	指向保存源信号向量的数组。	非空	输入
window	指向保存窗向量的数组。	可以为空。为空时默认为全位1，长度为windowLen。	输入
dst	功率谱密度的估计数组。	非空	输出
overlap	相邻两段数据之间的重叠部分占window长度的比例。	(0, 0.95]	输入
nfft	FFT点数。	[INT_MIN, windowLen]，小于零时默认为windowLen。	输入
len	源信号向量长度。	(0, INT_MAX]	输入
windowLen	窗向量长度。	(0, len]	输入

返回值

成功：返回HMPP_STS_NO_ERR。
失败：返回错误码。

错误码

错误码	描述
HMPP_STS_NULL_PTR_ERR	src、dst这两个入参中存在空指针。
HMPP_STS_SIZE_ERR	len小于或等于0。
HMPP_SYS_MALLOC_FAILED	函数执行中出现malloc失败。

示例

#define  BUFFER_SIZE_T 5
#define  NFFT 2
#define  WINDOWLEN 2
#define  OVERLAP 0.5
void PwelchExample()
{
    float src[BUFFER_SIZE_T] = {1, 2.5, 3.3, 1, 5};
    float window[2] = {1,1};
    float dst[BUFFER_SIZE_T];
    HMPPS_Zero_32f(dst, BUFFER_SIZE_T); // 数组初始化，将dst所有元素初始化为0.
    HmppResult result = HMPPS_Pwelch_32f(src, window, dst, BUFFER_SIZE_T, OVERLAP, WINDOWLEN, NFFT);
    if (result == HMPP_STS_NO_ERR) {
        printf("dst = ");
        for (int32_t i = 0; i < (NFFT + 1) / 2; i++) {
            printf("%.2f ", dst[i]);
        }
        printf("\n");
    }
}

运行结果：

dst = 0.97 0.18

父主题： 基础与通用运算