关于光盘的使用说明.docVIP

关于光盘的使用说明 本光盘共包含三个子目录，分别是DSP_FORTRAN, DSP_C和DSP_MATLAB。DSP_FORTRAN和DSP_C各含有约40个信号处理的子程序，概括了书中所涉及到的绝大部分算法。程序分别由FORTRAN语言和C语言编写（MA模型、ARMA模型及最小方差谱估计三个算法只给出了用C语言编写的程序, 没有给出相应的FORTRAN子程序），并在PC机上调试通过。编译环境是FORTRAN77 V5. 10和TURBO C2. 0。DSP_MATLAB含有近100个用MATLAB编写的信号处理程序，它们是本书各个章节的大部分例题，使用的是MATLAB6.1。 FORTRAN子程序名称的长度全都是6位，扩展名为 .for，C语言子程序的名称全部是7位，由相应的FORTRAN子程序在其名称前加字母m而形成，并将扩展名改为 .c。为了方便读者的使用，光盘中还给出了调用FORTRAN子程序的简单主程序。读者只需将此主程序和主程序指定的子程序作编译、连接和运行，即可得出相应的结果。FORTRAN主程序的名称为7位或8位，它是在原FORTRAN子程序前加字母h所构成的，扩展名仍是 .for。h后面的一个数（如果有的话）表示该程序是相应子程序的第几个主程序。例如，子程序desiir.for是用来设计IIR滤波器的FORTRAN子程序，对应的C程序是mdesiir.c，调用desiir.for的第一个主程序是h1desiir.for（设计低通IIR DF），依此类推。 用MATLAB编写的程序的名称由“exa”开头，接下来是所在的章、节及例题的序号，如exa010101,指的是第1章第1节（即1.1节）的第1个例题，即例1.1.1。如果该程序是为了说明某一个m文件的应用，则在上述名称的后面跟一个下划线，再在后面加上所说明的MATLAB文件的名称，如exa011001_rand，即是例1.10.1,该例用来说明rand.m文件的应用。应该说明的是，这些MATLAB程序不是像所附的FORTRAN和C程序那样作为一个个子程序应用，而是用来说明书上的例题及各个m文件的应用。 用FORTRAN和C语言编写的每一个子程序的功能及调用时各个参数的含义已在程序的开头作了较为详细的说明，此处不再赘述。所附程序中，绝大多数都是作者和其研究生编写的，也有少量是参考国外已公开发表的杂志和教科书，如经典的FFT和REMEZ算法等。 下面给出的是用FORTRAN语言和C语言编写的程序的名称、功能以及有关问题的说明，程序按字母顺序排列。 1．aftodf.for，maftodf.c 将模拟滤波器转变为数字滤波器。 2．ampres.for，mampres.c 由得到幅频响应。系统函数由下式给出： （A.1） 若，则为一全极点IIR系统，若，则为FIR系统。在以下程序中，凡涉及到系统，系统函数都按（A.1）式给出，不再一一说明。 3．arburg.for（harburg.for），marburg.c 用Burg算法求AR模型的参数。 4．ar1psd.for，mar1psd.c 由AR模型参数得到功率谱。 5．aryuwa.for（haryuwa.for），maryuwa.c 用Levinson算法求解Yule-Walker方程以得到阶AR模型的参数。 6．biline.for，mbiline.c 实现双线性Z变换。 7．butwcf.for，mbutwcf.c 设计巴特沃斯模拟低通滤波器，求出转移函数。 8．chebsh.for，mchebsh.c 设计切比雪夫I型模拟低通滤波器，求出转移函数。 9．cmpdft.for（hcmpdft.for），mcmpdft.c 直接由定义求点复序列的DFT。 10．cmpfft.for（hcmpfft.for），mcmpfft.c 利用经典的Cooley-Tukey基2算法求复序列的DFT。 11．convo1.for（hconvo1.for），mconvo1.c 直接按定义求序列的线性卷积。 12．convo2.for（hconvo2.for），mcnvo2.c 利用DFT的卷积性质求两个复序列的线性卷积。 13．corpsd.for（hcorpsd.for），mcorpsd.c 用Blackman-Tukey法（自相关法）对信号作功率谱估计。 14．corre1.for（hcorre1.for），mcorre1.c 直接按定义计算的互相关函数，若，则求出的是自相关，为复序列。 15．corre2.for（hcorre2.for），mcorre2.c 用FFT实现相关函数快速估计。 16．cztfft.for（hcztfft.for），mcztfft.c 实现Chirp