---
title: PowerSpectr
description: "计算复数向量的功率谱，计算公式如下：。"
url: https://www.hikunpeng.com/document/detail/zh/kunpengboostkithistory/2200/accel/kunpengaccel_hmpp_06_0123.html
sourcePath: /source/zh/kunpengboostkithistory/2200/accel/kunpengaccel_hmpp_06_0123.html
indexId: c80b6c7e8c64a38c50aa274d06c6ce4d26cfe203fdc68f920a5e373163c5511175
---
# PowerSpectr

计算复数向量的功率谱，计算公式如下：。

函数接口声明如下：

- 整型数的计算：
  HmppResult HMPPS_PowerSpectr_16s32f(const int16_t* srcRe, const int16_t* srcIm, float* dst, int32_t len);

  HmppResult HMPPS_PowerSpectr_16sc32f(const Hmpp16sc* src, float* dst, int32_t len);

- 浮点数的运算：
  HmppResult HMPPS_PowerSpectr_32f(const float* srcRe, const float* srcIm, float* dst, int32_t len);

  HmppResult HMPPS_PowerSpectr_64f(const double* srcRe, const double* srcIm, double* dst, int32_t len);

  HmppResult HMPPS_PowerSpectr_32fc(const Hmpp32fc* src, float* dst, int32_t len);

  HmppResult HMPPS_PowerSpectr_64fc(const Hmpp64fc* src, double* dst, int32_tl en);

- 带缩放的整型数计算：
  HmppResult HMPPS_PowerSpectr_16s_S(const int16_t* srcRe, cons tint16_t* srcIm, int16_t* dst, int32_t len, double scale);

  HmppResult HMPPS_PowerSpectr_16sc_S(const Hmpp16sc* src, int16_t* dst, int32_t len, double scale);


#### 参数

| 参数名 | 描述 | 取值范围 | 输入/输出 |
| --- | --- | --- | --- |
| srcRe | 指向输入实部向量的指针。 | 非空 | 输入 |
| srcIm | 指向输入虚部向量的指针。 | 非空 | 输入 |
| src | 指向输入复数向量的指针。 | 非空 | 输入 |
| dst | 指向输出功率谱密度向量的指针。 | 非空 | 输出 |
| scale | 缩放因数。 | (0,INF)且输入为2^n | 输入 |
| len | 向量长度。 | (0, INT\_MAX] | 输入 |


#### 返回值

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


#### 错误码

| 错误码 | 描述 |
| --- | --- |
| HMPP\_STS\_NULL\_PTR\_ERR | srcRe，srcIm或src为空指针。 |
| HMPP\_STS\_SIZE\_ERR | len小于或等于0。 |
| HMPP\_STS\_SCALE\_ERR | scale不在(0,INF)范围内或输入为nan。 |


#### 示例

```
#define BUFFER_SIZE_T 5
void PowerSpectrExample(void)
{
float srcRe[BUFFER_SIZE_T] = {0.0, 1.0, 2.0, 4.5, 150.0};
float srcIm[BUFFER_SIZE_T] = {0.5, 1.5, 2.5, 4.0, 150.0};
float dst[BUFFER_SIZE_T];
int32_t i;
HmppResult result = HMPPS_PowerSpectr_32f(srcRe, srcIm, dst, BUFFER_SIZE_T);
printf("result = %d \n", result);
if (result != HMPP_STS_NO_ERR) {
return;
}
printf("dst =");
for (int i = 0; i < BUFFER_SIZE_T; i++) {
printf(" %f  ", dst[i]);
}
printf("\n");
}
```

运行结果：

```
result = 0
dst = 0.250000  3.250000  10.250000  36.250000  45000.000000
```