dsp输出正弦波 c语言,快速计算正弦波 - DSP 与 ARM 处理器 - 电子工程世界-论坛...

快速计算正弦波

[复制链接]

在DSP运用中，经常需要产生正弦波。如果直接用c的数学函数sin，当然可以产生正弦波，但是由于sin函数本身的效率很低，产生正弦波所需要的MIPS就会占去DSP处理能力的相当大的一部分。本章介用递推数列算正弦波的方法，先介绍原理，推导出递推公式，然后用浮点小数实现计算，再用定点小数进一步优化算法，最后进行误差分析，并提出更精确的定点小数算法。

先来看看如何推导出递推数列的公式。

我们所要产生的正弦波，其实是一系列的整数，把这些整数按照一定的取样频率发送给数模转换器，就可以变成真正的正弦波了。假设取样周期是Ts，产生的正弦波的圆频率为w，那么我们需要产生的数列就是：

sin(0), sin(w*Ts), sin(2*w*Ts), ... sin(n*w*Ts)

假设f(n)= sin(n*w*Ts)，则问题就变成，从f(n-1), f(n-2), f(n-3),..., 如何计算f(n)了。解决了这个问题，也就找到了递推公式。

下面是这个递推公式的求解过程，假设x=w*Ts：

公式：sin( a + b) = sin(a)*cos(b) + cos(a)*sin(b)

sin( x + (n-1)x ) = sin(x)*cos( (n-1)x ) + cos(x)*sin( (n-1)x )

公式：Sin(a)*cos(b) = 1/2 * [ sin( a+b ) + sin( a-b ) ]

sin(x)*cos( (n-1)x ) = 1/2 * [ sin