[小学教育]数字信号处理MATLAB上机.pptVIP

Matlab 的特点与主要功能: Matlab 的特点与主要功能 Matlab 语言规则 1. 直接输入向量 x1=[1 2 4], x2=[1,2,1], x3=x1 3.生成线性等分向量 指令x=linspace(a,b,n) 在[a,b]区间产生 n 个等分点(包括端点) 例如：x=linspace(0,1,5) 结果 x = 0 0.2500 0.5000 0.7500 1.0000 矩阵： A=[1 4 3 2;2 1 4 3;3 2 1 4;4 3 2 1] B=[1;1;1;1] C=A*B 线性卷积计算循环卷积 xn=[1,2,3,4]; hn=[1,1,1]; yln=conv(xn,hn); ny1=(0:1:length(yln)-1); subplot(3,1,1); stem(ny1,yln);ylabel(‘线性卷积); axis([0 10 0 10]) xn=[1,2,3,4,0,0]; hn=[1,1,1,0,0,0]; ycn=conv(xn,hn); ny2=(0:1:length(ycn)-1); subplot(3,1,3);stem(ny2,ycn);ylabel(‘6 点循环卷积’)； axis([0 10 0 10]) 实验内容： 调试参考资料所给程序； 用matlab实现教材P109页题24（1）及题25 用FFT对信号作频谱分析 FFT只是DFT的一种快速算法，所以FFT的运算结果必然满足DFT的基本性质； y=fft(x)； y=fft(x,n); y=ifft(x,n) abs( ); ones(1,4); zeros(1行,4列); DFT-FFT DFT-FFT IIR数字滤波器设计 设数字滤波器系统函数为 模拟滤波器的系统函数为 函数butter和cheby1可以确定Butterworth和Chebyshev I型滤波器的系统函数。 例：设计butterworth低通滤波器 例7.8 程序运行结果 脉冲响应不变法 脉冲响应不变法: 用代换Ha(s)中的(s-sk)即可得到H(z)，从而将模拟滤波器转换为数字滤波器格式。 可用函数impinvar实现，调用格式为 [bz, az]= impinvar(b, a, fs) 其中，fs为取样频率。 例：模拟滤波器转换数字滤波器 程序段为 [b,a]= butter(4,.3,’s’); [bz,az]= impinvar(b,a,10); 程序运行结果为 bz=1.0e-006 * -0.0000 0.1324 0.5192 0.1273 0 az=1.0000 -3.9216 5.7679 -3.7709 0.9246 双线性变换法 双线性变换法: 用代换Ha(s)中的s即可得到H(z)，从而将模拟滤波器转换为数字滤波器格式。 可用函数bilinear实现，调用格式为 [zd, pd, kd]= bilinear(z, p, k, fs) 其中，z，p，k和zd，pd，kd分别为s域和z域系统函数的零点、极点和增益。 窗函数法设计FIR滤波器 窗函数法是通过对理想滤波器的单位取样响应加窗来逼近理想滤波器的。 函数fir1用于设计标准的低通、带通、高通和带阻滤波器。 函数fir1的调用格式 函数fir1的调用格式为 b= fir1(n, Wc, ftype, Windows) 其中，n为滤波器阶数，Wc为截止频率 ftype决定滤波器类型，ftype= high，设计高通FIR滤波器，ftype= stop，设计带阻FIR滤波器。 Windows指定窗函数类型，默认为Hamming窗；可选Hanning、Hamming、Blackman、triangle、bartlett和boxcar窗，每种窗都可以由Matlab的相应函数生成。 例：设计butterworth低通滤波器 例8.3 程序运行结果 FIR滤波器的优化设计 前面介绍了FIR滤波器的优化设计方法，通过迭代的方法求解FIR滤波器，过程十分复杂。 在Matlab中，可以调用函数remez实现滤波器的设计。 b= remez(n, f, m) 函数remez采用Parks-McClellan算法设计线性相位FIR滤波器，n为滤波器阶数，其幅频特性由f和m指定。 例：设计带通滤波器 程序段为 f= [0 0.3 0.4 0.6 0.7 1]; m= [0 0 1 1