数字信号插值

1、 时域测试技术综合实验数字示波器中的信号插值数字示波器中的信号插值电子测试技术及仪器教育部工程研究中心电子测试技术及仪器教育部工程研究中心潘卉青潘卉青信号插值信号插值v为什么插值： 显示所需采样率大于显示所需采样率大于ADC最高采样率，显示所最高采样率，显示所需的两个点的时间间隔小于需的两个点的时间间隔小于ADC输出的采样点的时输出的采样点的时间间隔间间隔时，需要插值。时基时基3实验内容实验内容v实验内容实验内容 使用Matlab设计一个插值倍数I=4倍的数字插值滤波器； 通过DSP实现对采样数据的4倍插值； 在Visual DSP+中完成调试，并送屏幕显示，确认运算结果的正确性。实验目的实验

2、目的v实验目的实验目的 掌握基于DSP的数字信号处理基本流程； 掌握信号插值数字滤波器的基本设计方法； 掌握数字信号插值在DSP中的实现流程。v实验设备与工具实验设备与工具 硬件：PG1000嵌入式开发平台，PC 机，信号发生器。 软件：VisualDSP+ Environment、Matlab5主要内容主要内容v信号插值的基本原理信号插值的基本原理v数字滤波器v基于Matlab的信号插值滤波器设计vDSP中数字滤波器的实现信号插值信号插值v为什么插值： 显示所需采样率大于ADC最高采样率，显示所需的两个点的时间间隔小于ADC输出采样点的时间间隔时，需要插值。v常用插值方式信号插值信号插值线性

3、内插线性内插v线性内插 已知直线上两点的坐标(x1, y1)和(x2, y2)，在由这两点决定的直线上插入若干个给定点的x的值y 由直线的两点式公式： 得线性插值公式： 121121yyyyxxxx211121()yyyyxxxxv正弦内插多抽样频率下信号的处理称为“多抽样率信号处理多抽样率信号处理” Multirate Signal Processing信号插值信号插值正弦内插正弦内插信号插值信号插值正弦内插正弦内插v正弦内插v将x(n)的采样频率增加I倍(整数倍)相邻两点之间等间隔插入 I-1个零点，I 称为插值因子v内插后的序列信号插值信号插值正弦内插正弦内插( )0, 2( )0enx

4、nIIx nIelseI( )ex n( )x nsfsIf信号插值信号插值正弦内插正弦内插v采样频率转换( / )0, 2( )0Ix n InIIx nelse()( )( / )( )jj nIInj nj kInnXex n ex n I ex k e()()jj IIXeX e信号抽取前后频域的关系X (e j )周期：- ; XI (e j)周期：- / I / I信号插值信号插值正弦内插正弦内插v低通滤波器1|()10|jjIeIH eI通带内幅度通带内幅度=I截止频率截止频率=1/I主要内容主要内容v信号插值的基本原理v数字滤波器数字滤波器v基于Matlab的信号插值滤波器设计

5、vDSP中数字滤波器的实现Digital Signal Processingv数字信号处理数字信号处理 信号用数字或符号表示的序列，通过计算机或通用(专用)信号处理设备，用数字的数值计算方法处理(例如等)，以达到提取有用信息便于应用的目的。v模拟滤波器：模拟滤波器： 用电阻、电容、电感及有源器件等构成滤波器对信号进行滤波v数字滤波器：数字滤波器： 通过对输入信号的进行来实现滤波v优点优点： 高精度、稳定、体积小、重量轻、灵活，不要求阻抗匹配，可实现特殊滤波功能数字滤波器数字滤波器( )x n1z0b1a( )y n1616DSP在数字存储示波器中的应用在数字存储示波器中的应用v信号滤波（减噪）

6、v频域补偿 (幅度平坦化、带宽增强、相位校正等)v信号重建 (非均匀校正等)数字滤波器的分类数字滤波器的分类vIIR数字滤波器：数字滤波器： 可以利用模拟滤波器设计 但相位非线性vFIR数字滤波器：数字滤波器： 可以严格线性相位，又可任意幅度特性 因果稳定系统 可用FFT计算 但阶次比IIR滤波器要高得多18信号插值信号插值低通滤波低通滤波v低通滤波器 窗函数设计法 频率抽样设计法 1|()10|jjIeIH eI通带内幅度通带内幅度=I截止频率截止频率=1/I通带内幅度通带内幅度=I截止频率截止频率=1/I19主要内容主要内容v数字滤波器v信号插值v基于基于Matlab的信号插值滤波器设计的

7、信号插值滤波器设计vDSP中数字滤波器的实现MATLAB启动与运行启动与运行菜单栏菜单栏工具栏工具栏文件编辑窗口文件编辑窗口工作空间工作空间历史命令历史命令窗口窗口Start Start 菜单菜单命令窗口命令窗口命令提示符命令提示符v具有很好的帮助功能。具有很好的帮助功能。 提供十分详细的帮助文件提供十分详细的帮助文件（PDF、HTML、demo文件）。 联机查询指令：联机查询指令： help指令指令（例：help elfun，help exp，help simulink） lookfor关键词关键词（例： lookfor fourier） doc指令指令（例：doc conv） 大量的网上资

8、源和用户社区。大量的网上资源和用户社区。MATLAB特点特点v具有若干功能强大的应用工具箱具有若干功能强大的应用工具箱。 MATLAB主工具箱主工具箱 符号数学工具箱符号数学工具箱 SIMULINK仿真工具箱仿真工具箱 控制系统工具箱控制系统工具箱 图象处理工具箱图象处理工具箱 通讯工具箱通讯工具箱 系统辨识工具箱系统辨识工具箱 神经元网络工具箱神经元网络工具箱 金融工具箱金融工具箱MATLAB特点特点Signal Processing ToolboxFilter AnalysisFIR Filter DesignFIR数字滤波器数字滤波器 窗函数设计法窗函数设计法v滤波器系数的确定滤波器系数

9、的确定 Window-based finite impulse response filter designb=fir1(n,wn,ftype,window) b = fir1(N,Wn); 若wn为标量，则设计低通滤波器， 若wn是12的向量，则用来设计带通滤波器，若wn是1L的向量，则可用来设计L带滤波器 b = fir1(N,Wn,high); b = fir1(N,Wn, stop); 在上述格式中，若不指定窗函数的类型，fir1自动选择Hamming窗25FIR数字滤波器数字滤波器 窗函数设计法窗函数设计法v滤波器系数的确定滤波器系数的确定 Window-based finite i

10、mpulse response filter designb=fir1(n,wn,ftype,window) b = fir1(15,1/I); b = I*b;FIR数字滤波器数字滤波器 频率抽样设计法频率抽样设计法v滤波器系数的确定滤波器系数的确定 Frequency sampling-based finite impulse response filter design b=fir2(n, f, m) f - 归一化频率向量，其值在01之间， m - 与f相对应的期望幅频响应。27FIR数字滤波器数字滤波器 频率抽样设计法频率抽样设计法v滤波器系数的确定滤波器系数的确定 Frequenc

11、y sampling-based finite impulse response filter design b=fir2(15,0 1/I 1/I 1,I I 0 0)v频率响应频率响应 freqz(num,den,n) freqz(num,den,n,Fs) 其中： num-系统传递函数分子系数组成的行向量 den-系统传递函数分母系数组成的行向量 n-是频率响应的点数，最好为2的幂缺省值512 Fs-采样频率FIR数字滤波器数字滤波器 Filter Analysisb = fir1(15,1/I);H,w = freqz(b,1,512);plot(w/pi,abs(H1);Filter

12、滤波滤波v因果因果LSI系统的的零状态响应，在系统的的零状态响应，在matlab中中可以用函数可以用函数 y=filter(num,den,x) 其中： num-系统传递函数分子系数组成的行向量 den-系统传递函数分母系数组成的行向量 x-输入的离散序列 y-输出的离散序列，y的长度与x的长度一样MATLAB - FDAToolvMATLAB中提供的滤波器设计辅助设计软件（在命令窗口中键入“fdatool”即可启动）定点化Simulink产生HDL幅度响应相位响应MATLAB - LowPassFilterv低通滤波器的设计规范主要内容主要内容v数字滤波器v信号插值v基于Matlab的信号插

13、值滤波器设计vDSP中数字滤波的实现中数字滤波的实现基本数据采集系统基本数据采集系统v本实验原理框图 ADC: AD9481, 250MSPS FPGA：Xilinx公司的Spartan-3A系列的XC3S400A-4，开发工具：ISE DSP: ADI公司的Blackfin系列BF531 开发工具：VisualDSP+数字滤波的实现数字滤波的实现卷积卷积v设两序列x(n)、 h(n)，则其卷积和定义为： 群延时 多相( )( ) ()( )( )my nx m h nmx nh n数字滤波的实现数字滤波的实现 convolve(卷积卷积)convolveconvolutionSynopsis

14、#include void convolve_fr16(const fract16 input_x,int length_x,const fract16 input_y,int length_y,fract16 output);Fract16 & FloatvFLOAT_TO_FR16Function:Convert a floating-point value to a fract16Synopsis: #include fract16 float_to_fr16(float x);vFR16_TO_FLOATFunction: Convert a fract16 to a floa

15、ting-point valueSynopsis: #include float fr16_to_float(fract16 x);基基本本流流程程图图38Matlab示例示例I = 4;N_filter = 16;h1 = fir1(N_filter-1,1/I,hamming(N_filter);h1 = I*h1;sum(h1)H1,w1 = freqz(h1,1,512);figure(1);plot(w1/pi,abs(H1);xlabel(Normalized Frequency (timespi rad/sample)ylabel(Magnitude)disp( sprintf(

16、%.4f,h1(1:N_filter-1) sprintf(%.4f,h1(N_filter) ;)39DSP示例示例#include #include static const float SinX_I4_16_float16=xxx,xxx,xxx;const short *InPtr; float Ibuffer1200=0; /输入fract16 Ibuffer_fr161200=0; /输入short Len=300; /插值前信号点数 fract16 SinX_I4_16_fr1616=0;/插值系数short Obuffer1215=0; /输出fract16 Obuffer_fr161215=0; /输出InPtr = (short *) gOriginBuffer;40DSP示例示例/对输入信号归一化，插入零值，完成数据类型的转换/完成插值系数的