DSP试验报告高通滤波器

1、高通滤波器设计报告学号：172030085:徐军一、实验目的：1、了解FIR滤波器的原理及使用方法；2、了解使用Matlab语言设计FIR滤波器的方法；3、了解用DSP来实现FIR滤波器的设计及编程方法；4、熟悉在CCS环境下对FIR滤波器的调试方法。二、实验要求：设计一个FIR高通滤波器，输入信号是频率为25Hz, 50Hz的合成等幅信号, 要求滤去25Hz信号成分，保留50Hz的信号成分。三、实验设计：本实验要求滤去25Hz的信号成分，保留50Hz的信号，根据耐特斯特准侧， 釆样频率需要大于最高频率的2倍，即设计的釆样频率至少为100Hz的高通滤波 器。先在MATLAB中验证设计思路，将得

2、到的滤波器参数调用到DSP程序中，这 样可以简化编译汇编连接的过程。然后通过图形仿真查看滤波前后的波形和频谱 图。四、实验步骤：1、滤波器的MATLAB语言设计2、在MATLAB中确定滤波器的各种参数3、滤波器的程序设计4、在CCS环境下调试程序5、比较滤波前后的效果、观测滤波前后的波形五、实验程序及结果：(1) MATLAB 程序：main.mf 1=50;咎信号频率Hzf2=25;%信号频率Hzfs=1000;%采样频率HzN=200;%采样点数t= (0 :N-l) /fs;乡采样时间ssignall=sin(2*pi*f1*t);signal2=sin(2*pi*f2*t);y=sig

3、nall+signal2;令原始信号时域波形图figure (1);plot(y)；%axis ( 0 100 一25 25);hold on;plot(signail,1r1);legend (被污染的信号*, 理想信号* );辛*原始信号频谱图fy=fftshift(fft(y);f=linspace(-fs/2f fs/2z N);figure (2);plot(f,abs(fy);t itle(*原始信号频谱J;xlabelCf/Hz1);ylabel (1 幅度);axis( 0 100 0 150);辛滤波后的时域波形图figure (3);Hd = high;output=fte

4、r(Hd,y);plot(output);t itleCilg波后的波形J;%*滤波后的频域波形图fy=fftshift(fft(output);f=linspace(-fs/2/fs/2z N);figure (4);plot(f,abs(fy);title (*滤波后信号频谱,);xlabelCf/Hz1);ylabel (1 幅度);axis( 0 100 0 150);high.mfunction Hd = high%HIGH Returns a disCrete-time filter object% MATLAB Code% Generated by MATLAB(R) 8.3 a

5、nd the Signal Processing Toolbox 621% Generated on: 12-Jul-2018 09:57:59% Equiripple Highpass filter designed using the FIRPM function % All frequency values are in HzFs = 1000;% Sampling FrequencyFstop = 25;% Stopband FrequencyFpass = 50;% Passband FrequencyDstop = 0.0001;% Stopband AttenuationDpas

6、s = 0.5;% Passband Rippledens = 20;% Density Factor Nz Foz Aoz W = f irpmord ( Fstop,Fpass/(Fs/2),0 1r Dstop, Dpass);b = firpm(Nz Fof Aof WF dens);Hd = dfiltdffir(b);实验结果：Matlab环境下的滤波前后的时域波形:(2) DSP在CCS下的程序: 源程序：ttinclude ttinclude ttinclude tiinclude ttinclude tiinclude fdacoefs. h #define pi 3.141

7、59 ftdefine N 30 ttdefine length 256 long yn; int inputlength; int outputlength; void main()int m, n;int *x;for(n=0;n=length-l;n+)inputn二0;outputn二0;for(n=0;n=length-l;n+)inputn二50*sin(2*pi*n*25/200)+50*sin(2*pi*n*50/200);for(n=0;n15;while (1);fdacoefes.h 如下:/* Filter Coefficients (C Source) generat

8、ed by the Filter Design and Analysis Tool* Generated by MATLAB(R) 9.0 and the Signal Processing Toolbox 7.2.* Generated on: 27-Jul-2016 11:32:00*/Discrete-Time FIR Filter (real)* Filter Structure* Filter Length* Stable* Linear Phase*/:Direct-Form FIR:43:Yes:Yes (Type 1)/* General type conversion for

9、 MATLAB generated C-code */#include “tmwtypes.h/* Expected path to tmwtypes.h* D:Program FilesMATLABR2016aexternincludetmwtypes.hT/* Warning - Filter coefficients were truncated to fit specified data type.* The resulti ng resp onse may not match gen erated theoretical response * Use the Filter Desig

10、n & Analysis Tool to design accurate* intl6 filter coefficients.*/const int BL = 29;const intl6_T B29 = -106,390,-296,33&135,619,281,-782, -1127z335,2277,1564,330&969&20109,969&330&1564,2277,335,-1127,-782,281,619,135,33&296,390,-106;（3）滤波器的仿真测试新建工程并编译成功后会在“工程所在目录/debeg文件夹下产生sheji2.out 文件，在CCS软件的Run

11、-*Load-*Load Program里打开这个.out文件，单击OK。 查看滤波效果：1）选择菜单栏上的ToolsGraph-Dual Time选项，弹出对话框，按照下图所示设 置时域波形参数。2）选择菜单栏上的ToolsGraph FFT Magnitude选项，弹出对话框，按照下图 所示设置频域波形参数。Graph Properties. 亠 jt 丄PropertyValue Dat PropertiesAcquisition Buffer Size256Dsp Dat3 Type16 bit signed integerIndex IncrementQ_ValucSampling Rate Hz101Signal TypeRealStan Addressinput Display Properties Axis Display Data Plot Style0 trueLineFrequency Display Uni Hz