第MATLAB在信号处理中的应用.ppt

* 第 4 章 MATLAB在信号处理中的应用 第4章 MATLAB在信号处理中的应用 4.1 信号及其表示 4.2 信号的基本运算 4.3 信号的能量和功率 4.4 线性时不变系统 4.5 线性时不变系统的响应 4.6 线性时不变系统的频率响应 4.7 傅里叶(Fourier)变换 4.8 IIR数字滤波器的设计方法 4.9 FIR数字滤波器设计 4.1 信号及其表示 离散时间信号(序列)：时间离散，如x(nT)=x(t)|t=nT. 4.1.1 连续时间信号的表示 连续时间信号：时间变化连续。如y=x(t) 例4-1：用MATLAB命令绘出连续时间信号 关于t的曲线，t的范围为0--30s，并以 0.1s递增。 t=0:0.1:30; x=exp(-0.707*t).*sin(2/3.*t); plot(t,x);grid; ylabel(x(t));xlabel(Time(sec)) 运行结果如图4.1所示。 图4.1 连续时间信号图形 4.1.2 工具箱中的信号产生函数 电压控制振荡器 vco 产生高斯单脉冲信号 gmonopuls 产生高斯正弦脉冲 gauspuls 产生Dirichlet或周期sinc函数 diric 产生调频余弦信号 chirp 产生非周期的三角波 tripuls 产生sinc函数波形 sinc 产生非周期的方波信号 rectpule 产生方波信号 square 产生冲激串 pulstran 产生锯齿波或三角波 sawtooth 功能 函数名 功能 函数名 1、产生锯齿波或三角波信号函数sawtooth 格式一：x=sawtooth(t) 功能：产生周期为 ，幅值从-1～+1的锯齿波。 格式二：x=sawtooth(t,width) 功能：产生三角波，width在0～1之间。 例4-2 产生周期为0.02的三角波。 Fs=10000;t=0:1/Fs:1; x1=sawtooth(2*pi*50*t,0); x2=sawtooth(2*pi*50*t,1) subplot(2,1,1),plot(t,x1),axis([0,0.2,-1,1]); subplot(2,1,2),plot(t,x2),axis([0,0.2,-1,1]); 图4.2 三角波波形 2、产生方波信号函数square 格式一：x=square(t) 功能：产生周期为 ，幅值从-1～+1的方波。 格式二：x=square(t,duty) 功能：产生指定周期的方波，duty为正半周期的比例。 例4-3 产生周期为0.02的方波。 Fs=10000;t=0:1/Fs:1; x1=square(2*pi*50*t,20); x2=square(2*pi*50*t,80) subplot(2,1,1),plot(t,x1),axis([0,0.2,-1.5,1.5]); subplot(2,1,2),plot(t,x2),axis([0,0.2,-1.5,1.5]); 图4.3 方波波形 4.1.3离散时间信号的表示 在MATLAB中，离散时间信号x(n)的表示：需用一个向量x表示序列幅值，用另一个等长的定位时间变量n，才能完整地表示一个序列。 [例4-10] 绘制离散时间信号的棒状图。其中x(-1)=-1, x(0)=1, x(1)=2, x(2)=1, x(3)=0, x(4)=-1。 n=-3:5; %定位时间变量 x=[0,0,-1,1,2,1,0,-1,0]; stem(n,x); grid; % 绘制棒状图 line([-3,5],[0,0]); %画x轴线 xlabel(n); ylabel(x[n]) 图 4.11 离散时间信号图形 4.1.4几种常用离散时间信号的表示 1．单位脉冲序列 直接实现：x=zeros(1,n); x(1,n0)=1; 或者： n=[ns:nf] x=[(n-n0)==0] 序列x的起始点用ns表示，终止点用nf表示。因此，序列的长度length(x)可写为： n=[ns:nf] 或者n=[ns:ns+length(x)-1] 直接实现：n=[ns:nf]; x=[(n-n0)=0]; 2．单位阶跃序列 3．实指数序列 直接实现：n=[ns:nf]; x=a.^n; 4．复指数序列 直接实现：n=[ns:nf]; x=exp((sigema+jw)*n); 5