DSP應(yīng)用技術(shù)課程設(shè)計(jì)

1、課 程 設(shè) 計(jì) 報(bào) 告（理工類）課程名稱：DSP應(yīng)用技術(shù)課程設(shè)計(jì)專業(yè)班級：08M通信工程專轉(zhuǎn)本1班 學(xué)生學(xué)號： 0821511027 學(xué)生姓名： 李惠敏 所屬院部： 龍?bào)磳W(xué)院 指導(dǎo)教師： 姜志鵬 2009 2010 學(xué)年 第 1 學(xué)期 金陵科技學(xué)院教務(wù)處制實(shí)驗(yàn)項(xiàng)目名稱： 基于DSP的數(shù)字濾波器設(shè)計(jì)與仿真 實(shí)驗(yàn)學(xué)時(shí)： 2周（分散40學(xué)時(shí)） 同組學(xué)生姓名： 實(shí)驗(yàn)地點(diǎn)： B402 實(shí)驗(yàn)日期： 2009.11.16-11.27 實(shí)驗(yàn)成績： 批改教師： 批改時(shí)間： 一、課程設(shè)計(jì)的目的和要求通過課程設(shè)計(jì)，加深對DSP芯片TMS320C54x的結(jié)構(gòu)、工作原理的理解，獲得DSP應(yīng)用技術(shù)的實(shí)際訓(xùn)練，掌握設(shè)計(jì)較復(fù)

2、雜DSP系統(tǒng)的基本方法。通過使用匯編語言編寫具有完整功能的圖形處理程序或信息系統(tǒng)，使學(xué)生加深對所學(xué)知識的理解，進(jìn)一步鞏固匯編語言講法規(guī)則。學(xué)會編制結(jié)構(gòu)清晰、風(fēng)格良好、數(shù)據(jù)結(jié)構(gòu)適當(dāng)?shù)膮R編語言程序，從而具備解決綜合性實(shí)際問題的能力。二、系統(tǒng)功能介紹及總體設(shè)計(jì)方案1、系統(tǒng)功能介紹一個(gè)實(shí)際的應(yīng)用系統(tǒng)中，總存在各種干擾。數(shù)字濾波器在語音信號處理、信號頻譜估計(jì)、信號去噪、無線通信中的數(shù)字變頻以及圖像信號等各種信號處理中都有廣泛的應(yīng)用，數(shù)字濾波器也是使用最為廣泛的信號處理算法之一。在本設(shè)計(jì)中，使用MATLAB模擬產(chǎn)生合成信號，然后利用CCS進(jìn)行濾波。設(shè)定模擬信號的采樣頻率為32000Hz，。設(shè)計(jì)一個(gè)FIR低

3、通濾波器，其參數(shù)為：濾波器名稱： FIR低通濾波器采樣頻率： Fs=40000Hz通帶： 4000Hz4500Hz過渡帶： 2500Hz3000Hz，3500Hz4000Hz帶內(nèi)波動： 0.5dB阻帶衰減： 50dB濾波器級數(shù)： N=154濾波器系數(shù)： 由MATLAB根據(jù)前述參數(shù)求得。2、總體設(shè)計(jì)方案流程圖 使用MATLAB編寫程序產(chǎn)生待濾波數(shù)據(jù)使用MATLAB濾波器的設(shè)計(jì)分析工具FDATOOL設(shè)計(jì)FIR濾波器使用CCS編寫程序?qū)崿F(xiàn)濾波功能查看程序運(yùn)行后的輸入輸出波形并比較濾波效果 圖1 總體設(shè)計(jì)方案三、主要內(nèi)容和步驟1、濾波器原理對于一個(gè)FIR濾波器系統(tǒng)，它的沖擊響應(yīng)總是又限長的，其系統(tǒng)函數(shù)

4、可記為：其中N-1是FIR的濾波器的階數(shù)，為延時(shí)結(jié)，為端口信號函數(shù)。最基本的FIR濾波器可用下式表示：其中輸入采樣序列，是濾波器系數(shù)，N是濾波器的階數(shù)Y（n）表示濾波器的輸出序列，也可以用卷積來表示輸出序列y(n)與x(n)、h(n)的關(guān)系，如下： y(n)=x(n)*h(n)2、DSP 實(shí)現(xiàn)FIR濾波的關(guān)鍵技術(shù)用線性緩沖區(qū)實(shí)現(xiàn)FIR濾波器 程序清單： .title “FIR1.ASM” .mmregs .def start x .usect “x”，8 PA0 .set 0 PA1 .set 1 .data COEF： .word 1*32768/10 .word 2*32768/10 .w

5、ord -4*32768/10 .word 3*32768/10 .word -4*32768/10 .word 2*32768/10 .word 1*32768/10 .textstart: SSBX FRCT STM #x+7,AR2 STM #6,AR0 LD #x+1,DP PORTR PA1,x+1FIR1: RPTZ A,#6 MACD *AR2-,COEF,A STH A,*AR2 PORTW *AR2+,PA0 BD FIR1 PORTR PA1,*AR2+0 .end 用循環(huán)緩沖區(qū)實(shí)現(xiàn)FIR濾波器 程序清單： .title “FIR2.ASM” .mmregs .def st

6、art .bss y,1xn .usect “xn”,7b0 .usect “b0”,7PA0 .set 0PA1 .set 1 .datatable: .word 1*32768/10 .word 2*32768/10 .word 3*32768/10 .word 4*32768/10 .word 5*32768/10 .word 6*32768/10 .word 7*32768/10 .text start: SSBX FRCT STM #b0,AR1 RPT #6 MVPD table,*AR1+ STM #xn+6,AR2 STM #b0+6,AR3 STM #7,BK STM #-1

7、,AR0 LD #xn,DP PORTR PA1,xnFIR2: RPTZ A,#6 MAC *AR2+0%,*AR3+0%,A STH A,y PORTW y,PA0 BD FIR2 PORTR PA1,*AR2+0% .end 3.操作步驟（1）打開FDATOOL,根據(jù)濾波要求設(shè)置濾波器類型、通帶截止頻率、指定階數(shù)、采樣頻率等。指定完設(shè)計(jì)參數(shù)后單擊按鈕Design Filter，生成濾波器系數(shù)。（2）把生成的濾波器系數(shù)傳到目標(biāo)DSP。選擇菜單Targets-Export to Code Composer Studio(tm)IDE，打開Export to C Header File對話框，

8、選擇C header file，指定變量名(濾波器階數(shù)和系數(shù)向量)，輸出數(shù)據(jù)類型可選浮點(diǎn)型或32 b，16 b整型等，根據(jù)自己安裝選擇目標(biāo)板板號和處理器號，單擊OK，保存該頭文件，需指定文件名(filtercoeffh)和路徑(保存在c：timyprojectsfir工程中)。（3）修改CCS匯編程序，刪掉數(shù)據(jù)前的所有文字，在開頭加上.data，第二行加coeff .word,在每行的前面加上.word，比且把每行的最后的逗號去掉。（4）編譯匯編程序，如果有錯(cuò)誤，按錯(cuò)誤進(jìn)行修改；沒錯(cuò)誤，則往下執(zhí)行。（5）加載初始化DATA數(shù)據(jù)。運(yùn)行程序，查看輸入輸出波形，修改相應(yīng)參數(shù)進(jìn)行調(diào)試四、詳細(xì)設(shè)計(jì)在本實(shí)

9、驗(yàn)中使用MATLAB模擬產(chǎn)生信號，觀察濾波前的時(shí)域波形和頻域波形。MATLAB仿真后，使用得到的濾波器參數(shù)，進(jìn)行DSP編程，在DSP中實(shí)現(xiàn)帶通濾波，并使用CCS的頻譜分析功能，查看DSP的濾波效果。1. MATLAB程序流程圖如圖1所示 初始化 產(chǎn)生需要慮波信號 畫出原始信號頻譜圖 生成fir濾波器 進(jìn)行濾波 畫出濾波后信號頻譜圖圖1 MATLAB程序流程圖2. CCS匯編程序流程圖如圖2所示 本次濾波結(jié)束待濾波系數(shù)減一為原始信號和系數(shù)指定寄存器 定義濾波器系數(shù) 輸入原始信號 初始化 做濾波算法是否有待濾波數(shù)據(jù)是 濾波結(jié)束是否還有待濾波數(shù)據(jù) 圖2 CCS匯編程序流程圖五、實(shí)驗(yàn)過程 1.匯編語言

10、1、實(shí)驗(yàn)步驟與內(nèi)容MATLAB輔助DSP 實(shí)現(xiàn)FIR ,其總體過程為在DSP 中編寫處理程序，在MATLAB中利用濾波器設(shè)計(jì)、分析工具( FDATOOL) ,根據(jù)指定的濾波器性能快速設(shè)計(jì)一個(gè)FIR ,然后把濾波器系數(shù)以頭文件形式導(dǎo)入CCS 中,頭文件中MATLAB 輔助DSP 實(shí)現(xiàn)FIR 數(shù)字濾波器含濾波器階數(shù)和系數(shù)數(shù)組,在MATLAB中調(diào)試、運(yùn)行DSP 程序并顯示、分析處理后的數(shù)據(jù)。使用該方法,便于采用匯編語言來實(shí)現(xiàn)程序。頭文件名不變,當(dāng)MATLAB中設(shè)計(jì)的濾波器系數(shù)改變時(shí),相應(yīng)頭文件中系數(shù)也改變,方便了程序調(diào)試、仿真。（一） 在MATLAB中先編寫程序，查看待濾波和已濾波的信號的時(shí)域和頻域

11、波形，如圖1和圖2所示。利用FDA TOOL 設(shè)計(jì)FIR 濾波器的參數(shù)，操作步驟如圖3所示。打開Export，把Numerator改為B,如圖4所示。圖6 EXPORT設(shè)置已濾波的信號的時(shí)域和頻域波形，如圖5和圖6所示。MATLAB濾波器設(shè)計(jì)工具在完成FIR設(shè)計(jì)后,如圖7所示，在Export as中選擇16bit符號整數(shù)輸出，然后單擊OK按鈕。生成fir.h C語言頭文件。將系數(shù)稍作修改，調(diào)整后拷貝到程序的系數(shù)初始化空間即可。圖7 MATLAB輸出數(shù)值類型轉(zhuǎn)換（二）在CCS中編寫匯編語言程序，進(jìn)行調(diào)試，實(shí)現(xiàn)帶通濾波的功能。在CCS IDE 中建立LHM.pjt工程,用匯編語言編寫處理主程序fi

12、r.asm。另外根據(jù)板上的存儲器配置方式,編寫存儲器配置文件fir.cmd文件,將MATLAB生成的LHM.h和input1.dat文件拷貝到LHM.pjt工程文件夾下，進(jìn)行編譯、鏈接,生成可執(zhí)行文件LHM.out。 （三）加載初始化data數(shù)據(jù)，圖8所示。圖8 data初始化運(yùn)行程序，查看輸入輸出波形，修改相應(yīng)參數(shù)進(jìn)行調(diào)試。2、實(shí)驗(yàn)過程中出現(xiàn)的錯(cuò)誤及解決的辦法（1）在MATLAB程序設(shè)計(jì)中，采樣頻率設(shè)置的過小，截止頻率大于采樣頻率的一半，運(yùn)行的時(shí)候圖形出現(xiàn)錯(cuò)誤。（2）FIR.m中的采樣頻率要和開始設(shè)置時(shí)的采樣頻率一致，否則結(jié)果會出現(xiàn)偏差。（3）在CCS中未定義標(biāo)號，程序運(yùn)行錯(cuò)誤。（4）在.h

13、文件中未把coeff頂格寫，編譯時(shí)出現(xiàn)錯(cuò)誤。（5）未將fir.h和input.dat文件考入工程文件中。（6）間接尋址過程中，程序中丟失。 (7) 在View的Graph中單擊Time/frequency出現(xiàn)graph property dialog框，未修改抽樣點(diǎn)數(shù)。顯示的圖形出現(xiàn)差異3、CCS程序運(yùn)行后的各種輸出結(jié)果在View的Graph中單擊Time/frequency出現(xiàn)graph property dialog框。將顯示類型，圖形名稱，起始地址，抽樣點(diǎn)數(shù)，數(shù)據(jù)類型等分別進(jìn)行設(shè)置，輸出各種波形。輸出濾波前的信號時(shí)域波形。首先，設(shè)置如圖9所示，輸出濾波前的信號時(shí)域波形。單擊OK后生成如下

14、圖10波形。圖10為濾波前的信號，波形很雜亂，從時(shí)域上很難看出信號的周期性。修改相應(yīng)設(shè)置，輸出如圖11所示濾波前信號頻譜。 經(jīng)過濾波后，時(shí)域波形如圖12所示。C語言1、實(shí)驗(yàn)步驟與內(nèi)容MATLAB輔助DSP 實(shí)現(xiàn)FIR ,其總體過程為在DSP 中編寫處理程序，在MATLAB中利用濾波器設(shè)計(jì)、分析工具( FDATOOL) ,根據(jù)指定的濾波器性能快速設(shè)計(jì)一個(gè)FIR ,然后把濾波器系數(shù)以頭文件形式導(dǎo)入CCS 中,頭文件中MATLAB 輔助DSP 實(shí)現(xiàn)FIR 數(shù)字濾波器含濾波器階數(shù)和系數(shù)數(shù)組,在MATLAB中調(diào)試、運(yùn)行DSP 程序并顯示、分析處理后的數(shù)據(jù)。使用該方法,便于采用C語言來實(shí)現(xiàn)程序。頭文件名不

15、變,當(dāng)MATLAB中設(shè)計(jì)的濾波器系數(shù)改變時(shí),相應(yīng)頭文件中系數(shù)也改變,方便了程序調(diào)試、仿真。（二） 在MATLAB中先編寫程序，查看待濾波和已濾波的信號的時(shí)域和頻域波形，如圖1和圖2所示。利用FDA TOOL 設(shè)計(jì)FIR 濾波器的參數(shù)，操作步驟如圖3所示。打開Export，把Numerator改為B,如圖4所示。已濾波的信號的時(shí)域和頻域波形，如圖5和圖6所示。MATLAB濾波器設(shè)計(jì)工具在完成FIR設(shè)計(jì)后,如圖7所示，在Export as中選擇16bit符號整數(shù)輸出，然后單擊OK按鈕。生成fir.h C語言頭文件。將系數(shù)稍作修改，調(diào)整后拷貝到程序的系數(shù)初始化空間即可。（二）在CCS中編寫匯編語言程

16、序，進(jìn)行調(diào)試，實(shí)現(xiàn)帶通濾波的功能。在CCS IDE 中建立LHM.pjt工程,用匯編語言編寫處理主程序fir.asm。另外根據(jù)板上的存儲器配置方式,編寫存儲器配置文件fir.cmd文件,將MATLAB生成的LHM.h和input1.dat文件拷貝到LHM.pjt工程文件夾下，進(jìn)行編譯、鏈接,生成可執(zhí)行文件LHM.out。 （三）加載初始化data數(shù)據(jù)，圖8所示。運(yùn)行程序，查看輸入輸出波形，修改相應(yīng)參數(shù)進(jìn)行調(diào)試。2、實(shí)驗(yàn)過程中出現(xiàn)的錯(cuò)誤及解決的辦法（1）在MATLAB程序設(shè)計(jì)中，采樣頻率設(shè)置的過小，截止頻率大于采樣頻率的一半，運(yùn)行的時(shí)候圖形出現(xiàn)錯(cuò)誤。（2）FIR.m中的采樣頻率要和開始設(shè)置時(shí)的采

17、樣頻率一致否則結(jié)果會出現(xiàn)偏差。（3）在CCS中未定義標(biāo)號，程序運(yùn)行錯(cuò)誤。（4）在.h文件中未把coeff頂格寫，編譯時(shí)出現(xiàn)錯(cuò)誤。（5）未將fir.h和input.dat文件考入工程文件中。（6）間接尋址過程中，程序中丟失。 (7) 在View的Graph中單擊Time/frequency出現(xiàn)graph property dialog框，未修改抽樣點(diǎn)數(shù)。顯示的圖形出現(xiàn)差異3、CCS程序運(yùn)行后的各種輸出結(jié)果在View的Graph中單擊Time/frequency出現(xiàn)graph property dialog框。將顯示類型，圖形名稱，起始地址，抽樣點(diǎn)數(shù)，數(shù)據(jù)類型等分別進(jìn)行設(shè)置，輸出各種波形。輸出濾波

18、前的信號時(shí)域波形，設(shè)置如圖所示，輸出濾波前的信號時(shí)域波形。單擊OK后生成如下圖10波形。圖10為濾波前的信號，波形很雜亂，從時(shí)域上很難看出信號的周期性。修改相應(yīng)設(shè)置，輸出如圖11所示濾波前信號頻譜。 經(jīng)過濾波后，時(shí)域波形如圖12所示。六、結(jié)論與體會通過本次課程設(shè)計(jì)，讓我了解了數(shù)字濾波是信號處理技術(shù)中的重要部分，研究了數(shù)字濾波器的基本理論知識以及它實(shí)現(xiàn)方法。學(xué)習(xí)了數(shù)字濾波器的結(jié)構(gòu)、設(shè)計(jì)理論，掌握了各種數(shù)字濾波器的原理和特性。并且復(fù)習(xí)了MATLAB的相關(guān)知識，并且利用MATLAB與CCS設(shè)計(jì)了低通濾波器以及高通濾波器，還研究了如何在定點(diǎn)DSP中實(shí)現(xiàn)數(shù)字濾波器的算法，學(xué)到了16位定點(diǎn)DSP芯片TMS

19、320C5402的硬件、軟件結(jié)構(gòu)和特性，掌握了CCS環(huán)境下的程序開發(fā)方法、調(diào)試工具的使用及優(yōu)化級別的選擇等。. .研究了MATLAB環(huán)境下FIR數(shù)字濾波器的設(shè)計(jì)方法，利用MATLAB軟件編程實(shí)現(xiàn)FIR濾波器設(shè)計(jì)。利用DSP來快速設(shè)計(jì)FIR數(shù)字濾波器的方法，尋找系數(shù)的快速傳遞，MATLAB中調(diào)試仿真DSP程序。DSP作為一門新興學(xué)科，越來越引起人們的關(guān)注，目前已廣泛應(yīng)用在各個(gè)領(lǐng)域。20世紀(jì)60年代以來，隨著計(jì)算機(jī)和信息技術(shù)的飛速發(fā)展，數(shù)字信號處理技術(shù)應(yīng)運(yùn)而生并得到迅速的發(fā)展。數(shù)字信號處理是(DSP)利用計(jì)算機(jī)或?qū)Ｓ锰幚碓O(shè)備，以數(shù)字形式對信號進(jìn)行采集、變換、濾波、估值、增強(qiáng)、壓縮、識別等處理，以得

20、到符合們需要的信號形式。數(shù)字信號處理(DSP)是圍繞著數(shù)字信號處理的理論、實(shí)現(xiàn)和應(yīng)用等幾個(gè)方面發(fā)展起來的。數(shù)字信號處理在理論上的發(fā)展推動了數(shù)字信號處理應(yīng)用的發(fā)展。反過來，數(shù)字信號處理的應(yīng)用又促進(jìn)了數(shù)字信號處理理論的提高。而數(shù)字信號處理的實(shí)現(xiàn)則是理論和應(yīng)用之間的橋梁。數(shù)字信號處理DSP是以眾多學(xué)科為理論基礎(chǔ)的，它所涉及的范圍極其廣泛。所以在今后的工作和學(xué)習(xí)生活中，綜合運(yùn)用模擬電子、數(shù)字電子和DSP基本原理等課程中所學(xué)的理論知識去獨(dú)立完成一個(gè)項(xiàng)目的設(shè)計(jì)。不僅要完善以上的MATLAB設(shè)計(jì)方法，更要學(xué)會利用多種途徑設(shè)計(jì)各種DSP最小系統(tǒng)，熟練C語言編程。七、參考文獻(xiàn)1 李真芳，蘇濤等.DSP程序開發(fā)一

21、MATLAB調(diào)試及直接目標(biāo)代碼生成M.西安:西安電子科技大學(xué)出版社，2003.87一163.2 樓順天，李博函.基于MATLAB的系統(tǒng)分析與設(shè)計(jì)一信號處理M.西安:西安電子科技大學(xué)出版社，1998.81一88.3王秀芳，關(guān)凌濤.基于MatLab與DSP的濾波器的快速設(shè)計(jì)方法期刊文章.現(xiàn)代計(jì)算機(jī)，2008，總第二七七期：97994趙秀玲.基于定點(diǎn)DSP的FIR數(shù)字濾波器的研究與設(shè)計(jì)碩士論文.曲阜師范大學(xué)，2007.15325羅軍輝，羅勇江等.MATLAB在數(shù)字信號處理中的應(yīng)用M.北京:機(jī)械工業(yè)出版社，2005.63一128.附件：源程序清單匯編程序清單：lhm.h文件內(nèi)容： .datacoeff

22、 .word -85, -64, -61, -36, 8, 62, 110, 136, 131 .word 96, 42, -11, -44, -44, -11, 39, 85, 106 .word 88, 37, -29, -83, -101, -73, -7, 70, 124 .word 131, 82, -5, -96, -151, -144, -72, 37, 140 .word 191, 162, 61, -76, -191, -231, -174, -35, 132 .word 256, 278, 180, -7, -208, -337, -330, -176, 72 .wor 3

23、16, 445, 392, 156, -178, -475, -599, -470, -108 .word 361, 745, 857, 594, -6, -748, -1336, -1456, -893 .word 386, 2192, 4154, 5816, 6768, 6768, 5816, 4154, 2192 .word 386, -893, -1456, -1336, -748, -6, 594, 857, 745 .word 361, -108, -470, -599, -475, -178, 156, 392, 445 .word 316, 72, -176, -330, -3

24、37, -208, -7, 180, 278 .word 256, 132, -35, -174, -231, -191, -76, 61, 162 .word 191, 140, 37, -72, -144, -151, -96, -5, 82 .word 131, 124, 70, -7, -73, -101, -83, -29, 37 .word 88, 106, 85, 39, -11, -44, -44, -11, 42 .word 96, 131, 136, 110, 62, 8, -36, -61, -64 .word -85Fir.asm內(nèi)容：.titlefir.asm.mmr

25、egs.global_c_int00ORDER .set154D_LEN.set1024; The Length of Input Dataxn.usectxn,(ORDER-1)a0.usecta0,(ORDER-1)input.usectinput,D_LENoutput.usectoutput,D_LEN.copy HM.h; Get coefs from coef.h.text.asgAR0, FIR_INDEX.asgAR2, FIR_DATA.asgAR3, FIR_COEF.asgAR5, DATA_IN.asgAR6, DATA_OUT_c_int00:SSBX FRCT ST

26、M #a0,FIR_COEF; Copy a0(coefs) to FIR_COEF(AR3) RPT #ORDER-1 MVPD #coeff,*FIR_COEF+ STM#1, FIR_INDEX STM#xn, FIR_DATA; Copy xn(data) to FIR_DATA(AR2) RPTZA,#ORDER-1 STLA, *FIR_DATA+STM #(xn+ORDER-1), FIR_DATASTM #(a0+ORDER-1), FIR_COEFSTM#input, DATA_IN; Get data from inputSTM#output, DATA_OUT; Writ

27、e data to outputSTM#D_LEN-1, BRCRPTBDnext-1STM#ORDER, BKLD*DATA_IN+, AFIR:STLA, *FIR_DATA+%RPTZA, (ORDER-1)MAC*FIR_DATA+0%, *FIR_COEF+0%, A; FIRSSTHA, *DATA_OUT+nextFIR_END:B FIR_END.endFIR.m的程序如下：fs=40000; %/采樣HzN=1024 %數(shù)據(jù)個(gè)數(shù)T=1/fs; %采樣周期n=0:N-1;df=n*(fs/N)%待濾波信號波形xin=randn(1,1024)figure(1)plot(xin)

28、%待濾波信號頻譜xinff=abs(fft(xin);figure(2)plot(df,xinff)%濾波后信號波形y_filter_out=filter(B,1,xin)% Y = FILTER(B,A,X)figure(3)plot(y_filter_out)%濾波后信號頻譜yff=fft(y_filter_out);figure(4)plot(df,yff)xin=xin/max(xin);%歸一化xto_ccs=round(32767*xin)%取整fid=fopen(input1.dat,w);%打開文件fprintf(fid,1651 1 0 0 0n);%輸出文件頭fprintf

29、(fid,%dn,xto_ccs);%輸出fclose(fid);C程序清單lhm.h文件內(nèi)容：/* * Filter Coefficients (C Source) generated by the Filter Design and Analysis Tool * * Generated by MATLAB(R) 7.0.1 and the Signal Processing Toolbox 6.2.1. * * Generated on: 11-Dec-2009 10:27:37 * */* * Discrete-Time FIR Filter (real) * - * Filter S

30、tructure : Direct-Form FIR * Filter Order : 110 * Stable : Yes * Linear Phase : Yes (Type 1) */* General type conversion for MATLAB generated C-code */#include C:MATLAB701externincludetmwtypes.h/* * Expected path to tmwtypes.h * C:MATLAB701externincludetmwtypes.h */* * Warning - Filter coefficients

31、were truncated to fit specified data type. * The resulting response may not match generated theoretical response. * Use the Filter Design & Analysis Tool to design accurate * int16 filter coefficients. */const int BL = 111;const int16_T B111 = -100, 654, -605, -216, 140, 234, 138, -25, -145, -158, -69, 63, 155, 145, 36, -104, -181, -138, 6, 159, 212, 121, -65, -223, -237, -85, 145, 295, 250, 21, -248, -370, -239, 79, 379, 443, 196, -227, -544, -511, -100, 452, 765, 569, -90, -823, -1105, -614, 502, 1602, 1844