dsp原理與應用課程設計信號采集和濾波dsp程序設計

1、<p>　　課 程 設 計 報 告</p><p>　　名 稱 DSP原理與應用課程設計 </p><p>　　題 目 信號采集和濾波DSP程序設計 </p><p>　　指導教師 潘建軍/沈冰夏 </p><p>　　設計起止日期 2015.11.16至2015.12.1

2、8</p><p>　　學 院 信息與通信工程學院 </p><p>　　專 業(yè) 電子信息工程 </p><p>　　學生姓名 李巖 </p><p>　　班級/學號 電信1301/20130109 </p>&

3、lt;p>　　成 績 </p><p><b>　　目 錄</b></p><p><b>　　第一章 概述1</b></p><p>　　1.1 課程設計目的1</p><p>　　1.2 課程設計所用設備1</p&

4、gt;<p>　　第二章 課程設計原理1</p><p>　　2.1 A/D原理1</p><p>　　2.1.1 模數(shù)轉(zhuǎn)換器結(jié)構(gòu)和時序1</p><p>　　2.1.2 模數(shù)轉(zhuǎn)換器的寄存器3</p><p>　　2.2 FIR 濾波器設計原理6</p><p>　　2.3 FIR濾波的DSP程

5、序設計6</p><p>　　第三章FIR數(shù)字低通濾波器的MATLAB設計8</p><p>　　第四章 ADC信號采集和濾波DSP程序設計9</p><p>　　4.1 采集信號頻率9</p><p>　　4.2 課設結(jié)果截圖9</p><p>　　4.2.1. AD采集兩路信號的時頻圖9</p&g

6、t;<p>　　4.2.2. 兩路信號混合后波形圖10</p><p>　　4.2.3. 濾波后波形圖11</p><p>　　第五章 結(jié)束語13</p><p><b>　　第一章 概述</b></p><p>　　1.1 課程設計目的 </p><p>　　掌握TMS320

7、VC5509A片內(nèi)模數(shù)轉(zhuǎn)換器（ADC）的控制方法。</p><p>　　掌握TI DSP 系統(tǒng)的軟件開發(fā)。</p><p>　　設計TMS320VC5509A DSP系統(tǒng)的AD采集程序。</p><p>　　設計FIR濾波器系數(shù)</p><p>　　對采集的信號進行FIR濾波</p><p>　　1.2 課程設計所用設

8、備</p><p>　　計算機：安裝軟件開發(fā)工具CCStudio v3.3，完成系統(tǒng)的軟件開發(fā)，進行硬件仿真調(diào)試，MATLAB。</p><p>　　ICETEK–VC5509-AG-EDU 實驗箱：包括</p><p><b>　　信號源</b></p><p>　　實現(xiàn)硬件仿真調(diào)試時與硬件系統(tǒng)的通信、控制和讀取硬件系

9、統(tǒng)狀態(tài)和數(shù)據(jù)的開發(fā)系統(tǒng)(ICETEK 5100-USB 仿真器模塊)</p><p>　　提供軟件運行和調(diào)試的平臺和用戶系統(tǒng)開發(fā)的參照的評估模塊(ICETEK VC5509-A)</p><p>　　示波器：觀察DSP輸入輸出信號的波形和頻率。</p><p>　　第二章 課程設計原理</p><p><b>　　2.1 A/D原理

10、</b></p><p>　　2.1.1 模數(shù)轉(zhuǎn)換器結(jié)構(gòu)和時序</p><p>　　在數(shù)字信號處理器的具體應用中往往需要采集一些模擬信號量，如電池電壓、面板旋鈕輸入值等，模數(shù)轉(zhuǎn)換器就是用來將這些模擬量轉(zhuǎn)化為數(shù)字量來供DSP使用。</p><p>　　TMS320VC5507/5509片內(nèi)模擬-數(shù)字轉(zhuǎn)換器（ADC）方框圖如圖1所示。TMS320VC5509

11、A有三種封裝，本次課程設計所用瑞泰創(chuàng)新公司生產(chǎn)的TMS320VC5509A DSP系統(tǒng)的評估板上焊接了144管腳薄型四方扁平封裝的TMS320VC5509A（PGE 后綴）。本課程設計所使用的TMS320VC5509A DSP只能同時支持兩路的信號采集。TMS320VC5509APGE的片內(nèi)模數(shù)轉(zhuǎn)換器支持兩路輸入，管腳分別為AIN0和AIN0。片內(nèi)ADC具有如下特性：</p><p>　　帶內(nèi)置采樣和保持的10

12、位模數(shù)轉(zhuǎn)換模塊ADC</p><p>　　最小轉(zhuǎn)換時間為500ns</p><p>　　最大采樣率為21.5kHz （最小采樣周期：103/21.5=46.51163μs） </p><p><b>　　圖1 ADC方框圖</b></p><p>　　模數(shù)轉(zhuǎn)換器采用連續(xù)逼近式結(jié)構(gòu)，在模數(shù)轉(zhuǎn)換器內(nèi)部用了三個可編程分頻器來靈

13、活地產(chǎn)生用戶需要的采樣率。</p><p>　　整個模數(shù)轉(zhuǎn)換過程分為兩個周期——采樣/保持周期及轉(zhuǎn)換周期，如圖2所示。</p><p>　?。?）采樣/保持周期是采樣/保持電路采集模擬信號的時間，這個周期大于或等于40μs；</p><p>　　（2）轉(zhuǎn)換周期是RC比較網(wǎng)絡在一次采樣中完成逼近處理并輸出A/D轉(zhuǎn)換結(jié)果的時間，這需要13個轉(zhuǎn)換時鐘周期。ADC轉(zhuǎn)換時鐘的

14、最大頻率為2MHz。</p><p><b>　　圖2 轉(zhuǎn)換時序圖</b></p><p>　　模數(shù)轉(zhuǎn)換模塊接到啟動轉(zhuǎn)換信號后，按照設置開始進行相應通道的數(shù)據(jù)采樣轉(zhuǎn)換。</p><p>　　經(jīng)過一個采樣時間的延遲后，將采樣結(jié)果放入AD 數(shù)據(jù)寄存器中保存。</p><p>　　轉(zhuǎn)換結(jié)束，設置標志。</p>&

15、lt;p>　　等待下一個啟動信號。</p><p>　　由于TMS320VC5509A DSP 片內(nèi)的A/D 轉(zhuǎn)換精度是10 位的，ADC數(shù)據(jù)寄存器(16 位)的最高位(第15 位)表示轉(zhuǎn)換值是否有效(0 有效)，第14-12 位表示轉(zhuǎn)換的通道號，第11-10位為保留字段，低10 位為轉(zhuǎn)換數(shù)值，所以在保留時應注意取出結(jié)果的低10 位，再根據(jù)高4 位進行相應保存。</p><p>　

16、　2.1.2 模數(shù)轉(zhuǎn)換器的寄存器</p><p>　　模數(shù)轉(zhuǎn)換器的寄存器包括ADC控制寄存器（ADCCTL）、ADC數(shù)據(jù)寄存器（ADCDATA）、ADC時鐘分頻寄存器（ADCCLKDIV）和ADC時鐘控制寄存器（ADCCLKCTL）。</p><p>　　表1 模數(shù)轉(zhuǎn)換器 (ADC) 的寄存器</p><p>　　? Hardware reset; x denote

17、s a “don’t care.”</p><p>　　ADC控制寄存器(ADCCTL)：發(fā)送轉(zhuǎn)換通道號和啟動命令，然后循環(huán)等待轉(zhuǎn)換結(jié)果，最后將結(jié)果保存。</p><p>　　表2 ADC控制寄存器ADCCTL</p><p>　　ADC數(shù)據(jù)寄存器(ADCDATA)：只讀寄存器，顯示AD轉(zhuǎn)換是否正在進行中以及通道號，存放采樣數(shù)據(jù)。</p><p&

18、gt;　　表3 ADC數(shù)據(jù)寄存器ADCDATA</p><p>　　ADC時鐘分頻寄存器(ADCCLKDIV)：顯示轉(zhuǎn)換時鐘的分頻值和采樣保持時間。</p><p>　　表4 ADC時鐘分頻寄存器ADCCLKDIV</p><p>　　ADC時鐘控制寄存器(ADCCLKCTL)：存放CPU時鐘分頻數(shù)，并確定是否將ADC置于低功耗狀態(tài)。</p><

19、p>　　表5 ADC時鐘控制寄存器ADCCLKCTL </p><p>　　本設計中DSP系統(tǒng)CPU時鐘為144MHz。</p><p>　　（1）首先對CPU時鐘分頻，產(chǎn)生模數(shù)轉(zhuǎn)換器時鐘(ADC Clock)，該時鐘應盡量運行在較低頻率下，以降低功率消耗，在本例中模數(shù)轉(zhuǎn)換器時鐘是通過對系統(tǒng)主時鐘36分頻產(chǎn)生的，則此時模數(shù)轉(zhuǎn)換器時鐘=144MHz/36=4MHz，根據(jù)公式</p

20、><p>　　ADC Clock = (CPU Clock) / (CPUCLKDIV + 1)</p><p>　　得出CPUCLKDIV =35；</p><p>　?。?）對模數(shù)轉(zhuǎn)換器時鐘分頻產(chǎn)生模數(shù)轉(zhuǎn)換器轉(zhuǎn)換時鐘（ADC Conversion Clock），該時鐘最大值為2MHz，為了獲得2MHz的模數(shù)轉(zhuǎn)換器轉(zhuǎn)換時鐘，則需要對模數(shù)轉(zhuǎn)換器時鐘2分頻，由</p

21、><p>　　ADC Conversion Clock = (ADC Clock) / (2 ×(CONVRATEDIV + 1))</p><p>　　得出CONVRATEDIV=0 </p><p>　　ADC 轉(zhuǎn)換時間 = 13 ×(1 / ADC Conversion Clock)</p><p>　　ADC 轉(zhuǎn)換時間

22、 = 13 × (1 / (2 MHz)) = 6.5μs；</p><p>　　（3）對采樣和保持周期進行設置，這個值必須大于或等于40μs </p><p>　　ADC Sample and Hold Period =</p><p>　　(1 / (ADC Clock)) × (2 ×(CONRATEDIV + 1 + SAMPT

23、IMEDIV))</p><p>　　= (1 / (4 MHz)) × (2 ×(0 + 1 + SAMPTIMEDIV))</p><p>　　= 250 ns× (2 ×(0 + 1 + 79)) = 40μs</p><p>　　SAMPTIMEDIV =79；</p><p>　?。?）ADC

24、整個轉(zhuǎn)換時間（ADC Total Conversion Time）為40μs的采樣和保持時間加上6.5μs的轉(zhuǎn)換時間等于46.5μs，采樣率=1/46.5μs =21.5 kHz。</p><p>　　2.2 FIR 濾波器設計原理</p><p>　　設計FIR低通。低通濾波器的幅度響應和性能如圖3所示。</p><p>　　圖3 低通濾波器的幅度響應和性能<

25、;/p><p>　　ωp為 通帶截止頻率，ωs為阻帶截止頻率，ωc為截止頻率。δp峰值通帶波紋。δs為最小阻帶衰減。</p><p>　　2.3 FIR濾波的DSP程序設計</p><p>　　FIR濾波器輸出y(n)是輸入數(shù)據(jù)和系數(shù)的卷積和。</p><p>　　任何時候進行計算，都需要把數(shù)據(jù)和相對應的系數(shù)相乘，再累加到輸出結(jié)果中。 </

26、p><p>　　FIR濾波器中，輸入信號在每個采樣時間T過后都發(fā)生改變。也就是說，x(n)將在T時候過后，成為x(n-1)，然后是x(n-2)，……，直到該數(shù)據(jù)被拋棄。</p><p>　　本程序設計將采用MATLAB設計的FIR濾波器系數(shù)導入到DSP程序中。對信號的FIR濾波函數(shù)可以采用設計FIR濾波的浮點實現(xiàn)C函數(shù)和兩種方法調(diào)用德州儀器的TMS320C55x DSPLIB的FIR濾波函數(shù)，

27、實現(xiàn)對A/D采集的輸入信號進行濾波處理。詳細FIR濾波程序設計實現(xiàn)方法見DSP原理與應用講義。</p><p>　　(1) 設計FIR 濾波的浮點實現(xiàn)C函數(shù)</p><p>　　(2) 調(diào)用德州儀器的TMS320C55x DSPLIB的FIR濾波函數(shù)</p><p>　　德州儀器的TMS320C55x DSPLIB是專門針對TMS320C55x系列處理器開發(fā)的數(shù)

28、字信號處理函數(shù)庫，它包括50個經(jīng)過匯編級優(yōu)化的數(shù)字信號處理函數(shù)，這些函數(shù)可由C調(diào)用。</p><p>　　數(shù)字信號處理庫由頭文件dsplib.h，目標庫55xdsp.lib（小模式）、 55xdspx.lib（大模式） ，源文件55xdsp.src等組成。</p><p>　　數(shù)字信號處理庫可以從www.ti.com網(wǎng)站下載，下載的文件是一個壓縮文件sprc100.zip （最新版本日期：

29、 Apr-10-2007   ），解壓之后生成一個dsplib_2.40.00目錄，用戶可以把這個目錄及目錄下的全部內(nèi)容復制到…\c5500\目錄下。例子程序在examples子目錄下。 </p><p>　　庫函數(shù)中主要的數(shù)據(jù)類型是Q.15（Q15, 16位字），其他數(shù)據(jù)類型還包括Q.31（Q31, 32位雙字）和Q.3.12（ Q12，3位整數(shù)，12位小數(shù)）。</p><

30、;p>　　函數(shù)中參數(shù)的傳遞大多采用數(shù)組的形式，其中數(shù)組成員都在內(nèi)存中連續(xù)保存。</p><p>　　當存放的是復數(shù)時，則以實部－虛部的格式存放。</p><p>　　為了節(jié)省存儲器的空間，源操作數(shù)和目的操作數(shù)可以使用相同的地址。</p><p>　　如果在函數(shù)運行時發(fā)生了溢出，函數(shù)可以返回溢出標志，但是由于C55x的寄存器為40位，其中有8位的保護位，因此在運

31、算中可能發(fā)生了32位溢出但結(jié)果仍然正確，在這種情況下，溢出標志代表報警而不是一個錯誤。</p><p>　　調(diào)用DSPLIB庫需要包含dsplib.h，并在工程中加入55xdsp.lib或55xdspx.lib庫或在工程中加入包含調(diào)用函數(shù)的匯編源代碼的匯編文件。TMS320C55x DSPLIB的FIR 濾波的函數(shù)說明如下。</p><p>　　fir2 FIR2 Filter FI

32、R direct form (Optimized to use DUAL?MAC)</p><p>　　Function </p><p>　　ushort oflag = fir2 (DATA *x, DATA *h, DATA *r, DATA *dbuffer, ushort nx, ushort nh)</p><p><b>　　Argume

33、nts</b></p><p>　　x[nx] Pointer to input vector of nx real elements.</p><p>　　h[nh] Pointer to coefficient vector of size nh in normal order.For example, if nh=6, then h[nh] = {h0, h1

34、, h2, h3, h4, h5} where h0 resides at the lowest memory address in the array.</p><p>　　r[nx] Pointer to output vector of nx real elements. In-place computation (r = x) is allowed.</p><p>　　dbu

35、ffer[nh+2] Pointer to delay buffer of length nh + 2</p><p>　　nx Number of input samples</p><p>　　nh The number of coefficients of the filter. For example, if the filter coefficients ar

36、e {h0, h1, h2, h3, h4, h5}, then nh= 6. Must be a minimum value of 3. For smaller filters, zero pad the coefficients to meet the minimum value.</p><p>　　oflag Overflow error flag (returned value)</p>

37、;<p>　　If oflag = 1, a 32-bit data overflow occurred in an intermediate or final result.</p><p>　　If oflag = 0, a 32-bit overflow has not occurred.</p><p>　　Benchmarks 基準測試</p>

38、<p>　　Cycles? Core: nx * (3 + nh/2)</p><p>　　Overhead: 25</p><p>　　Code size 107</p><p>　　(in bytes)</p><p>　　Assumes all data is in on-chip dual-acce

39、ss RAM (provided linker command file reflects those</p><p>　　conditions).</p><p>　　第三章FIR數(shù)字低通濾波器的MATLAB設計</p><p>　　本設計選用Kaiser窗設計FIR低通濾波器，采樣頻率Fs=21500Hz, 與AD采集部分計算出的采樣頻率一致，通帶截止頻率Fpa

40、ss=1000Hz, 阻帶截止頻率Fstop=2550Hz，阻帶最小衰減Astop=40dB，濾波器如下圖：</p><p>　　圖3.1 低通濾波器圖</p><p>　　第四章 ADC信號采集和濾波DSP程序設計</p><p>　　4.1 采集信號頻率 </p><p>　　fs=21.5KHz</p><p>

41、;　　4.2 課設結(jié)果截圖</p><p>　　4.2.1. AD采集兩路信號的時頻圖</p><p>　　圖4.2.1 AD采集兩路信號時頻圖</p><p>　　4.2.2. 兩路信號混合后波形圖</p><p>　　圖4.2.2 Q15定點混合時域圖</p><p>　　圖4.2.3 Q15定點混合頻域圖&l

42、t;/p><p>　　圖4.2.4 浮點混合時域圖</p><p>　　圖4.2.5 浮點混合頻域圖</p><p>　　4.2.3. 濾波后波形圖</p><p>　　圖4.3.1 Q15定點濾波后時域圖</p><p>　　圖4.3.2 Q15定點濾波后頻域圖</p><p>　　圖

43、4.3.3 浮點濾波后時域圖</p><p>　　圖4.3.4 浮點濾波后頻域圖</p><p><b>　　結(jié)束語</b></p><p><b>　　感悟：</b></p><p>　　通過本次課程設計，我學到了很多知識，同時也認識到自己編程技術的薄弱，我會在以后的學習中加強編程知識的學習。本

44、次試驗中我對CSS軟件開發(fā)環(huán)境有了較深的認識和理解。通過完整的實驗我對本專業(yè)的知識體系有了系統(tǒng)性的理解，使我受益匪淺。濾波器設計中我理解了濾波器參數(shù)的實際意義，讓我對課程中的知識有了更深入的理解。</p><p>　　試驗中遇到不少問題都是在老師和同學的幫助下解決的，謝謝老師的幫助。</p><p><b>　　參考文獻</b></p><p>

45、;　　汪春梅, 孫洪波. TMS320C55x DSP原理及應用(第 3 版). 北京: 電子工業(yè)出版社, 2011.</p><p>　　Texas Instruments Incorporated. TMS320VC5507/5509 DSP Analog-to-Digital Converter (ADC) Reference Guide[OL]. Texas Instruments Incorporate

46、d, SPRU586B, 2004. http://www.ti.com.</p><p>　　Texas Instruments Incorporated. TMS320VC5509A Fixed-Point Digital Signal Processor. Texas Instruments Incorporated, SPRS205K, 2008. http://www.ti.com.</p>

47、<p>　　DSP原理與應用課程講義.</p><p>　　瑞泰創(chuàng)新公司. ICETEK–VC5509-A評估板硬件使用指導. 北京:瑞泰創(chuàng)新公司, 2006.</p><p>　　瑞泰創(chuàng)新公司. ICETEK-VC5509-A-USB-EDU教學實驗系統(tǒng)使用指導. 北京:瑞泰創(chuàng)新公司, 2006.</p><p>　　瑞泰創(chuàng)新公司. ICETEK–V

48、C5509-A-USB-EDU教學實驗系統(tǒng)軟件實驗指導. 北京:瑞泰創(chuàng)新公司, 2006.</p><p>　　譚浩強. C程序設計(第三版)[M]. 北京: 清華大學出版社, 2005.</p><p>　　Texas Instruments Incorporated. Code Composer Studio Help. </p><p>　　Texas Ins

49、truments Incorporated. TMS320C55x DSP Programmer’s Guide[OL]. Texas Instruments Incorporated, SPRU376A, 2001. http://www.ti.com. </p><p>　　Texas Instruments Incorporated. TMS320C55x Optimizing C/C++ Compiler

50、 User’s Guide [OL]. Texas Instruments Incorporated, SPRU281F, 2003. http://www.ti.com.</p><p>　　Texas Instruments Incorporated. TMS320C55x DSP Peripherals Overview Reference Guide. Texas Instruments Incorpor

51、ated. SPRU317H，2006. http://www.ti.com.</p><p>　　Sen M. Kuo, Bob H. Lee. Real-Time Digital Signal Processing Implementations, </p><p>　　Applications and Experiments With the TMS320C55x [M]. 北京:

52、 清華大學出版社, 2003.</p><p>　　焦瑞莉等編著. 數(shù)字信號處理[M]. 北京: 機械工業(yè)出版社. 2011.</p><p><b>　　附錄：程序源代碼</b></p><p><b>　　main.c</b></p><p>　　#include "myapp.h&q

53、uot;</p><p>　　#include "ICETEK-VC5509-EDU.h"</p><p>　　#include "scancode.h"</p><p>　　void InitADC();</p><p>　　void wait( unsigned int cycles );</

54、p><p>　　void EnableAPLL( );</p><p>　　int nADC0[256],nADC1[256];</p><p>　　int mixQ15[256];</p><p>　　float mixfloat[256];</p><p>　　int FIROUTQ15[256];</p>

55、<p>　　int DBBufferQ15[34];</p><p>　　float Bufferfloat[32];</p><p>　　float FIROUTfloat[256];</p><p>　　const int BhQ15[32] = {</p><p>　　98, 141, 138, 57,

56、 -117, -360, -604, -747, -671,</p><p>　　-284, 451, 1490, 2700, 3888, 4838, 5366, 5366, 4838,</p><p>　　3888, 2700, 1490, 451, -284, -671, -747, -604, -

57、360,</p><p>　　-117, 57, 138, 141, 98</p><p><b>　　};</b></p><p>　　const float Bhfloat[32] = {</p><p>　　0.002999524353, 0.004298123531, 0.004219

58、199065, 0.001731093973,-0.003574917093,</p><p>　　-0.01098363195, -0.01843946427, -0.02278480493, -0.02048166841,-0.008664992638,</p><p>　　0.01377640851, 0.04545684904, 0.08239003271, 0.11864

59、56978, 0.1476465911,</p><p>　　0.163765952, 0.163765952, 0.1476465911, 0.1186456978, 0.08239003271,</p><p>　　0.04545684904, 0.01377640851,-0.008664992638, -0.02048166841, -0.0227848049

60、3,</p><p>　　-0.01843946427, -0.01098363195,-0.003574917093, 0.001731093973, 0.004219199065,</p><p>　　0.004298123531, 0.002999524353</p><p><b>　　};</b></p><p&g

61、t;<b>　　main()</b></p><p><b>　　{</b></p><p><b>　　int i;</b></p><p>　　unsigned int uWork;</p><p>　　EnableAPLL();</p><p>　

62、　SDRAM_init();</p><p>　　InitADC();</p><p>　　PLL_Init(144);</p><p>　　for (i=0;i<34;i++) </p><p>　　DBBufferQ15[i]=0;</p><p>　　for (i=0;i<32;i++) </p

63、><p>　　Bufferfloat[i]=0;</p><p>　　while ( 1 )</p><p><b>　　{ </b></p><p>　　for ( i=0;i<256;i++ )</p><p><b>　　{</b></p><p&

64、gt;　　ADCCTL=0x8000;// 啟動AD轉(zhuǎn)換，通道0</p><p><b>　　do</b></p><p><b>　　{</b></p><p>　　uWork=ADCDATA;</p><p>　　} while ( uWork&0x8000 );</p>

65、<p>　　nADC0[i]=uWork&0x0fff;</p><p><b>　　}</b></p><p>　　for ( i=0;i<256;i++ )</p><p><b>　　{</b></p><p>　　ADCCTL=0x9000;// 啟動AD轉(zhuǎn)換，通

66、道1</p><p><b>　　do</b></p><p><b>　　{</b></p><p>　　uWork=ADCDATA;</p><p>　　} while ( uWork&0x8000 );</p><p>　　nADC1[i]=uWork&

67、0x0fff;</p><p>　　mixQ15[i]=(nADC0[i]+nADC1[i]-1024)*32;</p><p>　　mixfloat[i]=(float)(nADC0[i]+nADC1[i]-1024)/1024;</p><p><b>　　}</b></p><p>　　fir2(mixQ15,Bh

68、Q15,FIROUTQ15,DBBufferQ15,256,32);</p><p>　　firc(mixfloat,256,Bhfloat,32,FIROUTfloat,Bufferfloat);</p><p>　　asm( " nop");// break </p><p><b>　　}</b></p>

69、;<p><b>　　}</b></p><p>　　void InitADC()</p><p><b>　　{</b></p><p>　　ADCCLKCTL=0x23; // 4MHz ADCLK</p><p>　　ADCCLKDIV=0x4f00;</p>&l

70、t;p><b>　　}</b></p><p>　　void wait( unsigned int cycles )</p><p><b>　　{</b></p><p><b>　　int i;</b></p><p>　　for ( i = 0 ; i < c

71、ycles ; i++ ){ }</p><p><b>　　}</b></p><p>　　void EnableAPLL( )</p><p><b>　　{</b></p><p>　　/* Enusre DPLL is running */</p><p>　　*(

72、ioport volatile unsigned short* )0x1f00 = 4;</p><p>　　wait( 25 );</p><p>　　*( ioport volatile unsigned short* )0x1f00 = 0;</p><p>　　// MULITPLY</p><p>　　*( ioport volat

73、ile unsigned short* )0x1f00 = 0x3000;</p><p><b>　　// COUNT</b></p><p>　　*( ioport volatile unsigned short* )0x1f00 |= 0x4F8;</p><p>　　wait( 25 );</p><p>　　/

74、/*( ioport volatile unsigned short* )0x1f00 |= 0x800</p><p><b>　　// MODE</b></p><p>　　*( ioport volatile unsigned short* )0x1f00 |= 2;</p><p>　　wait( 30000 );</p>

75、<p>　　// APLL Select</p><p>　　*( ioport volatile unsigned short* )0x1e80 = 1;</p><p><b>　　// DELAY</b></p><p>　　wait( 60000 );</p><p><b>　　}<

76、/b></p><p><b>　　Fir.c</b></p><p>　　void firc(float *in, int M, float *h, int L, float *out, float *x)</p><p>　　{ </p><p><

77、b>　　float yn;</b></p><p>　　int i,j;</p><p>　　for (j=0; j<M; j++)</p><p><b>　　{ </b></p><p>　　x[0] = in[j];</p><p>　　for(yn=0.0,

78、i=0;i<L;i++) </p><p><b>　　{</b></p><p>　　yn+=h[i]*x[i];</p><p><b>　　}</b></p><p>　　out[j]=yn;</p><p>　　for (i=L-1; i>0;

79、 i--)</p><p><b>　　{ </b></p><p>　　x[i] = x[i-1]; //x[L-1]=x[L-2], … , x[1]=x[0]</p><p><b>　　} </b></p><p><b>　　}</b></p>