试验一信号时域分析.DOC

实验一 信号时域分析 1 实验目的 (1) 熟悉和掌握常用信号仿真分析的MATLAB函数。 (2) 掌握连续时间和离散时间信号的MATLAB产生。 (3) 掌握将非周期信号经周期延拓形成周期信号的MATLAB编程方法。 (4) 掌握用MATLAB实现信号的基本运算。 2 实验原理 信号连续时间信号模拟信号指数信号Aeat y = A*exp(a*t) 指数序列ak 幂运算a.^k实现 正弦型信号 内部函数cos( )和sin( ) 抽样函数Sa(t) sinc(t/π) 矩形脉冲信号 y=rectpuls(t,width) 2.1.2用信号的波形图来描述信号 用函数曲线表示一个信号，如图1-1所示。 图1-1 连续、离散时间信号的波形 2.1.3将信号用一个数据序列来表示 对于离散时间信号，还可以表示成一个数的序列，例如： x[n]={...., 0.1, 1.1, -1.2, 0, 1.3, ….} ↑n=0 在《信号与系统》中上述三种信号的描述方法经常使用。 2.2 信号的MATLAB仿真 在MATLAB中，无论是连续时间信号还是离散时间信号，MATLAB都是用一个数字序列来表示信号，这个数字序列在MATLAB中叫做向量。通常的情况下，需要与时间变量相对应。 如前所述，MATLAB有很多内部数学函数可以用来产生这样的数字序列，例如sin()、cos()、exp()等函数可以直接产生一个按照正弦、余弦或指数规律变化的数字序列。 在绘制信号的波形图时，有时我们需要将若干个图形绘制在同一个图形窗口中，这就需要使用MATLAB的图形分割函数subplot()，其用法是在绘图函数stem或plot之前，使用图形分割函数subplot(n1,n2,n3)，其中的参数n1，n2和n3的含义是，该函数将把一个图形窗口分割成n1* n2个子图，即将绘制的图形绘制在第n3个子图中。 2.2.1连续时间信号的仿真 (1) 正弦信号 % Program1_1 dt = 0.001;t = -1:dt:1; w=pi;ft= sin(2*w*t); plot(t,ft),title(f(t)=sin(2πt)) xlabel(t),grid on axis([-1,1,-1.5,1.5]) (2) 抽样信号f(t)=Sa(t)： % Program1_2 t=-10:0.03:10; plot(t,sinc(t/pi)) (3) 矩形脉冲信号 % Program1_3 T=2; %脉冲宽度 A=2; %脉冲高度 t=-3*T:0.001: 3*T; ft=A*rectpuls(t-2*T,T); %脉冲中心在2*T位置 plot(t,ft),axis([2 6 -1 3]),box off 运行结果如图1-2所示。 图1-2 矩形脉冲信号 (4) 实指数信号 % Program1_4 t=-2:0.01:2; ft1=exp(t);ft2=exp(-2*t); subplot(211),plot(t,ft1) subplot(212),plot(t,ft2) (5) 虚指数信号 % Program1_5 t=0:0.01:15;w=pi/4;A=2; X=A*exp(j*w*t); Xr=real(X);Xi=imag(X); Xa=abs(X);Xn=angle(X); subplot(221),plot(t, Xr) axis([0 15 -(max(Xa)+0.5) max(Xa)+0.5]),title(实部) subplot(223),plot(t, Xi) axis([0 15 -(max(Xa)+0.5) max(Xa)+0.5]),title(虚部) subplot(222),plot(t, Xa) axis([0 15 0 max(Xa)+1]),title(模) subplot(224),plot(t, Xn) axis([0 15 -(max(Xn)+1) max(Xn)+1]),title(相角) (6) 复指数信号 % Program1_6 t=0:0.01:3;a=-1;b=10; z=exp((a+j*b)*t); subplot(221),plot(t,real(z)),title(实部) subplot(223),plot(t,imag(z)),title(虚部) subplot(222),plot(t,abs(z)), title(模) subplot(224),plot(t,angle(z)),title(相角) 在《信号与系统》课程中，单位阶跃信