首页 > 其他分享 >MCU上三角函数运算速度对比

MCU上三角函数运算速度对比

时间：2024-03-07 17:25:31浏览次数：21  
标签：运算 三角函数 CORDIC float us time output MCU type

前言

使用170M主频M4f核芯片(STM32G431),启用浮点加速，启用arm_dsp库。
对比条件：

方式 输入 输出
cordic q15 q15
arm_dsp_q15 q15 q15
arm_dsp_q31 q31 q31
arm_dsp_f32 float float
math float float
math double double
taylor6 float float

测试代码

硬件cordic

/**
  * @brief  Trigonometrical functions type definition
  */
typedef struct
{
  int16_t hCos;
  int16_t hSin;
} Trig_Components;
/* CORDIC FUNCTION: COSINE q1.15 */
#define CORDIC_CONFIG_COSINE    (LL_CORDIC_FUNCTION_COSINE | LL_CORDIC_PRECISION_6CYCLES | LL_CORDIC_SCALE_0 |\
         LL_CORDIC_NBWRITE_1 | LL_CORDIC_NBREAD_1 |\
         LL_CORDIC_INSIZE_16BITS | LL_CORDIC_OUTSIZE_16BITS)

__weak Trig_Components MCM_Trig_Functions(int16_t hAngle)
{
  union u32toi16x2 {
    uint32_t CordicRdata;
    Trig_Components Components;
  } CosSin;

  WRITE_REG(CORDIC->CSR, CORDIC_CONFIG_COSINE);
  LL_CORDIC_WriteData(CORDIC, ((uint32_t)0x7FFF0000) + ((uint32_t)hAngle));
  CosSin.CordicRdata = LL_CORDIC_ReadData(CORDIC);
  return (CosSin.Components);
 
}

泰勒展开函数

float factorial(int n) {
    if (n == 0) {
        return 1.0f;
    } else {
        return n * factorial(n - 1);
    }
}

float taylor_sin(float x, int n) {
    float result = 0.0f;
    float sign = 1.0f;
    for (int i = 0; i < n; i++) {
        result += (sign * powf(x, 2 * i + 1)) / (factorial(2 * i + 1));
        sign = -sign;
    }
    return result;
}


float taylor_cos(float x, int n) {
    float result = 0.0f;
    float sign = 1.0f;
    for (int i = 0; i < n; i++) {
        result += (sign * powf(x, 2 * i)) / (factorial(2 * i));
        sign = -sign;
  }
    return result;
}

测试流程

#include "arm_math.h"
#include "math.h"

void math_test(void)
{
  usb_printf("\r\n\r\n######################\r\n");
  HAL_Delay(0);
  int16_t angle = rand() % 0xffff - 0x7fff;
  float anglef = (float)angle / 0x7fff * PI;
  usb_printf("[input]:%.6f\r\n",anglef * 360 / 2 / PI);
  HAL_Delay(0);
  
  
  int time_base = 0,time_cost = 0;
  float sindata = 0;
  float cosdata = 0;
  Trig_Components sincos = {0};

  time_base = time.time();
  sincos = MCM_Trig_Functions(angle);
  sindata = (float)sincos.hSin / 0x7fff;
  cosdata = (float)sincos.hCos / 0x7fff;
  time_cost = time.time() - time_base;
  usb_printf("[type]:cordq15\t[time]:%fus [output]:%.6f,%.6f\r\n",(float)time_cost / 170,sindata,cosdata);
  HAL_Delay(10);
  
  time_base = time.time();
  sindata = (float)arm_sin_q15(angle >> 1) / 0x7fff;
  cosdata = (float)arm_cos_q15(angle >> 1) / 0x7fff;
  time_cost = time.time() - time_base;
  usb_printf("[type]:dspq15\t[time]:%fus [output]:%.6f,%.6f\r\n",(float)time_cost / 170,sindata,cosdata);
  HAL_Delay(10);
  
  time_base = time.time();
  sindata = (float)arm_sin_q31(angle << 15) / 0x7fffffff;
  cosdata = (float)arm_cos_q31(angle << 15) / 0x7fffffff;
  time_cost = time.time() - time_base;
  usb_printf("[type]:dspq31\t[time]:%fus [output]:%.6f,%.6f\r\n",(float)time_cost / 170,sindata,cosdata);
  HAL_Delay(10);
  
  time_base = time.time();
  sindata = arm_sin_f32(anglef);
  cosdata = arm_cos_f32(anglef);
  time_cost = time.time() - time_base;
  usb_printf("[type]:dspf32\t[time]:%fus [output]:%.6f,%.6f\r\n",(float)time_cost / 170,sindata,cosdata);
  HAL_Delay(10);
  
  time_base = time.time();
  sindata = sinf(anglef);
  cosdata = cosf(anglef);
  time_cost = time.time() - time_base;
  usb_printf("[type]:mathf\t[time]:%fus [output]:%.6f,%.6f\r\n",(float)time_cost / 170,sindata,cosdata);
  HAL_Delay(10);

  time_base = time.time();
  sindata = sin(anglef);
  cosdata = cos(anglef);
  time_cost = time.time() - time_base;
  usb_printf("[type]:mathd\t[time]:%fus [output]:%.6f,%.6f\r\n",(float)time_cost / 170,sindata,cosdata);
  HAL_Delay(10);
  
  time_base = time.time();
  sindata = taylor_sin(anglef,6);
  cosdata = taylor_cos(anglef,6);
  time_cost = time.time() - time_base;
  usb_printf("[type]:taylor6\t[time]:%fus [output]:%.6f,%.6f\r\n",(float)time_cost / 170,sindata,cosdata);
  HAL_Delay(10);
}

输出结果

######################
[input]:140.343643
[type]:cordq15	[time]:0.923529us [output]:0.638234,-0.769829
[type]:dspq15	[time]:1.441176us [output]:0.638203,-0.769860
[type]:dspq31	[time]:1.447059us [output]:0.638228,-0.769824
[type]:dspf32	[time]:1.682353us [output]:0.638170,-0.769872
[type]:mathf	[time]:2.458823us [output]:0.638182,-0.769886
[type]:mathd	[time]:36.694118us [output]:0.638182,-0.769886
[type]:taylor6	[time]:38.758823us [output]:0.638164,-0.769980


######################
[input]:67.578964
[type]:cordq15	[time]:0.923529us [output]:0.924406,0.381451
[type]:dspq15	[time]:1.382353us [output]:0.924345,0.381420
[type]:dspq31	[time]:1.552941us [output]:0.924386,0.381440
[type]:dspf32	[time]:1.629412us [output]:0.924399,0.381407
[type]:mathf	[time]:2.458823us [output]:0.924406,0.381410
[type]:mathd	[time]:30.770588us [output]:0.924406,0.381410
[type]:taylor6	[time]:39.188236us [output]:0.924406,0.381410


######################
[input]:51.236301
[type]:cordq15	[time]:0.923529us [output]:0.779717,0.626118
[type]:dspq15	[time]:1.382353us [output]:0.779626,0.626118
[type]:dspq31	[time]:1.835294us [output]:0.779711,0.626126
[type]:dspf32	[time]:1.858824us [output]:0.779728,0.626104
[type]:mathf	[time]:2.458823us [output]:0.779735,0.626110
[type]:mathd	[time]:29.305882us [output]:0.779735,0.626110
[type]:taylor6	[time]:39.223530us [output]:0.779735,0.626110


######################
[input]:-95.369728
[type]:cordq15	[time]:0.923529us [output]:-0.995575,-0.093539
[type]:dspq15	[time]:1.388235us [output]:-0.995697,-0.093692
[type]:dspq31	[time]:1.558824us [output]:-0.995599,-0.093530
[type]:dspf32	[time]:1.623529us [output]:-0.995594,-0.093581
[type]:mathf	[time]:2.505882us [output]:-0.995612,-0.093582
[type]:mathd	[time]:27.211765us [output]:-0.995612,-0.093582
[type]:taylor6	[time]:41.788235us [output]:-0.995612,-0.093583

标签：运算,三角函数,CORDIC,float,us,time,output,MCU,type
From： https://www.cnblogs.com/USTHzhanglu/p/18059299

相关文章

  • 编程语言中，差、交、并、自然连接、选择、投影、笛卡尔积分别都是什么运算...
    原文:https://blog.csdn.net/muzihuaner/article/details/119529646交(Intersection)：关系R与关系S的交由既属于R又属于S的元组组成，即R与S中相同的元组，组成一个新关系，其结果仍为n目关系。记作：R∩S={t|t∈R∧t∈S}简单来说，运算结果就是两或多个实体集所共有的部分 并(Union)：......
  • c语言四则运算小程序
    本文源程序代码来源于csdn一位博主：文章链接：http://t.csdnimg.cn/L29fs原程序简洁凝练，以简短的代码写出了一个加法运算器。以下是运行结果（本文全程用visualstudio2022进行代码的编写与运算）：通过对该代码的运算，我个人认为有以下问题：（1）尽管该代码简洁明了，但只有加法运算，缺少减......
  • php7中的三元运算符的区别
    <?php$tmparr=['cover'=>'http://img.immomo.com.cn'];echoisset($tmparr['cover'])."\n";echo!empty($tmparr['cover'])."\n";$rep=isset($tmparr['cover'])&&!empty......
  • C语言基础-1、逻辑类型和运算
    一、逻辑类型和运算#include<stdbool.h>之后就可以使用bool和true、falseex1:#include<stdio.h>#include<stdbool.h>intmain(){ boolb=6>5; boolt=true; printf("%d\n",t); t=2; printf("%d\n",t); printf("%d\n&q......
  • numpy基础运算
    numpy基础运算importnumpyasnpt1=np.array([1,2,3,4,5])#numpy数组类型为numpy.ndarrayprint("type(np.array)=",type(t1))t2=np.array(range(6))print("t1:",t1)print("t2:",t2)#np.arange([start,]stop[,stop,],dtype=......
  • 基本运算符
    Java语言支持如下运算符重点掌握的运算符算数运算符：+、-、*、/、%、++、--%(取余、模运算)：例如11/5=2..1（余1）+、-、*、/运算代码举例publicclassDemo01{publicstaticvoidmain(String[]args){inta=10;intb=20;in......
  • 4.Python3 运算符
    Python3运算符1.Python算术运算符以下假设变量a=10，变量b=21：运算符描述实例+加-两个对象相加a+b输出结果31-减-得到负数或是一个数减去另一个数a-b输出结果-11*乘-两个数相乘或是返回一个被重复若干次的字符串a*b输出结果210/除......
  • JAVA学习笔记--运算符
    运算符注意：()的优先级最高，因此可以多打一些()提高代码的可读性！！算术运算符：+、-、*、/、%(模：取余)、++(自增)、--(自减）publicclassDemo1{publicstaticvoidmain(String[]args){inta=10;intb=20;System.out.println(a+b);......
  • 基本操作之——位运算
    dev_clear_window()dev_disp_text('数128位左移两位为'+lsh(128,2),'window','center','center','black','box_color','#00ffffc0') dev_clear_window()dev_disp_text('数128位右移两位为'+rsh(12......
  • 学java11运算符
    运算符java支持以下运算符优先级()多用括号1.算术运算符:+,-,*,/,%(取余运算，模运算),++,--//++a先自增再给c赋值a=a+1//a++先给b赋值再自增a=a+1a--和--a原理同上2.赋值运算符：=3.关系运算符:>,<,>=,<=,==,(不等于)!=,instanceof4.逻辑运算符:&&（与）,||（或）,“!（）”（非）&&逻......

赞助商

阅读排行

网站主页关于我们联系我们网站地图

本网站内容转载自其他媒体，侵权联系[admin##ips99.com]。

Copyright © 2020-2023 IPS99 版权所有 IPS99