M0生成正弦波，波形很差怎么办？

只看该作者 · 2017-2-19 00:51

本帖最后由 yjmwxwx 于 2017-2-19 00:58 编辑

用的是STM32F030F4P6,在定时器中断里面实时计算正弦，频率高了分辨率不够了，波形很难看。

只看该作者 · 2017-2-19 07:40

是你的正弦波数据有误引起的。直接套用官方例程就会出现这个问题。

只看该作者 · 2017-2-19 14:34

sky.sun.zz 发表于 2017-2-19 07:40
是你的正弦波数据有误引起的。直接套用官方例程就会出现这个问题。

经过你的提醒我才想起来还有电机库，于是下载了个叫STSW-STM32100 4.3.0的文件，安装以后在文件夹里找到一个名叫STO_CORDIC_SpeednPosFdbkClass.o 的文件，凭感觉这里面有正弦波的核心算法，没学过电机库。。

这个文件有没有C语言的？

只看该作者 · 2017-2-19 16:50

直接计算三角函数自然是慢，何不查表？。如果非要计算三角函数，可以用arm提供的dsp库，比编译器自带的math库快些。

1万 · 2017-2-19 17:04

查表就行了.

1万 · 2017-2-19 17:48

用的DAC吗？分辨率是多少？你计算时候保留了多少位？

只看该作者 · 2017-2-19 18:01

icecut 发表于 2017-2-19 17:04
查表就行了.

最近在学这方面的东西，查表法固定频率的学会了，但是怎么改变频率啊？比如步长0.1HZ ？

只看该作者 · 2017-2-19 18:13

本帖最后由 yjmwxwx 于 2017-2-19 18:18 编辑

1223657347 发表于 2017-2-19 16:50
直接计算三角函数自然是慢，何不查表？。如果非要计算三角函数，可以用arm提供的dsp库，比编译器自带的math ...

在网上抄了个CORDIC算法，但是算出来数很大放不到定时器里，非要用乘除法，用移位法计算乘除频率步长太长了。
M0没硬件除法，算一个除法时间比算正弦还慢。没办法降低精度结果波形很难看。还有个毛病定时器中断里面计算，如果小于10跟本来不及计算，中断速度太快了

只看该作者 · 2017-2-19 18:17

734774645 发表于 2017-2-19 17:48
用的DAC吗？分辨率是多少？你计算时候保留了多少位？

用的定时器，一个周期最少50个点

只看该作者 · 2017-2-19 18:33

复制

//Cordic in 16 bit signed fixed point math

//Function is valid for arguments in range -pi/2 -- pi/2

//for values pi/2--pi: value = half_pi-(theta-half_pi) and similarly for values -pi---pi/2

//

// 1.0 = 16384

// 1/k = 0.6072529350088812561694

// pi = 3.1415926536897932384626

//Constants

#define cordic_1K 0x000026DD

#define half_pi 0x00006487

#define MUL 16384.000000

#define CORDIC_NTAB 16

int cordic_ctab [] = {0x00003243, 0x00001DAC, 0x00000FAD, 0x000007F5, 0x000003FE, 0x000001FF, 0x000000FF, 0x0000007F, 0x0000003F, 0x0000001F, 0x0000000F, 0x00000007, 0x00000003, 0x00000001, 0x00000000, 0x00000000, };



void cordic(int theta, int *s, int *c, int n)

{

  int k, d, tx, ty, tz;

  int x=cordic_1K,y=0,z=theta;

  n = (n>CORDIC_NTAB) ? CORDIC_NTAB : n;

  for (k=0; k<n; ++k)

  {

    d = z>>15;

    //get sign. for other architectures, you might want to use the more portable version

    //d = z>=0 ? 0 : -1;

    tx = x - (((y>>k) ^ d) - d);

    ty = y + (((x>>k) ^ d) - d);

    tz = z - ((cordic_ctab[k] ^ d) - d);

    x = tx; y = ty; z = tz;

  }  

 *c = x; *s = y;

}