MATLAB应用技术 作者 于润伟 MATLAB应用技术6.pptVIP

MATLAB应用技术 主编：于润伟 第6章 信号处理 本章要点 连续系统分析 离散系统分析 数字滤波器的设计 6.1 时间信号 信号是信息的载体，信息是通过各种各样的信号来进行传递的，携带着不同信息的不同信号之间可能相互叠加，而且信号在传递过程中也不可避免地受到噪声的干扰。因此，如何从所收到的信号中提取所需要的信息，是信号处理的基本任务。 时间信号是指将信号表示成时间的函数，信号的时间特性是指信号的波形出现的先后、持续时间的长短、随时间变化的快慢等。 一、连续时间信号 从严格意义上讲，MATLAB并不能处理连续信号。在MATLAB中是用连续信号在等间隔时间的采样值来近似表示连续信号的。当取样时间间隔足够小时，这些离散的样点就能够很好地近似出连续信号。 1．单位阶跃信号 程序如下： clear clc t=-3:0.01:3; u=t=0; plot(t,u) axis([-3,3,0,1.5]) xlabel(时间（t）);ylabel(幅度（f）);title(单位阶跃信号) 2．单位冲击信号 程序如下： clear clc t0=0; t1=-2; dt=0.001; t2=2; t=t1:dt:t2; n=length(t); ％时间样本的长度 x=zeros(1,n); ％各样本点信号的赋值为零 k=floor((t0-t1)/dt); ％向无穷小方向取整 x(k)=1/dt; ％在第k点给出的样本值赋为1/dt，近似于无穷大的表示 stairs(t,x) axis([-2,2,0,20]); xlabel(时间（t）);ylabel(幅度（f）);title(单位冲击信号) 3．单边衰减指数信号 程序如下： clear clc t1=-1; dt=0.1; t2=10; t=t1:dt:t2; k=1.5; a=0.5; x=k*exp(-a*t); plot(t,x) axis([t1,t2,-0.5,2]); xlabel(时间（t）);ylabel(幅度（f）);title(单边衰减指数信号) 4．复指数信号 clear clc t1=0; dt=0.01; t2=3; t=t1:dt:t2; a=3;b=4; z=exp(-(a+j*b)*t); subplot(221);plot(t,imag(z)); xlabel(时间);ylabel(幅度);title(实部); subplot(222);plot(t,real(z)); xlabel(时间);ylabel(幅度);title(虚部); subplot(223);plot(t,abs(z)); xlabel(时间);ylabel(幅度);title(模); subplot(224);plot(t,angle(z)); xlabel(时间);ylabel(幅度);title(相角) 二、离散时间信号 离散时间信号是仅在某些不连续的时间上有数值，而在其他时间点上没有定义的一类信号。离散信号一般可以利用模数转换由连续信号采样得到。 1．单位脉冲信号 clear clc for n=-5:5 y=0; if n==0 y=1; end hold on stem(n,y,r*) axis([-5 5 0 2]) xlabel(n) ylabel(x(n)) end hold off 2．单位阶跃信号 clear clc t=0; for n=-5:5 if n=t y=1; end hold on stem(n,y,r*) axis([-5 5 0 2]) xlabel(n) ylabel(x(n)) end hold off 3．实数指数信号 clear clc n=-5:5; y=2.^n; stem(n,y,r*) axis([-5 5 0 40]) xlabel(n) ylabel(x(n)) 4．复指数信号 clear clc n=-10:10; a=0.3; w=0.4; M_y=exp(a*n); A_y=(180/pi)*atan(w/a)*n; subplot(121) stem(n,M_y,r*);title(模); subplot(122) stem(n,A_y,r*);title(幅角); 6.2 连续系统分析 连续系统是指输入、输出信号在时间上均为连续变化的系统。通常使用微分方程对连续系统进行描述。 一、LTI连续系统时域分析 LTI（线性时不变）连续系统是最为常见的一类系统，线性包含可