特别说明：完整45期数字信号处理教程，原创高性能示波器代码全开源地址：链接
, z7 s2 |3 S6 O9 r, w) T# f

第20章 MatrixFunctions的使用（二）

5 B( C" k; ]: y8 I6 N

    本期教程主要讲解矩阵运算中的放缩，乘法和转置。

    20.1 矩阵放缩 MatScale

    20.2 矩阵乘法 MatMult

    20.3 转置矩阵 MatTrans

    20.4 总结

20.1 矩阵放缩 MatScale

20.1.1 arm_mat_scale_f32

公式描述：

函数定义如下：

    arm_status arm_mat_scale_f32(

        const arm_matrix_instance_f32 * pSrc,

        float32_t scale,

        arm_matrix_instance_f32 * pDst)

参数定义：

    [in]       *pSrc   points to input matrix structure        

    [in]       scale   scale factor to be applied         

    [out]      *pDst  points to output matrix structure        

    return    The function returns either <code>ARM_MATH_SIZE_MISMATCH</code>   

; I, x" h4 }! j* d* G* o20.1.2 arm_mat_scale_q31

函数定义如下：

    arm_status arm_mat_scale_q31(

        const arm_matrix_instance_q31 * pSrc,

        q31_t scaleFract,

        int32_t shift,

        arm_matrix_instance_q31 * pDst)

参数定义：

    [in]      *pSrc      points to input matrix        

    [in]      scaleFract  fractional portion of the scale factor        

    [in]      shift       number of bits to shift the result by        

    [out]    *pDst      points to output matrix structure        

    return     The function returns either   

注意事项：

    1. 两个1.31格式的数据相乘产生2.62格式的数据，最终结果要做偏移和饱和运算产生1.31格式数据。

    2. 定点数的最终放缩比例计算是：scale = scaleFract * 2^shift.

& J  i3 B: \( Q+ _20.1.3 arm_mat_scale_q15

函数定义如下：

    arm_status arm_mat_scale_q15(

        const arm_matrix_instance_q15 * pSrc,

        q15_t scaleFract,

        int32_t shift,

        arm_matrix_instance_q15 * pDst)

参数定义：

    [in,out] *S         points to an instance of the floating-point matrix structure.   

    [in]     nRows    number of rows in the matrix.   

    [in]     nColumns number of columns in the matrix.   

    [in]     *pData           points to the matrix data array.  

注意事项：

    1. 两个1.15格式的数据相乘产生2.30格式的数据，最终结果要做偏移和饱和运算产生1.15格式数据。

    2. 定点数的最终放缩比例计算是：scale = scaleFract * 2^shift.

20.1.4 实例讲解

实验目的：

    1. 学习MatrixFunctions中矩阵的放缩

实验内容：

    1. 按下按键K1, 串口打印函数DSP_MatScale的输出结果

实验现象：

    通过窗口上位机软件SecureCRT（V5光盘里面有此软件）查看打印信息现象如下：

程序设计：

1. /*
   
2. *********************************************************************************************************
     F- _' N" x) A  D7 K5 k9 g
3. *        函 数 名: DSP_MatScale
   
4. *        功能说明: 矩阵放缩
   
5. *        形    参：无
   
6. *        返 回 值: 无
   ; v, F8 E5 B- A/ t
7. *********************************************************************************************************
   
8. */
   
9. static void DSP_MatScale(void)
   0 {, w3 F/ @: \9 |' F
10. {
    * l0 x7 Z8 d$ o- q, F7 E$ P2 Z
11. uint8_t i;
    " q8 a1 R' Y- h7 j9 {7 i
12. /****浮点数数组******************************************************************/
    " o- \  j) s& c- z- r  i
13. float32_t pDataA[9] = {1.1f, 1.1f, 2.1f, 2.1f, 3.1f, 3.1f, 4.1f, 4.1f, 5.1f};
    
14. float32_t scale = 1.1f;
    
15. float32_t pDataDst[9];
    
16. arm_matrix_instance_f32 pSrcA; //3行3列数据
      x7 F9 E4 c, H2 K
17. arm_matrix_instance_f32 pDst;
    
18. /****定点数Q31数组******************************************************************/
    
19. q31_t pDataA1[9] = {1, 1, 2, 2, 3, 3, 4, 4, 5};
    
20. q31_t scaleFract = 10;
    " F2 F6 v0 i& |* V1 ~/ v
21.      int32_t shift = 0;
    " l* e" z4 h( s4 \7 X* H1 g
22. q31_t pDataDst1[9];
    " l0 q8 R* O5 s( ^2 V
23. arm_matrix_instance_q31 pSrcA1; //3行3列数据
    
24. arm_matrix_instance_q31 pDst1;
    
25. /****定点数Q15数组******************************************************************/
    ! R: Q+ _; @2 g3 ]) ?9 G/ v0 b1 G# n
26. q15_t pDataA2[9] = {1, 1, 2, 2, 3, 3, 4, 4, 5};
    ! y5 y+ @$ f' B  w4 }5 J
27. q15_t scaleFract1 = 10;
    : J( Q2 e. F9 S- s( F
28.      int32_t shift1 = 0;
    
29. q15_t pDataDst2[9];
    8 f5 @: E# C9 B
30. arm_matrix_instance_q15 pSrcA2; //3行3列数据
    
31. arm_matrix_instance_q15 pDst2;
    
32. /****浮点数***********************************************************************/
    
33. pSrcA.numCols = 3;
    7 `( G! R) ~6 i. X7 G
34. pSrcA.numRows = 3;
    + u; |& F  r1 P0 F7 N! U% a; B+ j
35. pSrcA.pData = pDataA;
      f$ y9 T6 c8 M
36. 
    
37. pDst.numCols = 3;
    # Q1 u# x$ }2 }( W1 i; |" [* e
38. pDst.numRows = 3;
    
39. pDst.pData = pDataDst;
    ( @$ ^1 n, |" D: g7 z
40. printf("****浮点数******************************************\r\n");
    
41. arm_mat_scale_f32(&pSrcA, scale, &pDst);
    
42. for(i = 0; i < 9; i++)
    
43. {
    ) c" O# a/ ~% p9 v8 u+ R7 D6 U. n
44. printf("pDataDst[%d] = %f\r\n", i, pDataDst[i]);
    
45. }
    : [( [8 K* y$ s
46. /****定点数Q31***********************************************************************/
    
47. pSrcA1.numCols = 3;
    
48. pSrcA1.numRows = 3;
    
49. pSrcA1.pData = pDataA1;
    
50. pDst1.numCols = 3;
    2 A: I: f/ U; n$ @& m
51. pDst1.numRows = 3;
    
52. pDst1.pData = pDataDst1;
    4 m% E; b- s8 l6 F( S! c' S. E$ Y
53. printf("****定点数Q31******************************************\r\n");
    
54. arm_mat_scale_q31(&pSrcA1, scaleFract, shift, &pDst1);
    
55. for(i = 0; i < 9; i++)
    ! a( S6 \( X- |
56. {
    
57. printf("pDataDst1[%d] = %d\r\n", i, pDataDst1[i]);
    
58. }
    ) `6 I8 I4 b/ `+ {# m! X! A/ U
59. /****定点数Q15***********************************************************************/
    
60. pSrcA2.numCols = 3;
    : b0 A2 L! v8 Z
61. pSrcA2.numRows = 3;
    
62. pSrcA2.pData = pDataA2;
    
63. pDst2.numCols = 3;
    
64. pDst2.numRows = 3;
    9 U! Z* o- l0 s/ n( g
65. pDst2.pData = pDataDst2;
    
66. printf("****定点数Q15******************************************\r\n");
    2 D0 m' E$ D' d- U1 m' v
67. arm_mat_scale_q15(&pSrcA2, scaleFract1, shift1, &pDst2);
    
68. for(i = 0; i < 9; i++)
    
69. {
    6 t# O6 j# Q- U7 z3 c; R
70. printf("pDataDst2[%d] = %d\r\n", i, pDataDst2[i]);
    2 ^8 C3 T4 ~4 F7 }6 j6 x
71. }
    7 q- y. |0 s& ?8 t- D. e
72. }

复制代码

1. 下面通过matlab来实现矩阵的放缩：

8 _4 ^. p) x9 G; h- I
, e& E' v- [! j% P$ y! p% o5 g

20.2 矩阵乘法 MatMult
20.2.1 arm_mat_mult_f32

公式描述：

函数定义如下：

    arm_status arm_mat_mult_f32(

        const arm_matrix_instance_f32 * pSrcA,

        const arm_matrix_instance_f32 * pSrcB,

        arm_matrix_instance_f32 * pDst)

参数定义：

    [in]       *pSrcA  points to the first input matrix structure   

    [in]       *pSrcB  points to the second input matrix structure   

    [out]      *pDst  points to output matrix structure   

    return     The function returns either   

注意事项：

    1. 两个矩阵M x N和N x P相乘的结果是M x P.（必须保证一个矩形的列数等于另一个矩阵的行数）。

0 K$ q+ \' n, I5 O" m7 i, r- a20.2.2 arm_mat_mult_q31

函数定义如下：

    arm_status arm_mat_mult_q31(

        const arm_matrix_instance_q31 * pSrcA,

        const arm_matrix_instance_q31 * pSrcB,

        arm_matrix_instance_q31 * pDst)

参数定义：

    [in]    *pSrcA  points to the first input matrix structure   

    [in]    *pSrcB  points to the second input matrix structure   

    [out]  *pDst   points to output matrix structure   

    return                     The function returns either

注意事项：

    1. 两个1.31格式的数据相乘产生2.62格式的数据，最终结果要做偏移和饱和运算产生1.31格式数据。

    2. 两个矩阵M x N和N x P相乘的结果是M x P.（必须保证一个矩形的列数等于另一个矩阵的行数）。

8 P; P% H! R+ ^* l20.2.3 arm_mat_mult_q15

函数定义如下：

    arm_status arm_mat_mult_q15(

        const arm_matrix_instance_q15 * pSrcA,

        const arm_matrix_instance_q15 * pSrcB,

       arm_matrix_instance_q15 * pDst,

        q15_t * pState CMSIS_UNUSED)

参数定义：

    [in]     *pSrcA   points to the first input matrix structure   

    [in]     *pSrcB   points to the second input matrix structure   

    [out]   *pDst    points to output matrix structure   

    [in]                *pState  points to the array for storing intermediate results   

    return                      The function returns either   

注意事项：

    1. 两个1.15格式数据相乘是2.30格式，函数的内部使用了64位的累加器，那个就是34.30格式，最终结果将低15位截取掉并作饱和处理为1.15格式。

    2. 两个矩阵M x N和N x P相乘的结果是M x P.（必须保证一个矩形的列数等于另一个矩阵的行数）。

" j* d* _7 j, T8 b, n2 Q5 S/ h20.2.4 arm_mat_mult_fast_q31

函数定义如下：

    arm_status arm_mat_mult_fast_q31(

        const arm_matrix_instance_q31 * pSrcA,

        const arm_matrix_instance_q31 * pSrcB,

        arm_matrix_instance_q31 * pDst)

参数定义：

    [in]    *pSrcA  points to the first input matrix structure   

    [in]    *pSrcB  points to the second input matrix structure   

    [out]  *pDst   points to output matrix structure   

    return                     The function returns either

注意事项：

    1. 两个1.31格式的数据相乘产生2.62格式的数据，最终结果要做偏移和饱和运算产生1.31格式数据。

    2. 两个矩阵M x N和N x P相乘的结果是M x P.（必须保证一个矩形的列数等于另一个矩阵的行数）。

    3. 函数arm_mat_mult_fast_q31是arm_mat_mult_q31的快速算法。

20.2.5 arm_mat_mult_fast_q15

函数定义如下：

     arm_status arm_mat_mult_fast_q15(

        const arm_matrix_instance_q15 * pSrcA,

        const arm_matrix_instance_q15 * pSrcB,

        arm_matrix_instance_q15 * pDst,

        q15_t * pState)

参数定义：

    [in]     *pSrcA   points to the first input matrix structure   

    [in]     *pSrcB   points to the second input matrix structure   

    [out]   *pDst    points to output matrix structure   

    [in]                *pState  points to the array for storing intermediate results   

    return                      The function returns either   

注意事项：

    1. 两个1.15格式数据相乘是2.30格式，函数的内部使用了64位的累加器，那个就是34.30格式，最终结果将低15位截取掉并作饱和处理为1.15格式。

    2. 两个矩阵M x N和N x P相乘的结果是M x P.（必须保证一个矩形的列数等于另一个矩阵的行数）。

    3. 函数arm_mat_mult_fast_q15是arm_mat_mult_q15的快速算法。

, \2 i* G9 H" P0 |6 k" Y9 e. y20.2.6 实例讲解

实验目的：

    1. 学习MatrixFunctions中矩阵乘法

实验内容：

    1. 按下按键K2, 串口打印函数DSP_MatMult的输出结果

实验现象：

     通过窗口上位机软件SecureCRT（V5光盘里面有此软件）查看打印信息现象如下：

程序设计：

1. /*
   
2. *********************************************************************************************************
   / r6 e! u) A8 D. s2 r* {6 A8 I
3. *        函 数 名: DSP_MatMult
   
4. *        功能说明: 矩阵乘法
   
5. *        形    参：无
   
6. *        返 回 值: 无
   
7. *********************************************************************************************************
   ( k  X+ A( w2 ?6 I' i( A& T
8. */
   - ~8 j8 D! Z) [/ X; v
9. static void DSP_MatMult(void)
   
10. {
    $ J- a* I. z: O7 o
11. uint8_t i;
    
12. /****浮点数数组******************************************************************/
    
13. float32_t pDataA[9] = {1.1f, 1.1f, 2.1f, 2.1f, 3.1f, 3.1f, 4.1f, 4.1f, 5.1f};
    
14. float32_t pDataB[9] = {1.1f, 1.1f, 2.1f, 2.1f, 3.1f, 3.1f, 4.1f, 4.1f, 5.1f};
    
15. float32_t pDataDst[9];
    
16. arm_matrix_instance_f32 pSrcA; //3行3列数据
    
17. arm_matrix_instance_f32 pSrcB; //3行3列数据
    
18. arm_matrix_instance_f32 pDst;
    
19. /****定点数Q31数组******************************************************************/
    
20. q31_t pDataA1[9] = {1, 1, 2, 2, 3, 3, 4, 4, 5};
    & a5 a6 Q3 x1 @9 m1 C2 z
21. q31_t pDataB1[9] = {1, 1, 2, 2, 3, 3, 4, 4, 5};
    ' P4 G/ C7 [! Z
22. q31_t pDataDst1[9];
    
23. arm_matrix_instance_q31 pSrcA1; //3行3列数据
    1 A4 \! o) r* H, k  [% L
24. arm_matrix_instance_q31 pSrcB1; //3行3列数据
    
25. arm_matrix_instance_q31 pDst1;
    
26. /****定点数Q15数组******************************************************************/
    ) Z/ S" x2 Z1 v, [. l4 i9 @
27. q15_t pDataA2[9] = {1, 1, 2, 2, 3, 3, 4, 4, 5};
    
28. q15_t pDataB2[9] = {1, 1, 2, 2, 3, 3, 4, 4, 5};
    & H' V9 a4 B" X5 |8 s0 p' @  d( _
29. q15_t pDataDst2[9];
    
30. arm_matrix_instance_q15 pSrcA2; //3行3列数据
    
31. arm_matrix_instance_q15 pSrcB2; //3行3列数据
    ; u) I% O! _( c  ]+ `! q
32. arm_matrix_instance_q15 pDst2;
    " @. h% I3 S5 m  F5 n0 u8 M) s
33. q15_t pState;
    
34. /****浮点数***********************************************************************/
    
35. pSrcA.numCols = 3;
    
36. pSrcA.numRows = 3;
    7 q* z0 x& f) z7 z9 h2 u: }
37. pSrcA.pData = pDataA;
    
38. pSrcB.numCols = 3;
    0 j, M8 l$ `' I2 |; G4 `' \
39. pSrcB.numRows = 3;
    
40. pSrcB.pData = pDataB;
      f6 [/ M# v9 f4 p- ~
41. pDst.numCols = 3;
    
42. pDst.numRows = 3;
    
43. pDst.pData = pDataDst;
    
44. printf("****浮点数******************************************\r\n");
    
45. arm_mat_mult_f32(&pSrcA, &pSrcB, &pDst);
    2 `: q5 T6 ~2 K% W( G! ]  I: [
46. for(i = 0; i < 9; i++)
    
47. {
    
48. printf("pDataDst[%d] = %f\r\n", i, pDataDst[i]);
    
49. }
    9 N9 h% R6 Y6 v% @. Z3 \
50. /****定点数Q31***********************************************************************/
    
51. pSrcA1.numCols = 3;
    % ]! n  ?4 g+ m$ p
52. pSrcA1.numRows = 3;
    
53. pSrcA1.pData = pDataA1;
    
54. pSrcB1.numCols = 3;
    ( D8 T" B) x: g, K
55. pSrcB1.numRows = 3;
    4 h) C# `, S# s' u0 l. s, H
56. pSrcB1.pData = pDataB1;
    
57. pDst1.numCols = 3;
    " Q: P. K' L! y
58. pDst1.numRows = 3;
    
59. pDst1.pData = pDataDst1;
    
60. printf("****定点数Q31******************************************\r\n");
    
61. arm_mat_mult_q31(&pSrcA1, &pSrcB1, &pDst1);
    
62. arm_mat_mult_fast_q31(&pSrcA1, &pSrcB1, &pDst1);
    + P( `9 J2 _9 N2 W
63. for(i = 0; i < 9; i++)
    & b( [  K( {$ o
64. {
    
65. printf("pDataDst1[%d] = %d\r\n", i, pDataDst1[i]);
    ! r9 \( _9 b1 h; o% ~; o/ w" @
66. }
    ! v$ z3 e0 E' A* B! A+ e
67. /****定点数Q15***********************************************************************/
    # ?0 F, ^3 Y; D# U+ W
68. pSrcA2.numCols = 3;
    . o6 T( p0 x! n+ O8 h- v6 j
69. pSrcA2.numRows = 3;
    " a1 z& c9 a5 |0 X( B
70. pSrcA2.pData = pDataA2;
    
71. pSrcB2.numCols = 3;
    2 ~$ S4 F2 _6 V6 O, A
72. pSrcB2.numRows = 3;
    
73. pSrcB2.pData = pDataB2;
    
74. pDst2.numCols = 3;
    
75. pDst2.numRows = 3;
    # {3 h4 o3 `! Z  m/ B* U$ K% J1 \
76. pDst2.pData = pDataDst2;
    
77. printf("****定点数Q15******************************************\r\n");
    
78. arm_mat_mult_q15(&pSrcA2, &pSrcB2, &pDst2, &pState);
    
79. arm_mat_mult_fast_q15(&pSrcA2, &pSrcB2, &pDst2, &pState);
    
80. for(i = 0; i < 9; i++)
    2 T: K2 y- K1 x
81. {
    
82. printf("pDataDst2[%d] = %d\r\n", i, pDataDst2[i]);
    
83. }
    ; r1 C1 {$ M5 f) k4 n5 ?) h
84. }

复制代码

1. 下面通过matlab实现矩阵的乘法：

0 |. u2 Q4 E' |1 z

20.3 转置矩阵 MatTrans
20.3.1 arm_mat_trans_f32

公式描述：

函数定义如下：

    arm_status arm_mat_trans_f32(

        const arm_matrix_instance_f32 * pSrc,

        arm_matrix_instance_f32 * pDst)

参数定义：

    [in]   *pSrc  points to the input matrix   

    [out]  *pDst  points to the output matrix   

    return         The function returns either  <code>ARM_MATH_SIZE_MISMATCH</code>  

注意事项：

    1. 矩阵M x N转置后是N x M。

20.3.2 arm_mat_trans_q31

函数定义如下：

    arm_status arm_mat_trans_q31(

        const arm_matrix_instance_q31 * pSrc,

        arm_matrix_instance_q31 * pDst)

参数定义：

    [in]  *pSrc points to the input matrix   

     [out] *pDst points to the output matrix   

     return         The function returns either  <code>ARM_MATH_SIZE_MISMATCH</code>  

注意事项：

    1. 矩阵M x N转置后是N x M。

20.3.3 arm_mat_trans_q15

函数定义如下：

    arm_status arm_mat_trans_q15(

        const arm_matrix_instance_q15 * pSrc,

        arm_matrix_instance_q15 * pDst)

参数定义：

    [in]  *pSrc points to the input matrix   

     [out] *pDst points to the output matrix   

     return         The function returns either  <code>ARM_MATH_SIZE_MISMATCH</code>  

注意事项：

    1. 矩阵M x N转置后是N x M。

20.3.4 实例讲解

实验目的：

    1. 学习MatrixFunctions中的转置矩阵

实验内容：

    1. 按下按键K3, 串口打印函数DSP_MatTrans的输出结果

实验现象：

    通过窗口上位机软件SecureCRT（V5光盘里面有此软件）查看打印信息现象如下：

程序设计：

1. /*
   3 y$ Y, r: K! ?& H" }; q# `* x
2. *********************************************************************************************************
   
3. *        函 数 名: DSP_MatTrans
   
4. *        功能说明: 求逆矩阵
   0 i! A% g2 `+ M" z( y% ]; p" i& [
5. *        形    参：无
   
6. *        返 回 值: 无
   $ E! }" a/ S, z( o9 n0 g
7. *********************************************************************************************************
   2 n% g/ s7 R6 q, A/ q$ M0 F. A
8. */
   7 i5 s- l6 b7 G8 _# a1 e+ Y& O& Q
9. static void DSP_MatTrans(void)
   2 r/ s9 c; r0 o3 G: U9 o9 J
10. {
    2 s- k3 P- ?6 b" H/ X+ c
11. uint8_t i;
    
12. /****浮点数数组******************************************************************/
    : R0 _  l( P( N
13. float32_t pDataA[9] = {1.1f, 1.1f, 2.1f, 2.1f, 3.1f, 3.1f, 4.1f, 4.1f, 5.1f};
    & G. A- ~. p/ Q+ U, T8 z* A  S) z
14. float32_t pDataDst[9];
    ) `% ~5 i  r. F! i, S- @
15. arm_matrix_instance_f32 pSrcA; //3行3列数据
    
16. arm_matrix_instance_f32 pDst;
    / b6 O7 A4 y" L: j1 S
17. /****定点数Q31数组******************************************************************/
    / q& j+ ~$ e- `) o
18. q31_t pDataA1[9] = {1, 1, 2, 2, 3, 3, 4, 4, 5};
    5 l  I2 U6 s! M  n
19. q31_t pDataDst1[9];
    ! s' H  i- S9 o! W' g2 w
20. arm_matrix_instance_q31 pSrcA1; //3行3列数据
    
21. arm_matrix_instance_q31 pDst1;
    
22. /****定点数Q15数组******************************************************************/
    
23. q15_t pDataA2[9] = {1, 1, 2, 2, 3, 3, 4, 4, 5};
    
24. q15_t pDataDst2[9];
    
25. arm_matrix_instance_q15 pSrcA2; //3行3列数据
    
26. arm_matrix_instance_q15 pDst2;
    
27. /****浮点数***********************************************************************/
    ' n" x3 J% u- C
28. pSrcA.numCols = 3;
    5 J7 k# d8 I- v. u, p
29. pSrcA.numRows = 3;
    
30. pSrcA.pData = pDataA;
    % N, F, U4 b9 b3 y3 B0 v
31. 
    
32. pDst.numCols = 3;
    / F) h+ w7 |: d' x8 \5 ?
33. pDst.numRows = 3;
    4 L1 W& O) H% O7 m/ V
34. pDst.pData = pDataDst;
    % d% p9 C3 K0 B' d
35. printf("****浮点数******************************************\r\n");
    
36. status = arm_mat_trans_f32(&pSrcA, &pDst);
    0 x3 ^& {, V5 J# N# e
37. for(i = 0; i < 9; i++)
    
38. {
    
39. printf("pDataDst[%d] = %f\r\n", i, pDataDst[i]);
    7 {& t4 }* U2 E
40. }
    8 [6 z1 \# T& N* e. o  R2 @
41. /****定点数Q31***********************************************************************/
    
42. pSrcA1.numCols = 3;
    
43. pSrcA1.numRows = 3;
    , w; U# M- e. A  W2 O  p
44. pSrcA1.pData = pDataA1;
    
45. pDst1.numCols = 3;
    
46. pDst1.numRows = 3;
    , v" j$ q6 c: Q. i
47. pDst1.pData = pDataDst1;
    " c4 u3 Z( x3 w- p( W" f1 w
48. printf("****定点数Q31******************************************\r\n");
    2 T5 Y2 e& a0 @( T, Y6 Q# F) {
49. status = arm_mat_trans_q31(&pSrcA1, &pDst1);
    
50. for(i = 0; i < 9; i++)
      R  a0 t/ e5 ~3 o; U
51. {
    
52. printf("pDataDst1[%d] = %d\r\n", i, pDataDst1[i]);
    
53. }
    
54. /****定点数Q15***********************************************************************/
    " D8 {0 i: A9 e( ^7 p# u
55. pSrcA2.numCols = 3;
    7 F7 W. R' |: c
56. pSrcA2.numRows = 3;
    
57. pSrcA2.pData = pDataA2;
    
58. pDst2.numCols = 3;
    
59. pDst2.numRows = 3;
    . T/ c0 x4 D# G8 k0 v: U
60. pDst2.pData = pDataDst2;
    ( i: X' i6 f0 F( T
61. printf("****定点数Q15******************************************\r\n");
    4 }* q, d# }) h8 l' \* ^
62. status = arm_mat_trans_q15(&pSrcA2, &pDst2);
    
63. for(i = 0; i < 9; i++)
    
64. {
    
65. printf("pDataDst2[%d] = %d\r\n", i, pDataDst2[i]);
    $ t0 ]- T0 E$ g7 P) {, O% y) i
66. }
    
67. }

复制代码

1. 下面通过matlab实现矩阵的转置：

20.4 总结

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

" U: l+ E/ f: {1 q" Z

谢谢分享啊