DSP C语言基础要点

年初到现在,接触DSP已经半年了,由于公司没有人指导,做了这么久都没出什么大的成果,也走了不少的弯路。以前做单片机用C编过一些程序,个人觉得还可以。现在看来也只是在吃以前在学校里面一些C语言皮毛的老本,加上一些网上说的所谓的技巧什么的结合以前汇编的经验做了一些针对处理器的优化而已,实在“鸡肋”得很。很多东西都没有接触应用到。
   直到接触到DSP,由于使用TI的库文件例程的缘故,对结构体,联合体等有了进一步的了解,也对一个工程变量的使用有了更深的认识。下面对谈谈对DSP C刚入门者的一些建议:
   1.DSP程序的定位配置---CMD文件,要求熟悉DSP的存储器结构。这一部分对接触硬件多一些的新手很容易,其实也就是配置一般的C语言编译器的一部分,比如AVR单片机的头文件,cpu选择等在DSP中就是这一部分完成。另外还要注意TI例程中对加密位,看门狗的处理。

   2.当然是C语言的一些基础了:运算符,优先级,几个基本流程控制语句。这些无论是什么C平台运行都一样的,基础中的基础。但是也很容易因为这些基础的简单的问题,导致一些很郁闷的结果。也是因为对这些基础的概念的理解深入程度影响你的软件思想。

   3.结构体,联合体,位域的运用,最好能有一个对比,了解运用的场合及优缺点。有了这些基础,就很容易明白TI例程的一些基本结构。

   4.变量的作用域以及生存期,这样才能从最初的单个的C文件过渡到C工程,才能更好地了解模块化编程地基本要素,实现基础。

   5.Q格式等的应用,因为大多数DSP都是定点的,而做运算处理的很多时候都不可避免要用到浮点数。使用Q格式才能更好的发挥dsp高速的特性。所以有人说“不会使用Q格式的DSP程序员不是优秀的程序员”!

   6.另外我觉得如果用DSP参与控制的话,最好能熟悉一些面向对象的程序语言,比如C++。要能了解用C实现面向对象的方法。为什么呢?一是面向对象我个人觉得在写程序方面模块化等要好一些,特别是控制领域,我看过一些高手的PLC程序也是采用的面向对象的思想构件模块的,从那以后我的单片机程序才开始这样靠拢,也确实发现了一些好处。比如对那些c流程中多次要操作的资源,有时候会出现一些共用操作的位置错误。二是TI的例程中有好多地方用到了这些思想,要明白这些东西,你才能更容易理解那些例程。这一部分我推荐一本书《DSP C2000程序员高手进阶》。

  7.数据结构知识,DSP要做运算,肯定要涉及到一些好的数据结构。才能更好的优化算法。

  8.matlab的运用。matlab提供了强大的数学运算能力,还能对DSP结合CCS进行仿真。可以把程序员从繁琐的底层编程中解放出来。个人感觉matlab给编程思想和底层程序提供了一个连接,一个平台。 

 

 

感觉不错的几本DSP书
TMS320C6000系列DSP编程工具与指南 {b9|r&</  
【作  者】 (美)Texas Instruments Incorporated  f ; OmR->  
【译  者】  田黎育 何佩琨 朱梦宇  $ ;(k35p  
【丛 书 名】  TI DSP系列中文手册 ,g?m?t KA  
【出 版 社】  清华大学出版社 @anxE%l?  
【书 号】 7302136157 {l|Hk![S  
【开 本】 185×260 dk*5u~]s  
【页 码】 676 -j/~sdrI  
【出版日期】 2006年9月 RgeNvY-  
gM$ r'4  
TI DSP/BIOS用户手册与驱动开发 y?C5'tR B  
【作  者】 (美)Texas Instruments Incorporated cns#n`  
【译  者】  王军宁;何迪;马娟;胡启龙 5aK||:Yp9  
【丛 书 名】  TI DSP系列中文手册 8/"-G@n  
【出 版 社】  清华大学出版社 a{bGgW<p3  
【书 号】 9787302141242 rkE34hRQ/  
【开 本】 16开 s za6:'zj  
【页 码】 255 J}IEMJI  
【出版日期】 2007年4月 hRxV$mLE  
slGtFbN2  
TMS320 DSP算法标准 iZtH_8A  
【作  者】 (美)Texas Instruments Incorporated D+ jc iu  
【译  者】  徐盛;胡剑凌 KfcD3Bv  
【丛 书 名】  TI DSP系列中文手册 rM8  
【出 版 社】  清华大学出版社 m0?p0 z  
【书 号】 9787302148722 SLDP+TNa  
【开 本】 16开 F@lA'QO8+  
【页 码】 130 mdCEN% P*8  
【出版日期】 2007年4月 +u8LQ G[  
%ljQz-[q  
TMS320C28x系列DSP指令和编程指南 [TtA/ef  
【作  者】 (美)Texas Instruments Incorporated ' ^, S  
【译  者】  刘和平 张卫宁 刘林等  :tj>~-`z  
【丛 书 名】  TI DSP系列中文手册 f"Em9*^B1  
【出 版 社】  清华大学出版社 {jHgdFr0  
【书 号】 7302104387 j^C OS  
【开 本】 185×260 `6R:^'x]B  
【页 码】 483 6%kb0QVh2T  
【出版日期】 2005年3月 o"X,#p<  
y|L.NSGDO  
TMS320C54x系列DSP的CPU与外 6C=7W /  
【作  者】 (美)Texas Instruments Incorporated  4}gb5E1  
【译  者】  梁晓雯 裴小平 李玉虎  vIJBv  
【丛 书 名】  TI DSP系列中文手册 3 z 1 >  
【出 版 社】  清华大学出版社 W yD?|#  
【书 号】 7302132216 mW3='e:Y  
【开 本】 185×260 |#?7,j5  
【页 码】 440 W925 @K~  
【出版日期】 2006年9月 ;L`Jes(j  
mB5e}{FF  
DSP系统的VLSI设计 49 j#Bg  
【作  者】 李哲英  54u-A  
【出 版 社】  机械工业出版社* #GclW-  
【书 号】 9787111204787 d!}ulTf  
【开 本】 16开 ;R&k|"]]  
【页 码】 296 1._@!U  
【出版日期】 2007年2月 )+<F[0 }7  
wTy'4 |2  
DSP C2000程序员高手进阶 .ww7<_$4?a  
【作  者】 王潞钢 陈林康 曾岳南 许贤昶 S"OEP;uA  
【出 版 社】  机械工业出版社 0W`lSM+  
【书 号】 7111152697 g''@5pru  
【开 本】 16开 43OJ|5C(<  
【页 码】 206 7Dr0PJj@G  
【出版日期】 2005年1月
  • 0
    点赞
  • 9
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
DSP算法大全C语言版本 完整版,共407页,审阅过的。包含多种数字信号产生、处理、分析方式,并附参考代码。第六章FIR数宇滤波器的设计………………… 227 §6.1窗函数方法……………… .227 §6.2频域最小误差平方设计…… “·自“238 §6.3切比雪夫逼近方法……………………………………………242 第三篇随机数字信号处理 第一章经典谱佔计……………………………………….264 的周期图方法………………·…264 12功率谱估计的相关方法 271 第二章现代谱估计……… 280 §2.1求解一般托布利兹方程组的莱文森算法……………0 82.2求解对称正定方程组的乔里斯基算法 §2.3求解尤利沃克方程的莱文森德宾算法…………灬….28 §24计算ARMA横型的功率谱密度18 §2.5尤利沃克谱佔计算法…………………*…………………22 §2.6协方差谱估计算法 ……ts29 §27Burg谱估计算法 ◆鲁b+吾·合品品‘山点亠4+·日叶中·中‘甲争早导 §2.8最大似然谱估计算法…… 308 第三章时频分析 甲甲甲手曾鲁卧鲁哲雪 §3.1堆格纳( wigner)分布 中中中“由节“昏音山曲画 32离散小波变换 委甲■即 318 第四章随机信号的数字滤波……… §41维纳( Wiener)数字滤波……… …………量·自330 §4-2卡尔曼( Kalman)数字滤波… 会血中自4B44品西4垂4+中如甲吾卧d古 §4-3最小均方(LMS)自适应数字滤波…… ●·中自·自·中中中平看 °·341 §4.4归一化LMS自适应数字滤波……34 §45递推最小二乘(RLs)自适应数字滤波 348 第四篇数字图像处理 第一章图像基本运算 ……·352 §1,1图像读取、存储与显示…… 81.2图像旋转……………………………………………………………366 1.3图像灰度级直方图的计算……………………………………*……3 §1.4图像二值化的固定阀值法… ma·也d■血dp §1.5图像二值化的自适应阀值法 导b血 第二章图像增颯… §2.1图像直方图均衡…………………… w"376 §2.2中值滤波……38 §23图像锐化…… 鲁平t自d “*………382 §2.4图像平滑 P····383 第三章图像边缘检测……………………….356 831 Roberts算子边缘检测…46 §32拉普拉斯算子边缘检测……… “………“敌…388 83.3 Sobel算子边缘检测……; 83.4 Robinson算子边缘检测 小392 §35 Kirsch算子边缘检测……………………………………………394 §3.6 Prewitt算子边缘检测 争昏平辛辛平中萨 396 第四章图像细化…… 喜即香看d画命合b分bb画品目如画如bL晶 品品4甲。自·。·中 399 §4.1 Hilditch细化算法… ●·命··“““…“399 §4.2 Pavlidis细化算法……404 §43 Rosenfeld细化算法………………… 0 第五篇人工神经网络 第一章神经网络模型 t……416 §1.1多层感知器神经网络 ·416 §1.2离散 Hopfield神经网络…………“………………………425 §13连续 Hopfield神经网络………………49434 §I4 Tank-Hopfield线性规划神经网络 ……437 参考文献 44命↓◆命啡4每◆普““女4古“4b中d●·4·面···4·= ,··442 第一篇常用数字信号的产生 第一章数字信号的产生 §1.1均匀分布的随机数 -、功能 产生(a,b)区间上均匀分布的随机数。 方法简介 均匀分布的概率密度函数为 ≤x≤b fCx) 其它 通常用U(b表示均匀分布的均值为“士方差为2 产生均匀分布随机数的方法如下 首先,由给定的初值x0,用混合同余法 T-1+c)(mod M) y; /M 产生(0,1)区间上的随机数y。其中a=2045c=1M=2然后,通过变换z=a+ (b…a)y产生(a,b)区间上的随杋数x 三、使用说明 1.子函数语句 double uniform (a, b, seed) 2形参说明 双精度实型变量。给定区间的下限 b—双精度实型变量。给定区间的上限。 seed——长整型指针变量。*seed为随机数的种子。 四、子函效程序(文件名: uniform,c) double uniform(a,bseed) long int seed double t; 2045爷( seed-#seed-(操Seed/1948576)*1048576; t=(“seed)/1048576.0; t=a+(b-a)头t; return(t)i 五、例题 产生50个0到1之间均匀分布的随机数。 主函数程序(文件名: uniform,m) 杜 include〃 stdio.h includ iform ain( doable a, b,x; int 1,J; g Int s double uniform (double, double, long int *) a-0.0;b=10;s=13579; for(i-0;i<10;i-+) for(j=0:j<5++) =unite printf (". 7fM,x>; intf("\n")y 运行结果 0.48263550.98959450.72067070.77158260.8864250 0.73916340.58915140.81457810.81212620.7979975 0.9483t 0.39095970.51266860.40730760.9440937 0.67162610.47535710.10517980.09266470.4993505 0.13187120.47657490.59566690,13878150,8082657 C.90332890.30759240.02644250.07497020.3141527 0.44230840.5207319 89678960.93463230.3230572 0.65192320.18290330.03722860.13245580.8721647 0,5768661 0.6912775 0.66249660.80548000.2066078 0.51299000.0645466099776740.43443300.4154520 §1.2正态分布的随机数 功能 产生正态分布N(,a2)的随机数。 二,方法简介 正态分布的概率密度函数为 √2πa 通常用、(p,d2)表示。式中p是均值,a2是方差。正态分布也称为高斯分布 产生正态分布随机数的方法如下 设r1 r为(0,1)上n个相互独立的均匀分布的随机数,由于E(n)=1, 根据中心极限定理可知,当n充分大时 r 的分布近似于正态分布N(,1)。通常取n=12,此时有 r一6 最后,再通过变换y=H+x,便可得到均值为八方差为m2的正态分布随机数y 三、使用说明 1.子函数语句 double gauss(mean, sigma, seed) 2.形参说明 mean—双精度实型变量。正态分布的均值k sigma—-双精度实型变量。正态分布的均方差o see 长整型指针变量。餐seed为随机数的种子。 四、子函数程序(文件名: gauss,c) include Uniform,ci double gauss(mean, sigma, s) double mean, sigma; long int *si i int i; double x, y; double uniform() for(x=0,=0;<12i++) x+= uniform(0,0,1.0,) X=x一6.0: y〓mean十x并sgm丑 return (y)f 五、例题 产生50个均值为0方差为1的正态分布的随机数。 主函数程序(文件名 tgauss.m) #include"stdio. hm #include"gauss, c main() t int i,j s long int s# double x, mean, sigma; double gauss(double, double, long int *) mean=0.0; sigma=1.05=13579 for(i=0<10i++) for(j=0<5i十+) x=gauss(mean, sigma, &s) printf(".7",x); printi("\n") 运行结果 2.8997211-0,90885730,2041950-0.2572155-0.8516827 0.79969980,9866190.04313851.919498702543507 0.36892511.2145863 ,05370901.70509531.6925945 0.49287221.9956684.-0.59806631.29232980.1707630 0.5213604-0.40513420.8358479-0.54450801.6452045 0.5338917-0.8120403-0.3886852-0.25463680.4690113 0.4013348-0.1117687-0.9708830.650224713179646 0.53624150.74646191.3275318-0.40414241.8053455 0.8525982-0,24906731.68234440.945543304819355 1.1704273-0.172575002068348-1.9993710.8360157 §1.3指数分布的随机数 功能 产生指数分布的随机数。 方法简介 1.产生随机变量的迆变换法 定理设F(x)是任一连续的分布函数如果~U(0,1)且=F-1(a),那么η~F(x)。 证明由于t~U(0,1),则有 P(≤x)=P(F-(x)≤x)=P(x≤F(x))=F(x) 所以,7~F(x),定理证毕 此定理给出了从均匀分布随机数到给定分布F(x)的随机数的变换。根据该变换可产 生分布函数为F(x)的机效x,其算法可用下列两个步骤实现 (1)产生均匀分布的随机数g,即4~U(0,1)g(2)计算x=F-(t) 2.产生數分布随机的方法 指数分布的概率密度函数为 x≥0 f(x) 0 其它 其分布函效为 FCr)= ,其它 指数分布的均值为,方为P2 根播上述的逆交换法,产生指数分布随机数的方法为 (1)产生均匀分布的随机数M,即w~U(0,1);(2)计算x=-ln(v) 、使用说明 1.子函微语句 double exponent(beta, s) 2.彩参说明 beta——双精皮实型变量。指数分布的均值。 s—长整型指针变量。¥s为随机数的种子。 四、子函数程序(文件名: exponent.[) include "math. h" Include uniform. c uble exponent(bcta, double beta; long int *s i double u,x double uniform) u=uniform(0. 0,1-0,s)+ ta i logt return(x) 五、例题 产生50个均值为2、方差为4的指数分布的随机数。 主函数程序(文件名; exponent n) toinclude stdio. h f include exponentc int i,j; long int s; ouble x 4a; double exponent()+ bea=2.仍;s=t3579; for(i=0;i<10;i++) {for(j=0;j5;j+) i x=exponent(beta, &s printf(".7f", x> 运行结果 45698710,02092010,6551459051862310,2411175 0.6044725 1.0581442 0.4101700 0.4161992 0.451299 0.20000171,87830141.33625131.79637310.1150597 0,7961070 .4873781450416854.75753501.3888938
FFT(快速傅里叶变换)是一种用于对信号进行频域分析和处理的算法,它在数字信号处理(DSP)中被广泛应用。DSP C语言编写FFT程序可以通过以下步骤实现: 1. 引入必要的头文件和定义变量:首先,需要引入相关的头文件,如stdio.h和math.h,并定义必要的变量,如采样点数N和复数结构体,用于存储实部和虚部。 2. 实现FFT算法:可以使用著名的Cooley-Tukey算法来实现FFT。该算法通过将DFT(离散傅里叶变换)分解为多个子问题的DFT来加速计算。具体实现步骤如下: a. 计算要进行FFT的信号的长度L,将其转换为二进制反转序列,以便进行迭代计算。 b. 对转换后的序列进行迭代计算,每一次迭代会将序列的长度减半,直到长度为1时计算结束。 c. 在每一次迭代中,进行蝶形运算(Butterfly Calculation)。蝶形运算将连续的两个输入值作为复数输入,并计算它们的DFT,然后将结果输出到不同的位置,以构建最终的输出序列。 3. 实现反变换(Inverse FFT):如果需要对FFT结果进行反变换,可以在FFT结束后,使用相同的步骤实现IFFT。反变换的结果将恢复原始信号。 4. 运行测试:编写过程中应该加入适当的数据输入输出和测试代码,以确保程序的正确性。 5. 优化和优化:对程序进行优化是重要的,可以通过选择最佳的FFT算法,提高代码的效率和性能。 编写DSP C语言的FFT程序可以通过上述步骤来实现,它可以用于处理各种类型的信号,如音频、图像和其他实时信号。通过使用FFT可以分析信号频谱、滤波、噪声消除等,是DSP领域中不可或缺的重要工具。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值