非常有用的关于C语言的数字信号处理程序集合,作为从Matlab上移植向C的良好参考35035数字信号处理C语言程序集殷福亮宋爱军主编辽宁科学技术出版社沈阳图书在版编目CI)数据数字信号处理C语,程序集/股楣亮,宋爱车主,一阳:江宁科7技术出社.197.2IsBN7.5381252131.数…,股…,C语言数学信号信号处应用程释.TN9172中国本图书C数据核字(97)第03286号辽宇科学技术出柱出暇沈阳市和平区北一马路1%号郎嫂编码11001)地方回营新民印崭总厂印刷新华书店北京发行所发行什不:787×1/16即条:28字数:7101,00019977月第11997年7月第1次印刷责任编辑:与东取式设计t干浪封画设计:雏君又责任酸对1东戈印数.150O0完价,出.∞0元目录第一篇常用数字信号的产生第一帝数字信号的产生…§1.1均匀分布的随机数………§1.2正态分布的随机数3§1.3指数分布的随机数§1.4拉普拉斯( Laplace)分布的随机数§1.5瑞利( Rayleigh)分布的随机数§1.6对数正态分布的随机数§1,7柯西( Cauchy)分布的随机数791381.8韦伯(Wabu)分布的随机数”1551,9爱尔朗( Erlang)分布的随机数s11U贝努里( bernoulli分布的随矾数19§1.11贝努旦高斯分布的随机数§1.12二项式分布的随机数…381.13泊松 (Poisson)分布的随机数51,14ARMA(p,q)模型数据的产生………21,1言有高斯白噪声的止弦组合信号的严生…305116解析信号的产生………………35第二篇数字信号处理第一章快速傅立叶变换……………………………§1.1离散傅立叶变换…9§1.2快速傅立叶变换§1,3基4快速傅立叶变换51.4分裂基快迷傅立叶变换*57§1,5实序列快速傅立叶变换(-)…………61§1,6实序列快速傅立叶变换(二)…56§1.7用一个N点复序列的FFT同时计算两个N点实序列离散傅立叶变换+70§18共轭对称序列的快速傅立叶反变换7381.9素因于快速傅立叶变换§1.10 Chirp Z变换算法第二章快速离散正交变換102§2.1快速哈特莱( Hartley)变换……102§2.2基4快速哈待莱( Hartley)变决a,aa44B、和和1062.3分裂基快速哈特菜( Hartley)变换…不110§24快速离散余弦变换……………………11525快谏离散余弦反变换…·,·,11R§2.6N一8点快速离散余弦变换”12132.7N一8点快速离散余弦反变换…125§2.8快速离散正弦变换…………119§2.9快速沃尔什(wash)变换“““““““……133§2.10快速希尔怡特变换(一)……137§2.11快速希尔伯特变换(二)141第三章快速卷积与相关……………………………143,1快速卷积144§3.2长序列的快速卷积……………147§3.3特别长序列的快速卷积…”15233,4饮速相关15悉第四章数字滤波器的时城和频城响应……………………16§4,1数字滤波器的频率响应163§4.2级联型数字滤波器的频率响应……166§4.3数字滤波器的时域响应………………………………171§44直接型IR数字滤波(一)………174§4.b直接型IK觌子据破〔二)177§46级联型IR数字滤波……81§4.7井联型IR数字滤波……………85第五章IR数字滤波器的设计…189§5.1巴特沃兹和切比雪夫数字滤波器的设计……………………………189§5.2任意幅度IR数字滤波器的优化设计…………208第六章FIR数字滤波器的设计………………….:27§61窗函数方法………………2§6z狽域敢小误差半方设计238§6.3切比雪夫逼近方法242第三篇随机数字信号处理第一章经典谱估计264§1.1功率谱估计的周期图方法…654§1.2功率谱估计的相关方法…某来某““271第二章现代谱估计…§2.1求解一般托布利兹方程组的莱文森算法………表280§22求解对称正定方程组的乔里斯基算法……283§2.3求解尤利沃克方程的莱文森德宾算法287§2.4计算ARMA模型的功率谱密度……289§2.5尤利沃克谱估计算法292§2.6协方差谱估计算法…297§z.了Burg谱估计算法303§2.8最大似然谐估计算法………第三章时频分析314§3,1维格纳( wigner)分布………“““31493,2离散小变换术318第四章随机信号的数宁滤波…41维纳( Wiener)数字滤波s4.2卡尔曼( Kalman)数字滤波…………§±0最小均方(LM3)白道应数宁滤被341s4.4归一化LMS自适应数字滤波………§45推最小二乘(RLS)自适应数字滤波………148第四篇数字图像处理第一章图像基本运算352§1,1图像读取、存储与显示…352§1,2图像旋转………36681.3图像灰度级直方图的计算………38§1,4图像二值化的固定阀值法…31§1,5图像二值化的自适应阀值法…372第二章图像增强§2.1图像直方图均衡…………………………………………………37§22中值滤波…““…““““水”…十””“”378§2.3图像锐化………………………………1382§2,4图像平滑383第三章图像边缘检测…………………38683.1 Roberts算子边缘检测…………6§3.2拉昔拉斯算子边缘检测388§3.3 Sobel算子边缘检测*390§3,4 Robinson算子边缘检测……12§3.5 Kirsch算子边缘检测“中394§3.6 Prewitt算子边缘检测396第四章图像细化…§41 Hilditch细化算法…39964.2 Pavlidis细化算法6404§4.3 Rosenfeld细化算法410第五篇人工神经网络第一章神经网络模型……………………………………………………4161.1多层感知器神经网络416§1,2离散 Hopfield神经网络……1125§1.3運绥 Hopfield柙经网络…434§1,4 Tank-Hopfield线性规划神经网络437参考文献…………………………………4第一篇常用数字信号的产生第一章数字信号的产生§1.1均匀分布的随机数功能产生(,b)区间上均匀分布的随机数,二、方法简介均匀分布的概率密度函数为≤≤b其它通常用U(a,b)表示。均匀分布的均值为十,方差为bt12产生均匀分布随机氨的方法如下:曾先,由给定的初值x,用混合同余法=(ax-1+ c)(od M)vi=x/M产生(0,1)区间上的随机数y其中:a=2045,c=1,M=2然后,通过变换24b-a)y产生(ab)区间上的随机数x三、使用说明1.子函数语句double uniform (a, b, seed)2.形参说明双精度实型变量,给定区间的下限b—双精度实强变量,给定区问的上限,seed长整型指针变量,·eed为随机数的种子四、子函数程序(文件名 uniform,c)double uniform (a, b seed)double a, biong tnt seed;k double t:seed=2045*(* seed)+l:eed="seed-(.sed/1048575)1048576;t=(.seed)/1048576.0;t=a+(b-a)*trreturn(t)五、例题严生50个0到1之同均匀分布的随机数。主函数程序(文件名; uniform.m):include"stdio. h'f inelude uniformmaino1 double a,b x: int i,rlong int s4double uniform(double, double, long int *)B=0.0;b=1.0;s=13579;for(i=0;i<10;i-+)forG-0;<5}+)uniform(a, b, &s)printf(%13.7f".x);printf")n”);运行结果0.48263550.98959450.72067070.77158260.88642500.73916340.58915140.81457810.81212620,79799750.90483670.39095970.51266860.40730760.94409370.67162610,47535710.10517980,09266470,49935050.17187120.47657490.595G6690.13878150.80820570.90332890.30759240.02644250.07497020.31415270.44230840.52073190.89678960,93463230.32305720.65192320.18290330.0372286.13245580.87216470.57686610.69127750.66249660.80548000.2066078加Nw