matlab第10章__MATLAB在信号探析.ppt

第10章 MATLAB在信号 处理中的应用　　　　;本章学习目标;主要内容;10.1 LTI系统的时域分析;②冲激函数 function y=chongji(t1,t2,t0) dt=0.01; t=t1:dt:t2; n=length(t); y=zeros(1,n); %先使所有的函数值为零 %找出t0时刻，使其对应的函数值为1/dt y(1,(-t0-t1)/dt+1)=1/dt; ;③门函数 可由系统内部函数rectpuls表示 y= rectpuls(t) 产生高度为1，宽度为1的门函数 y= rectpuls(t,W) 产生高度为1，宽度为W的门函数 ④抽样函数 可以用下面的自定义函数表示 function f=Sa(t) f=sinc(t./pi); ⑤符号函数 可由MATLAB内部函数sign表示 ; （4）;解： (1) t=-6:0.01:6 f =2*t.*exp(-3*abs(t)); plot(t,f,linewidth,2),grid on line([-6 6],[0 0],color,r‘) line([0 0],[min(f)-0.1 max(f)+0.1],color,r‘) axis([-6 6 min(f)-0.1 max(f)+0.1]) title(函数2te^{-3|t|}的波形‘) xlabel(t‘) ylabel(f(t)‘) ;(2)syms t f=sym(sin(pi*(t-1))*(u(2-t)-u(-t))); ezplot(f) grid on line([-6 6],[0 0],‘color’,‘r’) line([0 0],[-1.1 1.1],color,r‘) ;(3)t1=-3,t2=3,t0=-2,dt=0.01; t=t1:dt:t2 f=chongji(t1,t2,t0); stairs(t,f) grid on axis([t1,t2,0,1.2/dt title(冲激函数的波形’) xlabel(‘t’) ylabel(f(t)‘) ;(4) t=-3:0.01:3; f=u(1-abs(t)).*sign(t); subplot(221) plot(t,f,linewidth,2) grid on line([-3 3],[0 0],color,r‘) line([0 0],[min(f)-0.1 max(f)+0.1],color,r) axis([-3 3 min(f)-0.1 max(f)+0.1]) title(函数u(1-|t|)sgn(t)的波形‘) xlabel(t‘) ylabel(f(t)) ;2.离散信号的表示 ①单位序列 function [x,n]=impseq(k1,k2,k0) k=k1:k2; x=[(k+k0)==0]; ②单位阶跃序列 function [x,k]=stepseq(k1,k2,k0) k=k1:k2; x=[(k+k0)=0]; ; 在用MATLAB表示离散序列并将其可视化时，我们要注意以下几点：第一，与连续时间信号不同，离散时间信号无法用符号运算来表示；第二，由于在MATLAB中，矩阵的元素个数是有限的，因此，MATLAB无法表示无限序列；第三，在绘制离散信号波形时，要使用专门的绘制离散数据的stem命令，而不是plot命令。 ;【例10.2】用MATLAB画出以下离散信号的波形 （1） （2） （3） （4） ;k1=-2:3; f1=k1+1; subplot(221) stem(k1,f1,’filled’) title(f(k)=k+1) ,xlabel(k),ylabel(f(k)) k2=0:5;f2=2.^(-k2).*stepseq(0,5,-2); subplot(222),stem(k2,f2,’filled’) title(‘f(k)=2^{-k}u(k-2)’),xlabel(‘k’),ylabel(‘f(k)’) k3=0:5