1、實(shí)驗(yàn)六 IIR濾波器一、 實(shí)驗(yàn)?zāi)康?. 熟悉IIR濾波器C54X實(shí)現(xiàn)的編程方法。2. 測(cè)試IIR濾波器的單位沖擊響應(yīng)曲線。3. 檢查IIR濾波器的頻率特性。二、 實(shí)驗(yàn)條件1. 已經(jīng)設(shè)計(jì)出四階IIR濾波器的參數(shù)如下：通帶：0-200Hz過渡帶寬：200Hz-500Hz通帶內(nèi)波動(dòng)：0.5dB阻帶衰減：-20dB采樣頻率Fs：3600Hz脈沖傳遞函數(shù)H（z）： Hz=B0+B1Z-1+B2Z-2+B3Z-3+B4Z-41-A1Z-1-A2Z-2-A3Z-3-A4Z-4差分方程式為： y(n)=B0*x（n）+B1*x（n-1）+B2*x（n-2）+B3*x（n-3）+B4*x（n-4）+A1*y（n

2、-1）+A2*y（n-2）+A3*y（n-3）+A4*y（n-4） 其中：A1=-3.4647 A2=4.4615 A3=-2.8518 A4=0.6739 B0=0.0951 B1=-0.3139 B2=0.4460 B3=-0.3139 B4=0.0951生成正弦數(shù)據(jù)文件的高級(jí)語(yǔ)言程序。程序名為sin_flt.exe。2. 直接形式二階IIR濾波器程序Lab6.asm以及鏈接命令文件Lab6.cmd。三、 實(shí)驗(yàn)內(nèi)容1. 消化直接形式二階濾波器程序iir3.asm以及鏈接命令文件iir3.cmd。2. 修改Lab6.asm。定義循環(huán)緩沖區(qū)：Y(n-4)X(n-4)X： Y: Y(n-3)X(

3、n-3)Y(n-2)X(n-2)Y(n-1)X(n-1)Y(n)X(n) AR2AR4A4B4 B: A:A3B3A2B2A1B1A0B0 AR5AR3相應(yīng)的匯編命令為：X .usect X,5Y .usect Y,5B .usect B,5A .usect A,5修改系數(shù)表：初始值y（n-4）- y（n-1）以及x（n-4）- x（n-1）均設(shè)置為0。凡絕對(duì)值1的系數(shù)需化成絕對(duì)值1的系數(shù)。這樣，系數(shù)表示為：.word 3116 ;B4=0.0951.word -10286 ;B3=-0.3139.word 14615 ;B2=0.4460.word -10286 ;B1=-0.3139.wo

4、rd 3116 ;B0=0.0951.word -22082 ;A4=-0.6739.word 31149 ;A3/3=2.8518/3.word -30484 ;A2/5=-4.6515/5.word 28383 ;A1/4=3.4647/4修改傳送x的初始值、y的初始值、系數(shù)B、系數(shù)A的重復(fù)次數(shù)，分別為3、3、4、3。修改傳送數(shù)據(jù)的起始地址：傳送x MVPD #table,*AR1+傳送y MVPD #table+4,*AR1+傳送B MVPD #table+8,*AR1+傳送A MVPD #table+13,*AR1+修改BK寄存器，BK=5。為循環(huán)緩沖區(qū)設(shè)初始指針：AR2 X+4 AR

5、4 Y+4AR5 B+4 AR3 A+4由于歸一化成絕對(duì)值小于1的系數(shù)，因此反饋通道運(yùn)算必須作相應(yīng)的修正。對(duì)輸入數(shù)據(jù)文件中的數(shù)據(jù)定標(biāo)，右移4位，以防止運(yùn)算中溢出。3. 修改Lab.cmd文件。將X、Y、B、A的循環(huán)緩沖區(qū)均改為align(8)。4. 匯編和鏈接5. 測(cè)試本實(shí)驗(yàn)IIR濾波器的單位沖擊響應(yīng)曲線6. 檢查本實(shí)驗(yàn)IIR濾波器的頻率特性曲線四、 實(shí)驗(yàn)程序;*;* lab6.ASM IIR Filter * ;* .title lab6.asm .mmregs .def start .bss in,1 ;給in分配一個(gè)存儲(chǔ)單元 .bss out,1X .usect X,5 ;定義X的緩沖區(qū)

6、域Y .usect Y,5B .usect B,5A .usect A,5PA0 .set 0 ;I/O口地址賦值PA1 .set 1 .datatable: .word 0 ;X(N-4) ;為標(biāo)號(hào)table開始的存儲(chǔ)單元賦初值 .word 0 ;X(N-3) .word 0 ;X(N-2) .word 0 ;X(N-1) .word 0 ;Y(N-4) .word 0 ;Y(N-3) .word 0 ;Y(N-2) .word 0 ;Y(N-1) .word 3116 ;B4=0.0951 .word -10286 ;B3=-0.3139 .word 14615 ;B2=0.4460 .w

7、ord -10286 ;B1=-0.3139 .word 3116 ;B0=0.0951 .word -22082 ;A4=-0.6739 .word 31149 ;A3/3=2.8518/3 .word -30484 ;A2/5=-4.6515/5 .word 28383 ;A1/4=3.4647/4 .textstart: SSBX FRCT STM #X,AR1 ;AR1指向X的首地址 RPT #3 ;重復(fù)執(zhí)行下條語(yǔ)句3次 MVPD #table,*AR1+ ;從table首址重復(fù)傳遞4個(gè)數(shù)據(jù) STM #Y,AR1 RPT #3 MVPD #table+4,*AR1+ STM #B,AR

8、1 RPT #4 MVPD #table+8,*AR1+ STM #A,AR1 RPT #3 MVPD #table+13,*AR1+ STM #X+4,AR2 ;為循環(huán)緩沖區(qū)設(shè)初始指針 STM #Y+3,AR4 STM #B+4,AR5 STM #A+3,AR3 STM #5,BK ;將5賦值給BK STM #-1,AR0 ;將-1賦值給AR0STM #1000h,AR6 ;輸出數(shù)據(jù)緩沖區(qū)首址為#1000hSTM #0200h-1,AR7 ;循環(huán)計(jì)算512個(gè)樣本點(diǎn)LOOP: PORTR PA1,*AR2 ;x(n)/2，防止溢出;MVKD in,*AR2 LD *AR2,A ;將AR2的地址

9、加載到A累加器 STL A,-1,*AR2 ;A累加器低位右移1位送往AR2寄存器所指地址 MPY *AR2+0%,*AR5+0%,A ;作B0*X（n）+A=A運(yùn)算，并將地址減一 MAC *AR2+0%,*AR5+0%,A MAC *AR2+0%,*AR5+0%,A MAC *AR2+0%,*AR5+0%,A MAC *AR2,*AR5+0%,A MAC *AR4,*AR3,A ;作A1/4*Y(n-1) MAC *AR4,*AR3,A MAC *AR4,*AR3,A MAC *AR4+0%,*AR3+0%,A ;運(yùn)算四次后將寄存器地址減一 MAC *AR4,*AR3,A ; 作A2/5*Y

10、(n-2) MAC *AR4,*AR3,A MAC *AR4,*AR3,A MAC *AR4,*AR3,A MAC *AR4+0%,*AR3+0%,A ;運(yùn)算五次后將寄存器地址減一 MAC *AR4,*AR3,A ;作 A3/3*Y（n-3） MAC *AR4,*AR3,A MAC *AR4+0%,*AR3+0%,A ;運(yùn)算三次后將寄存器地址減一 MAC *AR4+0%,*AR3+0%,A ;作 A4*Y（n-4） MAR *AR3+0% STH A,*AR4 ;A累加器高位送往AR4寄存器所指地址 ;BD IIR;MVDK *AR4,out PORTW *AR4,PA0 ;輸出y（n）STH A,*AR6+ ;保存y（n） BANZ LOOP,*AR7 ;循環(huán)