DSP实验报告高通滤波器

高通滤波器设计报告学生编号：172030085姓名：徐军首先，实验的目的：1.了解FIR滤波器的原理和用法；2.理解用Matlab语言设计FIR滤波器的方法；3.了解基于DSP的FIR滤波器的设计和编程方法。4.熟悉CCS环境下FIR滤波器的调试方法。二、实验要求：设计了一种FIR高通滤波器。输入信号是频率为25Hz和50Hz的复合恒定幅度信号。要求滤除25Hz信号分量并保留50Hz信号分量。三、实验设计：这个实验需要滤除25Hz的信号成分并保留50Hz的信号。根据Nexter的准侧，采样频率需要大于最高频率的2倍，即设计了采样频率至少为100赫兹的高通滤波器。首先，在MATLAB中验证了设计思想，并将得到的滤波器参数传递给DSP程序，简化了编译和组装连接的过程。然后通过图形仿真，检查滤波前后的波形和频谱图。四、实验步骤：1.滤波器的MATLAB语言设计2.在MATLAB中确定滤波器的各种参数3.滤波器的程序设计4.在CCS环境下调试程序5.比较滤波前后的效果，观察滤波前后的波形。V.实验程序和结果：(1)MATLAB程序：main.mf1=50%信号频率赫兹f2=25%信号频率赫兹fs=1000%采样频率赫兹N=200。%采样点t=(0:N-1)/fs。%采样时间s信号1=sin(2 * pi * f1 * t)；信号2=sin(2 * pi * F2 * t)；y=信号1信号2；%原始信号的时域波形图(1)；图(y)；% %轴(0 100-2.5 2.5)；坚持住。绘图(信号1，r)；图例(污染信号、理想信号)；%原始信号谱图fy=FFT shift(FFT(y)；f=linspace(-fs/2，fs/2，N)；图(2)；曲线图(f，ABS(fy)；标题(原始信号频谱)；xlabel(频率/赫兹)；伊拉贝尔(振幅)；轴线(0 100 0 150 );%滤波时域波形图(3)；Hd=高电平。输出=滤波器(Hd，y)。绘图(输出)；标题(过滤波形)；%滤波频域波形fy=fftshift(fft(输出)；f=linspace(-fs/2，fs/2，N)；图(4)；曲线图(f，ABS(fy)；标题(滤波信号频谱)；xlabel(频率/赫兹)；伊拉贝尔(振幅)；轴线(0 100 0 150 );高m功能Hd=高%HIGH返回离散时间过滤器对象。% MATLAB代码%由MATLAB(R) 8.3和信号处理工具箱6.21生成。2018年7月12日：生成的百分比09:57:59使用FIRPM功能设计的%等波纹高通滤波器。%所有频率值都以赫兹为单位。Fs=1000。%采样频率Fstop=25。%阻带频率Fpass=50。%通带频率Dstop=0.0001。%阻带衰减Dpass=0.057501127785。%通带纹波密度=20；%密度因子N，Fo，Ao，W=firp word(Fstop，Fpass/(Fs/2)，0 1，Dstop，Dpass)；b=第一(N，Fo，Ao，W， dens )；HD=d filt . d fir(b)。实验结果：Matlab环境下滤波前后的时域波形；滤波前的时域信号(左)，滤波后的时域信号(右)Matlab环境下滤波前后的频谱图；滤波前的频域信号(左)和滤波后的频域信号(右)(CCS下的数字信号处理器程序：源程序：#包括#包括#包括#包括#包括#包括“fdacoefs.h”#定义pi 3.14159#定义N 30#定义长度256龙恩；int输入长度；int输出长度；void main()int m，n；int * x；对于(n=0；n=长度-1；n)输入n=0；输出n=0；对于(n=0；n=长度-1；n)输入n=50 * sin(2 * pi * n * 25/200)50 * sin(2 * pi * n * 50/200)；对于(n=0；n=长度-1；n)x=输入n；yn=0；对于(m=0；m=N-1；m)yn=Bm*(* x)* x；输出n=yn15；而(1)；Fdacoefes.h如下：/*滤波器设计和分析工具生成的滤波器系数(C源)*由MATLAB(R) 9.0和信号处理工具箱7.2生成。*生成于2016年7月27日： 11:32:00*/*离散时间FIR滤波器(实数)* -*滤波器结构：直接形式FIR*过滤器长度： 43*稳定：是*线性相位：是(类型1)*/* MATLAB生成的C代码的通用类型转换*/#包括“tmwtypes.h”/*预期的类型路径* D: 程序文件 MATLAB r 2016 a extern include tmw type . h*/*警告-过滤器系数被截断以适合指定的数据类型。*产生的响应可能与产生的理论响应不匹配。*使用过滤器设计分析工具进行精确设计* int16滤波器系数。*/常量整数BL=2929号国际机场16号线=-106，390，-296，-338，135，619，281，-782，-1127，335，2277，1564，-3308，-9698，20109，-9698，-3308，1564，2277，335，-1127，-782，281，619，135，-338，-296，390，-106；(3)过滤器模拟试验在新项目被成功构建和编译后，sheji2.out文件将在文件夹“项目目录/debeg”下生成。在CCS软件的运行加载加载程序中打开此. out文件，然后单击确定。要查看过滤效果：1)选择菜单栏上的工具图形双时间选项，弹出对话框，设置时域波