数字信号插值详解.ppt

Matlab示例 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 (\times\pi rad/sample)) ylabel(Magnitude) disp([{ sprintf(%.4f,,h1(1:N_filter-1)) … sprintf(%.4f,h1(N_filter)) };]) DSP示例 #include Filter.h #include fract2float_conv.h static const float SinX_I4_16_float[16]= xxx,xxx,…,xxx}; const short *InPtr; float Ibuffer[1200]=0; //输入 fract16 Ibuffer_fr16[1200]=0; //输入 short Len=300; //插值前信号点数 fract16 SinX_I4_16_fr16[16]=0; //插值系数 short Obuffer[1215]=0; //输出 fract16 Obuffer_fr16[1215]=0; //输出 InPtr = (short *) gOriginBuffer; DSP示例 //对输入信号归一化，插入零值，完成数据类型的转换 //完成插值系数的类型转换 //对Ibuffer_fr16低通滤波 //void convolve_fr16(const fract16 input_x[],int length_x, const fract16 input_y[],int length_y, fract16 output[]); //对Obuffer去归一化类型转换 //信号插值结果显示 显示结果示例 * * * * * * * * * 具有若干功能强大的应用工具箱。 MATLAB主工具箱 符号数学工具箱 SIMULINK仿真工具箱 控制系统工具箱 信号处理工具箱 图象处理工具箱 通讯工具箱 系统辨识工具箱 神经元网络工具箱 金融工具箱… MATLAB特点 Signal Processing Toolbox Filter Analysis FIR Filter Design FIR数字滤波器 — 窗函数设计法 滤波器系数的确定 Window-based finite impulse response filter design b=fir1(n,wn,’ftype’,window) b = fir1(N,Wn); 若wn为标量，则设计低通滤波器， 若wn是1×2的向量，则用来设计带通滤波器，若wn是1×L的向量，则可用来设计L带滤波器 b = fir1(N,Wn,‘high’); b = fir1(N,Wn, ‘stop’); 在上述格式中，若不指定窗函数的类型，fir1自动选择Hamming窗 FIR数字滤波器 — 窗函数设计法 滤波器系数的确定 Window-based finite impulse response filter design b=fir1(n,wn,’ftype’,window) b = fir1(15,1/I); b = I*b; FIR数字滤波器 — 频率抽样设计法 滤波器系数的确定 Frequency sampling-based finite impulse response filter design b=fir2(n, f, m) f ---- 归一化频率向量，其值在0～1之间， m ---- 与f相对应的期望幅频响应。 FIR数字滤波器 — 频率抽样设计法 滤波器系数的确定 Frequency sampling-based finite impulse response filter design b=fir2(15,[0 1/I 1/I 1],[I I 0 0]) 频率响应 freqz(num,den,n) freqz(num,den,n,Fs) 其中： num-----系统传递函数分子系数组