基于dsp的數(shù)字濾波器設(shè)計與仿真-閱讀頁
2024-12-02 15:33本頁面
【正文】 , (ORDER1) a0 .usect a0, (ORDER1) input .usect input,D_LEN output .usect output, D_LEN .copy 。 Copy a0(coefs) to FIR_COEF(AR3) RPT ORDER1 MVPD coeff,*FIR_COEF+ STM 1, FIR_INDEX STM xn, FIR_DATA 。 Get data from input STM output, DATA_OUT 。 FIRS STH A, *DATA_OUT+ next FIR_END: B FIR_END .end 的程序如下： fs=40000。 %采樣周期 n=0:N1。 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)。%歸一化 xto_ccs=round(32767*xin)%取整 fid=fopen(39。,39。)。1651 1 0 0 0\n39。%輸出文件頭 fprintf(fid,39。,xto_ccs)。 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。 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 }。 %/Hz f12=5500。 fs=32020。 %采樣周期 n=0:N。 x12=*sin(2*pi*f12*n*T)。 x_base=(x11+x12+x13)。 y_f_fft=abs(fft(y_f_out))。%歸一化 xto_ccs=round(32767*xout) fid=fopen(39。,39。)。1651 1 0 0 0\n39。%輸出文件頭 fprintf(fid,39。,xto_ccs)。 ： include include define N 111 define LEN 501 long yn。 int output[LEN]。 int *x。jLEN1。input[j]。 for(i=0。 i++) yn += B[i]*(*x++)。 } while
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1