1、FIR 高通數字濾波器的設計及 DSP 實現班級:學號:姓名:現場操作（30書面報告（70分）成績分）1.Matlab設計：1設計目的（FIR 高通濾波器設計）2. 設計要求：在 MATLA 課程中,熟悉軟件的使用以及 FIR 濾波器的設計原理,鍛煉程序的設計能 力和整體思想3. 設計內容：設計 FIR 高通濾波器：通帶截止頻率 15KHZ 采樣頻率 100KHZ4. 設計原理：數字濾波器的設計問題就是尋找一組系 ai 和 bi，使得其性能在某種意義上逼近所要求的特性。如果在s 平面上去逼近，就得到模擬濾波器，如果在z 平面上去逼Matlab設計（20分）DSP程序（20分）輸入輸出分析（20

2、分）格式規范（10總分 分）等級近，得到數字濾波器。數字濾波器是將輸入的信號序列，按規定的算法進行處理從而得到所期望的輸出序列。一個線性位移不變系統的輸出序列y(n 和輸入序列 x(n 之間的關系，應該滿足常系數線性差分方程x(n 為輸入序列，y( n 為輸出序列，ai、bi 為濾波器系數，N 為濾波器的階數。5.程序代碼：wc=0.3*pi;b=fir1(80,wc,high;freqz(b;figure;stem(b,.;lin e(0,25,0,0;xlabel( n;ylabel(h( n;M=ro un d(100000*b;6實現數據：51,-46,39,-31,20,-4,-17

3、,45,-80,123,-172,228,-287,347,-404,453,-490,507,-500,463,-390,277,-120,-84,335,-633,975,-1357,1772,-2213,2670,-3133,3591,-4031,4442,-4812,5131,-5390,5580,-5697,-5697,5580,-5390,5131,-4812,4442,-4031,3591,-3133,2670,-2213,1772,-1357,975,-633,335,-84,-120,277,-390,463,-500,507,-490,453,-404,347,-287,2

4、28,-172,123,-80,45,-17,-4,20,-31,39,-46,512.輸入信號產生：1信號函數：Y二si n(0.03*pi* n+si n(0.2*pi* n+si n(0.4*pi* n2.matlab程序：n=0:1:200;y=ro und (1000*(s in (0.03*pi* n+si n(0.2*pi* n+si n(0.4*pi*n +1000;plot(n,y;fid二fope n(s in .txt,w：fprin tf(fid,%dn,y;fclose(fid;3.產生的離散數據：500,1316,1363,821,502,727,950,625,7

5、3,106,905,1700,1722,1152,803,994,1181,818,227,219,976,17281708,1095,704,854,1000,599,-29,-71,655,1378,1332,697,287,422,557,149,-482,-524,206,939,905,287,-104,54,217,-162,-761,-767,0,772,778,201,-147,54,259,-77,-634,-600,206,1015,1057,512,194,422,650,334,-207,-160,655,1468,1510,963,637,854,1067,732,1

6、69,189,976,1757,1765,1181,817,994,1166,789,183,161,905,1644,1612,988,586,727,866,458,-176,-222,500,1222,1176,542,134,273,414,12,-612,-644,95,839,817,211,-166,6,183,-181,-765,-757,24,811,831,268,-67,146,363,37,-510,-468,345,1160,1207,666,350,578,806,488,-57,-15,794,1600,1634,1077,741,946,1147,799,222

7、,228,1000,1767,1761,1162,783,946,1104,713,95,61,794,1524,1482,851,443,578,713,303,-332,-378,345,1071,1029,401,0,146,296,-95,-708,-728,24,781,773,182,-181,6,197,-152,-722,-700,95,894,927,375,50,273,498,179,-363,-316三.DSP程序：ASM:.title fir.asm.mmregs.def start.bss y,1 保存濾波結果in data .usect buffer,1 輸入信號

8、緩沖區outdata .usect buffer輸出信號緩沖區xn .usect x n,80 ;xnh .usect h,80 ; h.datatable:.word 51,-46,39,-31,20,-4,-17,45.word -80,123,-172,228,-287,347,-404,453.word -490,507,-500,463,-390,277,-120,-84.word 335,-633,975,-1357,1772,-2213,2670,-3133.word 3591,-4031,4442,-4812,5131,-5390,5580,-5697.word -5697,5

9、580,-5390,5131,-4812,4442,-4031,3591.word -3133,2670,-2213,1772,-1357,975,-633,335.word -84,-120,277,-390,463,-500,507,-490.word 453,-404,347,-287,228,-172,123, -80.word 45,-17,-4,20,-31,39,-46,51.textstart:SSBX FRCT ;小數乘法STM #xn ,AR1RPT #79ST #0,*AR1+ ;把 x(n-x(n-9 賦初值 0STM #h,AR1RPT #79MVPD #table,

10、*AR1+ ;把參數表復制到數據存儲區STM #xn+79,AR3 ;AR3 的值給 x(n-79STM #h+79,AR4 ;AR4 的值給 h(n-79STM #80,BK ;循環緩沖區大小 80STM #-1,AR0 ;指針調整值-1STM #i ndata,AR5STM #outdata,AR2;LD #xn,DP ;DP 指向 xn 所在頁;PORTR PA1,xn ;輸入數據;LD #y,DP ;DP 指向 y 所在頁FIR:NOPMVDD *AR5,*AR3+0% ;設置斷點和探針從 indat 里面讀入數據RPTZ A,#79MAC *AR3+0%,*AR4+0%,ASTH A

11、,*AR1MVDD *AR1,*AR2 ;設置斷點和探針從 outdat 里面寫入數據NOP;STH A,y;PORTW y,PA0B FIR;B FIR;PORTR PA1,*AR3+0%.endCMD:fir.obj -o fir.out-m fir.map-estartMEMORYPAGE 0:EPROM: org=0190H,le n=0280HPAGE 1:SPRAM: org=1000H,le n=1000HDARAM: org=2000H,le n=2000H SECTIONS.text :EPROM PAGE 0.data : EPROM PAGE 0.bss : SPRAM PAGE 1xn : alig n(128 DARAM PAGE 1h : align(128 DARAM PAGE 1buffer :DARAM PAGE 14.DSP仿真運行圖：5.輸入和輸出結果分析:1. 輸入數據是通過，matlab 產生的離散數據產生時是三個頻率的正弦波疊加的結 果，通過我們的仿真軟件我們可以發現三種波形的疊加效果如圖中上半部分所示，由于 matlab 產生的離散的點，所以這些波形的疊加就會出現這樣的情況。2. 本實