频谱分析和滤波器相关设计.pptVIP

滤波器浏览器的概貌: 使用滤波器浏览器可以分析所选滤波器下列的响应特征 幅度响应,相位响应 脉冲响应,阶跃响应 群延迟 极点和零点位置 滤波器浏览器能在任何时候显示出选定滤波器的多达六种不同响应特征曲线. 滤波器浏览器可提供以下功能 缩放 测量滤波器响应 分析和修改滤波器响应 修改显示方法：例如频率范围或者重幅度单位等参数 使信号通过滤波器得到输出 频谱浏览器的概貌:谱分析 估计和分析一信号的功率谱密度(PSD). 用PSD估计值得知信号的频率内容. 使用谱浏览器可以 分析和比较谱的密度曲线； 评价不同谱分析的方法; 比较窗函数,FFT长度,采样频率等的影响; 声音信号的采集 用话筒及windows的录音机录下声音文件： 后缀.wav,存储采样频率默认22500Hz，可选 用MATLAB中的wavread命令读入.wav文件，参数为样本长度； 再用save命令将它改为.mat数据格式存储； 以后即可在sptool中直接调用。 * 此处是以Ingle书为蓝本时才有sigadd, sigmult, sigshift(移位）等函数 归纳：频谱分析和滤波器设计 一。MATLAB表述的信号和系统 一。离散信号的表述与运算 1。信号的表述： 用 n 及 x 两个数组，如书 impseq,stepseq,... 2。信号的初等运算；相加，相乘，移位等,...(书 sigadd, sigmult, sigshift） 3。信号的卷积、相关、抽选等；conv, xcorr, xcov, decimate, ... 二。系统的表述 (连续系统） 状态方程模型 （参数矩阵A,B,C,D） 传递函数模型 （参数数组 a,b） 零极增益模型 （参数数组 z,p,k） 极点留数模型 （参数数组 r,p） 二阶串接模型 （参数矩阵 [b0k,b1k,b2k,1,a1k,a2k]） 二。系统的表述 (离散系统） 状态方程模型 （参数矩阵A,B,C,D） 传递函数模型 （参数数组 a,b） 零极增益模型 （参数数组 z,p,k） 极点留数模型 （参数数组 r,p） 二阶串接模型 （参数矩阵 [b0k,b1k,b2k,1,a1k,a2k]） 三。系统的表述 (离散系统） 状态方程模型（不用） 传递函数模型（参数数组 ad,bd） ad =1 FIR滤波器（系统） ad ~=1 IIR滤波器（系统） 零极增益模型（参数数组 z,p,k） 极点留数模型（参数数组 r,p,C） 二阶串接模型 连续系统和离散系统之间的转换函数模拟/数字滤波器转换 1。双线性变换函数 bilinear [NUMd,DENd] = bilinear(NUM,DEN,Fs) H(z) = H(s) | | s = 2*Fs*(z-1)/(z+1) 2。脉冲响应不变法impinvar [BZ,AZ] = impinvar(B,A,Fs) 信号处理工具箱提供的线性模型转换函数 用MATLAB与信号流图求系统 信号流图方程列写 设Ki个输入节点ui K个中间和输出节点xj 方程为 写成矩阵形式： X= [x1;x2;…xK]为状态变量列向量； Q为K乘K维的状态传送矩阵； U = [u1;u2;…uL]为输入信号列向量 P为K乘L维的为信号输入矩阵 四。信号通过线性系统（滤波器实现） 信号与系统脉冲过渡函数的卷积 y=conv(u,h),... 差分方程解法 y=filter(bd,ad,u) 滤波器实行中的其他函数: fftfilt; latcfilt，sosfilt…等 Z-变换及反变换 部分分式展开 长除 二．频谱分析（第三、五章） 傅立叶变换（包括连续和离散形式 ） 正变换 反变换  离散时间傅立叶正反变换 离散傅立叶正反变换 连续和离散傅立叶变换的对比 在MATLAB中计算傅立叶变换 在MATLAB中的计算相当方便，只要下列两行八个语句：xn=input(x = ); N= (length(xn)); n = [0:1:N-1]; k = nWN = exp(-j*2*pi/N); nk = n*k; WNnk = WN .^ nk; Xk = xn * WNnk 在《数字信号处理及其MATLAB实现》书pp119-120中，定义它为dft函数，相应也定义了idft函数，即X=dft(x,N)x=idft(X,N) 但其运算速度低。 在MATLAB中设有基本函数fft，它是快速傅立叶变换的缩写，调用形式为：X=fft(x,N) x=idft(X,N) N取为2的幂次时，速度最快。（令N=2^nextpow2(N);) 它的功能和dft相同，不过速度快得多，而且有许多出错提示