你的浏览器版本过低,可能导致网站不能正常访问!
为了你能正常使用网站功能,请使用这些浏览器。

【经验分享】STM32H7的样条插补实现,波形拟合丝滑顺畅

[复制链接]
STMCU小助手 发布时间:2022-1-1 21:00
50.1 初学者重要提示
DSP库支持了样条插补,双线性插补和线性插补,我们这里主要介绍样条插补的实现。

50.2 样条插补介绍
在数学学科数值分析中,样条是一种特殊的函数,由多项式分段定义。样条的英语单词spline来源于可变形的样条工具,那是一种在造船和工程制图时用来画出光滑形状的工具。在中国大陆,早期曾经被称做“齿函数”。后来因为工程学术语中“放样”一词而得名。在插值问题中,样条插值通常比多项式插值好用。用低阶的样条插值能产生和高阶的多项式插值类似的效果,并且可以避免被称为龙格现象的数值不稳定的出现。并且低阶的样条插值还具有“保凸”的重要性质。在计算机科学的计算机辅助设计和计算机图形学中,样条通常是指分段定义的多项式参数曲线。由于样条构造简单,使用方便,拟合准确,并能近似曲线拟合和交互式曲线设计中复杂的形状,样条是这些领域中曲线的常用表示方法

50.3 样条插补实现
样条插补主要通过下面两个函数实现。

50.3.1 函数arm_spline_init_f32
函数原型:

  1. void arm_spline_init_f32(
  2.         arm_spline_instance_f32 * S,
  3.         arm_spline_type type,
  4.   const float32_t * x,
  5.   const float32_t * y,
  6.         uint32_t n,
  7.         float32_t * coeffs,
  8.         float32_t * tempBuffer)
复制代码

函数描述:

此函数用于样条函数初始化。

函数参数:

  第1个参数是arm_spline_instance_f32类型结构体变量。
  第2个参数是样条类型选择:
  ARM_SPLINE_NATURAL 表自然样条。
  ARM_SPLINE_PARABOLIC_RUNOUT 表示抛物线样条。
  第3个参数是原始数据x轴坐标值。
  第4个参数是原始数据y轴坐标值。
  第5个参数是原始数据个数。
  第6个参数是插补因数缓存。
  第7个参数是临时缓冲。
注意事项:

  x轴坐标数据必须是递增方式。
  第6个参数插补因数缓存大小问题,如果原始数据个数是n,那么插补因数个数必须要大于等于3*(n-1)。
  第7个参数临时缓冲大小问题,如果原始数据个数是n,那么临时缓冲大小必须大于等于2*n - 1

50.3.2 函数arm_spline_f32
函数原型:

  1. void arm_spline_f32(
  2.         arm_spline_instance_f32 * S,
  3.   const float32_t * xq,
  4.         float32_t * pDst,
  5.         uint32_t blockSize)
复制代码

函数描述:

此函数用于样条插补实现。

函数参数:

  第1个参数是arm_spline_instance_f32类型结构体变量。
  第2个参数是插补后的x轴坐标值,需要用户指定,注意坐标值一定是递增的。
  第3个参数是经过插补计算后输出的y轴数值
  第4个参数是数据输出个数
50.3.3 使用样条插补函数的关键点
样条插补的主要作用是使得波形更加平滑。比如一帧是128点,步大小是8个像素,我们可以通过插补实现步长为1, 1024点的波形,本质是你的总步长大小不能变,我们这里都是1024,这个不能变,在这个基础上做插补,效果就出来了。

这个认识非常重要,否则无法正常使用插补算法。

50.3.4 自然样条插补测试
样条测试代码的实现如下:

  1. #define INPUT_TEST_LENGTH_SAMPLES         128  /* 输入数据个数 */
  2. #define OUT_TEST_LENGTH_SAMPLES           1024   /* 输出数据个数 */

  3. #define SpineTab OUT_TEST_LENGTH_SAMPLES/INPUT_TEST_LENGTH_SAMPLES  /* 插补末尾的8个坐标值不使用 */


  4. float32_t xn[INPUT_TEST_LENGTH_SAMPLES];   /* 输入数据x轴坐标 */
  5. float32_t yn[INPUT_TEST_LENGTH_SAMPLES];   /* 输入数据y轴坐标 */

  6. float32_t coeffs[3*(INPUT_TEST_LENGTH_SAMPLES - 1)];     /* 插补系数缓冲 */  
  7. float32_t tempBuffer[2 * INPUT_TEST_LENGTH_SAMPLES - 1]; /* 插补临时缓冲 */  

  8. float32_t xnpos[OUT_TEST_LENGTH_SAMPLES];  /* 插补计算后X轴坐标值 */
  9. float32_t ynpos[OUT_TEST_LENGTH_SAMPLES];  /* 插补计算后Y轴数值 */

  10. /*
  11. *********************************************************************************************************
  12. *        函 数 名: main
  13. *        功能说明: c程序入口
  14. *        形    参: 无
  15. *        返 回 值: 错误代码(无需处理)
  16. *********************************************************************************************************
  17. */
  18. int main(void)
  19. {
  20.         uint32_t i;
  21.         uint32_t idx2;
  22.         uint8_t ucKeyCode;        
  23.         arm_spline_instance_f32 S;
  24.         

  25.         bsp_Init();                /* 硬件初始化 */
  26.         PrintfLogo();        /* 打印例程名称和版本等信息 */
  27.         PrintfHelp();        /* 打印操作提示 */
  28.         
  29.         bsp_StartAutoTimer(0, 100);        /* 启动1个100ms的自动重装的定时器 */
  30.         
  31.         
  32.         /* 原始x轴数值和y轴数值 */
  33.         for(i=0; i<INPUT_TEST_LENGTH_SAMPLES; i++)
  34.         {
  35.                 xn<i> = i*SpineTab;
  36.              </i>   yn = 1 + cos(2*3.1415926*50*i/256 + 3.1415926/3);
  37.         }
  38.         
  39.         /* 插补后X轴坐标值,这个是需要用户设置的 */
  40.         for(i=0; i<OUT_TEST_LENGTH_SAMPLES; i++)
  41.         {
  42.                 xnpos = i;
  43.         }
  44.         
  45.         while (1)
  46.         {
  47.                 bsp_Idle();                /* 这个函数在bsp.c文件。用户可以修改这个函数实现CPU休眠和喂狗 */

  48.                 /* 判断定时器超时时间 */
  49.                 if (bsp_CheckTimer(0))        
  50.                 {
  51.                         /* 每隔100ms 进来一次 */  
  52.                         bsp_LedToggle(2);
  53.                 }

  54.                 ucKeyCode = bsp_GetKey();        /* 读取键值, 无键按下时返回 KEY_NONE = 0 */
  55.                 if (ucKeyCode != KEY_NONE)
  56.                 {
  57.                         switch (ucKeyCode)
  58.                         {
  59.                                 case KEY_DOWN_K1:  /* K1键按下,自然样条插补 */
  60.                                         /* 样条初始化 */
  61.                                         arm_spline_init_f32(&S,
  62.                                                                                 ARM_SPLINE_NATURAL ,
  63.                                                                                 xn,
  64.                                                                                 yn,
  65.                                                                                 INPUT_TEST_LENGTH_SAMPLES,
  66.                                                                                 coeffs,
  67.                                                                                 tempBuffer);
  68.                                         /* 样条计算 */
  69.                                         arm_spline_f32        (&S,
  70.                                                                          xnpos,
  71.                                                                          ynpos,
  72.                                                                          OUT_TEST_LENGTH_SAMPLES);

  73.                                 
  74.                                         /* 打印输出输出 */
  75.                                         idx2 = 0;
  76.                                         for (i = 0; i < OUT_TEST_LENGTH_SAMPLES-SpineTab; i++)
  77.                                         {        
  78.                                                         if ((i % SpineTab) == 0)
  79.                                                         {
  80.                                                                         printf("%f,%f\r\n", ynpos, yn[idx2++]);
  81.                                                         }
  82.                                                         else
  83.                                                         {
  84.                                                                         printf("%f,\r\n", ynpos);
  85.                                                         }
  86.                                         }
  87.                                         break;
  88.                                 
  89.                                 default:
  90.                                         /* 其它的键值不处理 */
  91.                                         break;
  92.                         }
  93.                 }
  94.         }
  95. }
复制代码

代码里面的几个关键地方:

  原始坐标数组xn和yn是128组,而我们通过插补生成的是1024组xnpos和ynpos,其中1024组的xnpos需要用户设置初值,这点不能忽略。
  函数arm_spline_init_f32用于样条函数初始化,这里特别注意,此函数主要是对原始数据的操作。自然样条插补用的ARM_SPLINE_NATURAL。
  函数arm_spline_f32用于样条函数计算。
实际输出效果如下:

5ad87cda57524d4c99fc092c22b611c2.png




50.3.5 抛物线样条插补测试
样条测试代码的实现如下:

  1. #define INPUT_TEST_LENGTH_SAMPLES         128  /* 输入数据个数 */
  2. #define OUT_TEST_LENGTH_SAMPLES           1024   /* 输出数据个数 */

  3. #define SpineTab OUT_TEST_LENGTH_SAMPLES/INPUT_TEST_LENGTH_SAMPLES  /* 插补末尾的8个坐标值不使用 */


  4. float32_t xn[INPUT_TEST_LENGTH_SAMPLES];   /* 输入数据x轴坐标 */
  5. float32_t yn[INPUT_TEST_LENGTH_SAMPLES];   /* 输入数据y轴坐标 */

  6. float32_t coeffs[3*(INPUT_TEST_LENGTH_SAMPLES - 1)];     /* 插补系数缓冲 */  
  7. float32_t tempBuffer[2 * INPUT_TEST_LENGTH_SAMPLES - 1]; /* 插补临时缓冲 */  

  8. float32_t xnpos[OUT_TEST_LENGTH_SAMPLES];  /* 插补计算后X轴坐标值 */
  9. float32_t ynpos[OUT_TEST_LENGTH_SAMPLES];  /* 插补计算后Y轴数值 */

  10. /*
  11. *********************************************************************************************************
  12. *        函 数 名: main
  13. *        功能说明: c程序入口
  14. *        形    参: 无
  15. *        返 回 值: 错误代码(无需处理)
  16. *********************************************************************************************************
  17. */
  18. int main(void)
  19. {
  20.         uint32_t i;
  21.         uint32_t idx2;
  22.         uint8_t ucKeyCode;        
  23.         arm_spline_instance_f32 S;
  24.         

  25.         bsp_Init();                /* 硬件初始化 */
  26.         PrintfLogo();        /* 打印例程名称和版本等信息 */
  27.         PrintfHelp();        /* 打印操作提示 */
  28.         
  29.         bsp_StartAutoTimer(0, 100);        /* 启动1个100ms的自动重装的定时器 */
  30.         
  31.         
  32.         /* 原始x轴数值和y轴数值 */
  33.         for(i=0; i<INPUT_TEST_LENGTH_SAMPLES; i++)
  34.         {
  35.                 xn<span style="font-style: italic;"><span style="font-style: normal;"> = i*SpineTab;
  36.                 yn</span><span style="font-style: normal;"> = 1 + cos(2*3.1415926*50*i/256 + 3.1415926/3);
  37.         }
  38.         
  39.         /* 插补后X轴坐标值,这个是需要用户设置的 */
  40.         for(i=0; i<OUT_TEST_LENGTH_SAMPLES; i++)
  41.         {
  42.                 xnpos</span><span style="font-style: normal;"> = i;
  43.         }
  44.         
  45.         while (1)
  46.         {
  47.                 bsp_Idle();                /* 这个函数在bsp.c文件。用户可以修改这个函数实现CPU休眠和喂狗 */

  48.                 /* 判断定时器超时时间 */
  49.                 if (bsp_CheckTimer(0))        
  50.                 {
  51.                         /* 每隔100ms 进来一次 */  
  52.                         bsp_LedToggle(2);
  53.                 }

  54.                 ucKeyCode = bsp_GetKey();        /* 读取键值, 无键按下时返回 KEY_NONE = 0 */
  55.                 if (ucKeyCode != KEY_NONE)
  56.                 {
  57.                         switch (ucKeyCode)
  58.                         {
  59.                                 case KEY_DOWN_K2:                        /* K2键按下,抛物线样条插补 */
  60.                                         /* 样条初始化 */
  61.                                         arm_spline_init_f32(&S,
  62.                                                                                 ARM_SPLINE_PARABOLIC_RUNOUT ,
  63.                                                                                 xn,
  64.                                                                                 yn,
  65.                                                                                 INPUT_TEST_LENGTH_SAMPLES,
  66.                                                                                 coeffs,
  67.                                                                                 tempBuffer);
  68.                                         /* 样条计算 */
  69.                                         arm_spline_f32        (&S,
  70.                                                                          xnpos,
  71.                                                                          ynpos,
  72.                                                                          OUT_TEST_LENGTH_SAMPLES);

  73.                                 
  74.                                         /* 打印输出输出 */
  75.                                         idx2 = 0;
  76.                                         for (i = 0; i < OUT_TEST_LENGTH_SAMPLES-SpineTab; i++)
  77.                                         {        
  78.                                                         if ((i % SpineTab) == 0)
  79.                                                         {
  80.                                                                         printf("%f,%f\r\n", ynpos</span><span style="font-style: normal;">, yn[idx2++]);
  81.                                                         }
  82.                                                         else
  83.                                                         {
  84.                                                                         printf("%f,\r\n", ynpos</span><span style="font-style: normal;">);
  85.                                                         }
  86.                                         }
  87.                                         break;
  88.                                 
  89.                                 default:
  90.                                         /* 其它的键值不处理 */
  91.                                         break;
  92.                         }
  93.                 }
  94.         }
  95. }</span></span>
复制代码

代码里面的几个关键地方:

  原始坐标数组xn和yn是128组,而我们通过插补生成的是1024组xnpos和ynpos,其中1024组的xnpos需要用户设置初值,这点不能忽略。
  函数arm_spline_init_f32用于样条函数初始化,这里特别注意,此函数主要是对原始数据的操作。抛物线样条插补用的ARM_SPLINE_PARABOLIC_RUNOUT。
  函数arm_spline_f32用于样条函数计算。


实际输出效果如下:

13a066ab945a455cbf93927d3e1bcfe7.png




50.4 实验例程说明(MDK)

配套例子:

V7-235_样条插补,波形拟合丝滑顺畅

实验目的:

学习样条插补的实现。

实验内容:

启动一个自动重装软件定时器,每100ms翻转一次LED2。
K1键按下,自然样条插补测试。
K2键按下,抛物线样插补测试。


使用AC6注意事项

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

上电后串口打印的信息:

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

5b12ec8275204957827f235f8dfa148d.png


RTT方式打印信息:

05ddd5a9bd0a424aab56d8613ae53d45.png



程序设计:

  系统栈大小分配:

3759d3e86623406aafb0011de1bae045.png



  RAM空间用的DTCM:

f3bab02ea37543ca8f62633409bf45e6.png



  硬件外设初始化

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

  1. /*
  2. *********************************************************************************************************
  3. *        函 数 名: bsp_Init
  4. *        功能说明: 初始化所有的硬件设备。该函数配置CPU寄存器和外设的寄存器并初始化一些全局变量。只需要调用一次
  5. *        形    参:无
  6. *        返 回 值: 无
  7. *********************************************************************************************************
  8. */
  9. void bsp_Init(void)
  10. {
  11.     /* 配置MPU */
  12.         MPU_Config();
  13.         
  14.         /* 使能L1 Cache */
  15.         CPU_CACHE_Enable();

  16.         /*
  17.        STM32H7xx HAL 库初始化,此时系统用的还是H7自带的64MHz,HSI时钟:
  18.            - 调用函数HAL_InitTick,初始化滴答时钟中断1ms。
  19.            - 设置NVIC优先级分组为4。
  20.          */
  21.         HAL_Init();

  22.         /*
  23.        配置系统时钟到400MHz
  24.        - 切换使用HSE。
  25.        - 此函数会更新全局变量SystemCoreClock,并重新配置HAL_InitTick。
  26.     */
  27.         SystemClock_Config();

  28.         /*
  29.            Event Recorder:
  30.            - 可用于代码执行时间测量,MDK5.25及其以上版本才支持,IAR不支持。
  31.            - 默认不开启,如果要使能此选项,务必看V7开发板用户手册第8章
  32.         */        
  33. #if Enable_EventRecorder == 1  
  34.         /* 初始化EventRecorder并开启 */
  35.         EventRecorderInitialize(EventRecordAll, 1U);
  36.         EventRecorderStart();
  37. #endif
  38.         
  39.         bsp_InitKey();            /* 按键初始化,要放在滴答定时器之前,因为按钮检测是通过滴答定时器扫描 */
  40.         bsp_InitTimer();          /* 初始化滴答定时器 */
  41.         bsp_InitUart();        /* 初始化串口 */
  42.         bsp_InitExtIO();        /* 初始化FMC总线74HC574扩展IO. 必须在 bsp_InitLed()前执行 */        
  43.         bsp_InitLed();            /* 初始化LED */        
  44. }
复制代码

  MPU配置和Cache配置:

数据Cache和指令Cache都开启。配置了AXI SRAM区(本例子未用到AXI SRAM),FMC的扩展IO区。

  1. /*
  2. *********************************************************************************************************
  3. *        函 数 名: MPU_Config
  4. *        功能说明: 配置MPU
  5. *        形    参: 无
  6. *        返 回 值: 无
  7. *********************************************************************************************************
  8. */
  9. static void MPU_Config( void )
  10. {
  11.         MPU_Region_InitTypeDef MPU_InitStruct;

  12.         /* 禁止 MPU */
  13.         HAL_MPU_Disable();

  14.         /* 配置AXI SRAM的MPU属性为关闭读Cache和写Cache */
  15.         MPU_InitStruct.Enable           = MPU_REGION_ENABLE;
  16.         MPU_InitStruct.BaseAddress      = 0x24000000;
  17.         MPU_InitStruct.Size             = MPU_REGION_SIZE_512KB;
  18.         MPU_InitStruct.AccessPermission = MPU_REGION_FULL_ACCESS;
  19.         MPU_InitStruct.IsBufferable     = MPU_ACCESS_NOT _BUFFERABLE;
  20.         MPU_InitStruct.IsCacheable      = MPU_ACCESS_NOT _CACHEABLE;
  21.         MPU_InitStruct.IsShareable      = MPU_ACCESS_NOT_SHAREABLE;
  22.         MPU_InitStruct.Number           = MPU_REGION_NUMBER0;
  23.         MPU_InitStruct.TypeExtField     = MPU_TEX_LEVEL0;
  24.         MPU_InitStruct.SubRegionDisable = 0x00;
  25.         MPU_InitStruct.DisableExec      = MPU_INSTRUCTION_ACCESS_ENABLE;

  26.         HAL_MPU_ConfigRegion(&MPU_InitStruct);
  27.         
  28.         
  29.         /* 配置FMC扩展IO的MPU属性为Device或者Strongly Ordered */
  30.         MPU_InitStruct.Enable           = MPU_REGION_ENABLE;
  31.         MPU_InitStruct.BaseAddress      = 0x60000000;
  32.         MPU_InitStruct.Size             = ARM_MPU_REGION_SIZE_64KB;        
  33.         MPU_InitStruct.AccessPermission = MPU_REGION_FULL_ACCESS;
  34.         MPU_InitStruct.IsBufferable     = MPU_ACCESS_BUFFERABLE;
  35.         MPU_InitStruct.IsCacheable      = MPU_ACCESS_NOT_CACHEABLE;        
  36.         MPU_InitStruct.IsShareable      = MPU_ACCESS_NOT_SHAREABLE;
  37.         MPU_InitStruct.Number           = MPU_REGION_NUMBER1;
  38.         MPU_InitStruct.TypeExtField     = MPU_TEX_LEVEL0;
  39.         MPU_InitStruct.SubRegionDisable = 0x00;
  40.         MPU_InitStruct.DisableExec      = MPU_INSTRUCTION_ACCESS_ENABLE;
  41.         
  42.         HAL_MPU_ConfigRegion(&MPU_InitStruct);

  43.         /*使能 MPU */
  44.         HAL_MPU_Enable(MPU_PRIVILEGED_DEFAULT);
  45. }

  46. /*
  47. *********************************************************************************************************
  48. *        函 数 名: CPU_CACHE_Enable
  49. *        功能说明: 使能L1 Cache
  50. *        形    参: 无
  51. *        返 回 值: 无
  52. *********************************************************************************************************
  53. */
  54. static void CPU_CACHE_Enable(void)
  55. {
  56.         /* 使能 I-Cache */
  57.         SCB_EnableICache();

  58.         /* 使能 D-Cache */
  59.         SCB_EnableDCache();
  60. }
复制代码

  主功能:

主程序实现如下操作:

  启动一个自动重装软件定时器,每100ms翻转一次LED2。
  K1键按下,自然样条插补测试。
  K2键按下,抛物线样插补测试。
  1. /*
  2. *********************************************************************************************************
  3. *        函 数 名: main
  4. *        功能说明: c程序入口
  5. *        形    参: 无
  6. *        返 回 值: 错误代码(无需处理)
  7. *********************************************************************************************************
  8. */
  9. int main(void)
  10. {
  11.         uint32_t i;
  12.         uint32_t idx2;
  13.         uint8_t ucKeyCode;        
  14.         arm_spline_instance_f32 S;
  15.         

  16.         bsp_Init();                /* 硬件初始化 */
  17.         PrintfLogo();        /* 打印例程名称和版本等信息 */
  18.         PrintfHelp();        /* 打印操作提示 */
  19.         
  20.         bsp_StartAutoTimer(0, 100);        /* 启动1个100ms的自动重装的定时器 */
  21.         
  22.         
  23.         /* 原始x轴数值和y轴数值 */
  24.         for(i=0; i<INPUT_TEST_LENGTH_SAMPLES; i++)
  25.         {
  26.                 xn<span style="font-style: italic;"><span style="font-style: normal;"> = i*SpineTab;
  27.                 yn</span><span style="font-style: normal;"> = 1 + cos(2*3.1415926*50*i/256 + 3.1415926/3);
  28.         }
  29.         
  30.         /* 插补后X轴坐标值,这个是需要用户设置的 */
  31.         for(i=0; i<OUT_TEST_LENGTH_SAMPLES; i++)
  32.         {
  33.                 xnpos</span><span style="font-style: normal;"> = i;
  34.         }
  35.         
  36.         while (1)
  37.         {
  38.                 bsp_Idle();                /* 这个函数在bsp.c文件。用户可以修改这个函数实现CPU休眠和喂狗 */

  39.                 /* 判断定时器超时时间 */
  40.                 if (bsp_CheckTimer(0))        
  41.                 {
  42.                         /* 每隔100ms 进来一次 */  
  43.                         bsp_LedToggle(2);
  44.                 }

  45.                 ucKeyCode = bsp_GetKey();        /* 读取键值, 无键按下时返回 KEY_NONE = 0 */
  46.                 if (ucKeyCode != KEY_NONE)
  47.                 {
  48.                         switch (ucKeyCode)
  49.                         {
  50.                                 case KEY_DOWN_K1:  /* K1键按下,自然样条插补 */
  51.                                         /* 样条初始化 */
  52.                                         arm_spline_init_f32(&S,
  53.                                                                                 ARM_SPLINE_NATURAL ,
  54.                                                                                 xn,
  55.                                                                                 yn,
  56.                                                                                 INPUT_TEST_LENGTH_SAMPLES,
  57.                                                                                 coeffs,
  58.                                                                                 tempBuffer);
  59.                                         /* 样条计算 */
  60.                                         arm_spline_f32        (&S,
  61.                                                                          xnpos,
  62.                                                                          ynpos,
  63.                                                                          OUT_TEST_LENGTH_SAMPLES);

  64.                                 
  65.                                         /* 打印输出输出 */
  66.                                         idx2 = 0;
  67.                                         for (i = 0; i < OUT_TEST_LENGTH_SAMPLES-SpineTab; i++)
  68.                                         {        
  69.                                                         if ((i % SpineTab) == 0)
  70.                                                         {
  71.                                                                         printf("%f,%f\r\n", ynpos</span><span style="font-style: normal;">, yn[idx2++]);
  72.                                                         }
  73.                                                         else
  74.                                                         {
  75.                                                                         printf("%f,\r\n", ynpos</span><span style="font-style: normal;">);
  76.                                                         }
  77.                                         }
  78.                                         break;

  79.                                 case KEY_DOWN_K2:                        /* K2键按下,抛物线样条插补 */
  80.                                         /* 样条初始化 */
  81.                                         arm_spline_init_f32(&S,
  82.                                                                                 ARM_SPLINE_PARABOLIC_RUNOUT ,
  83.                                                                                 xn,
  84.                                                                                 yn,
  85.                                                                                 INPUT_TEST_LENGTH_SAMPLES,
  86.                                                                                 coeffs,
  87.                                                                                 tempBuffer);
  88.                                         /* 样条计算 */
  89.                                         arm_spline_f32        (&S,
  90.                                                                          xnpos,
  91.                                                                          ynpos,
  92.                                                                          OUT_TEST_LENGTH_SAMPLES);

  93.                                 
  94.                                         /* 打印输出输出 */
  95.                                         idx2 = 0;
  96.                                         for (i = 0; i < OUT_TEST_LENGTH_SAMPLES-SpineTab; i++)
  97.                                         {        
  98.                                                         if ((i % SpineTab) == 0)
  99.                                                         {
  100.                                                                         printf("%f,%f\r\n", ynpos</span><span style="font-style: normal;">, yn[idx2++]);
  101.                                                         }
  102.                                                         else
  103.                                                         {
  104.                                                                         printf("%f,\r\n", ynpos</span><span style="font-style: normal;">);
  105.                                                         }
  106.                                         }
  107.                                         break;
  108.                                 
  109.                                 default:
  110.                                         /* 其它的键值不处理 */
  111.                                         break;
  112.                         }
  113.                 }
  114.         }
  115. }</span></span>
复制代码

50.5 实验例程说明(IAR)
配套例子:

V7-235_样条插补,波形拟合丝滑顺畅

实验目的:

学习样条插补的实现。


实验内容:

启动一个自动重装软件定时器,每100ms翻转一次LED2。
K1键按下,自然样条插补测试。
K2键按下,抛物线样插补测试。


上电后串口打印的信息:

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

3c1c8f8b6298498f93b65e9bd514b096.png


RTT方式打印信息:

1bb1f5ed7beb49abbddd27af92c6f5c9.png




程序设计:

  系统栈大小分配:

90d4d24f3be74d308bff1f31cb282b7f.png




  RAM空间用的DTCM:

3f73ce7e46c941d4a7737bccd79524ce.png




  硬件外设初始化

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

  1. /*
  2. *********************************************************************************************************
  3. *        函 数 名: bsp_Init
  4. *        功能说明: 初始化所有的硬件设备。该函数配置CPU寄存器和外设的寄存器并初始化一些全局变量。只需要调用一次
  5. *        形    参:无
  6. *        返 回 值: 无
  7. *********************************************************************************************************
  8. */
  9. void bsp_Init(void)
  10. {
  11.     /* 配置MPU */
  12.         MPU_Config();
  13.         
  14.         /* 使能L1 Cache */
  15.         CPU_CACHE_Enable();

  16.         /*
  17.        STM32H7xx HAL 库初始化,此时系统用的还是H7自带的64MHz,HSI时钟:
  18.            - 调用函数HAL_InitTick,初始化滴答时钟中断1ms。
  19.            - 设置NVIC优先级分组为4。
  20.          */
  21.         HAL_Init();

  22.         /*
  23.        配置系统时钟到400MHz
  24.        - 切换使用HSE。
  25.        - 此函数会更新全局变量SystemCoreClock,并重新配置HAL_InitTick。
  26.     */
  27.         SystemClock_Config();

  28.         /*
  29.            Event Recorder:
  30.            - 可用于代码执行时间测量,MDK5.25及其以上版本才支持,IAR不支持。
  31.            - 默认不开启,如果要使能此选项,务必看V7开发板用户手册第8章
  32.         */        
  33. #if Enable_EventRecorder == 1  
  34.         /* 初始化EventRecorder并开启 */
  35.         EventRecorderInitialize(EventRecordAll, 1U);
  36.         EventRecorderStart();
  37. #endif
  38.         
  39.         bsp_InitKey();            /* 按键初始化,要放在滴答定时器之前,因为按钮检测是通过滴答定时器扫描 */
  40.         bsp_InitTimer();          /* 初始化滴答定时器 */
  41.         bsp_InitUart();        /* 初始化串口 */
  42.         bsp_InitExtIO();        /* 初始化FMC总线74HC574扩展IO. 必须在 bsp_InitLed()前执行 */        
  43.         bsp_InitLed();            /* 初始化LED */        
  44. }
复制代码

  MPU配置和Cache配置:

数据Cache和指令Cache都开启。配置了AXI SRAM区(本例子未用到AXI SRAM),FMC的扩展IO区。

  1. /*
  2. *********************************************************************************************************
  3. *        函 数 名: MPU_Config
  4. *        功能说明: 配置MPU
  5. *        形    参: 无
  6. *        返 回 值: 无
  7. *********************************************************************************************************
  8. */
  9. static void MPU_Config( void )
  10. {
  11.         MPU_Region_InitTypeDef MPU_InitStruct;

  12.         /* 禁止 MPU */
  13.         HAL_MPU_Disable();

  14.         /* 配置AXI SRAM的MPU属性为Write back, Read allocate,Write allocate */
  15.         MPU_InitStruct.Enable           = MPU_REGION_ENABLE;
  16.         MPU_InitStruct.BaseAddress      = 0x24000000;
  17.         MPU_InitStruct.Size             = MPU_REGION_SIZE_512KB;
  18.         MPU_InitStruct.AccessPermission = MPU_REGION_FULL_ACCESS;
  19.         MPU_InitStruct.IsBufferable     = MPU_ACCESS_BUFFERABLE;
  20.         MPU_InitStruct.IsCacheable      = MPU_ACCESS_CACHEABLE;
  21.         MPU_InitStruct.IsShareable      = MPU_ACCESS_NOT_SHAREABLE;
  22.         MPU_InitStruct.Number           = MPU_REGION_NUMBER0;
  23.         MPU_InitStruct.TypeExtField     = MPU_TEX_LEVEL1;
  24.         MPU_InitStruct.SubRegionDisable = 0x00;
  25.         MPU_InitStruct.DisableExec      = MPU_INSTRUCTION_ACCESS_ENABLE;

  26.         HAL_MPU_ConfigRegion(&MPU_InitStruct);
  27.         
  28.         
  29.         /* 配置FMC扩展IO的MPU属性为Device或者Strongly Ordered */
  30.         MPU_InitStruct.Enable           = MPU_REGION_ENABLE;
  31.         MPU_InitStruct.BaseAddress      = 0x60000000;
  32.         MPU_InitStruct.Size             = ARM_MPU_REGION_SIZE_64KB;        
  33.         MPU_InitStruct.AccessPermission = MPU_REGION_FULL_ACCESS;
  34.         MPU_InitStruct.IsBufferable     = MPU_ACCESS_BUFFERABLE;
  35.         MPU_InitStruct.IsCacheable      = MPU_ACCESS_NOT_CACHEABLE;        
  36.         MPU_InitStruct.IsShareable      = MPU_ACCESS_NOT_SHAREABLE;
  37.         MPU_InitStruct.Number           = MPU_REGION_NUMBER1;
  38.         MPU_InitStruct.TypeExtField     = MPU_TEX_LEVEL0;
  39.         MPU_InitStruct.SubRegionDisable = 0x00;
  40.         MPU_InitStruct.DisableExec      = MPU_INSTRUCTION_ACCESS_ENABLE;
  41.         
  42.         HAL_MPU_ConfigRegion(&MPU_InitStruct);

  43.         /*使能 MPU */
  44.         HAL_MPU_Enable(MPU_PRIVILEGED_DEFAULT);
  45. }

  46. /*
  47. *********************************************************************************************************
  48. *        函 数 名: CPU_CACHE_Enable
  49. *        功能说明: 使能L1 Cache
  50. *        形    参: 无
  51. *        返 回 值: 无
  52. *********************************************************************************************************
  53. */
  54. static void CPU_CACHE_Enable(void)
  55. {
  56.         /* 使能 I-Cache */
  57.         SCB_EnableICache();

  58.         /* 使能 D-Cache */
  59.         SCB_EnableDCache();
  60. }
复制代码

  主功能:

主程序实现如下操作:

  启动一个自动重装软件定时器,每100ms翻转一次LED2。
  K1键按下,自然样条插补测试。
  K2键按下,抛物线样插补测试。
  1. /*
  2. *********************************************************************************************************
  3. *        函 数 名: main
  4. *        功能说明: c程序入口
  5. *        形    参: 无
  6. *        返 回 值: 错误代码(无需处理)
  7. *********************************************************************************************************
  8. */
  9. int main(void)
  10. {
  11.         uint32_t i;
  12.         uint32_t idx2;
  13.         uint8_t ucKeyCode;        
  14.         arm_spline_instance_f32 S;
  15.         

  16.         bsp_Init();                /* 硬件初始化 */
  17.         PrintfLogo();        /* 打印例程名称和版本等信息 */
  18.         PrintfHelp();        /* 打印操作提示 */
  19.         
  20.         bsp_StartAutoTimer(0, 100);        /* 启动1个100ms的自动重装的定时器 */
  21.         
  22.         
  23.         /* 原始x轴数值和y轴数值 */
  24.         for(i=0; i<INPUT_TEST_LENGTH_SAMPLES; i++)
  25.         {
  26.                 xn<span style="font-style: italic;"><span style="font-style: normal;"> = i*SpineTab;
  27.                 yn</span><span style="font-style: normal;"> = 1 + cos(2*3.1415926*50*i/256 + 3.1415926/3);
  28.         }
  29.         
  30.         /* 插补后X轴坐标值,这个是需要用户设置的 */
  31.         for(i=0; i<OUT_TEST_LENGTH_SAMPLES; i++)
  32.         {
  33.                 xnpos</span><span style="font-style: normal;"> = i;
  34.         }
  35.         
  36.         while (1)
  37.         {
  38.                 bsp_Idle();                /* 这个函数在bsp.c文件。用户可以修改这个函数实现CPU休眠和喂狗 */

  39.                 /* 判断定时器超时时间 */
  40.                 if (bsp_CheckTimer(0))        
  41.                 {
  42.                         /* 每隔100ms 进来一次 */  
  43.                         bsp_LedToggle(2);
  44.                 }

  45.                 ucKeyCode = bsp_GetKey();        /* 读取键值, 无键按下时返回 KEY_NONE = 0 */
  46.                 if (ucKeyCode != KEY_NONE)
  47.                 {
  48.                         switch (ucKeyCode)
  49.                         {
  50.                                 case KEY_DOWN_K1:  /* K1键按下,自然样条插补 */
  51.                                         /* 样条初始化 */
  52.                                         arm_spline_init_f32(&S,
  53.                                                                                 ARM_SPLINE_NATURAL ,
  54.                                                                                 xn,
  55.                                                                                 yn,
  56.                                                                                 INPUT_TEST_LENGTH_SAMPLES,
  57.                                                                                 coeffs,
  58.                                                                                 tempBuffer);
  59.                                         /* 样条计算 */
  60.                                         arm_spline_f32        (&S,
  61.                                                                          xnpos,
  62.                                                                          ynpos,
  63.                                                                          OUT_TEST_LENGTH_SAMPLES);

  64.                                 
  65.                                         /* 打印输出输出 */
  66.                                         idx2 = 0;
  67.                                         for (i = 0; i < OUT_TEST_LENGTH_SAMPLES-SpineTab; i++)
  68.                                         {        
  69.                                                         if ((i % SpineTab) == 0)
  70.                                                         {
  71.                                                                         printf("%f,%f\r\n", ynpos</span><span style="font-style: normal;">, yn[idx2++]);
  72.                                                         }
  73.                                                         else
  74.                                                         {
  75.                                                                         printf("%f,\r\n", ynpos</span><span style="font-style: normal;">);
  76.                                                         }
  77.                                         }
  78.                                         break;

  79.                                 case KEY_DOWN_K2:                        /* K2键按下,抛物线样条插补 */
  80.                                         /* 样条初始化 */
  81.                                         arm_spline_init_f32(&S,
  82.                                                                                 ARM_SPLINE_PARABOLIC_RUNOUT ,
  83.                                                                                 xn,
  84.                                                                                 yn,
  85.                                                                                 INPUT_TEST_LENGTH_SAMPLES,
  86.                                                                                 coeffs,
  87.                                                                                 tempBuffer);
  88.                                         /* 样条计算 */
  89.                                         arm_spline_f32        (&S,
  90.                                                                          xnpos,
  91.                                                                          ynpos,
  92.                                                                          OUT_TEST_LENGTH_SAMPLES);

  93.                                 
  94.                                         /* 打印输出输出 */
  95.                                         idx2 = 0;
  96.                                         for (i = 0; i < OUT_TEST_LENGTH_SAMPLES-SpineTab; i++)
  97.                                         {        
  98.                                                         if ((i % SpineTab) == 0)
  99.                                                         {
  100.                                                                         printf("%f,%f\r\n", ynpos</span><span style="font-style: normal;">, yn[idx2++]);
  101.                                                         }
  102.                                                         else
  103.                                                         {
  104.                                                                         printf("%f,\r\n", ynpos</span><span style="font-style: normal;">);
  105.                                                         }
  106.                                         }
  107.                                         break;
  108.                                 
  109.                                 default:
  110.                                         /* 其它的键值不处理 */
  111.                                         break;
  112.                         }
  113.                 }
  114.         }
  115. }</span></span>
复制代码
50.6 总结
本章节主要讲解了样条插补的实现,实际项目比较实用,有兴趣可以深入源码了解。



收藏 评论0 发布时间:2022-1-1 21:00

举报

0个回答

所属标签

相似分享

官网相关资源

关于
我们是谁
投资者关系
意法半导体可持续发展举措
创新与技术
意法半导体官网
联系我们
联系ST分支机构
寻找销售人员和分销渠道
社区
媒体中心
活动与培训
隐私策略
隐私策略
Cookies管理
行使您的权利
官方最新发布
STM32N6 AI生态系统
STM32MCU,MPU高性能GUI
ST ACEPACK电源模块
意法半导体生物传感器
STM32Cube扩展软件包
关注我们
st-img 微信公众号
st-img 手机版