【STM32F407的DSP教程】第21章 DSP矩阵运算-加法，减法和逆矩阵

article/2025/10/7 13:02:07

完整版教程下载地址：http://www.armbbs.cn/forum.php?mod=viewthread&tid=94547

第21章       DSP矩阵运算-加法，减法和逆矩阵

本期教程主要讲解矩阵运算中的初始化，加法，逆矩阵和减法。

目录

第21章       DSP矩阵运算-加法，减法和逆矩阵

21.1 初学者重要提示

21.2 DSP基础运算指令

21.3 矩阵初始化（MatInit）

21.3.1 函数arm_mat_init_f32

21.3.2 函数arm_mat_init_q31

21.3.3 函数arm_mat_init_q15

21.3.4 使用举例

21.4 矩阵加法（MatAdd）

21.4.1 函数arm_mat_add_f32

21.4.2 函数arm_mat_add_q31

21.4.3 函数arm_mat_add_q15

21.4.4 用举例（含Matlab实现）

21.5 矩阵减法（MatSub）

21.5.1 函数arm_mat_sub_f32

21.5.2 函数arm_mat_sub_q31

21.5.3 函数arm_mat_sub_q15

21.5.4 使用举例（含Matlab实现）

21.6 逆矩阵（MatInverse）

21.6.1 函数arm_mat_inverse_f64

21.6.2 函数arm_mat_inverse_f32

21.6.3 使用举例（含Matlab实现）

21.7 实验例程说明（MDK）

21.8 实验例程说明（IAR）

21.9 总结

 

21.1 初学者重要提示

  1.   复数运算比较重要，后面FFT章节要用到，如果印象不深的话，需要温习下高数知识了。
  2.   ARM提供的DSP库逆矩阵求法有局限性，通过Matlab验证是可以求逆矩阵的，而DSP库却不能正确求解。

21.2 DSP基础运算指令

本章用到的DSP指令在前面章节都已经讲解过。

21.3 矩阵初始化（MatInit）

主要用于矩阵结构体成员的初始化，浮点格式矩阵结构体定义如下：

typedef struct
{
uint16_t numRows;     // 矩阵行数.
uint16_t numCols;     // 矩阵列数
float32_t *pData;     // 矩阵地址
} arm_matrix_instance_f32

定点数Q31格式矩阵结构体定义如下：

typedef struct
{uint16_t numRows;     //矩阵行数uint16_t numCols;     //矩阵列数q31_t *pData;         //矩阵地址
} arm_matrix_instance_q31;

定点数Q15格式矩阵结构体定义如下：

typedef struct
{uint16_t numRows;     //矩阵行数uint16_t numCols;     //矩阵列数q15_t *pData;         //矩阵地址
} arm_matrix_instance_q15;

 

21.3.1 函数arm_mat_init_f32

函数原型：

void arm_mat_init_f32(

  arm_matrix_instance_f32 * S,

  uint16_t nRows,

  uint16_t nColumns,

  float32_t * pData)

函数描述：

这个函数用于浮点格式的矩阵数据初始化。

函数参数：

  •   第1个参数是arm_matrix_instance_f32类型矩阵结构体指针变量。
  •  第2个参数是矩阵行数。
  •   第3个参数是矩阵列数。
  •   第4个参数是矩阵数据地址。

 

21.3.2 函数arm_mat_init_q31

函数原型：

void arm_mat_init_f32(

  arm_matrix_instance_f32 * S,

  uint16_t nRows,

  uint16_t nColumns,

  float32_t * pData)

函数描述：

这个函数用于定点数Q31格式的矩阵数据初始化。

函数参数：

  •   第1个参数是arm_matrix_instance_q31类型矩阵结构体指针变量。
  •   第2个参数是矩阵行数。
  •   第3个参数是矩阵列数。
  •   第4个参数是矩阵数据地址。

 

21.3.3 函数arm_mat_init_q15

函数原型：

void arm_mat_init_f32(

  arm_matrix_instance_f32 * S,

  uint16_t nRows,

  uint16_t nColumns,

  float32_t * pData)

函数描述：

这个函数用于定点数Q15格式的矩阵数据初始化。

函数参数：

  •   第1个参数是arm_matrix_instance_q15类型矩阵结构体指针变量。
  •   第2个参数是矩阵行数。
  •   第3个参数是矩阵列数。
  •   第4个参数是矩阵数据地址。

21.3.4 使用举例

程序设计：

/*
*********************************************************************************************************
*    函 数 名: DSP_MatInit
*    功能说明: 矩阵数据初始化
*    形    参：无
*    返 回 值: 无
*********************************************************************************************************
*/
static void DSP_MatInit(void)
{uint8_t i;/****浮点数数组******************************************************************/float32_t pDataA[9] = {1.1f, 1.1f, 2.1f, 2.1f, 3.1f, 3.1f, 4.1f, 4.1f, 5.1f};arm_matrix_instance_f32 pSrcA; //3行3列数据/****定点数Q31数组******************************************************************/q31_t pDataA1[9] = {1, 1, 2, 2, 3, 3, 4, 4, 5};arm_matrix_instance_q31 pSrcA1; //3行3列数据/****定点数Q15数组******************************************************************/q15_t pDataA2[9] = {1, 1, 2, 2, 3, 3, 4, 4, 5};arm_matrix_instance_q15 pSrcA2; //3行3列数据/****浮点数***********************************************************************/    printf("****浮点数******************************************\r\n");arm_mat_init_f32(&pSrcA, 3,3, pDataA);for(i = 0; i < 9; i++){printf("pDataA[%d] = %f\r\n", i, pDataA[i]);}/****定点数Q31***********************************************************************/printf("****浮点数******************************************\r\n");arm_mat_init_q31(&pSrcA1, 3,3, pDataA1);for(i = 0; i < 9; i++){printf("pDataA1[%d] = %d\r\n", i, pDataA1[i]);}/****定点数Q15***********************************************************************/printf("****浮点数******************************************\r\n");arm_mat_init_q15(&pSrcA2, 3,3, pDataA2);for(i = 0; i < 9; i++){printf("pDataA2[%d] = %d\r\n", i, pDataA2[i]);}
}

 

实验现象（按下K1按键后串口打印模平方）：

21.4 矩阵加法（MatAdd）

以3*3矩阵为例，矩阵加法的实现公式如下：

21.4.1 函数arm_mat_add_f32

函数原型：

arm_status arm_mat_add_f32(

  const arm_matrix_instance_f32 * pSrcA,

  const arm_matrix_instance_f32 * pSrcB,

  arm_matrix_instance_f32 * pDst)

函数描述：

这个函数用于浮点数的矩阵加法。

函数参数：

  •   第1个参数是矩阵A的源地址。
  •   第2个参数是矩阵B的源地址。
  •   第3个参数是矩阵A + 矩阵B计算结果存储的地址。
  •   返回值，ARM_MATH_SUCCESS表示成功，ARM_MATH_SIZE_MISMATCH表示矩阵大小不一致。

注意事项：

  1. pSrcA，pSrcB，pDst的行数和列数必须是相同的，否则没有办法使用加法运算。
  2. 矩阵在数组中的存储是从左到右，再从上到下。

21.4.2 函数arm_mat_add_q31

函数原型：

arm_status arm_mat_add_f32(

  const arm_matrix_instance_f32 * pSrcA,

  const arm_matrix_instance_f32 * pSrcB,

  arm_matrix_instance_f32 * pDst)

函数描述：

这个函数用于定点数Q31的矩阵加法。

函数参数：

  •   第1个参数是矩阵A的源地址。
  •   第2个参数是矩阵B的源地址。
  •   第3个参数是矩阵A + 矩阵B计算结果存储的地址。
  •   返回值，ARM_MATH_SUCCESS表示成功，ARM_MATH_SIZE_MISMATCH表示矩阵大小不一致。

注意事项：

  1. 使用了饱和运算，输出结果范围[0x80000000 0x7FFFFFFF]。
  2. pSrcA，pSrcB，pDst的行数和列数必须是相同的，否则没有办法使用加法运算。
  3. 矩阵在数组中的存储是从左到右，再从上到下。

 

21.4.3 函数arm_mat_add_q15

函数原型：

arm_status arm_mat_add_f32(

  const arm_matrix_instance_f32 * pSrcA,

  const arm_matrix_instance_f32 * pSrcB,

  arm_matrix_instance_f32 * pDst)

函数描述：

这个函数用于定点数Q15的矩阵加法。

函数参数：

  •   第1个参数是矩阵A的源地址。
  •   第2个参数是矩阵B的源地址。
  •   第3个参数是矩阵A + 矩阵B计算结果存储的地址。
  •   返回值，ARM_MATH_SUCCESS表示成功，ARM_MATH_SIZE_MISMATCH表示矩阵大小不一致。

注意事项：

  1. 使用了饱和运算，输出结果范围[0x8000 0x7FFF]。
  2. pSrcA，pSrcB，pDst的行数和列数必须是相同的，否则没有办法使用加法运算。
  3. 矩阵在数组中的存储是从左到右，再从上到。

21.4.4 用举例（含Matlab实现）

程序设计：

/*
*********************************************************************************************************
*    函 数 名: DSP_MatAdd
*    功能说明: 矩阵求和
*    形    参：无
*    返 回 值: 无
*********************************************************************************************************
*/
static void DSP_MatAdd(void)
{uint8_t i;/****浮点数数组******************************************************************/float32_t pDataA[9] = {1.1f, 1.1f, 2.1f, 2.1f, 3.1f, 3.1f, 4.1f, 4.1f, 5.1f};float32_t pDataB[9] = {1.1f, 1.1f, 2.1f, 2.1f, 3.1f, 3.1f, 4.1f, 4.1f, 5.1f};float32_t pDataDst[9];arm_matrix_instance_f32 pSrcA; //3行3列数据arm_matrix_instance_f32 pSrcB; //3行3列数据arm_matrix_instance_f32 pDst;/****定点数Q31数组******************************************************************/q31_t pDataA1[9] = {1, 1, 2, 2, 3, 3, 4, 4, 5};q31_t pDataB1[9] = {1, 1, 2, 2, 3, 3, 4, 4, 5};q31_t pDataDst1[9];arm_matrix_instance_q31 pSrcA1; //3行3列数据arm_matrix_instance_q31 pSrcB1; //3行3列数据arm_matrix_instance_q31 pDst1;/****定点数Q15数组******************************************************************/q15_t pDataA2[9] = {1, 1, 2, 2, 3, 3, 4, 4, 5};q15_t pDataB2[9] = {1, 1, 2, 2, 3, 3, 4, 4, 5};q15_t pDataDst2[9];arm_matrix_instance_q15 pSrcA2; //3行3列数据arm_matrix_instance_q15 pSrcB2; //3行3列数据arm_matrix_instance_q15 pDst2;/****浮点数***********************************************************************/pSrcA.numCols = 3;pSrcA.numRows = 3;pSrcA.pData = pDataA;pSrcB.numCols = 3;pSrcB.numRows = 3;pSrcB.pData = pDataB;pDst.numCols = 3;pDst.numRows = 3;pDst.pData = pDataDst;printf("****浮点数******************************************\r\n");arm_mat_add_f32(&pSrcA, &pSrcB, &pDst);for(i = 0; i < 9; i++){printf("pDataDst[%d] = %f\r\n", i, pDataDst[i]);}/****定点数Q31***********************************************************************/pSrcA1.numCols = 3;pSrcA1.numRows = 3;pSrcA1.pData = pDataA1;pSrcB1.numCols = 3;pSrcB1.numRows = 3;pSrcB1.pData = pDataB1;pDst1.numCols = 3;pDst1.numRows = 3;pDst1.pData = pDataDst1;printf("****定点数Q31******************************************\r\n");arm_mat_add_q31(&pSrcA1, &pSrcB1, &pDst1);for(i = 0; i < 9; i++){printf("pDataDst1[%d] = %d\r\n", i, pDataDst1[i]);}/****定点数Q15***********************************************************************/pSrcA2.numCols = 3;pSrcA2.numRows = 3;pSrcA2.pData = pDataA2;pSrcB2.numCols = 3;pSrcB2.numRows = 3;pSrcB2.pData = pDataB2;pDst2.numCols = 3;pDst2.numRows = 3;pDst2.pData = pDataDst2;printf("****定点数Q15******************************************\r\n");arm_mat_add_q15(&pSrcA2, &pSrcB2, &pDst2);for(i = 0; i < 9; i++){printf("pDataDst2[%d] = %d\r\n", i, pDataDst2[i]);}
}

 

实验现象（按下K2按键后串口打印矩阵加法）：

下面通过Matlab来求解矩阵和（在命令窗口输入）：

 

21.5 矩阵减法（MatSub）

以3*3矩阵为例，矩阵减法的实现公式如下：

21.5.1 函数arm_mat_sub_f32

函数原型：

arm_status arm_mat_sub_f32(

  const arm_matrix_instance_f32 * pSrcA,

  const arm_matrix_instance_f32 * pSrcB,

  arm_matrix_instance_f32 * pDst)

函数描述：

这个函数用于浮点数的矩阵减法。

函数参数：

  •   第1个参数是矩阵A的源地址。
  •   第2个参数是矩阵B的源地址。
  •   第3个参数是矩阵A减去矩阵B计算结果存储的地址。
  •   返回值，ARM_MATH_SUCCESS表示成功，ARM_MATH_SIZE_MISMATCH表示矩阵大小不一致。

注意事项：

  1. pSrcA，pSrcB，pDst的行数和列数必须是相同的，否则没有办法使用加法运算。
  2. 矩阵在数组中的存储是从左到右，再从上到下。

21.5.2 函数arm_mat_sub_q31

函数原型：

arm_status arm_mat_sub_q31(

  const arm_matrix_instance_q31 * pSrcA,

  const arm_matrix_instance_q31 * pSrcB,

        arm_matrix_instance_q31 * pDst)

函数描述：

这个函数用于定点数Q31的矩阵减法。

函数参数：

  •   第1个参数是矩阵A的源地址。
  •   第2个参数是矩阵B的源地址。
  •   第3个参数是矩阵A减去矩阵B计算结果存储的地址。
  •   返回值，ARM_MATH_SUCCESS表示成功，ARM_MATH_SIZE_MISMATCH表示矩阵大小不一致。

注意事项：

  1. 使用了饱和运算，输出结果范围[0x80000000 0x7FFFFFFF]。
  2. pSrcA，pSrcB，pDst的行数和列数必须是相同的，否则没有办法使用加法运算。
  3. 矩阵在数组中的存储是从左到右，再从上到下。

21.5.3 函数arm_mat_sub_q15

函数原型：

arm_status arm_mat_sub_q15(

  const arm_matrix_instance_q15 * pSrcA,

  const arm_matrix_instance_q15 * pSrcB,

        arm_matrix_instance_q15 * pDst)

函数描述：

这个函数用于定点数Q15的矩阵减法。

函数参数：

  •   第1个参数是矩阵A的源地址。
  •   第2个参数是矩阵B的源地址。
  •   第3个参数是矩阵A减去矩阵B计算结果存储的地址。
  •   返回值，ARM_MATH_SUCCESS表示成功，ARM_MATH_SIZE_MISMATCH表示矩阵大小不一致。

注意事项：

  1. 使用了饱和运算，输出结果范围[0x8000 0x7FFF]。
  2. pSrcA，pSrcB，pDst的行数和列数必须是相同的，否则没有办法使用加法运算。
  3. 矩阵在数组中的存储是从左到右，再从上。

 

21.5.4 使用举例（含Matlab实现）

程序设计：

/*
*********************************************************************************************************
*    函 数 名: DSP_MatSub
*    功能说明: 矩阵减法
*    形    参：无
*    返 回 值: 无
*********************************************************************************************************
*/
static void DSP_MatSub(void)
{uint8_t i;/****浮点数数组******************************************************************/float32_t pDataA[9] = {1.1f, 1.1f, 2.1f, 2.1f, 3.1f, 3.1f, 4.1f, 4.1f, 5.1f};float32_t pDataB[9] = {1.1f, 1.1f, 2.1f, 2.1f, 3.1f, 3.1f, 4.1f, 4.1f, 5.1f};float32_t pDataDst[9];arm_matrix_instance_f32 pSrcA; //3行3列数据arm_matrix_instance_f32 pSrcB; //3行3列数据arm_matrix_instance_f32 pDst;/****定点数Q31数组******************************************************************/q31_t pDataA1[9] = {1, 1, 2, 2, 3, 3, 4, 4, 5};q31_t pDataB1[9] = {2, 2, 2, 2, 2, 2, 2, 2, 2};q31_t pDataDst1[9];arm_matrix_instance_q31 pSrcA1; //3行3列数据arm_matrix_instance_q31 pSrcB1; //3行3列数据arm_matrix_instance_q31 pDst1;/****定点数Q15数组******************************************************************/q15_t pDataA2[9] = {1, 1, 2, 2, 3, 3, 4, 4, 5};q15_t pDataB2[9] = {2, 2, 2, 2, 23, 2, 2, 2, 2};q15_t pDataDst2[9];arm_matrix_instance_q15 pSrcA2; //3行3列数据arm_matrix_instance_q15 pSrcB2; //3行3列数据arm_matrix_instance_q15 pDst2;/****浮点数***********************************************************************/pSrcA.numCols = 3;pSrcA.numRows = 3;pSrcA.pData = pDataA;pSrcB.numCols = 3;pSrcB.numRows = 3;pSrcB.pData = pDataB;pDst.numCols = 3;pDst.numRows = 3;pDst.pData = pDataDst;printf("****浮点数******************************************\r\n");arm_mat_sub_f32(&pSrcA, &pSrcB, &pDst);for(i = 0; i < 9; i++){printf("pDataDst[%d] = %f\r\n", i, pDataDst[i]);}/****定点数Q31***********************************************************************/pSrcA1.numCols = 3;pSrcA1.numRows = 3;pSrcA1.pData = pDataA1;pSrcB1.numCols = 3;pSrcB1.numRows = 3;pSrcB1.pData = pDataB1;pDst1.numCols = 3;pDst1.numRows = 3;pDst1.pData = pDataDst1;printf("****定点数Q31******************************************\r\n");arm_mat_sub_q31(&pSrcA1, &pSrcB1, &pDst1);for(i = 0; i < 9; i++){printf("pDataDst1[%d] = %d\r\n", i, pDataDst1[i]);}/****定点数Q15***********************************************************************/pSrcA2.numCols = 3;pSrcA2.numRows = 3;pSrcA2.pData = pDataA2;pSrcB2.numCols = 3;pSrcB2.numRows = 3;pSrcB2.pData = pDataB2;pDst2.numCols = 3;pDst2.numRows = 3;pDst2.pData = pDataDst2;printf("****定点数Q15******************************************\r\n");arm_mat_sub_q15(&pSrcA2, &pSrcB2, &pDst2);for(i = 0; i < 9; i++){printf("pDataDst2[%d] = %d\r\n", i, pDataDst2[i]);}
}

 

实验现象（按下OK按键后串口打印矩阵减法）：

下面通过Matlab来求解矩阵减法（在命令窗口输入）。

21.6 逆矩阵（MatInverse）

以3*3矩阵为例，逆矩阵的实现公式如下（Gauss-Jordan法求逆矩阵）：

21.6.1 函数arm_mat_inverse_f64

函数原型：

arm_status arm_mat_inverse_f64(

  const arm_matrix_instance_f64 * pSrc,

        arm_matrix_instance_f64 * pDst)

函数描述：

这个函数用于64bit浮点数的逆矩阵求解。

函数参数：

  •   第1个参数是矩阵源地址。
  •   第2个参数是求逆后的矩阵地址。
  •   返回值，ARM_MATH_SUCCESS表示成功，ARM_MATH_SIZE_MISMATCH表示矩阵大小不一致。ARM_MATH_SINGULAR表示矩阵不可逆。

注意事项：

  1. pSrc必须得是方阵（行数和列数相同）。
  2. pSrc和pDst必须是相同的方阵。
  3. 输入的矩阵可逆，函数会返回ARM_MATH_SUCCESS，如果不可逆，返回ARM_MATH_SINGULAR。
  4. ARM官方库只提供了浮点数矩阵求逆矩阵。

21.6.2 函数arm_mat_inverse_f32

函数原型：

arm_status arm_mat_inverse_f32(

  const arm_matrix_instance_f32 * pSrc,

  arm_matrix_instance_f32 * pDst)

函数描述：

这个函数用于32bit浮点数的逆矩阵求解。

函数参数：

  •   第1个参数是矩阵源地址。
  •   第2个参数是求逆后的矩阵地址。
  •   返回值，ARM_MATH_SUCCESS表示成功，ARM_MATH_SIZE_MISMATCH表示矩阵大小不一致。ARM_MATH_SINGULAR表示矩阵不可逆。

注意事项：

  1. pSrc必须得是方阵（行数和列数相同）。
  2. pSrc和pDst必须是相同的方阵。
  3. 输入的矩阵可逆，函数会返回ARM_MATH_SUCCESS，如果不可逆，返回ARM_MATH_SINGULAR。
  4. ARM官方库只提供了浮点数矩阵求逆矩阵。   

21.6.3 使用举例（含Matlab实现）

程序设计：

/*
*********************************************************************************************************
*    函 数 名: DSP_MatInverse
*    功能说明: 求逆矩阵
*    形    参: 无
*    返 回 值: 无
*********************************************************************************************************
*/
static void DSP_MatInverse(void)
{uint8_t i;arm_status sta;/****浮点数数组******************************************************************/float32_t pDataB[36];float32_t pDataA[36] = {1.0f,   0.0f,   0.0f,  0.0f,  0.0f,  0.0f,0.0f,   1.0f,   0.0f,  0.0f,  1.0f,  0.0f,0.0f,   0.0f,   2.0f,  0.0f,  0.0f,  0.0f,0.0f,   0.0f,   0.0f,  2.0f,  0.0f,  1.0f,0.0f,   0.0f,   0.0f,  0.0f,  3.0f,  0.0f,0.0f,   0.0f,   0.0f,  0.0f,  0.0f,  4.0f};arm_matrix_instance_f32 pSrcA; //6行6列数据arm_matrix_instance_f32 pSrcB; //6行6列数据;/****浮点数***********************************************************************/pSrcA.numCols = 6;pSrcA.numRows = 6;pSrcA.pData = pDataA;pSrcB.numCols = 6;pSrcB.numRows = 6;pSrcB.pData = pDataB;sta = arm_mat_inverse_f32(&pSrcA, &pSrcB);/*sta = ARM_MATH_SUCCESS, 即返回0，表示求逆矩阵成功。sta = ARM_MATH_SINGULAR, 即返回-5，表示求逆矩阵失败，也表示不可逆。注意，ARM提供的DSP库逆矩阵求发有局限性，通过Matlab验证是可以求逆矩阵的，而DSP库却不能正确求解。*/printf("----sta %d\r\n", sta);for(i = 0; i < 36; i++){printf("pDataB[%d] = %f\r\n", i, pDataB[i]);}    
}

 

实验现象（按下K3按键后串口打印逆矩阵）：

下面我们通过Matlab来实现求逆矩阵（在命令窗口输入）：

21.7 实验例程说明（MDK）

配套例子：

V6-216_DSP矩阵运算（加法，减法和逆矩阵）

实验目的：

  1. 学习DSP复数运算（加法，减法和逆矩阵）

实验内容：

  1. 启动一个自动重装软件定时器，每100ms翻转一次LED2。
  2. 按下按键K1，串口打函数DSP_MatInit的输出数据。
  3. 按下按键K2，串口打函数DSP_MatAdd的输出数据。
  4. 按下按键K3，串口打函数DSP_MatInverse的输出数据。
  5. 按下摇杆OK键，串口打函数DSP_MatSub的输出数据。

使用AC6注意事项

特别注意附件章节C的问题

上电后串口打印的信息：

波特率 115200，数据位 8，奇偶校验位无，停止位 1。

详见本章的3.4 ，4.4，5.4和6.3小节。

程序设计：

  系统栈大小分配：

  硬件外设初始化

硬件外设的初始化是在 bsp.c 文件实现：

/*
*********************************************************************************************************
*    函 数 名: bsp_Init
*    功能说明: 初始化所有的硬件设备。该函数配置CPU寄存器和外设的寄存器并初始化一些全局变量。只需要调用一次
*    形    参：无
*    返 回 值: 无
*********************************************************************************************************
*/
void bsp_Init(void)
{/* STM32F407 HAL 库初始化，此时系统用的还是F407自带的16MHz，HSI时钟:- 调用函数HAL_InitTick，初始化滴答时钟中断1ms。- 设置NVIV优先级分组为4。*/HAL_Init();/* 配置系统时钟到168MHz- 切换使用HSE。- 此函数会更新全局变量SystemCoreClock，并重新配置HAL_InitTick。*/SystemClock_Config();/* Event Recorder：- 可用于代码执行时间测量，MDK5.25及其以上版本才支持，IAR不支持。- 默认不开启，如果要使能此选项，务必看V5开发板用户手册第8章*/    
#if Enable_EventRecorder == 1  /* 初始化EventRecorder并开启 */EventRecorderInitialize(EventRecordAll, 1U);EventRecorderStart();
#endifbsp_InitKey();        /* 按键初始化，要放在滴答定时器之前，因为按钮检测是通过滴答定时器扫描 */bsp_InitTimer();      /* 初始化滴答定时器 */bsp_InitUart();    /* 初始化串口 */bsp_InitExtIO();   /* 初始化扩展IO */bsp_InitLed();        /* 初始化LED */        
}

 

  主功能：

主程序实现如下操作：

  •   启动一个自动重装软件定时器，每100ms翻转一次LED2。
  •   按下按键K1，串口打函数DSP_MatInit的输出数据。
  •   按下按键K2，串口打函数DSP_MatAdd的输出数据。
  •   按下按键K3，串口打函数DSP_MatInverse的输出数据。
  •   按下摇杆OK键，串口打函数DSP_MatSub的输出数据。
/*
*********************************************************************************************************
*    函 数 名: main
*    功能说明: c程序入口
*    形    参: 无
*    返 回 值: 错误代码(无需处理)
*********************************************************************************************************
*/
int main(void)
{uint8_t ucKeyCode;        /* 按键代码 */bsp_Init();        /* 硬件初始化 */PrintfLogo();    /* 打印例程信息到串口1 */PrintfHelp();    /* 打印操作提示信息 */bsp_StartAutoTimer(0, 100);    /* 启动1个100ms的自动重装的定时器 *//* 进入主程序循环体 */while (1){bsp_Idle();        /* 这个函数在bsp.c文件。用户可以修改这个函数实现CPU休眠和喂狗 *//* 判断定时器超时时间 */if (bsp_CheckTimer(0))    {/* 每隔100ms 进来一次 */  bsp_LedToggle(2);}ucKeyCode = bsp_GetKey();    /* 读取键值, 无键按下时返回 KEY_NONE = 0 */if (ucKeyCode != KEY_NONE){switch (ucKeyCode){case KEY_DOWN_K1:            /* 按下按键K1，串口打函数DSP_MatInit的输出数据 */DSP_MatInit();break;case KEY_DOWN_K2:            /* 按下按键K2，串口打函数DSP_MatAdd的输出数据 */DSP_MatAdd();break;case KEY_DOWN_K3:            /* 按下按键K3，串口打函数DSP_MatInverse的输出数据 */DSP_MatInverse();break;case JOY_DOWN_OK:          /* 按下摇杆OK键，串口打函数DSP_MatSub的输出数据 */DSP_MatSub();break;default:/* 其他的键值不处理 */break;}}}
}

 

21.8 实验例程说明（IAR）

配套例子：

V6-216_DSP矩阵运算（加法，减法和逆矩阵）

实验目的：

  1. 学习DSP复数运算（加法，减法和逆矩阵）

实验内容：

  1. 启动一个自动重装软件定时器，每100ms翻转一次LED2。
  2. 按下按键K1，串口打函数DSP_MatInit的输出数据。
  3. 按下按键K2，串口打函数DSP_MatAdd的输出数据。
  4. 按下按键K3，串口打函数DSP_MatInverse的输出数据。
  5. 按下摇杆OK键，串口打函数DSP_MatSub的输出数据。

上电后串口打印的信息：

波特率 115200，数据位 8，奇偶校验位无，停止位 1。

详见本章的3.4 ，4.4，5.4和6.3小节。

程序设计：

  系统栈大小分配：

 

  硬件外设初始化

硬件外设的初始化是在 bsp.c 文件实现：

/*
*********************************************************************************************************
*    函 数 名: bsp_Init
*    功能说明: 初始化所有的硬件设备。该函数配置CPU寄存器和外设的寄存器并初始化一些全局变量。只需要调用一次
*    形    参：无
*    返 回 值: 无
*********************************************************************************************************
*/
void bsp_Init(void)
{/* STM32F407 HAL 库初始化，此时系统用的还是F407自带的16MHz，HSI时钟:- 调用函数HAL_InitTick，初始化滴答时钟中断1ms。- 设置NVIV优先级分组为4。*/HAL_Init();/* 配置系统时钟到168MHz- 切换使用HSE。- 此函数会更新全局变量SystemCoreClock，并重新配置HAL_InitTick。*/SystemClock_Config();/* Event Recorder：- 可用于代码执行时间测量，MDK5.25及其以上版本才支持，IAR不支持。- 默认不开启，如果要使能此选项，务必看V5开发板用户手册第8章*/    
#if Enable_EventRecorder == 1  /* 初始化EventRecorder并开启 */EventRecorderInitialize(EventRecordAll, 1U);EventRecorderStart();
#endifbsp_InitKey();        /* 按键初始化，要放在滴答定时器之前，因为按钮检测是通过滴答定时器扫描 */bsp_InitTimer();      /* 初始化滴答定时器 */bsp_InitUart();    /* 初始化串口 */bsp_InitExtIO();   /* 初始化扩展IO */bsp_InitLed();        /* 初始化LED */        
}

 

  主功能：

主程序实现如下操作：

  •   启动一个自动重装软件定时器，每100ms翻转一次LED2。
  •   按下按键K1，串口打函数DSP_MatInit的输出数据。
  •   按下按键K2，串口打函数DSP_MatAdd的输出数据。
  •   按下按键K3，串口打函数DSP_MatInverse的输出数据。
  •   按下摇杆OK键，串口打函数DSP_MatSub的输出数据。
/*
*********************************************************************************************************
*    函 数 名: main
*    功能说明: c程序入口
*    形    参: 无
*    返 回 值: 错误代码(无需处理)
*********************************************************************************************************
*/
int main(void)
{uint8_t ucKeyCode;        /* 按键代码 */bsp_Init();        /* 硬件初始化 */PrintfLogo();    /* 打印例程信息到串口1 */PrintfHelp();    /* 打印操作提示信息 */bsp_StartAutoTimer(0, 100);    /* 启动1个100ms的自动重装的定时器 *//* 进入主程序循环体 */while (1){bsp_Idle();        /* 这个函数在bsp.c文件。用户可以修改这个函数实现CPU休眠和喂狗 *//* 判断定时器超时时间 */if (bsp_CheckTimer(0))    {/* 每隔100ms 进来一次 */  bsp_LedToggle(2);}ucKeyCode = bsp_GetKey();    /* 读取键值, 无键按下时返回 KEY_NONE = 0 */if (ucKeyCode != KEY_NONE){switch (ucKeyCode){case KEY_DOWN_K1:            /* 按下按键K1，串口打函数DSP_MatInit的输出数据 */DSP_MatInit();break;case KEY_DOWN_K2:            /* 按下按键K2，串口打函数DSP_MatAdd的输出数据 */DSP_MatAdd();break;case KEY_DOWN_K3:            /* 按下按键K3，串口打函数DSP_MatInverse的输出数据 */DSP_MatInverse();break;case JOY_DOWN_OK:          /* 按下摇杆OK键，串口打函数DSP_MatSub的输出数据 */DSP_MatSub();break;default:/* 其他的键值不处理 */break;}}}
}

 

21.9 总结

本期教程就跟大家讲这么多，有兴趣的可以深入研究下算法的具体实现。

 


http://chatgpt.dhexx.cn/article/Ml9ajb4j.shtml

相关文章

【STM32H7的DSP教程】第17章 DSP功能函数-定点数互转

【STM32H7的DSP教程】第17章 DSP功能函数-定点数互转

完整版教程下载地址&#xff1a;http://www.armbbs.cn/forum.php?modviewthread&tid94547 第17章 DSP功能函数-定点数互转 本期教程主要讲解功能函数中的Q7&#xff0c;Q15和Q31分别向其它类型数据转换。 目录 第17章 DSP功能函数-定点数互转 17.1 初学者…
阅读更多...
【STM32F407的DSP教程】第19章 DSP复数运算-共轭，点乘和求模

【STM32F407的DSP教程】第19章 DSP复数运算-共轭，点乘和求模

完整版教程下载地址&#xff1a;http://www.armbbs.cn/forum.php?modviewthread&tid94547 第19章 DSP复数运算-共轭&#xff0c;点乘和求模 本期教程主要讲解复数运算中的共轭&#xff0c;点乘和模的求解。 目录 第19章 DSP复数运算-共轭&#xff0c;点乘…
阅读更多...
【STM32F407的DSP教程】第1章 初学数字信号处理准备工作

【STM32F407的DSP教程】第1章 初学数字信号处理准备工作

完整版教程下载地址&#xff1a;http://www.armbbs.cn/forum.php?modviewthread&tid94547 第1章 初学数字信号处理准备工作 本期教程开始带领大家学习DSP教程&#xff0c;学习前首先要搞明白一个概念&#xff0c;DSP有两层含义&#xff0c;一个是DSP芯片也就是Digital…
阅读更多...
【STM32H7的DSP教程】第11章 DSP基础函数-绝对值，求和，乘法和点乘

【STM32H7的DSP教程】第11章 DSP基础函数-绝对值，求和，乘法和点乘

完整版教程下载地址&#xff1a;http://www.armbbs.cn/forum.php?modviewthread&tid94547 第11章 DSP基础函数-绝对值&#xff0c;求和&#xff0c;乘法和点乘 本期教程开始学习ARM官方的DSP库&#xff0c;这里我们先从基本数学函数开始。本期教程主要讲绝对值&am…
阅读更多...
c++中ifstream和ofstream的详细用法

c++中ifstream和ofstream的详细用法

TITLE&#xff1a;因为不经常用这两个操作&#xff0c;每次用都得找&#xff0c;网上的文章良莠不齐&#xff0c;有时候不容易找到合适的&#xff0c;故将其记录下来&#xff0c;供以后使用方便。 1. fstream fstream提供了三个类&#xff0c;用来实现c对文件的操作。&#x…
阅读更多...
C++文件读取ifstream

C++文件读取ifstream

使用ifstream流来读取文件 说明&#xff1a; 1.ifstream类的对象创建成功的时候会返回非空值&#xff0c;借此判断是否创建文件对象成功 2.ifstream有个函数eof()用来判断文件是否读到尾部,没读到尾部返回false&#xff0c;否则返回true。 若尾部有回车&#xff0c;那么最后一…
阅读更多...
C++中读取文件ifstream、ofstream

C++中读取文件ifstream、ofstream

现在&#xff0c;我想对一个utf-8的json文件进行读写操作。在此过程中&#xff0c;遇到了一些问题。 头文件的添加&#xff1a; #include <iostream> #include <cstring> #include <fstream> #include <sstream> &#xff08;一&#xff09;ifstrea…
阅读更多...
C++ 文件读写操作std::ofstream和std::ifstream

C++ 文件读写操作std::ofstream和std::ifstream

文章目录 前言一、文件操作的介绍二、使用步骤1.读文件操作&#xff08;std::ifstream&#xff09;2.写文件操作&#xff08;std::ofstream&#xff09; 前言 例如&#xff1a;随着人工智能的不断发展&#xff0c;机器学习这门技术也越来越重要&#xff0c;很多人都开启了学习…
阅读更多...
C++ 中 ifstream读取txt文件内容

C++ 中 ifstream读取txt文件内容

文章目录 头文件把txt文件放在当前目录下1、逐行读入文件2、逐个字符读入&#xff08;忽略空格与回车&#xff09;3、逐个字符读入&#xff08;包括空格与回车&#xff09;示例代码 头文件 #include <iostream> #include <fstream> #include <cassert> #inc…
阅读更多...
【C++11】文件操作ifstreamofstream

【C++11】文件操作ifstreamofstream

文章目录 文件输入流创建ifstream对象读取文件数据 文件输出流文件流的打开模式创建ofstream对象写入文件数据 文件输入输出示例 文件输入流 创建ifstream对象 使用ifstream类创建ifstream对象&#xff0c;所获取的对象能够像cin一样使用>>运算符从所绑定文件中读取数据…
阅读更多...
C++ ifstream eof()的使用

C++ ifstream eof()的使用

C ifstream eof() 的使用 eof() 的使用方法1 ifstream.eof() 读到文件结束符时返回true。 大家可能有一个误区&#xff0c;认为读到文件结束符就是读到文件的最后一个字符。 其实不然&#xff0c;文件结束符是文件最后一个字符的下一个字符0xFF&#xff0c;eof() 读到文件结束…
阅读更多...
C++文件读写详解（ofstream,ifstream,fstream）

C++文件读写详解（ofstream,ifstream,fstream）

在看C编程思想中&#xff0c;每个练习基本都是使用ofstream,ifstream,fstream&#xff0c;以前粗略知道其用法和含义&#xff0c;在看了几位大牛的博文后&#xff0c;进行整理和总结&#xff1a; 这里主要是讨论fstream的内容&#xff1a; #include <fstream>ofstream …
阅读更多...
c++输入文件流ifstream用法详解

c++输入文件流ifstream用法详解

目录 文章目录 输入流的继承关系&#xff1a;成员函数Public member functions 1&#xff0c; (constructor) 2&#xff0c;ifstream::open 3&#xff0c;ifstream:: is_open 4&#xff0c;ifstream:: close 5&#xff0c;ifstream:: rdbuf 6&#xff0c;ifstream:: operator …
阅读更多...
ifstream的使用

ifstream的使用

fstream提供了三个类&#xff0c;用来实现c对文件的操作。&#xff08;文件的创建、读、写&#xff09;。 ifstream – 从已有的文件读入 ofstream – 向文件写内容 fstream - 打开文件供读写 文件打开模式&#xff1a; ios::in 只读 ios::out 只写 ios::app 从文件末尾开始写&…
阅读更多...
c++ 输入文件流ifstream用法详解

c++ 输入文件流ifstream用法详解

文章目录 c 输入文件流ifstream用法详解输入流的继承关系&#xff1a;C 使用标准库类来处理面向流的输入和输出&#xff1a;成员函数Public member functions1. **(constructor)**2. **ifstream::open**3. **ifstream:: is_open**4. **ifstream:: close**5. **ifstream:: rdbuf…
阅读更多...
计算机工业控制高职教材,计算机控制技术(21世纪高职高专系列规划教材)

计算机工业控制高职教材,计算机控制技术(21世纪高职高专系列规划教材)

导语 本书以工业控制计算机(IPC)为主线&#xff0c;理论联系实际&#xff0c;突出工程应用&#xff0c;阐述了计算机控制技术及其工程实现方法。全书分为8章&#xff0c;内容包括&#xff1a;计算机控制系统概述&#xff0c;计算机控制过程通道&#xff0c;数字控制技术&#x…
阅读更多...
高职高专信息工程学院专业设置

高职高专信息工程学院专业设置

学院全面落实立德树人根本任务，注重实习实训，着力培养德才兼备的技能型信息技术人才。与华为、腾讯、百度、阿里巴巴、新浪、搜狐、网易等多家知名IT企业合作，建成多所紧贴行业前沿的实习实训基地，保证人才培养与企业需求无缝对接。学院现有计算机应用技术、云计算技术应用…
阅读更多...
湖北省高职计算机本科学校有哪些,盘点最新湖北十大高职高专院校排名，湖北最好的高职院校有哪些？...

湖北省高职计算机本科学校有哪些,盘点最新湖北十大高职高专院校排名，湖北最好的高职院校有哪些？...

高职高专就是高等职业学院和高等专科学校的简称&#xff0c;是专科(大专)层次的普通高等学校。简单点来说&#xff0c;高职高专院校就是职业技术教育&#xff0c;是职业技术教育的高等阶段。今天小编就来给大家盘点下最新湖北十大高职高专院校排名&#xff0c;湖北最好的高职院…
阅读更多...
江西省计算机学会高职高专,我校应邀出席江西省计算机学会高职高专工作委员会成立大会...

江西省计算机学会高职高专,我校应邀出席江西省计算机学会高职高专工作委员会成立大会...

8月21日下午&#xff0c;江西省计算机学会高职高专工作委员会成立大会在南昌召开&#xff0c;中国计算机学会职业教育发展委员会、江西省计算机学会、广东省计算机学会高职高专分会、省内相关高职院校领导和企业代表等90余人参加成立大会。我校作为主任委员单位&#xff0c;副校…
阅读更多...
武汉高职高专计算机专业分数线,武汉高职高专学校有哪些及分数线

武汉高职高专计算机专业分数线,武汉高职高专学校有哪些及分数线

武汉市高职高专众多&#xff0c;其中也包含不少野鸡大学&#xff0c;哪些高职高专是值得2020年高考生选择的正规高校&#xff0c;7月9日教育部官网已公布全国高校名单&#xff0c;其中武汉市有37所公办或民办的高职学校&#xff0c;已整理各高职院校2019年名单及最低录取分数线…
阅读更多...
推荐文章

Copyright @ 2022~2023