资源描述：
《实验一、MATLAB编程基础及典型实例》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、实验九连续信号时域分析的MATLAB实现实验九连续信号时域分析的MATLAB实现一.实验目的1.熟悉MATLAB软件平台；2.掌握MATLAB编程方法、常用语句和可视化绘图技术；3.编程实现常用信号及其运算MATLAB实现方法。二.实验原理连续信号是指自变量的取值范围是连续的，且对于一切自变量的取值，除了有若干个不连续点之外，信号都有确定的值与之对应。严格来说，MATLAB并不能处理连续信号，而是用等时间间隔点的样值来近似表示连续信号。当取样时间间隔足够小时，这些离散的样值就能较好地近似连续信号。矩阵是MATLAB进行数据处理的基本单元，矩阵运算是MATLAB最重要的运算。通常意义上的数
2、量（也称为标量）在MATLAB系统中是作为1×1的矩阵来处理的，而向量实际上是仅有一行或者一列的矩阵。通常用向量表示信号的时间取值范围，如n=-5:5，但信号x(n)、向量n本身的下标都是从1开始的，因此必须用一个与向量x等长的定位时间变量n，以及向量x，才能完整地表示序列x(n)。在MATLAB可视化绘图中，对于以t为自变量的连续信号，在绘图时统一用plot函数；而对n为自变量的离散序列，在绘图时统一用stem函数。1.单位阶跃信号e(t)和单位冲激信号d(t)的m文件在自己的工作目录中创建如下两个m文件。（1）单位阶跃信号e(t)functionf=heaviside(t)%阶跃信号
3、f=(t>0);（2）单位冲激信号d(t)functionchongji(t1,t2,t0)%冲激信号d(t-t0)，t1和t2分为起始时间和终止时间dt=0.01;t=t1:dt:t2;52实验九连续信号时域分析的MATLAB实现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,
4、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(
5、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：信号f1(t)、f2(t)如图9-1(a)、(b)所示，编制一个m文件，绘出f1(t)=e-0.5t[e(t)-e(t-2)]、f2(t)=sin(5pt)[e(t)-e(t-2)]、f1(-2t+1)、f1(t)f2(t-1)的波形。解：实现所要求运算的m文件如下，symstf1=sym('exp(-0.5*t)*(heaviside(t)-heaviside(t-2))');f2=sym('sin(5*p
6、i*t)*(heaviside(t)-heaviside(t-2))');52实验九连续信号时域分析的MATLAB实现图9-1subplot(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(
7、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(kD)、f2(kD)；（2）构造f1(kD)、f2(kD)与相对应的时间向量k1和k2；（3）调用conv()函数计算卷积积分f(t)的近
