Matlab数字信号处理基础..ppt

第2章 Matlab数字信号处理基础 主 要 内 容 小结 ? Tripuls函数：调用方式如下 x=tripuls(t)：产生非周期单位高度的三角波，三角波的中心位置在t=0处； x=tripuls(t,width,s)：产生倾斜度为s的三角波。 x=tripuls(t,width)：产生width宽度的三角波； 【例2-8】用tripuls函数产生长度为1s，宽度为0.6s，倾斜度分别为0和0.9的非周期三角波。 ? Pulstran 函数：调用方式如下 x=pulstran(t,d,’func’)：其中参数func取值为gauspuls(高斯调制正弦信号)；rectpuls(非周期方波)；tripuls（非周期三角波）。该函数产生以d为采样间隔的func指定形状的冲激串； x=pulstran(t,d,’func’,p1,p2)：将参数p1和p2传递给func函数； x= pulstran(t,d,p,Fs)：向量p表示原始序列，Fs为采样率。对原始序列多次延迟相加得到输出序列。 【例2-9】用pulstran函数产生三角波冲激串。 ? chirp函数 chirp(t,f0,t1,f1)：产生线性调频余弦信号。f0和f1分别是0时刻和t1时刻的瞬时频率。 chirp(t,f0,t1,f1,method)：参数method指定不同的扫频方式,取值方式三种包括: Linear线性的 quadratic二次的 logarithmic对数的 【例2-10】用chirp函数产生二次扫频信号，绘出时域波形和时频图。 2、傅里叶变换函数 常用的FFT及反变换函数如下表所示。 计算离散傅立叶反变换 ifft 调整fft函数的输出顺序，将零频位置移到频谱的中心 fftshift 计算快速离散傅立叶变换 fft 说明 函数 ? fft函数：调用方式如下 y=fft(x)：计算信号x的快速傅立叶变换y。当x的长度为2的幂时，用基2算法，否则采用较慢的分裂基算法。 y=fft(x,n)：计算n点FFT。当length(x)n时，截断x，否则补零。 【例2-11】产生一个正弦信号频率为60Hz，并用fft函数计算并绘出其幅度谱。 ? fftshift函数：调用方式如下 y=fftshift(x)：如果x为向量，fftshift(x)直接将x的左右两部分交换；如果x为矩阵（多通道信号），将x的左上、右下和右上、左下四个部分两两交换。 【例2-12】产生一个正弦信号频率为60Hz,采样率为1000Hz,用fftshift将其零频位置搬到频谱中心。 例2-12 60Hz正弦信号的幅度谱及搬移后的频谱 ? ifft函数：调用方式如下 y=ifft(x)：计算信号x的傅立叶反变换。 【例2-13】计算方波信号的傅立叶反变换。 例2-13 矩形信号及其IFFT变换 y=ifft(x,n)：计算n点IFFT。如果length(x)n，以n为长度截短x，否则补零。 3、滤波器分析与实现函数 直接Ⅱ型滤波器 filter 离散系统的零极点图 zplane 求相角 angle 数字滤波器的冲激响应 impz 求幅值 Abs 求卷积 conv 函数功能 函数名 函数功能 函数名 常用滤波器分析与实现函数 ? conv：调用方式如下 【例2-14】用conv求两个向量的卷积。 ? impz：调用方式如下 [h,t]=impz(b,a,n)：返回n点冲激响应。 c=conv(a,b)：返回向量a、b的卷积c。 [h,t]=impz(b,a)：b、a分别为系统传递函数的分子和分母的系数向量。返回系统(b,a)的冲激响应h和相应的时间轴向量t。 [h,t]=impz(b,a，n,Fs)：指定冲激响应采样点间隔1/Fs。Fs为相对频率，缺省值为1。 【例2-15】计算线性系统(b,a)的冲激响应： 图2-15 系统的冲激响应 zplane：调用方式如下 zplane(z,p)：绘制系统零极点图，“o”表示零点，“x” 表示极点。z，p分别为零点和极点向量。 zplane(b,a)：b、a分别为系统传递函数的分子和分母系 数向量。 * ? MATLAB7.0语言特点及工作环境 ? M文件介绍 ? MATLAB基本控制语句 ? 数字信号处理常用MATLAB函数 2.1 MATLAB 7.0语言特点 1)开发环境 ? 新的桌面提供多文档管理、个人定制桌面以及常用  命令定义快捷键； ? 功能更强的数组编辑器和工作空间浏览器，用户可 以更方便地浏览、编辑和图形化变量； ? M_Lint代码分析器可以方便用户修改代码，不仅提 高了程序性能，而且增强了可维护性； ? 改进后的编辑器可以方便用户选择执行M文件中的