21ic电子技术开发论坛 › 单片机与嵌入式系统 › 新唐MCU › M4的DSP内核数学计算
返回列表 发新帖我要提问

[DemoCode下载] M4的DSP内核数学计算

[复制链接]
1086|4
 楼主| xinpian101 发表于 2018-3-8 11:03 | 显示全部楼层 |阅读模式
DSP, 内核, 向量, define, se
展示使用DSP 基本数**算，其中包含
1. 向量绝对值
2. 向量加法
3. 向量减法
4. 向量乘法
5. 向量内积
6. 向量反向
7. 向量缩放
8. 向量偏移
用户可以直接使用这些函式，来实现自己的数学方程式运算。程序内也比较了有无使用DSP计算时间的差异，使用者把#define USE_DSP给取消即可不使用DSP来运算这些方程式。

回复

收藏 举报

相关帖子
 楼主| xinpian101 发表于 2018-3-8 11:03 | 显示全部楼层
复制
  1. /****************************************************************************
    

  2.  * [url=home.php?mod=space&uid=288409]@file[/url]     main.c
    

  3.  * [url=home.php?mod=space&uid=895143]@version[/url]  V1.0
    

  4.  * $Date: 15/09/02 10:04a $
    

  5.  * @brief
    

  6.  *           Display how to use DSP Basic math functions.
    

  7.  *                      and compare with calcultion without DSP
    

  8.  * @note
    

  9.  * Copyright (C) 2016 Nuvoton Technology Corp. All rights reserved.
    

  10.  *
    

  11.  ******************************************************************************/
    

  12. #include <stdio.h>
    

  13. #include <stdlib.h>
    

  14. #include "M451Series.h"
    

  15. #include "arm_math.h" 
    

  16. #define blockSize 32
    

  17. #define USE_DSP
    


  18. 

  19. float32_t srcA_buf_f32[blockSize] =   
    

  20. {   
    

  21. -0.4325648115282207,        -1.6655843782380970,        0.1253323064748307,         
    

  22.  0.2876764203585489,    -1.1464713506814637,        1.1909154656429988,         
    

  23.  1.1891642016521031,    -0.0376332765933176,        0.3272923614086541,         
    

  24.  0.1746391428209245,    -0.1867085776814394,        0.7257905482933027,         
    

  25. -0.5883165430141887,         2.1831858181971011,   -0.1363958830865957,         
    

  26.  0.1139313135208096,         1.0667682113591888,        0.0592814605236053,         
    

  27. -0.0956484054836690,        -0.8323494636500225,        0.2944108163926404,         
    

  28. -1.3361818579378040,         0.7143245518189522,        1.6235620644462707,         
    

  29. -0.6917757017022868,         0.8579966728282626,        1.2540014216025324,         
    

  30. -1.5937295764474768,        -1.4409644319010200,        0.5711476236581780,         
    

  31. -0.3998855777153632,         0.6899973754643451 
    

  32. };   
    

  33.   
    

  34. /* ----------------------------------------------------------------------  
    

  35. ** Test input data of srcB for blockSize 32   
    

  36. ** ------------------------------------------------------------------- */   
    

  37. float32_t srcB_buf_f32[blockSize] =   
    

  38. {   
    

  39.  1.7491401329284098,        0.1325982188803279,         0.3252281811989881,         
    

  40. -0.7938091410349637,        0.3149236145048914,        -0.5272704888029532,         
    

  41.  0.9322666565031119,        1.1646643544607362,        -2.0456694357357357,         
    

  42. -0.6443728590041911,        1.7410657940825480,         0.4867684246821860,         
    

  43.  1.0488288293660140,        1.4885752747099299,         1.2705014969484090,         
    

  44. -1.8561241921210170,        2.1343209047321410,  1.4358467535865909,         
    

  45. -0.9173023332875400,   -1.1060770780029008,         0.8105708062681296,         
    

  46.  0.6985430696369063,   -0.4015827425012831,         1.2687512030669628,         
    

  47. -0.7836083053674872,        0.2132664971465569,         0.7878984786088954,         
    

  48.  0.8966819356782295,   -0.1869172943544062,         1.0131816724341454,         
    

  49.  0.2484350696132857,        0.0596083377937976 
    

  50. };
    

  51. float32_t dotoutput_f32,suboutput[blockSize],multoutput[blockSize],absoutput[blockSize];
    

  52. float32_t scaleoutput[blockSize],offsetoutput[blockSize],addoutput[blockSize],negoutput[blockSize];
    

  53. int CalTime,i;
    

  54. float sum,subout[blockSize],multout[blockSize],absout[blockSize];
    

  55. float scaleout[blockSize],offsetout[blockSize],addout[blockSize],negout[blockSize];
    

  56. void DOT(int size)
    

  57. {
    

  58.         sum=0;
    

  59.         for(i=0;i<size;i++)
    

  60.         {
    

  61.                 sum=sum+(srcA_buf_f32[i]*srcB_buf_f32[i]);
    

  62.         }
    

  63. }
    

  64. void SYS_Init(void)
    

  65. {
    

  66.    /*---------------------------------------------------------------------------------------------------------*/
    

  67.     /* Init System Clock                                                                                       */
    

  68.     /*---------------------------------------------------------------------------------------------------------*/
    

  69.     /* Enable HIRC clock */
    

  70.     CLK_EnableXtalRC(CLK_PWRCTL_HIRCEN_Msk);
    


  71. 

  72.     /* Waiting for HIRC clock ready */
    

  73.     CLK_WaitClockReady(CLK_STATUS_HIRCSTB_Msk);
    


  74. 

  75.     /* Switch HCLK clock source to HIRC */
    

  76.     CLK_SetHCLK(CLK_CLKSEL0_HCLKSEL_HIRC, CLK_CLKDIV0_HCLK(1));
    


  77. 

  78.     /* Enable HXT */
    

  79.     CLK_EnableXtalRC(CLK_PWRCTL_HXTEN_Msk);
    


  80. 

  81.     /* Waiting for clock ready */
    

  82.     CLK_WaitClockReady(CLK_STATUS_HXTSTB_Msk);
    


  83. 

  84.     /* Set core clock as PLL_CLOCK from PLL and SysTick source to HCLK/2*/
    

  85.     CLK_SetCoreClock(72000000);
    

  86.     CLK_SetSysTickClockSrc(CLK_CLKSEL0_STCLKSEL_HCLK_DIV2);
    


  87. 

  88.     /* Enable peripheral clock */
    

  89.     CLK_EnableModuleClock(UART0_MODULE);
    

  90.     CLK_EnableModuleClock(TMR0_MODULE);
    


  91. 

  92.     /* Peripheral clock source */
    

  93.     CLK_SetModuleClock(UART0_MODULE, CLK_CLKSEL1_TMR0SEL_HXT, CLK_CLKDIV0_UART(1));
    

  94.     CLK_SetModuleClock(TMR0_MODULE, CLK_CLKSEL1_TMR0SEL_HXT, 0);
    


  95. 

  96. }
    


  97. 

  98. void UART_Init(void)
    

  99. {
    

  100.     /* Set PD multi-function pins for UART0 RXD and TXD */
    

  101.     SYS->GPD_MFPL = SYS_GPD_MFPL_PD0MFP_UART0_RXD | SYS_GPD_MFPL_PD1MFP_UART0_TXD;
    

  102.         
    

  103.     /* Reset UART module */
    

  104.     SYS_ResetModule(UART0_RST);
    


  105. 

  106.     /* Configure UART0 and set UART0 Baudrate */
    

  107.     UART_Open(UART0, 115200);
    

  108. }
    

  109. /*---------------------------------------------------------------------------------------------------------*/
    

  110. /*  Main Function                                                                                          */
    

  111. /*---------------------------------------------------------------------------------------------------------*/
    

  112. int32_t main(void)
    

  113. {
    

  114.                 /* Unlock protected registers */
    

  115.     SYS_UnlockReg();
    

  116.     /* Init System, peripheral clock and multi-function I/O */
    

  117.     SYS_Init();
    

  118.     /* Lock protected registers */
    

  119.     SYS_LockReg();
    

  120.     /* Init UART for printf */
    

  121.     UART_Init();
    


  122. 

  123.            TIMER_Open(TIMER0, TIMER_CONTINUOUS_MODE, 1);
    

  124.                 TIMER_Start(TIMER0);
    


  125. 

  126. #ifdef USE_DSP        
    

  127.                 /*Calculate dot product (32 sample) with DSP*/
    

  128.                 arm_dot_prod_f32(srcA_buf_f32, srcB_buf_f32, blockSize, &dotoutput_f32);
    

  129.                 
    

  130.                 /*Calculate Absolute value (32 sample) with DSP*/
    

  131.                 arm_abs_f32(srcA_buf_f32, absoutput, blockSize );
    

  132.                 
    

  133.                 /*Calculate Addition (32 sample) with DSP*/
    

  134.                 arm_add_f32(srcA_buf_f32, srcB_buf_f32, addoutput, blockSize);
    

  135.                 
    

  136.                 /*Calculate Subtraction (32 sample) with DSP*/
    

  137.                 arm_sub_f32 (srcA_buf_f32, srcB_buf_f32, suboutput, blockSize);
    

  138.                 
    

  139.                 /*Calculate Multiplication (32 sample) with DSP*/
    

  140.                 arm_mult_f32 (srcA_buf_f32, srcB_buf_f32, multoutput, blockSize);
    

  141.                 
    

  142.                 /*Calculate Scale (32 sample) with DSP*/
    

  143.                 arm_scale_f32 (srcA_buf_f32, 10, scaleoutput, blockSize);        
    

  144.                 
    

  145.                 /*Calculate Offest (32 sample) with DSP*/
    

  146.                 arm_offset_f32 (srcA_buf_f32, 10, offsetoutput, blockSize);
    

  147.                 
    

  148.                 /*Calculate Negate (32 sample) with DSP*/
    

  149.                 arm_negate_f32 (srcA_buf_f32, negoutput, blockSize);
    


  150. 

  151. #else 
    

  152.                 /*Calculate Dot Product (32 sample) without DSP*/                
    

  153.                 DOT(blockSize);
    


  154. 

  155.                 /*Calculate Absolute value (32 sample) without DSP*/
    

  156.                 for(i=0;i<blockSize;i++){
    

  157.                 absout[i]=fabs(srcA_buf_f32[i]);
    

  158.                 }                
    

  159.                 /*Calculate Addition (32 sample) without DSP*/                
    

  160.                 for(i=0;i<blockSize;i++){
    

  161.                 addout[i]=srcA_buf_f32[i]+srcB_buf_f32[i];
    

  162.                 }
    

  163.                 /*Calculate Subtraction (32 sample) without DSP*/        
    

  164.                 for(i=0;i<blockSize;i++){
    

  165.                 subout[i]=srcA_buf_f32[i]-srcB_buf_f32[i];
    

  166.                 }
    

  167.                 /*Calculate Multiplication (32 sample) without DSP*/        
    

  168.                 for(i=0;i<blockSize;i++){
    

  169.                 multout[i]=srcA_buf_f32[i]*srcB_buf_f32[i];
    

  170.                 }
    

  171.                 /*Calculate Scale (32 sample) without DSP*/
    

  172.                 for(i=0;i<blockSize;i++){
    

  173.                 scaleout[i]=10*srcA_buf_f32[i];
    

  174.                 }                
    

  175.                 /*Calculate Offest (32 sample) without DSP*/
    

  176.                 for(i=0;i<blockSize;i++){
    

  177.                 offsetout[i]=10+srcA_buf_f32[i];
    

  178.                 }
    

  179.                 /*Calculate Negate (32 sample) without DSP*/
    

  180.                 for(i=0;i<blockSize;i++){
    

  181.                 negout[i]=-srcA_buf_f32[i];
    

  182.                 }
    

  183.                 
    

  184. #endif                
    

  185.                 TIMER_Close(TIMER0);
    

  186.                 CalTime=TIMER_GetCounter(TIMER0);
    

  187.                 /* Display function calculation time */
    

  188.           printf("time is %d \n",CalTime);                
    

  189.     while(1);   
    


  190. 

  191. }
回复

收藏 举报

dongnanxibei 发表于 2018-3-8 15:22 | 显示全部楼层
噢，用了两组对比，非常赞。
回复

收藏 举报

doit888 发表于 2018-3-9 13:25 | 显示全部楼层
请问一下 这个历程 是在那里 找到的？
回复

收藏 举报

734774645 发表于 2018-3-10 11:21 | 显示全部楼层
doit888 发表于 2018-3-9 13:25
请问一下 这个历程 是在那里 找到的？

肯定是这个M451的BSP里面。
回复

收藏 举报

返回列表 发新帖
高级模式
B Color Image Link Quote Code Smilies
您需要登录后才可以回帖 登录 | 注册

本版积分规则

xinpian101

129

主题

1650

帖子

1

粉丝
快速回复 在线客服 返回列表 返回顶部