[理学]第7章 MATLAB在信号处理中的应用.pptVIP

第7章 MATLAB在数字信号处理中的应用 ——本章将介绍数字信号处理的相关知识、傅里叶变换、IIR数字滤波器、FIR数字滤波器的基本理论和MATLAB实现。 （一）典型信号及其表示 1. 单位抽样信号（又称Kronecker函数） 在MATLAB中用zeros函数实现，但需要注意的是，MATLAB只能实现有限长的序列，选择产生N点的单位抽样序列： x=zeros(1, N); x(1)=1; 2.单位阶跃序列 在MATLAB中，可采用ones函数产生N点单位阶跃序列： x=ones(1,N); 3.正弦序列 其中，f是正弦波的频率， 是抽样周期， 为正弦波初相， 为正弦波幅度。在MATLAB中定义变量后，下列MATLAB语句给出正弦波序列： n=0:N-1; x=A*sin(2*pi*f*n*TS+pha); %TS为抽样周期，pha为正弦波初相 4. sinc函数 连续sinc函数定义为： MATLAB函数提供了sinc命令来产生N点采样序列： x=-2*pi:4*pi/N:2*pi; y=sinc(x); 5. 随机序列 MATLAB提供了多种命令产生不同概率密度分布的随机变量序列，其中产生均匀分布与高斯分布随机序列命令： （1）rand：产生均值为0.5，幅度在[0,1]间均匀分布的伪随机数，调用方式为u=rand(N)或u=rand(M,N)。 （2）randn：产生均值为0，方差为1，服从高斯（正态）分布的随机序列，调用方法与rand相同。 除此之外，MATLAB中还包含其它随机序列命令，如chirp、gauspuls、triplus等，读者若用到，可再进一步查找。 （二）离散时间系统 离散时间系统的特点 离散时间系统，可以抽象为一种变换，或者一种映射，即把输入序列 变换为输出序列 ： 其中T表示一种变换。 在数字信号处理中，通常研究的都是线性时不变（LTI）系统（亦可称作线性移不变LSI系统）。 1.线性：指该系统的输入、输出之间满足叠加原理。 2.时不变性：指对给定的输入，系统的输出和输入施加的时间无关。 离散时间系统的表示方法 离散时间系统可以用时域和频域两种方法描述。在时域，LTI系统的输入输出关系可以用差分方程或单位抽样响应 ；在频域还可以采用转移函数、二次分式和零极点增益等形式。 单位抽样响应 频域响应 零极点增益形式 1. 单位抽样响应 单位抽样响应是当输入为单位抽样信号时系统的输出响应，因此利用单位抽样信号来求单位抽样响应。在数字信号处理中，如果一个离散时间系统是用作对输入信号做滤波处理，离散时间系统又可称为数字滤波器。可以利用MATLAB中的filter函数和单位抽样信号来获得离散时间的单位抽样响应。 例：利用 求出系统的单位抽样响应。 2.频率响应  对于离散系统，频域响应由MATLAB函数freqz得到，函数用法为： [h,f]=freqz (b, a, n, ’whole’, fs)：b, a同例11.1由离散系统决定；fs为抽样频率，可自行定义，若fs=1，频率轴给出归一化频率；n是频率轴的分点数，建议取2的整次幂；whole指定计算的频率范围是从0~fs，缺省时是0~fs/2。 例：分区间显示函数值 function y=f(x) if x0 y=0; elseif x1 f=x; elseif x2 f=2-x; else f=0; end 例：利用 求出系统的频域响应。 3.零极点增益形式 零极点增益形式——利用MATLAB中的求根函数由离散系统差分方程得到系统的零极点。 例：根据上例差分方程求系统零极点 b=[0.6 -0.3 0.5]; a=[1 -0.7 -0.2]; zp=roots(b) pp=roots(a) 4.离散系统转换 同一个离散系统，可以用多种表示方式，如转移函数、零极点增益、状态空间等等，MATLAB工具箱中提供了实现不同表示方法之间进行转换的函数，如前一章所述。 （三）傅里叶分析 离散时间信号的Z变换和傅里叶变换 离散傅里叶变换 快速傅里叶变换 1.离散时间信号的Z变换和傅里叶变换 Z变换是对离散序列进行的一种数学变换