归纳频谱分析滤波器设计.pptVIP

归纳：频谱分析和滤波器设计 一。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） 二阶串接模型 二。系统的表述 (离散系统） 状态方程模型 （参数矩阵A,B,C,D） 传递函数模型 （参数数组 a,b） 零极增益模型 （参数数组 z,p,k） 极点留数模型 （参数数组 r,p） 二阶串接模型 三。系统的表述 (离散系统） 状态方程模型（不用） 传递函数模型（参数数组 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相同，不过速度快得多，而且有许多出错提示能力。 Matlab程序 ftdemo 离散频谱与连续频谱的转换 MATLAB只能处理离散量，把离散量取得很密时，就接近于连续量。在时间轴上取得很密，对应的频谱应该取得很窄，离散傅立叶变换的最高频率通常取为2pi，故只能取到0.2*pi以下的部分，才代表连续系统的频谱。 另外它抽去了时间和频率，应当补上。时间应乘以采样周期Ts，即X=Ts*fft(x) fft的频率轴应乘以2π/(N*Ts),单位：弧度/秒 或1 /(N*Ts)=Fs/N,单位：赫兹 作傅立叶反变换ifft时注意事项 1。x=ifft(X)中的X必须取全频段(-π~π]的值。通常给出的X为[0~ π],要扩展。 2。实序列x的X为复数序列，其实部（幅度）偶对称，虚部（相位）奇对称，这可作为扩展的规则。 3。反变换x应取实部x=real(ifft(X)) 4. 如X已含有采样周期T,反变换应除以T. 求连续信号频谱的程序例 x=input(x(1:N)序列); Nx=length(x); K=nextpow2(Nx);N=2^K; n=0:N-1; Fs=input(采样频率Fs= );T=1/Fs; X=T*fft(x); m1=0:N/2-1;D=2*pi/(T*N); subplot(2,1,1),plot(m1*D