实验九连续信号时域分析的MATLAB实现.docVIP

实验九 连续信号时域分析的MATLAB实现 一. 实验目的 熟悉MATLAB软件平台； 掌握MATLAB编程方法常用语句可视化绘图技术； 编程实现常用信号及其运算MATLAB方法。 二. 实验原理 连续信号是指自变量的取值范围是连续的，且对于一切自变量的取值，除了有若干个不连续点外，信号都有确定的值与之对应。严格来说，MATLAB并不能处理连续信号，而是用等时间间隔点的样值来近似表示连续信号。当取样时间间隔足够小时，这些离散的样值就能较好地近似连续信号。矩阵是MATLAB进行数据处理的基本单元，矩阵运算是MATLAB最重要的运算。通常意义上的数量（也称为标量）在MATLAB系统中是作为1×1的矩阵来处理的，而向量实际上是仅有一行或者一列的矩阵。通常用向量表示信号的时间取值范围，如n = -5:5，但信号x(n)、向量n本身的下标都是从1开始的，因此必须用一个与向量x等长的定位时间变量n，以及向量x，才能完整地表示序列x(n)。在MATLAB可视化绘图中，对于以t为自变量的连续信号，在绘图时统一用plot函数；而对n为自变量的离散序列，在绘图时统一用stem函数。 ?(t)和单位冲激信号?(t)的m文件 在自己的工作目录中创建如下两个m文件。 （1）单位阶跃信号?(t) function f=heaviside(t) %阶跃信号 f=(t0); （2）单位冲激信号?(t) function chongji(t1,t2,t0) %冲激信号?(t- t0)，t1和t2分为起始时间和终止时间 dt=0.01; t=t1:dt:t2; n=length(t); x=zeros(1,n); x(1,(t0-t1)/dt+1)=1/dt; stairs(t,x); %以阶梯方式绘画 axis([t1,t2,0,1.1/dt]) 2.连续信号的相加、相乘、时移、反转和尺度变换等基本运算 （1）两个连续信号的相加 在MATLAB中要实现两个连续信号f1(t)、f2(t)的相加，可用如下语句： x=symadd(f1,f2) 或 x=f1+f2 % x(t)= f1(t)+f2(t) （2）两个连续信号的相乘 在MATLAB中要实现两个连续信号f1(t)、f2(t)的相乘，可用如下语句： x=symmul(f1,f2) 或 x=f1*f2 % x(t)= f1(t) f2(t) （3）连续信号的平移 要实现连续信号f(t)向右平移t0，MATLAB语句格式为： x=subs(f,t,t-t0) % x(t)= f(t-t0) （4）连续信号的反转 要实现连续信号f(t)的反转，MATLAB语句格式为： x=subs(f,t,-t) % x(t)= f(-t) （5）连续信号的尺度变换 要实现连续信号f(t)的尺度变换，MATLAB语句格式为： x=subs(f,t,a*t) % x(t)= f(at) 要实现连续信号f(t)的平移、尺度变换的综合运算，MATLAB语句格式为： x=subs(f,t,a*t-b) % x(t)= f(at-b) 例9-1：信号编制一个m文件，?(t)-?(t-2)]、 f2(t)=sin(5?t)[?(t)-?(t-2)]、f1(-2t+1)、f1(t) f2(t-1)的波形。 解：实现所要求运算的m文件syms t f1=sym(exp(-0.5*t)*(heaviside(t)-heaviside(t-2))); f2=sym(sin(5*pi*t)*(heaviside(t)-heaviside(t-2))); subplot(2,2,1),ezplot(f1);title(f1); subplot(2,2,2), ezplot(f2);title(f2); f3=subs(f1,t,-2*t+1); subplot(2,2,3);ezplot(f3);title(f3); f4=f1*subs(f2,t,t-1); subplot(2,2,4);ezplot(f4);title(f4); 运算结果如图9-1所示。其中，图9-1(a)、(b)、(c)、(d)分别为f1(t)、f2(t) 、f1(-t+1) 、f1(t) f2(t-1)的波形。 3.连续信号的卷积运算 由于MATLAB运算的特点，两个连续信号f1(t)、f2(t)的卷积f(t)=f1(t)*f2(t)，用MATLAB实现的过程应为： （1）将连续信号f1(t)、f2(t)以时间间隔Δ进行取样，得离散序列f1(k?)、f2(k?)； （2）构造f1(k?)、f2(k?)与相对应的时间向量k1和k2； （3）调用con