DSP算法大全C语言版本
一个简单的dspC语言例子
一个简单的dsp C语言例子开发平台: CCS集成开发环境通过这个简单的例子, 可以大致了解用C语言开发dsp程序的原理。
程序要求: 用C语言编写产生正弦调幅波信号的源程序；正弦调幅波的公式在离散域中的表示:y(n) = (1 + M*sin(2 * PI * fb / fs * n)) * sin(2 * PI * fa / fs * n);编写文件1.sin_am.c#include<stdio.h>#include<math.h>#define TRUE 1#define pi 3.1415926536int y[500],i;float M;void main(){puts("amplitude modulation sinewave example started.\n");M = 50;for(i = 0; i < 500; i++)y[i]= 0;while(TRUE){for(i = 0; i < 500; i++)y[i]=(int)((1 + M / 100 * sin(i * 2 * pi * 20 / 4000))* sin(i * 2 * pi * 200 / 4000)* 16384);puts("program end");}}2.sin_am_v.asm (reset vector file).title "sin_am_v.asm".sect ".vectors".ref _c_int00RESET:B _c_int00.end..3.sin_am.cmdsin_am.objsin_am_v.obj-m sin_am.map-o sin_am.outMEMORY{PAGE 0:EPROG: origin = 0x1400, len = 0x7c00 VECT: origin = 0xff80, len = 0x80PAGE 1:USERREGS: origin = 0x60, len = 0x1c IDATA: origin = 0x80, len = 0x3000 }SECTIONS{.vectors:>VECT PAGE 0.text:>EPROG PAGE 0.cinit:>EPROG PAGE 0.bss:>IDATA PAGE 1.const:>IDATA PAGE 1.switch:>IDATA PAGE 1.system:>IDATA PAGE 1.stack:>IDATA PAGE 1}"*.cmd"文件说明:链接命令文件是实现对段的存储空间位置的定位, C语言程序中常用已初始化和未初始化段如下:已初始化段包括:.init 存放C程序中的变量的初值和常量, 放在ROM和RAM 中均可, 一般属于PAGE 0.const 存放C程序中的字符常量、浮点常量和用const声明的常量, 放在ROM和RAM中均可, 一般属于PAGE 1.text 存放C程序代码, 放在ROM和RAM中均可, 一般属于PAGE 0.switch 存放C程序中的语句的跳针表, 放在ROM和RAM中均可, 一般属于PAGE 0未初始化段包括:.bss 为C程序中的全局和静态变量保留存储空间, 一般存放于RAM中, 属于PAGE 1.stack 为C程序系统堆栈保留存储空间, 用于保存返回地址、函数间的参数传递、存储局部变量和保存中间结果, 一般存放于RAM中, 属于PAGE 1.sysmem 用于C程序中malloc、calloc和realloc函数动态分配存储空间, 一般存放于RAM中, 属于PAGE 14.vary_M.gelmenuitem "Myfunctions"slider vary_M(0, 100, 10, 1, Amount_of_modulation){M = Amount_of_modulation;}该文件用于调试的时候可随意改变变量M的值, 该文件通过file->load GEL File添加到工程中, 调试的时候可选择GEL->My Functions->vary_M来打开vary_M滑动条组件。
嵌入式C语言开发ADSP21系列DSP
嵌入式C语言开发ADSP21系列DSP嵌入式C语言开发ADSP21XX系列DSP摘要详细介绍使用VisualDSP开发工具进行ADSP21XX的C语言编程的方法；分析其C语言运行库的结构，并且结合实例介绍C语言工具的使用方法，包括C语言与汇编语言混合编程的方法，从C运行库提取代码用于自己的汇编语言程序的方法、修改运行库的源代码以适应自己开发需要的方法等。
关键词DSPVisualDSP嵌入式C语言汇编语言引言长期以来，在DSP系统开发中，一直把汇编语言作为主要的开发工具；但汇编语言与自然语言差距很大，不易常，而且汇编语言是依赖于处理器的，不利于软件的可重复利用和系统的稳定性，程序不易移植，给开发工作带来了很大的困难。
随着嵌入式系统复杂程度的不断提高，用汇编语言编写一个巨大的程度将是困难，甚至是不可能的。
为此，AD公司推出了针对ADSP21XX系列DSP的嵌入式C和C++语言集成开发工具，分别是VisualDSP和VisualDSP++系列，这些开发工具提供了C语言和C++语音的开发功能。
以下就以笔者在实际开发中的一些经验，结合VisualDSP6.1版本，介绍用C语言开发VisualDSP6.1版本，介绍用C语言开发ADSP21XX的方法。
VisualDSP提供了一个开放源码软件组织GNU的C编译器，和一套成熟稳定的C运行时间库（CRuntimeLibrary）等。
GNU 的编译器一向以编译效率高著称，在编译后的代码长度和运行速度方面非常优秀；C运行时间库则把很多重复性的工作，如浮点运行、三角函数、FFT等作为C语言的库函数，提供给用户，大大提高了用户的开发效率和程序的稳定性，降低了开发难度，另外，由于把这些库函数的源代码提供给了用户，还提高了C语言与汇编语言之间的透明性，使用户开发的程序兼具两者的优点。
1VisualDSP简介VisualDSP是AD公司的DSP开发工具，主要由可执行文件、库文件和各种帮助文档组成。
DSP_c2000_C语言编程要点
例 3、DSP 程序中比较通用的命名方法
1） 变量、常量和对象命名多用名词或形容词+名词,比如 ADC 的寄存器就可以命名为 “AdcRegs”。
例 2、Q 格式
由于 C2000 系列是定点 DSP，遇到小数时不能直接处理，所以就要想一些办法，通常 DSP 用 Q 格式来表示。当假定小数点位于第 0 位的右侧时，为 Q0；当小数电位于第 15 位的 右侧时为 Q15；若是 32 位处理器则最大可为 Q31。
浮点数和定点数转换的公示：
X 为浮点数， X 为定点数。
WDCR=0x0068; WDKEY=0X0055; WDKEY=0X00aa; printf(“hello,DSP world\n”); for(;;) ; } 首先，要禁止看门狗 watchdog。DSP 默认上电后看门狗会自动启动，时间一到就会令 DSP 复位。 其次，DSP 的 C 需要绝对定位，而主机的 C 的代码有操作系统定位。在“printf（“hello， DSP world\n”）；”语句后，需要加设死循环语句“for（；；） ；”。标准 C 运行在操作系统下， 退出 main（）函数后，控制权会交给操作系统。而一般 DSP 没有操作系统，退出 main（） 函数即意味着程序跑飞了，所以要在主函数里加一个死循环。 在程序开头是通过宏定义去表示 watchdog 的寄存器，因为例子中用到的 DSP 芯片是 F240 系列的，而其 WDCR 和 WDKEYD7 寄存器的地址分别是 0x7029h 和 0x7025h，（volitial unsigned int*）表明 0x7029 是一个可被修改的无符号整型指针变量。
DSP的C语言编程
72 71 70 69 68 67 66 65 64 63 62 61 60 59 58 57 56 55 54 53 52 51 50 49 48 47 46 45 44 43 42 41 40 39 38 37
DG ND B C L KX 1 DG ND HD 1 C V DD IN T3 IN T2 IN T1 IN T0 NM I HB IL IA C K B D X2 B D X0 HD 0 DG ND 3 .3 V HR DY B F SX2 B F SX0 C V DD HP INT DG ND B C LKX 2 B C LKX 0 BD R2 HC T1 BD R0 B F SR 2 B F SR 0 B C LKR 2 B C LKR 0 DG ND HC T0 B C L KR 1 DG ND
DSP的C语言程序入门
• 1、以最小系统板的LED灯闪烁程 序为例： 分别利用汇编语言、C语言、C语 言和汇编语言混合编程三种方法 实现DSP最小系统板LED灯D1闪烁。
C L KOU T TO UT
1 2
J1 C O N2 R2 4 .7 K DG ND
C LK OU T X1 H D3
D2 R ES ET D1 C LD 0IN K
H D7 A 11
A 10
A 12
A 13
A 14
A 15
NC
NC 36
1
2
3
4
5
6
7
8
9 10
11
12 13
14 15
16
17
18
19
20
21
22 23
24 25
26
27
28
29
第9讲 DSP的C及汇编语言
OnFileLoaded()函数
GEL文件—回调函数
该回调函数在加载program/symbol(.out)文件之后执行。 可以建立调试源搜索路径（在没有CCS工程文件的时候）， 设定断点和探针，完成软件的复位和重启。
GEL文件—回调函数
OnReset()函数
当目标处理器复位后该函数被调用。如果你需要每次重
无论用哪种方法定义中断服务函数，都须注意以下问题：
（1）中断处理函数必须是void类型，而且不能有任何
输入参数。
（2）进入中断服务函数，编译器将自动产生程序保护 所有必要的寄存器，并在中断服务函数结束时恢复运行 环境。 （3）进入中断服务函数，编译器只保护与运行上下文 相关的寄存器，而不是保护所有的寄存器。中断服务函 数可以任意修改不被保护的寄存器，如外设控制寄存器 等。 （4）要注意IMR、INTM等中断控制量的设置。通常 进入中断服务程序要设置相应寄存器将中断屏蔽，退出 中断服务程序时再打开，避免中断嵌套。
需要注意的问题
(2)指针类型转换是必须的，关系到接收分配 好的内存块的地址可以向前看多少字节。 如果不做指针类型转换会怎么样？ void * 指针存在的意义。
需要注意的问题
(3)内存块大小的可移植性问题
分配一个整型变量数组应使用： 数组元素个数 * sizeof(int) 确定内存块的大小 问题：sizeof和strlen函数的区别
/*进行CCS存储器映射，告知目标处理器哪些空间可以访问，哪些不可以访问。*/
Setup_Memory_Map( ); } OnTargetConnect( )//对处理器进行最小初始化 { Setup_Cache( ); //设置缓存L1P,L1D,L2 Setup_Pin_Mux( ); // 设置管脚 Setup_Psc_All_On( ); // 设置psc Setup_PLL0_594_MHz_OscIn( ); // 设置dsp主频[DSP @ 594 MHz][Core 1.20V] Setup_PLL1_DDR_135_MHz_OscIn();//设置ddr时钟频率 Setup_Aemif_8Bit_Bus( ); // 设置Async-EMIF[8-bit bus] }
DSP-常用算法代码
一、32 位加法//链接命令文件MEMORY{PAGE 0: ROM :origin=0080h,length=1000hROM :origin=0060h,length=10hPAGE 1: OTHER :origin=0400h,length=40h}SECTIONS{.text : { }>ROM PAGE 0.data :{ }>ROM PAGE 0.bss : { }>OTHER PAGE 1.stack : { }>OTHER PAGE 1}//汇编语言源代码.title "ADD32".mmregs.def start,_c_int00.bss xhi , 2,1,1.bss yhi , 2,1,1.bss zhi , 2,1, 1table .long 13578468H.long 1020B30AH.text_c_int00b startnopnopstart: LD #xhi ,DPSTM #xhi,AR1RPT #3MVPD table,*AR1+DLD xhi , ADADD yhi , ADST A ,zhiEND: B END.end二、64位加法/减法//链接命令文件/*ADDSUB64.cmd*/ADDSUB64.obj-m ADDSUB64.map-o ADDSUB64.outMEMORY{PAGE 0: ROM :origin=0080h,length=1000hROM :origin=0060h,length=10hPAGE 1: OTHER :origin=0400h,length=40h}SECTIONS{.text : { }>ROM PAGE 0.data :{ }>ROM PAGE 0.bss : { }>OTHER PAGE 1.stack : { }>OTHER PAGE 1}//汇编语言源代码.title "ADDSUB64".mmregs.def start,_c_int00.bss x1 , 2 ,1,1.bss x3 , 2 ,1,1.bss y1 , 2 ,1,1.bss y3 , 1.bss y2 , 1.bss z1 , 2 ,1,1.bss z3 , 1.bss z2 , 1.bss w1 , 2 ,1,1.bss w3 , 2 ,1,1table .long 12345678H ;x1x0.long 02468ACEH ;x3x2.long 22222222H ;y1y0.word 1357H,2468H ;y3,y2.long 44444444H ;z1z0.word 1020H,0B30AH ;z3,z2.text_c_int00b startnopnopstart: LD #x1 ,DPSTM #x1,AR1RPT #11MVPD table,*AR1+DLD x1,A ; A = X1 X0DADD y1,A ; A= X1 X0 + Y1 Y0,产生进位CDLD x3,B ; B = X3 X2ADDC y2,B ; B = X3 X2 + 00 Y2 + CADD y3,16,B ; B = X3 X2 + Y3 Y2 + CDSUB z1,A ; A= X1 X0 + Y1 Y0-Z1Z0，产生借位C'DST A,w1 ; W1W0= X1 X0 + Y1 Y0-Z1Z0SUBB z2,B ; B = X3 X2 + Y3 Y2 + C - 00 Z2-C'SUB z3,16,B ; B = X3 X2 + Y3 Y2 + C -Z3 Z2-C'DST B,w3 ; W3 W2= X3 X2 + Y3 Y2 + C - Z3 Z2-C' END: B END.end三、FIR滤波//链接命令文件fir.obj-m fir.map-o fir.outMEMORY{PAGE 0: ROM1(RIX) :ORIGIN=0080H,LENGTH=100HPAGE 1: INTRAM1(RW) :ORIGIN=2400H,LENGTH=0200HINTRAM2(RW) :ORIGIN=2600H,LENGTH=0100HINTRAM3(RW) :ORIGIN=2700H,LENGTH=0100HB2B(RW) :ORIGIN=0070H,LENGTH=10H}SECTIONS{.text : {}>ROM1 PAGE 0.data : {}>INTRAM1 PAGE 1FIR_COFF: {}>INTRAM2 PAGE 1FIR_BFR : {}>INTRAM3 PAGE 1.stack : {}>B2B PAGE 1}//汇编语言源代码;一个FIR滤波器源程序fir.asm.mmregs.global start.def start,_c_int00INDEX .set 1KS .set 256 ;模拟输入数据缓冲区大小N .set 17COFF_FIR .sect "COFF_FIR" ;FIR滤波器系数.word 0.word 158.word 264.word -290.word -1406.word -951.word 3187.word 9287.word 12272.word 9287.word 3187.word -951.word -1406.word -290.word 260.word 158.word 0.dataINPUT .copy "firin.inc" ;模拟输入在数据存储区0x2400 OUTPUT .space 1024 ;输出数据在数据区0x2500 ;FIR_DP .usect "FIR_V ARS",0;D_FIN .usect "FIR_V ARS",1;D_FOUT .usect "FIR_V ARS",1COFFTAB .usect "FIR_COFF",NDA TABUF .usect "FIR_BFR",NBOS .usect "STACK",0FhTOS .usect "STACK",1.text.asg AR0,INDEX_P.asg AR4,DA TA_P ;输入数据x(n)循环缓冲区指针.asg AR5,COFF_P ;FIR系数表指针.asg AR6,INBUF_P ;模拟输入数据指针.asg AR7,OUTBUF_P;FIR滤波器输出数据指针_c_int00b startnopnopstart: ssbx FRCTSTM #COFFTAB,COFF_PRPT #N-1 ;将FIR系数从程序存储器移动MVPD #COFF_FIR,*COFF_P+ ;到数据存储器STM #INDEX,INDEX_PSTM #DATABUF,DATA_PRPTZ A,#N-1STL A,*DATA_P+ ;将数据循环缓冲区清零STM #(DA TABUF+N-1),DATA_P ;数据缓冲区指针指向x[n-(N-1)]STM #COFFTAB,COFF_P ;FIR_TASK:STM #INPUT,INBUF_PSTM #OUTPUT,OUTBUF_PSTM #KS-1,BRCRPTBD LOOP-1STM #N,BK ;FIR循环缓冲区大小LD *INBUF_P+,A ;装载输入数据FIR_FILTER:STL A,*DATA_P+%RPTZ A,N-1MAC *DATA_P+0%,*COFF_P+0%,ASTH A,*OUTBUF_P+LOOP:EEND B EEND.end四、IIR滤波//链接命令文件MEMORY{PAGE 0:EPROM: org=0E00H, len=1000HVECS0: org=0FF80H, len=0080hPAGE 1:SPRAM: org=0060H, len=0020hDARAM: org=0080H, len=1380H}SECTIONS{.text :>EPROM PAGE 0.data :>EPROM PAGE 0X: align(8) {}>DARAM PAGE 1Y: align(8) {}>DARAM PAGE 1B: align(8) {}>DARAM PAGE 1A: align(8) {}>DARAM PAGE 1.vecs0:>VECS0, PAGE 0}//输入数据序列(512点)0x7fff0x00x0…0x0//汇编语言源代码.title "IIR.ASM".mmregs.def _c_int00X .usect "X", 5Y .usect "Y", 5B .usect "B", 5A .usect "A", 5PA0 .set 0PA1 .set 1.datatable: .word 0.word 0.word 0.word 0.word 0.word 0.word 0.word 0.word 3116.word -10286.word 14615.word -10286.word 3116.word -22082.word 31149.word -30484.word 28383.text_c_int00:SSBX FRCTSTM #X, AR1RPT #3MVPD #table, *AR1+STM #Y, AR1RPT #3MVPD #table+4, *AR1+STM #B, AR1RPT #4MVPD #table+8, *AR1+STM #A, AR1RPT #3MVPD #table+13, *AR1+STM #X+4, AR2STM #A+3, AR3STM #Y+3, AR4STM #B+4, AR5STM #5, BKSTM #-1, AR0STM #1000H, AR6STM #0200H-1, AR7loop:PORTR PA1, *AR2LD *AR2, ASTL A, -1, *AR2MPY *AR2+0%, *AR5+0%, AMAC *AR2+0%, *AR5+0%, AMAC *AR2+0%, *AR5+0%, AMAC *AR2+0%, *AR5+0%, AMAC *AR2, *AR5+0%, AMAC *AR4, *AR3, AMAC *AR4, *AR3, AMAC *AR4, *AR3, AMAC *AR4+0%, *AR3+0%, AMAC *AR4, *AR3, AMAC *AR4, *AR3, AMAC *AR4, *AR3, AMAC *AR4, *AR3, AMAC *AR4+0%, *AR3+0%, AMAC *AR4, *AR3, AMAC *AR4, *AR3, AMAC *AR4+0%, *AR3+0%, AMAC *AR4+0%, *AR3+0%, AMAR *AR3+0%STH A, *AR4PORTW *AR4, PA0STH A, *AR6+BANZ loop, *AR7-end: B end.end五、除法//链接命令文件/*chuf.cmd*/chuf.obj-m chuf.map-o chuf.outMEMORY{PAGE 0: ROM :origin=0080h,length=1000hROM :origin=0060h,length=10hPAGE 1: OTHER :origin=0400h,length=40h}SECTIONS{.text : {}>ROM PAGE 0.data :{}>ROM PAGE 0.stack : {}>OTHER PAGE 1.bss : {}>OTHER PAGE 1}//汇编语言源代码;*** 编制计算除法运算的程序段。
DSP的C语言编程
DSP的C语言开发一、流程步骤：main() {}下面是vectors.asm函数，该文件在每个DSP的project中（需手工加入），其中有对_c_int00的调用，而_c_int00在rts.lib中，在开发时要手工加入。
在reset 后，rom等外存中的程序已经转移到了L2 cache中，并且程序从0x0000 0000处开始执行，而0x0000 0000处的程序正是vectors.asm，以下就开始层层调用，进入main函数。
======== vectors.asm ========; Plug in the entry point at RESET in the interrupt vector table;;; ======== unused ========; plug inifinite loop -- with nested branches to; disable interrupts -- for all undefined vectors;unused .macro id.global unused:id:unused:id:b unused:id: ; nested branches to block interruptsnop 4b unused:id:nopnopnopnopnop.endm.sect ".vectors".ref _c_int00 ; C entry point.align 32*8*4 ; must be aligned on 256 word boundaryRESET: ; reset vectormvkl _c_int00,b0 ; load destination function address to b0mvkh _c_int00,b0b b0 ; start branch to destination functionmvc PCE1,b0 ; address of interrupt vectorsmvc b0,ISTP ; set table to point herenop 3 ; fill delay slotnopnop;; plug unused interrupts with infinite loops to; catch stray interrupts;unused 1unused 2unused 3unused 4unused 5unused 6unused 7unused 8unused 9unused 10unused 11unused 12unused 13unused 14unused 15Rts6000.lib来自于rts6000.src，该原文件是由多个.c和.cpp以及.asm组成的，其中关于int _args_main()的函数：这个函数就是void __interrupt c_int00()在初始化完成后调用的函数，int _args_main()函数中调用了main(argc, argv)从而正式转入main函数。
dspC54x的C语言程序设计
2020年2月22日
DSP原理及应用
11
第6章 TMS320C54xC语言程序设计
2.数据出输入输出设备不同
面向PC机的ANSI C所处理的数据来自硬盘、内存等 存储介质，或者是键盘的输入。其输出是显示器、通信口 或者实时输出。对人机界面提出了较高的要求。
本章重点介绍DSP的C语言程序设计的特点、DSP C语 言编程的设计技巧、C与汇编的混合编程实现等
2020年2月22日
DSP原理及应用
1
第6章 TMS320C54xC语言程序设计
6.1 C54x C语言介绍 6.2 C54x C语言编程 6.3 C54x C代码优化
2020年2月22日
DSP原理及应用
21
第6章 TMS320C54xC语言程序设计
3）寄存器变量 DSP的C语言扩展了寄存器变量。在一个C函数内部最
多可以使用两个寄存器变量，编译器在编译这两个寄存器 变量时用AR6和AR7这两个辅助寄存器分别做了对应：
★AR6被赋给第一个寄存器变量 ★AR7被赋给了第二个寄存器变量
寄存器变量的定义格式： register type AR6 或 register type AR7
23
第6章 TMS320C54xC语言程序设计
例如： int sum; int a[8] ={1,2,3,4,5,6,7,8}; for(i=0;i<8;i++) sum=sum+a[i];
在PC机上运行的结果为sum=36。
注意：在DSP上运行，结果可能不是sum=36。 因为变量sum没有初始化，初始值是未知的。