dsp算法大全c语言版本

‘壹’ DSP学习入门

DSP学习入门：
DSP的特点 ：
对于没有使用过DSP的初学者来说,第一个困惑就是DSP其他的嵌入式处理器究竟有什么不同,它和单片机,ARM有什么区别.事实上,DSP也是一种嵌入式处理器,它完全可以完成单片机的功能.。
唯一的重要的区别在于DSP支持单时钟周期的"乘-加"运算.这几乎是所有厂家的DSP芯片的一个共有特征.几乎所有的DSP处理器的指令集中都会有一条MAC指令,这条指令可以把两个操作数从RAM中取出相乘,然后加到一个累加器中,所有这些操作都在一个时钟周期内完成.拥有这样一条指令的处理器就具备了
DSP功能：
具有这条指令就称之为数字信号处理器的原因在于,所有的数字信号处理算法中最为常见的算术操作就是"乘-加".这是因为数字信号处理中大量使用了内积,或称"点积"的运算.无论是FIR滤波,FFT,信号相关,数字混频,下变频.所有这些数字信号处理的运算经常是将输入信号与一个系数表或者与一个本地参考信号相乘然后积分(累加),这就表现为将两个向量(或称序列)进行点积,在编程上就变成将输入的采样放在一个循环buffer里,本地的系数表或参考信号也放在一个buffer里,然后使用两个指针指向这两个buffer.这样就可以在一个loop里面使用一个MAC指令将二者进行点积运算.这样的点积运算对与处理器来说是最快的,因为仅需一个始终周期就可以完成一次乘加.
了解DSP的这一特点后,当我们设计一个嵌入式系统时,首先要考虑处理器所实现的算法中是否有点积运算,即是否要经常进行两个数组的乘加,(记住数字滤波,相关等都表现为两个数组的点积)如果有的话,每秒要做多少次,这样就能够决定是否采用DSP,采用多高性能的DSP了.
浮点与定点 ：
浮点与定点也是经常是初学者困惑的问题,在选择DSP器件的时候,是采用浮点还是采用定点,如果用定点是16位还是32位?其实这个问题和你的算法所要求的信号的动态范围有关.
定点的计算不过是把一个数据当作整数来处理,通常AD采样来的都是整数,这个数相对于真实的模拟信号有一个刻度因子,大家都知道用一个16位的AD去采样一个0到5V的信号,那么AD输出的整数除以2^16再乘以5V就是对应的电压.在定点DSP中是直接对这个16位的采样进行处理,并不将它转换成以小数表示的电压,因为定点DSP无法以足够的精度表示一个小数,它只能对整数进行计算.
而浮点DSP的优势在于它可以把这个采样得到的整数转换成小数表示的电压,并不损失精度(这个小数用科学记数法来表示),原因在于科学记数法可以表示很大的动态范围的一个信号,以IEEE754浮点数为例,
单精度浮点格式: [31] 1位符号 [30-23]8位指数 [22-00]23位小数
这样的能表示的最小的数是+-2^-149,最大的数是+-(2-2^23)*2^127.动态范围为20*log(最大的数/最小的数)=1667.6dB 这样大的动态范围使得我们在编程的时候几乎不必考虑乘法和累加的溢出,而如果使用定点处理器编程,对计算结果进行舍入和移位则是家常便饭,这在一定程度上会损失是精度.原因在于定点处理处理的信号的动态范围有限,比如16位定点DSP,可以表示整数范围为1-65536,其动态范围为20*log(65536/1)=96dB.对于32定点DSP,动态范围为20*log(2^32/1)=192dB,远小于32位ieee浮点数的1667.6dB,但是,实际上192dB对绝大多数应用所处理的信号已经足够了。
由于AD转换器的位数限制,一般输入信号的动态范围都比较小,但在DSP的信号处理中,由于点积运算会使中间节点信号的动态范围增加,所以主要考虑信号处理流程中中间结果的动态范围,以及算法对中间结果的精度要求,来选择相应的DSP.另外就是浮点的DSP更易于编程,定点DSP编程中程序员要不断调整中间结果的P,Q值,实际就是不断对中间结果进行移位调整和舍入.。
DSP与RTOS：
TI的CCS提供BIOS,ADI的VDSP提供VDK,都是基于各自DSP的嵌入式多任务内核.DSP编程可以用单用C,也可以用汇编,或者二者结合,一般软件编译工具都提供了很好的支持.我不想在这里多说BIOS,VDK怎么用这在相应的文档里说的很详细.我想给初学者说说DSP的RTOS原理.用短短几段话说这个复杂的东西也是挑战!
其实DSP的RTOS和基于其他处理器的通用RTOS没什么大的区别,现在几乎人人皆知的uCOSii也很容易移植到DSP上来,只要把寄存器保存与恢复部分和堆栈部分改改就可以.一般在用BIOS和VDK之前,先看看操作系统原理的书比较好.uCOS那本书也不错.
BIOS和VDK其实是一个RTOS内核函数集,DSP的应用程序会和这些函数连接成一个可执行文件.其实实现一个简单的多任务内核并不复杂,首先定义好内核的各种数据结构,然后写一个scheler函数,功能是从所有就绪任务中(通过查找就绪任务队列或就绪任务表)找出优先级最高的任务,并恢复其执行.然后在此基础上写几个用于任务间通信的函数就可以了,比如event,message box,等等.
RTOS一般采用抢先式的任务调度方式,举例说当任务A等待的资源available的时候,DSP会执行一个任务调度函数scheler,这个函数会检查当前任务是否比任务A优先级低,如果是的话,就会把它当前挂起,然后把任务A保存在堆栈里寄存器值全部pop到DSP处理器中(这就是所谓的任务现场恢复).接着scheler还会把从堆栈中取出任务A挂起时的程序执行的地址,pop到PC,使任务A继续执行.这样当前任务就被任务A抢先了.
使用RTOS之后,每个任务都会有一个主函数,这个函数的起始地址就是该任务的入口.一般每个任务的主函数里有一个死循环,这个循环使该任务周期地执行,完成一部分算法模块的功能,其实这个函数跟普通函数没任何区别,类似于C语言中的main函数.一个任务创建的时候,RTOS会把这个函数入口地址压入任务的堆栈中,好象这个函数(任务)刚发生过一次中断一样.一旦这个新创建任务的优先级在就绪队列中是最高的,RTOS就会从其堆栈中弹出其入口地址开始执行.
有一个疑问是,不使用RTOS,而是简单使用一个主循环在程序中调用各个函数模块,一样可以实现软件的调度执行.那么,这种常用的方法与使用RTOS相比有什么区别呢?其实,使用主循环的方法不过是一种没有优先级的顺序执行的调度策略而已.这种方法的缺点在于,主循环中调用的各个函数是顺序执行的,那么,即使是一个无关紧要的函数(比如闪烁一个LED),只要他不主动返回,也会一直执行直到结束,这时,如果发生一个重要的事件(比如DMA buffer full 中断),就会得不到及时的响应和处理,只能等到那个闪烁LED的函数执行完毕.这样就使整个DSP处理的优先次序十分不合理.而在使用了RTOS之后,当一个重要的事件发生时,中断处理会进入RTOS,并调用scheler,这时scheler 会让处理这一事件的任务抢占DSP处理器(因为它的优先级高).而哪个闪烁LED任务即使晚执行几毫秒都没任何影响.这样整个DSP的调度策略就十分合理。

‘贰’ DSP的算法移植问题

移植：
1）如果你的算法是基本opencv这样的基本上开发的，你需要脱离opencv的环境。

2）如果你的算法是C++语言，请你改成标准的C语言。虽然DSP的开发环境是支持C++的，但是不建议你这么做。

3）修改你算法的内存分配，尽量内存一次分配好，DSP在算法不断的申请和释放时会有隐患。优先使用静态数组，会减轻很多工作量。

4）在CCS下建立工程，来调试你的算法，内存分配函数需要使用TI提供的函数。如果你的算法能够长期稳定的运行，那么恭喜你，你的算法移植就完成了。

优化：
算法优化，需要你能懂算法，也懂DSP。如果你只会写DSP程序，而不会算法，这对整个产品来说，是不能达到最优的。有些公司怕算法泄密，给优化人员一段或几段程序让其优化。我觉得这样做是很不合理的。除非你自己能控制大局，精通优化，这样才可行。

1）你需要对算法原理做一个深刻苦的理解，阅读相关的文章。

2）对你拿到的算法做全方位的熟悉。

3）做好上面的准备工作后，你要对算法的结构做重新的整理。依据DSP的特点，比如内存的分布。算法结构调整完成后，你的算法在DSP上速度应该有一个明显的提高了。

4）结构调整完成后，找到算法中比较费时的部分。确定我们需要优化的重点，这部分内容多是每张图像都要处理一次或多次的部分。对于算法启动时初始化部分的内容，一般不需要优化。

5）确定优化内容后，你首先考虑从语言结构上去做优化，这个时候应该还是C语言的。我不建议大家用TI提供的在C语言中使用优化嵌入的C库函数。

6）你把需要优化的函数改写为线性汇编或汇编函数。不断的调整软件流水，提高速率。

‘叁’ DSP实验 让我用C语言编写程序完成计算sin(2.3π)+cos(1.7π)的值

sin(2.3*pi)+cos(1.7*pi)
娘的，楼下的你当我白痴啊，你不会用快速傅里叶变换啊，计算点数越大就越省时间，1024点计算，比你得到的那个时间起码要缩短一半以上，还有很多改进型的FFT，更省时间

‘肆’ DSP库函数cfft32_SCALE、cifft32_NOSCALE、cbrev32分别是干什么的，最好能给出C或C++语言实现

以上DSP的函数虽然是汇编的，但在库函数的目录中有相关文档，其中详细解释了C语言的调用规则，还有算法原理，大部分会有相应C代码。你可以去目录中查一下。
汇编是针对DSP进行过优化的，目的是为了实现效率最大化。

‘伍’ 请问DSP编程软件是什么，，应该用什么编程语言

DSP的编程软件是CCS开发平台，编程语言一般是C语言。

CCS有两种工作模式：

1、软件仿真器模式：可与DSP芯片分离，并在PC上模拟DSP指令集和工作机制，主要用于早期算法的实现和调试。

2、硬件在线编程模式：可以在DSP芯片上实时运行，结合硬件开发板的在线编程和调试应用程序。

CCS的开发系统主要由以下组件构成：

1、TMS320C54x集成代码生成工具。

2、CCS集成开发环境。

3、DSP / BIOS实时内核插件及其应用程序接口API。

4、RTDX插件，用于实时数据交换和相应的程序接口API。

5、TI以外的第三方提供的各种应用程序模块插件。

(5)dsp算法大全c语言版本扩展阅读：

dsp编程中双重循环、多重循环的优化：

1、将多个循环分成单层循环以减少循环数；

例如，在双循环中，一个周期仅使用一个乘法器。 分成单层环路后，两个乘法器可以使用一个周期，充分利用DSP乘法器资源，运算速度也将提高。

2、循环次数少的放在外层循环，循环次数多的放在内存循环；

3、二维数组的双循环：二维数组的行循环置于外循环，列循环置于内循环；

4、避免在循环内进行乘法和除法运算：将循环内的乘法和除法运算尽可能移至循环外，并用加法代替。