基于dsp的數(shù)字濾波器設(shè)計(jì)與仿真-資料下載頁(yè)
【導(dǎo)讀】得DSP應(yīng)用技術(shù)的實(shí)際訓(xùn)練，掌握設(shè)計(jì)較復(fù)雜DSP系統(tǒng)的基本方法。理解，進(jìn)一步鞏固匯編語(yǔ)言講法規(guī)則。學(xué)會(huì)編制結(jié)構(gòu)清晰、風(fēng)格良好、數(shù)據(jù)結(jié)構(gòu)。適當(dāng)?shù)膮R編語(yǔ)言程序，從而具備解決綜合性實(shí)際問(wèn)題的能力。數(shù)字濾波器在語(yǔ)音信號(hào)處理、信。都有廣泛的應(yīng)用，數(shù)字濾波器也是使用最為廣泛的信號(hào)處理算法之一。在本設(shè)計(jì)中，使用MATLAB模擬產(chǎn)生合成信號(hào)，然后利用CCS進(jìn)行濾波。定模擬信號(hào)的采樣頻率為32020Hz，。其中N-1是FIR的濾波器的階數(shù)，為延時(shí)結(jié)，為端口信號(hào)函數(shù)。數(shù)、采樣頻率等。指定完設(shè)計(jì)參數(shù)后單擊按鈕DesignFilter，生成濾波器系數(shù)。指定文件名和路徑。行加coeff.word,在每行的前面加上.word，比且把每行的最后的逗號(hào)去掉。編譯匯編程序，如果有錯(cuò)誤，按錯(cuò)誤進(jìn)行修改；沒(méi)錯(cuò)誤，則往下執(zhí)行。加載初始化DATA數(shù)據(jù)。運(yùn)行程序，查看輸入輸出波形，修改相應(yīng)參。在本實(shí)驗(yàn)中使用MATLAB模擬產(chǎn)生信號(hào)，觀察濾波前的時(shí)域波形和頻域波形。并使用CCS的頻譜分析功能，查看DSP的濾波效果。
【正文】 頭 fprintf(fid,39。%d\n39。,xto_ccs)。%輸出 fclose(fid)。 C 程序清單 文件內(nèi)容： /* * Filter Coefficients (C Source) generated by the Filter Design and Analysis Tool * * Generated by MATLAB(R) and the Signal Processing Toolbox . * * Generated on: 11Dec2020 10:27:37 * */ /* * DiscreteTime FIR Filter (real) * * Filter Structure : DirectForm FIR * Filter Order : 110 * Stable : Yes * Linear Phase : Yes (Type 1) */ /* General type conversion for MATLAB generated Ccode */ include C:\MATLAB701\extern\include\ /* * Expected path to * C:\MATLAB701\extern\include\ */ /* * Warning Filter coefficients were truncated to fit specified data type. * The resulting response may not match generated theoretical response. * Use the Filter Design amp。 Analysis Tool to design accurate * int16 filter coefficients. */ const int BL = 111。 const int16_T B[111] = { 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, 641, 1953, 5169, 7819, 23925, 7819, 5169, 1953, 641, 1844, 1602, 502, 614, 1105, 823, 90, 569, 765, 452, 100, 511, 544, 227, 196, 443, 379, 79, 239, 370, 248, 21, 250, 295, 145, 85, 237, 223, 65, 121, 212, 159, 6, 138, 181, 104, 36, 145, 155, 63, 69, 158, 145, 25, 138, 234, 140, 216, 605, 654, 100 }。 內(nèi)容： f11=2500。 %/Hz f12=5500。 %/Hz f13=9000。 fs=32020。 %/采樣 Hz N=500 %數(shù)據(jù)個(gè)數(shù) T=1/fs。 %采樣周期 n=0:N。 x11=sin(2*pi*f11*n*T)。 x12=*sin(2*pi*f12*n*T)。 x13=*sin(2*pi*f13*n*T)。 x_base=(x11+x12+x13)。 %待濾波信號(hào)波形 figure(1) plot(x_base) %待濾波信號(hào)頻譜 figure(2) yff=abs(fft(x_base)) df=n*(fs/N) plot(df,yff) y_f_out=filter(B,1,x_base)。 y_f_fft=abs(fft(y_f_out))。 figure(3) plot(y_f_out) figure(4) plot(df,y_f_fft) xout=x_base/max(x_base)。%歸一化 xto_ccs=round(32767*xout) fid=fopen(39。39。,39。w39。)。%打開(kāi)文件 fprintf(fid,39。1651 1 0 0 0\n39。)。%輸出文件頭 fprintf(fid,39。%d\n39。,xto_ccs)。%輸出 fclose(fid)。 ： include include define N 111 define LEN 501 long yn。 int input[LEN]。 int output[LEN]。 void main() { int i,j。 int *x。 for(j=0。jLEN1。j++) { x=amp。input[j]。 yn = 0。 for(i=0。 iN1。 i++) yn += B[i]*(*x++)。 output[j]=yn15。 } while(1)。 }
鄂ICP備17016276號(hào)-1