第九章 时域信号描述及运算课件.ppt

2）离散周期方波 Square函数 例：n=-10:10; x=square(pi/4*n,50); stem(n,x) * * 第二篇 MATLAB在信号处理及通信中应用 主要讲述： MATLAB与信号系统 MATLAB在数字信号处理中的应用 MATLAB与通信系统 第九章 时域信号描述及运算 9.1 连续时间信号 一、连续时间信号的表示方法 1、向量表示法 可以用两个行向量表示（即一个是信号值，另一个是时间轴的值）（其绘图命令为plot） 2、符号运算表示方法 定义符号函数的方法：f=sym(‘expression’) 符号函数绘图命令为ezplot 例：对于函数 t=-10:1.5:10 f=sin(t)./t plot(t,f) title(f(t)=Sa(t)) xlabel(t) axis([-10,10,-0.4,1.1]) 采样间隔加大 t=-10:0.02:10; f=sin(t)./t plot(t,f) title(f(t)=Sa(t)) xlabel(t) axis([-10,10,-0.4,1.1]) 例：f=sin(pi/4*t) f=sym(sin(pi/4*t)) ezplot(f,[-16,16]) 二、常用连续时间信号 1、周期信号 1）三角信号 2）周期方波信号 例：产生一个幅度为1，基频为3Hz，占空比为20%的周期方波 A=1; f0=3; t=0:0.0001:2.5; w0=2*pi*f0; y=A*square(w0*t,20); plot(t,y) axis([0 2.5 -2 2]) 3）锯齿波信号 例： A=1; f0=3; t=0:0.0001:2.5; w0=2*pi*f0; y=A*sawtooth(w0*t,0.8); plot(t,y) axis([0 2.5 -2 2]) 2、非周期信号 1）单位阶跃信号u(t) a、t=-2:0.02:6; y=t=0; b、f=[0 0 0…1 1 1…] t=[-t……..0 +…..] 2）实指数信号 例： A=1;a=-1; t=0:0.001:2.5; y=A*exp(a*t); plot(t,y) 3）矩形信号 rectpuls(t)：产生高度为1，宽度为1，关于t=0对称的矩形脉冲信号。 rectpuls(t,w)：W为宽度 例如： t=-5:0.01:5; plot(t,rectpuls(t,2)) 4）三角形脉冲信号 tripuls(t)：宽度为1，波峰在t=0时刻 tripuls(t,w)：W为宽度 tripuls(t,w,s)：s（-1，1），波峰向左右偏转的相对大小值。 例： t=-4:0.001:4; Plot(t,tripuls(t)); Plot(t,tripuls(t,2)); Plot(t,tripuls(t,2,0.5)); 5）取样函数 由sinc(x)得到sa(x)： T=-10:1/500:10; X=sinc(t/pi); Plot(t,x); 笛拉克函数： diric(x,N) 例： t=linspace(-pi,pi,1000); plot(t,diric(t,100)) 6）线性调频信号（chirp信号） y=chirp(t,f0,t1,f1) y=chirp(t,f0,t1,f1,’method’) y=chirp(t,f0,t1,f1,’method’,phi) method的选项： 例如： t=-5:0.001:5; y=chirp(t,1,10,10); plot(t,y) 3、复指数信号 t=0:0.01:15; x=2*exp(j*pi/4*t); xr=real(x); xi=imag(x); xa=abs(x); xn=angle(x); subplot(2,2,1) plot(t,xr) subplot(2,2,3) plot(t,xi) subplot(2,2,2) plot(t,xa) subplot(2,2,4) plot(t,xn) 例：求 的图形 t=0:0.01:3; a=-1;b=10; z=exp((a+j*b)*t); aa=real(z); bb=imag(z); cc=abs(z); dd=angle(z); subplot(2,2,1) plot(t,aa) subplot(2,2,3) plot(t,bb) subplot(2,2,2) plot(t,cc) subplot(2,2,4)，plot(t,dd) 例：求 的图形 t=0:0.01: