【正文】 ** Allocate 100 words in .bss for the **19 ** symbol named ARRAY。? z 運行連接器的選擇項。產生一個與中間文件格式相同的高效率版本。187 188源程序分析器 輸入是 C源文件。164Goertzel算法的特點v Goertzel算法的 IIR濾波器結構涉及兩個復共軛極點，但對于實際的 DTMF音頻檢測，只要有幅度信息就夠了（實際上使用幅度平方），從而簡化為只需要一個實系數參與計算。156 157CCITT關于 DTMF的指標v Low band： 69 770、 85 941（ Hz）RBW% v High band： 120 133 147 1633（ Hz） RBW% v TWIST ： STD 4dB REV 8dBv DYN RANGE： 25dB v GUARD TIME： MIN TONE TIME 45msv SNR： NOISE 24dBV 158算法設計 算法設計是非常關鍵的步驟，因為算法決定運算量和存儲器的大小，從而決定 DSP的選擇，進而決定系統(tǒng)的設計，包括硬件和軟件的設計。上面用多種不同的方法進行了濾波器的設計，在實現濾波器功能的前提下，程序的繁簡和對存儲器的使用情況是不同的，應用中應根據具體情況進行選擇。 此外，由于字長有限，每個二階基本節(jié)運算后都會帶來一定的誤差，合理安排各二階基本節(jié)的前后次序，將使系統(tǒng)的精度得到優(yōu)化。IIR濾波器系數 .word 63， 0， 188， 0， 188， 0， 63 ； b系數 .word 0， 11675， 0， 9663， 0， 2948 ； a系數 .data 143BN .usect BN,N+1AN .usect AN,N+1INBUF .usect INBUF,256 。 } fclose(fp)。for(i=0。39。 if((fp=fopen(,wt))==NULL) { printf(can39。rs=20 [n,wn]=cheb1ord(wp/8,ws/8,rp,rs) [b,a]=cheby1(n,rp,wn) 138設計結果為： N=3 wn = b0= a0= b1=0 a1=0 b2= a2= b3=0 a3=0 b4= a4= b5=0 a5=0 b6= a6=139（ 2）產生濾波器輸入信號的文件使用 CCS的 Simulator進 行 濾 波器特性的 測試時 ，需要 輸 入時間 信號 x(n)。這 4個循環(huán)緩沖區(qū)的結構如圖所示：137【例】 設計一個三階的切比雪夫 Ⅰ 型帶通數字濾波器，其采樣頻率 Fs=16kHz，其通頻帶 f，內損耗不大于 1dB； f f，其衰減大于20dB。但是，在有限精度的運算中，可能出現不穩(wěn)定現象。 新緩沖區(qū)向老緩沖區(qū)傳送一個數 。 重復執(zhí)行 (4)～ (9)步。(7)修正數據指針，讓 AR2和 AR3分別指向 New循環(huán)緩沖區(qū)最新的數據和 Old循環(huán)緩沖區(qū)中最老的數據。(3)在程序存儲器中設置系數表。如果利用對稱性，可將其改寫成： Y(n)=h0[x(n)+x(n7)]+h1[x(n1)+x(n6)]+h2[x(n2) +x(n5)]+h3[x(n3)+x(n4)] 變成 4次乘法和 7次加法，乘法運算的次數少了一半。120 Graph 屬性設置窗口121 輸入信號的時域波形 輸入信號的頻譜圖122 濾波器輸出信號時域波形 濾波器輸出信號頻譜圖123 FIR濾波器設計 系數 對 稱的 FIR濾 波器具有 線 性相位特性， 這 種 濾波器是用得最多的 FIR濾 波器，特 別 是 對 相位失真要求很高的 場 合。FIR濾波器階數COEF_FIR .sect COEF_FIR 。使用 CCS的 Simulator進行濾波器特性測試時，需要輸入時間信號 x(n)。乘法累加利用 MACD指令， 該 指令完成數據存 儲 器 單 元與程序存 儲 器 單 元相乘，并累加、移位的功能。 將系數 a0～ a4存放在數據存儲器中，并設置線性緩沖區(qū)存放輸入數據。它不斷地從輸入端讀入樣本值 x[n]，經延時 (z1)后做乘法累加，輸出濾波結果 y[n]。有兩條實現輸入和輸出的指令：PORTR PA, Smem ；將 PA端口的內容送 ；數據存儲器 SmemPORTW Smem, PA ；將地址為 Smem的數據 ；存儲器內容送端口 PA數據的輸入 /輸出 94 使用存儲器延時指令 DELAY，可以將數據存儲單元中的內容向較高地址的下一單元傳送。顯然，最主要的問題是延遲環(huán)節(jié)如何實現。在下一個循環(huán)中，pmad加 1。 filteringSTH A, *AR5+ 。 FIR circular buffer sizeLD *AR7+, A 。 load the input valuefir_filter:STL A, *FIR_DATA_P+% 。fs/2的區(qū)間66雙線性變換67有限沖擊響應（ FIR）濾波器= a0 + a1z1 + … + a Nz(N1)68FIR濾波器的結構69頻率響應70FIR濾波器的主要特點v單位沖擊響應只有有限多項v可以設計成線性相位系統(tǒng)v只在零點處有極點，因此系統(tǒng)總是穩(wěn)定的v便于 DSP實現并可用立即數乘加指令編程，（節(jié)約存儲器） 71系統(tǒng)非線性相移造成輸出信號失真v 系統(tǒng)相位特性決定了信號不同頻率的時延72 輸入波形 DFT變換忽略相位信息IDFT變換輸出波形忽略相位信息的后果73相位特性74要求線性相位的例子v通信系統(tǒng)：數據通信、調制解調器v希爾伯特變換器：要求輸入輸出信號正交。d(n2)*A2+d(n1)*A1DELAY *AR3 。input+d(n2)*B2MAC *AR4+, *AR3, A 。 IIR系統(tǒng)因為與過去某些輸出有關，所以網絡結構有反饋支路稱為遞歸結構， FIR系統(tǒng)只與激勵有關沒有反饋支路，稱為非遞歸結構。而濾波在數字信號處理中占有極其重要的地位，是語音處理、圖像處理、模式識別、頻譜分析等應用的基本處理算法。定義堆棧TOS .usect stack,1 .copy “” 。蝶形級數sav_grp .usect tempv,3 。283. 實序列 FFT匯編源程序及鏈接命令文件【例】 256點實序列 FFT的 DSP實現。26第三步，分離復數 FFT的輸出為奇部分和偶部分分離 FFT輸出為 RP、 RM、 IP和 IM四個序列，即偶實數、奇實數、偶虛數和奇虛數四部分。 DIT的復數乘法出現在加減之前； DIF的復數乘法出現在加減之后。傅立葉變換2 時域 —— 頻域 離散信號 周期譜 周期信號 離散譜 離散周期信號 離散周期譜34DFT的定義56DFT和 FFT計算量對比7 N點 蝶 2FFTx1[n] = x[2n]x2[n] = x[2n+1]將 N點的序列分為兩個 N/2點的序列89上式的最后一步是因為10v 我們已經將一個 N點的 DFT分解成為兩個 N/2點的 DFT。 作為一個非常重要的工具經常使用，甚至成為DSP運算能力 的一個考核標準。13將 N點 DFT分為兩個 N/2點 DFT14將 N/2點 DFT分為兩個 N/4點 DFT152點 DFT168點蝶 2FFT17按頻率抽?。?DIF） x1[n] = x[n] x2[n] = x[n+N/2] n=0， 1， …… ， N/211819DIT與 DIF對比 DIT輸入是混序的，頻域的輸出是順序的； DIF輸入是順序的，頻域的輸出是混序的。第二步， N點復數 FFT在 DATA數據處理緩沖器里進行 N點復數 FFT運算。第二種方法是建立與輸入信號對應的數據流文件，具體操作步驟見教材內容。復數點數 7 。信號功率譜BOS .usect stack,0Fh 。 余弦表系數由 ；出，起始地址標號 TWI2.text_c_int00b start nopnopstart:STM TOS,SPLD 0,DPSSBX FRCT …………31輸入信號時域波形、頻譜圖及其功率譜4. 觀察信號時域波形及其頻譜32二、數字濾波器的 DSP實現模擬濾波器回顧 數字信號處理主要是面向實時的密集型數據運算。54數字濾波器的差分方程分為無限沖激響應系統(tǒng) (IIR)和有限沖激響應系統(tǒng) (FIR)。AR4:coeff of IIR filter B2,B1,A2,A1,A0……….* Feedback pathMAC *AR4+, *AR3, A 。d(n2)*A2MAC *AR4+, *AR3, A 。d(n1)=d(n)………. 60級 連61多級二階節(jié)級連成多級 IIR62并 聯(lián)63IIR濾波器的設計(1)設計一個符合要求的模擬濾波器(2)從 S平面映射到 Z平面，得到相應的數字濾波器64沖激響應不變法 h(nts)= ha[k]65Z平面的單位圓映射為 S平面的 177。 FIR circular bffr sizeLD *INBUF_P+, A 。 Repeat 256 timesRPTBD fir_filter_loopSTM 16, BK 。 newest sampleRPTZ A, 15MAC *AR4+0%, *AR3+0%, A 。同時，存儲器操作數 Xmem和 Ymem相加，結果左移 16bits，然后裝入累加器 A。 FIR濾波器的 DSP實現90FIR濾波器的結構圖：由圖看出，實際算法由一組基本運算單元：乘法器、加法器、延遲器組成。(1)用 線 性 緩 沖區(qū)法 實現 z1 92N=6的線性緩沖區(qū)存儲器圖93C54x片內沒有 I/O資源， CPU通過外部譯碼可以尋址64K的 I/O空間。 (2)用循 環(huán)緩 沖區(qū)法 實現 z196N=6的循環(huán)緩沖區(qū)存儲器圖97 3． FIR濾波器的 DSP實現方法 FIR濾波器的輸出表達式寫為 y(n)=b0x(n)+b1x(n1)+……+ bn1x(nN+1)其結構圖FIR濾波算法實際上是一種乘累加運算。(1)用 線 性 緩 沖區(qū)和直接 尋 址方法 實現 FIR線性緩沖區(qū)安排：99 .title “” ；定義源程序名 .mmregs ；定義存儲器映象寄存器 .def start ；定義語句標號 start .bss y,1 ；為結果 y預留 1個單元的空間XN .usect “XN”,1 ；在自定義的未初始化段 “XN”XNM1 .usect “XN”,1 ；中保留 5個單元的空間XNM2 .usect “XN”,1XNM3 .usect “XN”,1 XNM4 .usect “XN”,1A0 .usect “A0”,1 ；在自定義的未初始化段 “A0”A1 .usect “A0”,1 ；中保留 5個單元的空間A2 .usect “A0”,1A3 .usect “A0”,1A4 .usect “A0”,1PA0 .set 0 ；定義 PA0為輸出端口PA1 .set 1 ；定義 PA1為輸入端口100 .datatable: .word 1*32768/1