IIR數(shù)字濾波器的DSP實現(xiàn)

1、 湖南科技大學(xué)信息與電氣工程學(xué)院課程設(shè)計報告題 目：IIR數(shù)字濾波器的DSP實現(xiàn)專 業(yè)： 電子信息工程 班 級： 電子二班 姓 名： 高二奎 學(xué) 號： 1104030205 指導(dǎo)教師： 尹艷群 2015年 1月 8 日信息與電氣工程學(xué)院課程設(shè)計任務(wù)書20142015學(xué)年第一學(xué)期專業(yè)： 電子信息工程 班級： 電子二班 學(xué)號： 1104030205姓名： 高二奎 課程設(shè)計名稱： DSP原理及應(yīng)用 設(shè)計題目： IIR數(shù)字濾波器的DSP實現(xiàn) 完成期限：自 2015 年 1 月 1 日至 2015 年 1 月 8 日共 1 周設(shè)計依據(jù)、要求及主要內(nèi)容（可另加附頁）：1、 設(shè)計目的：通過課程設(shè)計，使學(xué)生綜

2、合運用DSP技術(shù)課程和其他有關(guān)先修課程的理論和生產(chǎn)實際知識去分析和解決具體問題的能力得到提高，并使其所學(xué)知識得到進一步鞏固、深化和發(fā)展。通過課程設(shè)計初步培養(yǎng)學(xué)生對工程設(shè)計的獨立工作能力，學(xué)習(xí)設(shè)計的一般方法。通過課程設(shè)計樹立正確的設(shè)計思想，提高學(xué)生分析問題、解決問題的能力。通過課程設(shè)計訓(xùn)練學(xué)生的設(shè)計基本技能，如計算、繪圖、查閱設(shè)計資料和手冊，熟悉標(biāo)準與規(guī)范等。2、 要求：1熟悉DSP處理器及其結(jié)構(gòu)性能，掌握DSP芯片配套開發(fā)工具的使用方法。2按要求設(shè)計出硬件電路。3畫出硬件連接原理圖，并對硬件工作原理進行說明。4給出軟件流程圖及編寫程序，每一條指令的后面附上相應(yīng)的注釋。5進行軟、硬件調(diào)試，檢查是

3、否達到相關(guān)的功能。6寫出調(diào)試方法。7設(shè)計報告結(jié)尾附上心得體會。3、 主要內(nèi)容：熟悉5410DSP的MCBSP的使用，了解AD50的結(jié)構(gòu)，掌握AD50各寄存器的意義及其設(shè)置，掌握AD50與DSP的接口，AD50的通訊格式及AD50的DA實驗。指導(dǎo)教師（簽字）： 批準日期： 年 月 日目 錄一、 摘要3二、 數(shù)字濾波器介紹和IIR數(shù)字濾波器的理論分析32.1數(shù)字濾波器介紹32.2 IIR濾波器的設(shè)計方法及原理4三、 DSP軟件的簡單介紹63.1DSP系統(tǒng)的特點63.2DSP系統(tǒng)的設(shè)計流程73.3DSP系統(tǒng)的開發(fā)工具CCS7四、 TMS320C5402的介紹9五、 基于DSP的IIR程序設(shè)計及軟件調(diào)

4、試115.1IIR程序設(shè)計115.2軟件調(diào)試結(jié)果13六、 總結(jié)和心得15七、參考文獻一、摘要隨著計算機和信息技術(shù)的飛速發(fā)展，數(shù)字信號處理已經(jīng)成為高速實時處理的一項關(guān)鍵技術(shù)，廣泛應(yīng)用在語音識別、智能檢測、工業(yè)控制等各個領(lǐng)域。數(shù)字濾波器是對數(shù)字信號實現(xiàn)濾波的線性時不變系統(tǒng)。數(shù)字濾波實質(zhì)上是一種運算過程，實現(xiàn)對信號的運算處理。DSP數(shù)字信號處理(Digital Signal Processing，簡稱DSP)是一門涉及許多學(xué)科而又廣泛應(yīng)用于許多領(lǐng)域的新興學(xué)科。20世紀60年代以來，隨著計算機和信息技術(shù)的飛速發(fā)展，數(shù)字信號處理技術(shù)應(yīng)運而生并得到迅速的發(fā)展。傳感器數(shù)字信號處理是利用傳感器對模擬信號或數(shù)字

5、信號進行采集并把其轉(zhuǎn)換成計算機可識別的電信號，并利用計算機對信號進行處理以達到計算機輔助控制或是計算機自動控制的目的。DSP 芯片是一種特別適合數(shù)字信號處理運算的微處理器，主要用來實時、快速地實現(xiàn)各種數(shù)字信號處理算法。用DSP 芯片實現(xiàn)IIR 數(shù)字濾波器，不僅具有精確度高、不受環(huán)境影響等優(yōu)點，而且因DSP 芯片的可編程性，可方便地修改濾波器參數(shù)，從而改變?yōu)V波器的特性，設(shè)計十分靈活。本課題主要應(yīng)用MATLAB軟件設(shè)計IIR數(shù)字濾波器，并對所設(shè)計的濾波器進行仿真；應(yīng)用DSP集成開發(fā)環(huán)境CCS調(diào)試匯編程序，文章結(jié)合TM320C5509的結(jié)構(gòu)特點，介紹了一種IIR濾波器在TM320C5509中的實現(xiàn)方

6、法。文中程序已經(jīng)過硬件驗證，仿真結(jié)果表明該設(shè)計符合要求。 關(guān)鍵詞 數(shù)字濾波；IIR；DSP；TM320C5509；MATLAB二、數(shù)字濾波器介紹和IIR數(shù)字濾波器的理論分析2.1數(shù)字濾波器介紹數(shù)字濾波器是對數(shù)字信號實現(xiàn)濾波的線性時不變系統(tǒng)。數(shù)字濾波實質(zhì)上是一種運算過程，實現(xiàn)對信號的運算處理。輸入數(shù)字信號（數(shù)字序列）通過特定的運算轉(zhuǎn)變?yōu)檩敵龅臄?shù)字序列，因此，數(shù)字濾波器本質(zhì)上是一個完成特定運算的數(shù)字計算過程，也可以理解為是一臺計算機。描述離散系統(tǒng)輸出與輸入關(guān)系的卷積和差分方程只是給數(shù)字信號濾波器提供運算規(guī)則，使其按照這個規(guī)則完成對輸入數(shù)據(jù)的處理。時域離散系統(tǒng)的頻域特性:,其中、分別是數(shù)字濾波器的輸

7、出序列和輸入序列的頻域特性（或稱為頻譜特性）,是數(shù)字濾波器的單位取樣響應(yīng)的頻譜，又稱為數(shù)字濾波器的頻域響應(yīng)。輸入序列的頻譜經(jīng)過濾波后,因此，只要按照輸入信號頻譜的特點和處理信號的目的， 適當(dāng)選擇，使得濾波后的滿足設(shè)計的要求，這就是數(shù)字濾波器的濾波原理。數(shù)字濾波器根據(jù)其沖激響應(yīng)函數(shù)的時域特性，可分為兩種，即無限長沖激響應(yīng)(IIR)數(shù)字濾波器和有限長沖激響應(yīng)(FIR)數(shù)字濾波器。IIR 數(shù)字濾波器的特征是，具有無限持續(xù)時間沖激響應(yīng)，需要用遞歸模型來實現(xiàn)，其差分方程為：系統(tǒng)函數(shù)為：設(shè)計IIR濾波器的任務(wù)就是尋求一個物理上可實現(xiàn)的系統(tǒng)函數(shù)H(z)，使其頻率響應(yīng)H(z)滿足所希望得到的頻域指標(biāo)，即符合給

8、定的通帶截止頻率、阻帶截止頻率、通帶衰減系數(shù)和阻帶衰減系數(shù)。2.2 IIR濾波器的設(shè)計方法及原理IIR濾波器差分方程的一般表達式為:yn=i=0Nbixn-i-i=1Maiy(n-i)式中x(n)為輸入序列;y(n)為輸出序列; ai和bi為濾波器系數(shù).若所有系數(shù)ai等于0,則為FIR濾波器.IIR濾波器具有無限長的單位脈沖響應(yīng),在結(jié)構(gòu)上存在反饋回路,具有遞歸性,即IIR濾波器的輸出不僅與輸入有關(guān),而且與過去的輸出有關(guān).將上式展開得出y(n)表達式為:yn=b0xn+b1xn-1+bNxn-N-a1yn-1-a2yn-2-aMyn-M在零初始條件下,對上式進行z變換,得到:Yz=b0Xz+b1

9、z-1Xz+bNz-NXz-a1z-1Yz-a2z-2Yz-aMz-MYz設(shè)N=M,則傳遞函數(shù)為:Hz=Y(z)X(z)=b0+b1z-1+bNz-N1+a1z-1+aNz-N上式可寫成:Hz=b0zN+b1zN-1+bNzN+a1z-1+aN=Ci=1Nz-ziz-pi上式具有N個零點zi和N個極點pi.若有極點位于單位圓外將導(dǎo)致系統(tǒng)不穩(wěn)定.由于FIR濾波器所有的系數(shù)ai均為0,不存在極點,不會造成系數(shù)的不穩(wěn)定.對于IIR濾波器,系統(tǒng)穩(wěn)定的條件如下:若|pi|1,當(dāng)n時,h(n),系統(tǒng)不穩(wěn)定.IIR濾波器具有多種形式,主要有:直接型(也稱直接I型)、標(biāo)準型(也稱直接II型)、變換型、級聯(lián)型和

10、并聯(lián)型.三、DSP軟件的簡單介紹3.1DSP系統(tǒng)的特點DSP系統(tǒng)是以數(shù)字信號處理為基礎(chǔ)的，因此不但具有數(shù)字處理的全部優(yōu)點而且還具有以下特點24:1接口方便：DSP應(yīng)用系統(tǒng)與其他以現(xiàn)代數(shù)字技術(shù)為基礎(chǔ)的系統(tǒng)或設(shè)備都是相互兼容的，這樣的系統(tǒng)接口以實現(xiàn)某種功能要比模擬系統(tǒng)與這些系統(tǒng)接口要容易得多。2編程方便：DSP應(yīng)用系統(tǒng)中的可編程DSP芯片，能靈活方便地進行修改和升級。3穩(wěn)定性好：DSP應(yīng)用系統(tǒng)以數(shù)字處理為基礎(chǔ)，受環(huán)境溫度及噪聲的影響較小、可靠性高，無器件老化現(xiàn)象。4精度高：16位數(shù)字系統(tǒng)可以達到10-5級的精度。5可重復(fù)性好：模擬系統(tǒng)的性能受元器件參數(shù)性能變化的影響比較大，而數(shù)字系統(tǒng)基本不受影響，

11、因此數(shù)字系統(tǒng)便于測試、調(diào)試和大規(guī)模生產(chǎn)。6集成方便：DSP應(yīng)用系統(tǒng)中的數(shù)字部件有高度的規(guī)范性，便于大規(guī)模集成。當(dāng)然，數(shù)字信號處理也存在一些缺點。例如，對于簡單信號處理任務(wù)，若采用DSP則使成本增加。DSP系統(tǒng)中的高速時鐘可能帶來高頻干擾和電磁泄漏等問題，而且DSP系統(tǒng)消耗的功率也較大。此外，DSP技術(shù)更新速度快，對于數(shù)學(xué)知識要求高，開發(fā)和測試工具還有待進一步完善。3.2DSP系統(tǒng)的設(shè)計流程一個DSP系統(tǒng)的設(shè)計過程大概要有以下幾個步驟。1根據(jù)系統(tǒng)的任務(wù)要求，確定系統(tǒng)處理精度要求、速度要求、實時性要求等性能指標(biāo)。2根據(jù)系統(tǒng)的要求進行高級語言的算法模擬，比如使用MATLAB等仿真工具，驗證算法的可行

12、性，得出最佳的處理方法。3DSP的系統(tǒng)設(shè)計，主要分為硬件設(shè)計和軟件設(shè)計。硬件設(shè)計是指根據(jù)系統(tǒng)要求選擇合適的DSP芯片，然后設(shè)計相應(yīng)的外圍電路。軟件設(shè)計主要是指根據(jù)系統(tǒng)的要求和選用的DSP芯片編寫相應(yīng)的程序。程序的編寫可以使用匯編語言，匯編語言編寫的程序效率高，但比較煩雜；也可采用C語言，DSP的C語言基本上是標(biāo)準C語言，編寫比較簡單，但效率低。在實際系統(tǒng)開發(fā)時往往是兩種語言結(jié)合編寫，在算法運算量大的地方使用匯編語言，在運算量小的地方使用C語言，這樣既能縮短軟件的開發(fā)周期，提高程序的可讀性和可移植性，又滿足了系統(tǒng)的實時性要求。本文的設(shè)計采用匯編語言編寫、設(shè)計軟件程序。3.3DSP系統(tǒng)的開發(fā)工具C

13、CSCCS是一種針對TMS320系列DSP的集成開發(fā)環(huán)境,在Windows操作系統(tǒng)下，采用圖形接口界面，提供有環(huán)境配置、源文件編輯、程序調(diào)試、跟蹤和分析等工具。CCS有兩種工作模式，即軟件仿真器模式：可以脫離DSP芯片，在PC機上模擬DSP的指令集和工作機制，主要用于前期算法實現(xiàn)和調(diào)試。硬件在線編程模式：可以實時運行在DSP芯片上,與硬件開發(fā)板相結(jié)合在線編程和調(diào)試應(yīng)用程序。 CCS的開發(fā)系統(tǒng)主要由以下組件構(gòu)成： 1. TMS320C54x集成代碼產(chǎn)生工具；2. CCS集成開發(fā)環(huán)境；3. DSP/BIOS實時內(nèi)核插件及其應(yīng)用程序接口API；4. 實時數(shù)據(jù)交換的RTDX插件以及相應(yīng)的程序接口API

14、；5. 由TI公司以外的第三方提供的各種應(yīng)用模塊插件。 CCS的功能十分強大，它集成了代碼的編輯、編譯、鏈接和調(diào)試等諸多功能，而且支持C/C+和匯編的混合編程，其主要功能如下： 1具有集成可視化代碼編輯界面，用戶可通過其界面直接編寫C、匯編、.cmd文件等；2含有集成代碼生成工具，包括匯編器、優(yōu)化C編譯器、鏈接器等，將代碼的編輯、編譯、鏈接和調(diào)試等諸多功能集成到一個軟件環(huán)境中；3高性能編輯器支持匯編文件的動態(tài)語法加亮顯示，使用戶很容易閱讀代碼，發(fā)現(xiàn)語法錯誤；4工程項目管理工具可對用戶程序?qū)嵭许椖抗芾?。在生成目?biāo)程序和程序庫的過程中，建立不同程序的跟蹤信息，通過跟蹤信息對不同的程序進行分類管理；

15、5基本調(diào)試工具具有裝入執(zhí)行代碼、查看寄存器、存儲器、反匯編、變量窗口等功能，并支持C源代碼級調(diào)試；6斷點工具，能在調(diào)試程序的過程中，完成硬件斷點、軟件斷點和條件斷點的設(shè)置；7探測點工具，可用于算法的仿真，數(shù)據(jù)的實時監(jiān)視等；8分析工具，包括模擬器和仿真器分析，可用于模擬和監(jiān)視硬件的功能、評價代碼執(zhí)行的時鐘；9數(shù)據(jù)的圖形顯示工具，可以將運算結(jié)果用圖形顯示,包括顯示時域/頻域波形、眼圖、星座圖、圖像等，并能進行自動刷新； 10提供GEL工具。利用GEL擴展語言，用戶可以編寫自己的控制面板/菜單，設(shè)置GEL菜單選項，方便直觀地修改變量，配置參數(shù)等；11支持多DSP的調(diào)試；12支持RTDX技術(shù)，可在不中

16、斷目標(biāo)系統(tǒng)運行的情況下，實現(xiàn)DSP與其他應(yīng)用程序的數(shù)據(jù)交換；13提供DSP/BIOS工具，增強對代碼的實時分析能力。四、TMS320C5402的介紹TMS320C5509采用雙乘累加單元(MAC)結(jié)構(gòu)。整個處理器內(nèi)部分為5個大的功能單元：存儲器緩沖單元(M)、指令緩沖單元(I)、程序控制單元(P)、地址生成單元(A)和數(shù)據(jù)計算單元(D)，各個功能單元之間通過總線連接。TMS320C5509中共有12條總線：1條32位程序數(shù)據(jù)總線(PB)，1條24位程序地址總線(PAB)，5條16位的數(shù)據(jù)總線(BB、CB、DB、EB、FB)和5條24位的數(shù)據(jù)地址總線(BAB、CAB、DAB、EAB、FAB)。作

17、為嵌入式芯片的一種，DSP芯片是一種非常適合于進行數(shù)字信號處理的微處理器芯片，已經(jīng)廣泛應(yīng)用于實現(xiàn)各種數(shù)字信號處理運算。其顯著特點可以歸納如下：1哈佛結(jié)構(gòu)哈佛結(jié)構(gòu)是不同于傳統(tǒng)的馮諾曼(Von Neuman)結(jié)構(gòu)的并行體系結(jié)構(gòu)，其主要特點是將程序和數(shù)據(jù)存儲在不同的存儲空間中，即程序存儲器和數(shù)據(jù)存儲器是兩個相互獨立的存儲器，每個存儲器獨立編址，獨立訪問。與兩個存儲器相對應(yīng)的是系統(tǒng)中設(shè)置了程序總線和數(shù)據(jù)總線兩條總線，從而使數(shù)據(jù)的吞吐率提高了一倍。2流水線與哈佛結(jié)構(gòu)相關(guān)，DSP芯片廣泛采用流水線以減少指令執(zhí)行時間，從而增強了處理器的處理能力。TMS320系列處理器的流水線深度從2-6級不等。第一代TMS

18、320處理器采用二級流水線，第二代采用三級流水線，而第三代則采用四級流水線。也就是說，處理器可以并行處理2-6條指令，每條指令處于流水線上的不同階段。3專用的硬件乘法器在一般形式的FIR濾波器中，乘法是DSP的重要組成部分。對每個濾波器抽頭，必須做一次乘法和一次加法。乘法速度越快，DSP處理器的性能就越高。在通用的微處理器中，乘法指令是由一系列加法來實現(xiàn)的，故需許多個指令周期來完成。相比而言，DSP芯片的特征就是有一個專用的硬件乘法器。4特殊的DSP指令DSP芯片的另一個特征是采用特殊的指令。例如TMS320C10中的LTD指令，可單周期完成加載寄存器、數(shù)據(jù)移動、同時累加操作。還有DMOV指令

19、，它完成數(shù)據(jù)移位功能。在數(shù)字信號處理中，延遲操作非常重要，這個延遲就是由DMOV指令來實現(xiàn)的9。5快速的指令周期哈佛結(jié)構(gòu)、流水線操作、專用的硬件乘法器、特殊的DSP指令再加上集成電路的優(yōu)化設(shè)計，可使DSP芯片的指令周期縮短到200ns以下?，F(xiàn)在，許多DSP處理器的指令周期已經(jīng)從第一代的200ns降低至現(xiàn)在的20ns以下，甚至在10ns以內(nèi)?？焖俚闹噶钪芷谑沟肈SP芯片能夠?qū)崟r實現(xiàn)許多DSP應(yīng)用。6面向寄存器和累加器DSP所使用的不是一般的寄存器文件，而是專用寄存器，較新的DSP產(chǎn)品都有類似于RISC的寄存器文件。許多DSP還有大的累加器，可以在異常情況下對數(shù)據(jù)溢出進行處理。7.支持前、后臺處理

20、DSP支持復(fù)雜的內(nèi)循環(huán)處理，包括建立起X、Y內(nèi)存和分址/循環(huán)計數(shù)器。一些DSP在做內(nèi)循環(huán)處理中把中斷屏蔽了，另一些則以類似后臺處理的方式支持快速中斷。許多DSP使用硬連線的堆棧來保存有限的上下文，而有些則用隱蔽的寄存器來加快上下文轉(zhuǎn)換時間。8.擁有簡便的單片內(nèi)存和內(nèi)存接口DSP設(shè)法避免了大型緩沖器或復(fù)雜的內(nèi)存接口，減少了內(nèi)存訪問。一些DSP的內(nèi)循環(huán)是在其單片內(nèi)存中重復(fù)執(zhí)行指令或循環(huán)操作部分代碼，它多采用SRAM而不是DRAM，因為前者接口更簡便。2000年3月，德州儀器（TI）公司推出了高性能低功耗的TMS320C55x芯片10。其性能可以達到400800MIPS，但功耗低到0.05 mW/M

21、IPS。TMS320系列包括定點、浮點和多處理器等三種類型的數(shù)字信號處理器。它的結(jié)構(gòu)是專門針對實時信號處理而設(shè)計的，具有指令靈活、可操作性強、速度快以及支持并行運算和C語言等特點，是性價比較高的一類DSP，在通信設(shè)備中得到了廣泛的應(yīng)用。C5000系列DSP是針對個人便攜設(shè)備而設(shè)計的，如音樂播放器、3G蜂窩電話、數(shù)碼相機、高速音頻設(shè)備、高精度的信號和多通道應(yīng)用。通常人們把C54x和C55x系列通稱C5000系列，主要用于功耗低、便攜式的無線通信終端產(chǎn)品。C5000系列包括舊有的C5x、當(dāng)前主流的C54x和較新的C55x系列。其中，C54x采用改進的哈佛結(jié)構(gòu)，并集成有豐富的硬件邏輯和外部接口資源，

22、具有較高的性能，及較低的成本和體積。C55x是在C54x的基礎(chǔ)上發(fā)展起來的，也采用改進型哈佛結(jié)構(gòu)，其器件功耗較C54x更低，性能更高。TMS320C55x采用多總線結(jié)構(gòu)，由五組內(nèi)部數(shù)據(jù)總線（3個用于讀，2個用于寫）和一組內(nèi)部程序存儲總線構(gòu)成。32bit的程序總線（PB）傳送從程序存儲器來的指令代碼和立即數(shù)。三組16bit數(shù)據(jù)讀總線（BB，CB和DB）連接數(shù)據(jù)讀地址產(chǎn)生邏輯。CB和DB總線傳送雙操作數(shù)。DB總線傳送單操作數(shù)。BB總線提供第三種讀路徑并且能夠提供雙乘操作數(shù)的系數(shù)。兩組16bit數(shù)據(jù)寫總線(EB，F(xiàn)B)連接數(shù)據(jù)寫地址產(chǎn)生邏輯。六組24bit數(shù)據(jù)總線（PAB，CAB，DAB，BAB，E

23、AB，F(xiàn)AB）傳送執(zhí)行指令所需要的地址。還有一條附加總線為DMA控制器和外設(shè)控制器提供服務(wù)。TMS320C55x的CPU結(jié)構(gòu)主要包括四種功能單元：即指令緩沖單元（I單元），程序流程單元（P單元），地址數(shù)據(jù)流程單元（A單元）和數(shù)據(jù)計算單元（D單元）。其數(shù)據(jù)計算單元包括一個40bit的可以提供-32到31移位范圍的桶形移位器，一個40bit的算術(shù)邏輯電路（ALU），兩個乘累加器（MAC）可以在一個周期中執(zhí)行兩個MAC操作，以及四個40bit的累加器。在其地址數(shù)據(jù)流程單元還包括一個16bit的算術(shù)邏輯電路（ALU），為主ALU提供簡單的算術(shù)運算。55x的指令集功能強大而且使用靈活，它同時具有易于使用

24、和程序效率高的特點。尋址方式包括絕對尋址、寄存器間接尋址、直接尋址，這些尋址方式降低了算法所需要的指令數(shù)量，減小了代碼量，也提高了芯片運算速度。特別要提到的是C55x指令集中所提供的三操作數(shù)指令，這種指令可以實現(xiàn)3個操作數(shù)的同時處理（如寫入寄存器或存儲器等操作），大大壓縮了代碼的指令數(shù)量。五、基于DSP的IIR程序設(shè)計及軟件調(diào)試5.1IIR程序設(shè)計IIR.c程序#include #define IIRNUMBER 2#define SIGNAL1F 1000#define SIGNAL2F 4500#define SAMPLEF 10000#define PI 3.1415926float I

25、nputWave();float IIR();float fBnIIRNUMBER= 0.0,0.7757 ;float fAnIIRNUMBER= 0.1122,0.1122 ;float fXnIIRNUMBER= 0.0 ;float fYnIIRNUMBER= 0.0 ;float fInput,fOutput;float fSignal1,fSignal2;float fStepSignal1,fStepSignal2;float f2PI;int i;float fIn256,fOut256;int nIn,nOut;main()nIn=0; nOut=0;fInput=fOutp

26、ut=0;f2PI=2*PI;fSignal1=0.0;fSignal2=PI*0.1;/fStepSignal1=2*PI/30;/fStepSignal2=2*PI*1.4;fStepSignal1=2*PI/50;fStepSignal2=2*PI/2.5;while ( 1 )fInput=InputWave();fInnIn=fInput;nIn+; nIn%=256;fOutput=IIR();fOutnOut=fOutput;nOut+;/ break pointif ( nOut=256 )nOut=0;float InputWave()for ( i=IIRNUMBER-1;

27、i0;i- )fXni=fXni-1;fYni=fYni-1;fXn0=sin(double)fSignal1)+cos(double)fSignal2)/6.0;fYn0=0.0;fSignal1+=fStepSignal1; if ( fSignal1=f2PI )fSignal1-=f2PI;fSignal2+=fStepSignal2;if ( fSignal2=f2PI )fSignal2-=f2PI;return(fXn0);float IIR()float fSum;fSum=0.0;for ( i=0;i PRAM PAGE 0.data : PRAM PAGE 0.cinit : PRAM PAGE 0.switch : P