资源描述:
《MATLAB70使用详解第16章信号处理工具箱》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、第16章信号处理工具箱MATLAB信号处理工具箱提供函数主要用于处理信号与系统问题,并可对数字或离散的信号进行变换和滤波。工具箱为滤波器设计和谱分析提供了丰富的支持,通过信号处理工具箱的有关函数可以直接设计数字滤波器,也可以建立模拟原型并离散化。在这章里将主要介绍采样函数、信号生成、滤波器的设计和随机信号功率谱分析。16.1MATLAB7.0中的采样函数数字信号处理的对象是在采样时钟的控制之下,通过A/D转换器以一定的采样率对模拟信号进行采样得到的。由采样定理可知,采样率必须大于模拟信号的最高频率
2、的2倍。然而在很多情况下,需要对信号进行不同频率的采样,这就需要对采样信号进行处理。这一节将主要介绍常用的改变采样频率的相关MATLAB函数。16.1.1upfirdn()函数upfirdn()函数用于改变信号的采样率,适用于FIR滤波器。upfirdn()函数的调用格式为:upfirdn(X,H,P,Q)返回信号X是通过上采样、滤波和下采样3个级联系统后的输出结果。上采样系统:通过插入零值,使得采样频率上升为原来的P倍。滤波系统:用H给定的单位冲激响应的滤波器滤波。下采样系统:通过抽取采样点,使
3、得采样率下降Q倍。16.1.2decimat()下采样函数decimat()函数主要用于低通滤波的下采样。decimat()函数的调用格式为Y=decimat(X,R),函数返回的是向量X的重采样序列,其采样率为原来的1/R倍。decimat()使用的滤波器是一个8阶的ChebyshevI型低通滤波器,其截止频率为0.8*(Fs/2)/R,其中Fs为采样频率。16.1.3interp()上采样函数interp()函数主要用于上采样。interp()函数的调用格式为Y=interp(X,R),返回的
4、是以原来采样率R倍的重采样序列。当函数调用格式为Y=interp(X,R,L,ALPHA)时,参数L(缺省值为4)和ALPHA(缺省值为0.5)可以指定。2*L是进行插值时用到的原始数据的长度。理想情况下,L必须小于等于10。当调用格式为[Y,B]=interp(X,R,L,ALPHA)时,返回插入所用滤波器的系数向量B。16.1.4resample()重采样函数resample()函数主要用于改变信号的采样率。resample()函数的调用格式为Y=resample(X,P,Q),返回的是向量X
5、经重采样后的结果,重采样后采样率为原来的P/Q倍。在重采样时,resample()采用的是一个抗混叠低通FIR滤波器对X进行滤波。滤波器使用的是Kaiser窗。当调用格式为Y=resample(X,P,Q,N)时,在重采样时,x(n)的两边都取N点。resample()使用的滤波器长度与N(缺省值为10)成正比;通过增大N值可获得较大的精度,但这会增加计算的时间。当调用格式为Y=resample(X,P,Q,N,BETA)时,BETA为设计用到滤波器使用Kaiser窗的参数,缺省值为5。当调用格式
6、为Y=resample(X,P,Q,B)时,如果B是滤波器系数向量,则用B对X进行滤波。16.2MATLAB中的信号发生函数MATLAB是一种程序化语言,可以处理各种信号变量,在信号处理工具箱中提供了许多信号生成函数,在实际运用中,可以通过这些基本函数生成其它一些需要的“信号”,用以模拟现实中的信号源。这些函数大多数是需要事先输入时间向量的。16.2.1线性调频信号发生器线性调频信号发生器用chirp()函数来实现,其调用格式为Y=chirp(T,R0,T1,F1),产生一个频率随时间线性变化的信
7、号的采样,其时间轴的设置由数组T定义。时刻0的瞬时频率为F0;时刻T1的瞬时频率为F1。缺省情况下,F0=0Hz,T1=1,F1=100Hz。当函数的调用格式为Y=chirp(T,F0,T1,F1,'method')时,method指定改变扫频的方法。可用的方法有linear((线性调频)、quadratic(二次调频)、logarithmic(对数调频)。缺省时为linear。当函数的调用格式为Y=chirp(T,F0,T1,F1,'method',PHI)时,PHI指定信号的初始相位,缺省时P
8、HI的值为0。16.2.2周期函数发生器diric()即周期sinc函数,diric()函数的调用格式为Y=diric(X,N),其返回的是一个大小与X相同的矩阵,其元素为Dirichlet函数。N必须为正整数,该函数将0~2π等间隔地分成N等份。Dirichlet函数的定义是:d(x)=sin(N*x/2)./(N*sin(x/2))16.2.3锯齿波、三角波和矩形波发生器1.sawtooth()锯齿波和三角波发生器sawtooth(T),产生周期为2π,幅值为1的锯齿波,采样时