用matlab实现常用的连续时间信号及其时域运算
精品课件-信号与系统实验(MATLAB版)-第3章
(a) f3(t)=f1(-t)+f1(t) )+f1(t)
(c) f5(t)=f2(t)×f3(t)
(b) f4(t)=-［f1((d) f6(t)=f1(t)×f2(t)
subplot(2，3，1); ezplot(f，［-3，3］);
实验3 连续时间信号的时域基本运算
y1=subs(f，t，t+2); subplot(2，3，2); ezplot(y1，［-5， 1］); y2=subs(f，t，t-2) subplot(2，3，3); ezplot(y2， ［-1，5］); y3=subs(f，t，-t); subplot(2，3，4); ezplot(y3，［-3， 3］); y4=subs(f，t，2*t); subplot(2，3，5); ezplot(y4，［-2， 2］); y5=-f; subplot(2，3，6) ezplot(y5，［-3，3］) ; 注：在运行以上程序时，需先建立Heaviside的M
实验3 连续时间信号的时域基本运算
Heaviside的函数M function ［x，n］=Heaviside(n0，n1，n2) n=［n1：n2］; x=［(n-n0)==0］;
各个信号的波形如图3.9
实验3 连续时间信号的时域基本运算 图 3.9 各个信号的波形
实验3 连续时间信号的时域基本运算
a＜0 时，还必须包含翻转;
实验3 连续时间信号的时域基本运算
(6) 标量乘法： f(t)→af(t)
(7) 倒相：f(t)→-f(t)
实验一 连续时间信号在MATLAB中的表示
k=2;
a=-1.5; b=10;
ft=k*exp((a+i*b)*t);
subplot(2,2,1);plot(t,real(ft));title('实部');axis([0,3,-2,2]);grid on; subplot(2,2,2);plot(t,imag(ft));title('虚部');axis([0,3,-2,2]);grid on; subplot(2,2,3);plot(t,abs(ft));title('模');axis([0,3,0,2]);grid on; subplot(2,2,4);plot(t,angle(ft));title('相角');axis([0,3,-4,4]);grid on;
1.1 实验目的
• 学会使用MATLAB表示常用连续时间信号的 方法
• 观察并熟悉这些信号的波形和特性
1.2 实验原理及实例分析
在某一时间区间内，除若干个不连续点外，如果任意 时刻都可给出确定的函数值，则称该信号为连续时间信号。 从严格意义上讲，MATLAB并不能处理连续时间信号，然而 可以利用连续信号在等时间间隔点的取样值来近似表示连续
出其波形图。
1.2.1 典型信号的MATLAB表示
1. 实指数信号
MATLAB中用exp函数来表示实指数信号，其语句格式为
y=K*exp(a*t) 例1-1 用MATLAB产生单边衰减指数信号 2e
1.5t
u(t ) , 并绘出时间
0 t 3 的波形。
解：MATLAB源程序为
K=2; a=-1.5; t=0:0.01:3; ft=K*exp(a*t); plot(t,ft),grid on axis([0,3,0,2.2]) title('单边指数衰减信号')
实验2 连续时间信号在MATLAB中运算
实验原理：
信号 f (t)的时移就是将信号数学表达式中的自变量 t
用t t0替换，其中 t0为整实数
信号 f (t)的反折就是将信号数学表达式中的自变量 t
用 t 替换
信号 f (t)的尺度变换就是将信号数学表达式中的自变
量 t 用 at 替t0换
例子1：
function f = unc(t) f = heaviside(t+2)- heaviside(t)+(-t+1).*(heaviside(t)- heaviside(t-1));
实验2 连续时间信号在 MATLAB中的运算
实验目的：
1、学会运用MATLAB进行连续时间信号时移，反 折和尺度变换 2、学会运用MATLAB进行连续信号微分、积分运 算 3、学会运用MATLAB进行连续信号相加、相乘运 算 4、学会运用MATLAB进行连续信号的奇偶分解
一、信号的时移、反折和尺度变换
二、连续时间信号的微分与积分运算
实验原理：
diff ( function,' var iable', n)
int( function,' var iable', a,b)
例子1：
(1) y1 sin(ax2 )
(2) y2 x sin x ln x
clear all syms a x y1 y2 ; y1 = sin(a*x^2) ; y2 = x*sin(x)*log(x); dy1 = diff(y1,'x') ; dy2 = diff(y2) ;
例子3：
(1) y1 (x5 ax2 x )dx
2
(2) y2
1 0
ቤተ መጻሕፍቲ ባይዱ
实验一 连续时间信号的Matlab表示与计算
实验一连续时间信号的Matlab表示与计算一、实验目的1、初步学习MATLAB语言，熟悉MATLAB软件的基本使用。
2、掌握用MA TLAB描述连续时间信号方法，能够编写MATLAB程序，实现各种信号的时域变换和运算，并且以图形的方式再现各种信号的波形。
二、实验原理连续信号是指自变量的取值范围是连续的，且对于一切自变量的取值，除了有若干个不连续点之外，信号都有确定的值与之对应。
严格来说，MATLAB并不能处理连续信号，而是用等时间间隔点的样值来近似表示连续信号。
当取样时间间隔足够小时，这些离散的样值就能较好地近似连续信号。
在MATLAB可视化绘图中，对于以t为自变量的连续信号，在绘图时统一用plot函数；而对n为自变量的离散序列，在绘图时统一用stem函数。
对于连续时间信号f（t），可用f、t两个行向量来表示。
例：t=-10：1.5：10；f=sin(t)./ t ;可以产生t= -10~10，间隔1.5的序列以及t tf)sin(=的值。
用命令：plot(t,f)可得如下图形，显然显示效果较差，这是因为t的间隔过大，只要改变为：t=-10：0.5：10；可得图1.2。
图1.1 图1.21. 信号的时域表示方法MATLAB提供了大量用以生成基本信号的函数，比如最常用的指数信号、正弦信号等就是MATLAB的内部函数，即不需要安装任何工具箱就可以调用的函数。
1.1单位阶跃信号u(t)function y=heaviside(t) %阶跃信号y = (t>=0); % y = 1 for t > 0, else y = 01.2单位冲激信号δ(t)function chongji(t1,t2,t0) %冲激信号δ(t- t 0)，t 1和t 2分为起始时间和终止时间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]) 或function y = delta(t)dt = 0.01;y = (u(t)-u(t-dt))/dt;1.3指数信号指数信号atAe 在MATLAB 中可以用exp 函数表示，其调用形式为：y=A*exp(a*t)例如图1-3所示指数衰减信号的MATLAB 源程序如下（取A=1，a=-0.4）：%program7_1 Decaying expponential signalA=1;a=-0.4;t=0:0.01:10;ft=A*exp(a*t);plot(t,ft);grid on;1.4正弦信号正弦信号)cos(ϕω+t A o 和)sin(ϕω+t A o 分别用MATLAB 的内部函数cos 和sin 表示，其调用形式为：)*cos(*phi t A o +ω)*sin(*phi t A o +ω 例如图1-4所示MATLAB 源程序如下（取A=1，πω20=，6/πϕ=）：%program7_2 Sinusoidal signalA=1;w0=2*pi;phi=pi/6;t=0:0.01:8;ft=A*sin(w0*t+phi);plot(t,ft);grid on;图1-3 单边指数衰减信号 图1-4 正弦信号 除了内部函数外，在信号处理工具箱（Signal Processing Toolbox ）中还提供了诸如抽样 函数、矩形波、三角波、周期性矩形波和周期性三角波等在信号处理中常用的信号。
2连续时间信号在MATLAB中的表示
2连续时间信号在MATLAB中的表示2-1.利用MATLAB命令画出下列连续信号的波形图（1）>> t=0:0.01:3;>> ft=2*cos(3*t+pi/4);>> plot(t,ft),grid on;>> axis([0 3 -2.2 2.2]);>> title('2cos(3t+pi/4)')（2）>> t=0:0.01:3;>> ft=2-exp(-t);>> plot(t,ft),grid on;>> title('(2-exp(-t))u(t)')（3）>> t=-1:0.01:1;>> ft=t.*(uCT(t)-uCT(t-1));>> plot(t,ft),grid on>> axis([-1 1 -0.2 1.2]);>> title('t[u(t)-u(t-1)]')（4）>> t=-1:0.01:3;>> ft=(1+cos(pi*t)).*(uCT(t)-uCT(t-2)); >> plot(t,ft),grid on>> axis([-1 3 -0.2 2.2]);>> title('[1+cos(pi*t)][u(t)-u(t-2)]')2-2.利用MATLAB命令画出下列复信号的实部、虚部、模和辐角（1）>> t=0:0.01:3;>> ft=2+exp(i*(pi/4)*t)+exp(i*(pi/2)*t);>> subplot(2,2,1);plot(t,real(ft));title('实部');axis([0 3 0 4]);grid on;>> subplot(2,2,2);plot(t,imag(ft));title('虚部');axis([0 3 0 2]);grid on;>> subplot(2,2,3);plot(t,abs(ft));title('模');axis([0 3 0 4]);grid on;>> subplot(2,2,4);plot(t,angle(ft));title('相角');axis([0 3 0 2]);grid on;（2）t=0:0.01:3;>> ft=2*exp(i*(t+pi/4));>> subplot(2,2,1);plot(t,real(ft));title('实部');axis([0 3 0 2]);grid on;>> subplot(2,2,2);plot(t,imag(ft));title('虚部');axis([0 3 0 2]);grid on;>> subplot(2,2,3);plot(t,abs(ft));title('模');axis([0 3 0 4]);grid on;>> subplot(2,2,4);plot(t,angle(ft));title('相角');axis([0 3 0 4]);grid on;2-3.利用MATLAB命令产生幅度为1、周期为1、占空比为0.5的一个周期矩形脉冲信号>> t=-0.5:0.01:3;>> ft=square(2*pi*t,50);>> plot(t,ft);grid on;axis([-0.5 3 -1.2 1.2]);>> title('幅度为1、周期为1、占空比0.5的周期举行脉冲信号')3连续时间信号在MATLAB中的运算3-1.试用MATLAB命令绘出以下信号的波形图（1）>> syms x t;>> t=-1:0.01:1;>> x=exp(-t).*sin(10*pi*t)+exp(-0.5*t).*sin(9*pi*t);>> plot(t,x)（2）>> syms x t;>> t=-1:0.01:1;>> x=sinc(t).*cos(10*pi*t);>> plot(t,x)3-2.已知连续时间信号f(t)的波形如图3-6所示，试用MATLAB 命令画出下列信号的波形图先画出图3-6：>> t=-2:0.01:2;>>f=(-t-1).*(-uCT(t+2)+uCT(t+1))+uCT(t+1)+uCT(t)-uCT(t-1)-(t-1).*(uCT(t-1)-uCT(t-2))-uC T(t-2);>> plot(t,f)>> axis([-4 4 -1 2])>> title('图3-6')>> t=-2:0.01:2;>> f1=funct2(t-1);>> f2=funct2(2-t);>> f3=funct2(2*t+1);>> f4=funct2(4-t/2);>> f5=(funct2(t)+funct2(-t)).*uCT(t);>> subplot(231);plot(t,f1);grid on;title('f(t-1)');axis([-3 3 -1 2]);>> subplot(232);plot(t,f2);grid on;title('f(2-t)');axis([-3 3 -1 2]);>> subplot(233);plot(t,f3);grid on;title('f(2t-1)');axis([-3 3 -1 2]);>> subplot(234);plot(t,f4);grid on;title('f(4-t/2)');axis([-3 3 -1 2]);>> subplot(235);plot(t,f5);grid on;title('(f(t)+f(-t))u(t)');axis([-3 3 -1 2]);3-3.试用MATLAB命令绘出如图3-7所示信号的偶分量和奇分量>> t=0:0.01:2;>> f=(uCT(t)-uCT(t-2)).*(-t+1);>> plot(t,f);title('图3-7')>> f1=fliplr(f);>> fe=(f+f1)/2;fo=(f-f1)/2;>> subplot(211),plot(t,fe);grid on>> title('fe')>> subplot(212),plot(t,fo);grid on;title('fo')4连续时间信号的卷积计算4-1用MATLAB命令绘出下列信号的卷积积分的时域波形图>> dt=0.001;t1=-0.5:dt:3.5;>> f1=uCT(t1)-uCT(t1-2);>> t2=t1;>> f2=uCT(t2)+uCT(t2-1)-uCT(t2-2)-uCT(t2-3);>> [t,f]=ctsconv(f1,f2,t1,t2,dt);6周期信号的傅里叶级数及频谱分析6-1已知周期三角信号如图6-5所示，试求出该信号的傅里叶级数，利用MATLAB编程实现其各次谐波的叠加，并验证其收敛性。
实验二 连续信号时域分析的MATLAB实现
实验二 连续信号时域分析的MATLAB 实现一. 实验目的1. 熟悉MATLAB 软件平台；2. 掌握MATLAB 编程方法、常用语句和可视化绘图技术；3. 编程实现常用信号及其运算MATLAB 实现方法。
二. 实验原理信号一般是随时间而变化的某些物理量。
按照自变量的取值是否连续，信号分为连续时间信号和离散时间信号，一般用()f t 和()f k 来表示。
若对信号进行时域分析，就需要绘制其波形，如果信号比较复杂，则手工绘制波形就变得很困难，且难以精确。
MATLAB 强大的图形处理功能及符号运算功能，为实现信号的可视化及其时域分析提供了强有力的工具。
根据MATLAB 的数值计算功能和符号运算功能，在MATLAB 中，信号有两种表示方法，一种是用向量来表示，另一种则是用符号运算的方法。
在采用适当的MATLAB 语句表示出信号后，就可以利用MATLAB 中的绘图命令绘制出直观的信号波形了。
下面分别介绍连续时间信号和离散时间信号的MATLAB 表示及其波形绘制方法。
1.连续时间信号所谓连续时间信号，是指其自变量的取值是连续的，并且除了若干不连续的点外，对于一切自变量的取值，信号都有确定的值与之对应。
从严格意义上讲，MATLAB 并不能处理连续信号。
在MATLAB 中，是用连续信号在等时间间隔点上的样值来近似表示的，当取样时间间隔足够小时，这些离散的样值就能较好地近似出连续信号。
在MATLAB 中连续信号可用向量或符号运算功能来表示。
⑴ 向量表示法对于连续时间信号()f t ，可以用两个行向量f 和t 来表示，其中向量t 是用形如12::t t p t =的命令定义的时间范围向量，其中，1t 为信号起始时间，2t 为终止时间，p 为时间间隔。
向量f 为连续信号()f t 在向量t 所定义的时间点上的样值。
例如：对于连续信号sin()()()t f t Sa t t== ，我们可以将它表示成行向量形式，同时用绘图命令plot()函数绘制其波形。
用MATLAB实现常用的连续时间信号及其时域运算
用MATLAB实现常用的连续时间信号及其时域运算信息与通信工程学院通信133班卢承慧一.引言1.要求：1.1用MATLAB语言产生连续时间信号1.2对连续时间信号进行时域运算2.任务：①绘制用于产生以下信号的通用程序，要求对于任意给定的参数都能实现所要求的信号。
调试并运行这些程序，具体产生由指导教师制定的信号并绘制信号波形。
a. f（t）=δ(t-t );b. f (t) = Au(t-t )。
②已知信号波形如图7.6所示，使用MATLAB语言求出下列信号的表达式并绘制出各信号波形。
a.f（-t）；b.f（t-2）；c.f（1-2t）。
图7.6任务②中的f（t）3.思考题编制一通用程序用于产生信号)()cos()(0t t u t Ae t f at -=-ω，要求对于任意给定的参数都能实现所要求的信号。
二.基本原理1.1连续时间信号如果在所讨论的时间间隔内，除若干个不连续点之外，对于任意时间值都可以给出确定的函数值，此信号就称为连续信号。
从严格意义上来讲，MATLAB 不能处理连续时间信号。
在MATLAB 中，使用连续时间信号在等时间间隔点的样值来近似表示连续时间信号的。
当取样时间间隔足够小时，这些离散的样值就能较好的近似出连续时间信号。
由于在MATLAB 中，矩阵的元素个数是有限的，因此MATLAB 无法表示无限序列。
MATLAB 的绘图命令有很多种，其中比较常用的绘制连续时间信号的绘图命令有“plot ”，“stairs ”，“ezplot ”等。
“plot ”适用于绘制平滑的曲线，而“stairs ”适合于绘制具有阶跃形式的图形，“ezplot ”只能用于符号函数的绘图。
1.2单位阶跃信号单位阶跃信号的波形图如图1所示，通常以符号u （t ）表示⎩⎨⎧><=)0( 1)0( 0)(t t t u在跳变点t=0处，函数未定义，或在t=0处规定函数值u （0）=21。
图1用MATLAB 实现单位阶跃信号%t1:起始时刻；t2:终止时刻；t0：跳变时刻function u(t1,t0,t2)t=t1:0.01:t2; %步长值越小，图形越精确 n=length(t); tt=t0:0.01:t2; n1=length(tt);x=[zeros(1,n-n1),ones(1,n1)]; %产生单位阶跃信号 stairs(t,x),grid on title('单位阶跃信号')axis([t1 t2 -0.2 1.1]) %为方便波形顶部避开图 框，改变图框坐标1.3单位冲激信号单位冲激信号是持续时间无穷小、瞬间幅度无穷大、涵盖面积恒1的理想信号。
信号与系统 连续时间信号的Matlab表示与计算
t)
sin(2
t)
2
(b) f (t) sin(t) 2 cos(4t) sin(5t)
(c) f (t) sin( t) 2 cos(2t)
三．实验结果
1 . (a) 单位阶跃信号 u(t) heaviside.m文件： function y=heaviside(t) %定义heaviside函数 y=(t>=0);
(b) 单位冲激信(t) chongji.m 文件： function chongji(t1,t2,t0) %文件名与函数名一致， 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])
函数 f3(t)是由两个三角波脉冲信号相减而成。运行程序，得如右图波形。
3．用 MATALB 画出图 1-15 所示信号。其中，K=10，A=5。
① 程序如下：
K=10; t=0:0.01:2*pi; ft=abs(K*sin(t)); plot(t,ft); set(gca,'xtick',[0 pi 2*pi]) set(gca,'xticklabel','0|π|2π|') axis([0 7 0 11]);
第 2页
信号与系统上机实验一
x1a.m文件：
t=-1:0.01:3; y=heaviside(t); %调用heaviside函数。 plot(t,y) axis([-1,2,-0.5,1.5])
建立heaviside.m文件，输入程序保存， 保存文件名与定义函数名一致。然后运行 x1a.m文件，得到图形如右图：
实验二连续时间信号在MATLAB中的表示和运算
实验二连续时间信号在MATLAB中的表示和运算[实验目的]1．掌握连续时间信号在MATLAB中的表示法。
2．掌握信号的时域运算和变换。
3．进一步掌握部分绘图函数的应用。
[实验原理]在MATLAB中通常用两种方法来表示信号，一种是向量来表示信号，另一种则是用符号运算的方法来表示信号。
一. 向量表示法MATLAB的信号处理工具箱有大量的函数可用于产生信号，这些函数中大部分都要求样间隔的时间值矢量t，是用命令t=0:0.005:1;这个矢量包含每秒200个时间抽样点，或者说抽样频率为200Hz。
要产生离散时间信号的时间值矢量k，例如，从k＝0到k＝1000，使用命令k=0 :1000 ;给定了t或k，就可以开始产生需要的信号。
在MATLAB中，可以精确地表示离散时间信号，因为信号的值对应于矢量的各元素。
但是，MATLAB对连续时间信号只能提供近似表示，由各元素分别代表连续时间信号抽样值的一个矢量来近似。
当使用这种近似时，应该选择足够小的抽样间隔，以保证那些样本值能反映信号的全部细节。
t=0:0.1:30;f=exp(-.1*t).*sin(2/3*t);plot(t,f);gridylabel('f(t)')xlabel('Time(sec)')axis([0 30 -1 1]);在这段程序中，绘制的曲线时，时间坐标值作为元素保存在矢量中。
表达式exp(-.1*t)和sin(2/3*t)分别产生一个矢量，各矢量中的元素等于对应不同时间点处表达式的值。
由这两个表达式生成的两个矢量的对应元素相乘得到矢量，然后用plot 命令绘出该信号的时域波形。
plot命令可以将点与点间用直线连接，当点与点间的距离很小时，绘出的图形就成了光滑的曲线，如图。
二. 符号运算的表示法如前说述，MATLAB可以有两种方法来表示连续时间信号。
用这两种方法均可实现连续信号的时域运算和变换，但用符号运算的方法则较为简便。
实验2 连续时间信号在MATLAB中的运算
信息工程学院实验报告学会运用MA TLAB 进行连续信号的时移、反折和尺度变换；学会运用MATLAB 进行连续信号的相加、相乘运算；学会运用MATLAB 数值计算方法求连续信号的卷积。
实 验 环 境:Windows 7 MATLAB7.1实 验 内 容 及 过 程:4.1 试用MATLAB 命令绘制信号/2()sin(10)sin(9)tt f t e t et ππ--=+的波形图。
4.2 已知信号()()(1)(1)[(1)()]f t u t u t t u t u t =--+-+-，画出()f t 、(2)f t +、()f t -、(21)f t -+的波形。
4.3 求信号1()()(2)f t u t u t =--与2()()(1)(2)(3)f t u t u t u t u t =+-----的卷积结果12()()*()f t f t f t =，并画出12(),()f t f t 和()f t 的波形。
4.4 求信号1()(0.5)(0.5)f t u t u t =+--与自身的卷积结果11()()*()f t f t f t =，并画出1()f t 和()f t 的波形。
4.1、MATLAB 源程序为：clear; clc;t=-2:0.001:3;ft=exp(-t).*sin(10*pi*t)+exp((-t)/2).*sin(9*pi*t); plot(t,ft);grid on axis([-2 3 -2 2]); title('f(t)波形图');波形图如图T2-1所示：图T2-1 4.2、MATLAB源程序为：clear;clc;t=-3:0.01:4;ft1=funct1(t);ft2=funct1(t+2);ft3=funct1(-t);ft4=funct1(-2*t+1);subplot(221)plot(t,ft1);grid ontitle('f(t)');axis([-2 3 -2.5 1.5]);subplot(222)plot(t,ft2);grid ontitle('f(t+2)');axis([-4 3 -2.5 1.5]);subplot(223)plot(t,ft3);grid ontitle('f(-t)');axis([-3 4 -2.5 1.5]);subplot(224)plot(t,ft4);grid ontitle('f(-2*t+1)');axis([-2 3 -2.5 1.5]); 波形图如图T2-2所示：图T2-24.3、MATLAB源程序为：clear;clc;dt = 0.01; t = -1:dt:3.5;f1 = uCT(t)- uCT(t-2);f2 = uCT(t)+uCT(t-1)-uCT(t-2)-uCT(t-3);f = conv(f1,f2)*dt;n =length(f);tt = (0:n-1)*dt-2;subplot(221), plot(t,f1),grid on;axis([-1, 3.5, -0.2,1.2]);title('f1(t)'); xlabel('t');subplot(222), plot(t,f2),grid on;axis([-1, 3.5, -0.2,2.2]);title('f2(t)'); xlabel('t');subplot(212), plot(tt,f),grid on;title('f(t)=f1(t)*f2(t)'); xlabel('t');波形图如图T2-3所示：图T2-34.4、MATLAB源程序为：clear;clc;dt = 0.01; t = -2:dt:2.5;f1 = uCT(t+0.5)- uCT(t-0.5);f = conv(f1,f1)*dt;n =length(f);tt = (0:n-1)*dt-2;subplot(211), plot(t,f1),grid on;axis([-2, 2.5, -0.2,1.2]);title('f1(t)'); xlabel('t');subplot(212), plot(tt,f),grid on;title('f(t)=f1(t)*f1(t)'); xlabel('t'); 波形图如图T2-4所示：图T2-4阅读与思考：MATLAB运算符中.*和* 的区别？可结合例子说明。
实验二连续时间信号在MATLAB中的运算
实验二 连续时间信号在MATLAB 中的运算2.1实验目的１．学会运用MATLAB 进行连续信号时移、反折和尺度变换；２．学会运用MATLAB 进行连续信号微分、积分运算；3．学会运用MATLAB 进行连续信号相加、相乘运算；4．学会运用MATLAB 进行连续信号的奇偶分解。
2.２实验原理及实例分析1 信号的时移、反折和尺度变换信号f(t)的时移就是将信号数学表达式中的自变量t 用t ±t 0替换，其中t 0为正实数。
信号f(t)的反折就是将表达式中的自变量t 用-t 替换。
信号f(t)的尺度变换就是将表达式中的自变量t 用at 替换，其中a 为正实数。
2 连续时间信号的微分及积分运算微分 diff(function ,’variable ’,n)函数 变量 阶数积分 int(function ,’variable ’,a,b)函数 变量 下、上限3 信号的相加与相乘运算信号的相加与相乘是指在同一时刻信号取值的相加与相乘。
因此，MATLAB 对于时间信号的相加与相乘都是基于向量的点运算。
故只需将信号表达式进行相加与相乘即可。
4 信号的奇偶分解从波形角度看，求信号的偶分量和奇分量时，首先是将信号进行反折，得到f(-t)，然后与原信号f(t)进行相加减，再除以2，即可分别得到偶分量f e (t)和奇分量f o (t)。
2.3 编程练习1.试用MA TLAB 命令绘出下列信号的波形图。
（1）)9sin()10sin()(211t e t e t x t t ππ--+=t=0:0.0001:3;x=exp(-t).*sin((10*pi)*t)+exp(-(1/2)*t).*sin(9*pi*t);plot(t,x) axis([-1,2,-2,2])（2） )10cos()(sin )(2t t c t x π=x=sinc(t).*cos(10*pi*t);t=-2:0.01:3;plot(t,x); axis([-2,4,-1,2])2. 已知连续信号f(t)的时域波形如图2-1所示，试用MA TLAB 命令绘出f(t)及其时域变换信号f(-t)、 f(t-1.5)、 f(t+1.5) 、f(0.5t)和 -f(t)的波形。
《MATLAB》连续时间信号的频域分析和连续时间系统的时域分析实验报告
《MATLAB 》连续时间信号的频域分析和连续时间系统的时域分析实验报告1、编写程序Q3_1，绘制下面的信号的波形图：其中，ω0 = 0.5π，要求将一个图形窗口分割成四个子图，分别绘制cos(ω0t)、cos(3ω0t)、cos(5ω0t) 和x(t) 的波形图，给图形加title ，网格线和x 坐标标签，并且程序能够接受从键盘输入式中的项数ｎ。
2、给程序例3_1增加适当的语句，并以Q3_2存盘，使之能够计算例题3-1中的周期方波信号的傅里叶级数的系数，并绘制出信号的幅度谱和相位谱的谱线图。
-+-=)5cos(51)3cos(31)cos()(000t t t t x ωωω∑∞==10)cos()2sin(1n t n n nωπ3.3反复执行程序例3_2，每次执行该程序时，输入不同的N值，并观察所合成的周期方波信号。
通过观察，你了解的吉布斯现象的特点是什么？3.4分别手工计算x1(t) 和x2(t) 的傅里叶级数的系数。
1．利用MATLAB 求齐次微分方程，，起始条件为，，时系统的零输入响应、零状态响应和全响应。
2． 已知某LTI 系统的方程为：其中，。
利用MATLAB 绘出范围内系统零状态响应的波形图。
3．已知系统的微分方程如下，利用MATLAB 求系统冲激响应和阶跃响应的数值解，并绘出其时域波形图。
（1）'''()2''()'()'()y t y t y t x t ++=()()t x t e u t -=(0)1y -='(0)1y -=''(0)2y -=''()5'()6()6()y t y t y t x t ++=()10sin(2)()x t t u t π=05t ≤≤''()3'()2()()y t y t y t x t ++=（2）''()2'()2()'()y t y t y t x t ++=。
连续时间信号和系统时域分析及MATLAB实现(推荐文档)
课程设计任务书学生姓名:专业班级:指导教师:工作单位:题目:连续时间信号和系统时域分析及MATLAB实现初始条件：MATLAB 6.5要求完成的主要任务：一、用MATLAB实现常用连续时间信号的时域波形（通过改变参数，分析其时域特性）。
1、单位阶跃信号，2、单位冲激信号，3、正弦信号，4、实指数信号，5、虚指数信号，6、复指数信号。
二、用MATLAB实现信号的时域运算1、相加，2、相乘，3、数乘，4、微分，5、积分三、用MATLAB实现信号的时域变换（参数变化，分析波形变化）1、反转，2、使移（超时，延时），3、展缩，4、倒相，5、综合变化四、用MATLAB实现信号简单的时域分解1、信号的交直流分解，2、信号的奇偶分解五、用MATLAB实现连续时间系统的卷积积分的仿真波形给出几个典型例子，对每个例子，要求画出对应波形。
六、用MATLAB实现连续时间系统的冲激响应、阶跃响应的仿真波形。
给出几个典型例子，四种调用格式。
七、利用MATLAB实现连续时间系统对正弦信号、实指数信号的零状态响应的仿真波形。
给出几个典型例子，要求可以改变激励的参数，分析波形的变化。
时间安排：学习MATLAB语言的概况第1天学习MATLAB语言的基本知识第2、3天学习MATLAB语言的应用环境，调试命令，绘图能力第4、5天课程设计第6-9天答辩第10天指导教师签名：年月日系主任（或责任教师）签名：年月日目录摘要 (I)ABSTRACT (II)绪论 (1)1 MATLAB简介 (2)1.1MATLAB语言功能 (2)1.2MATLAB语言特点 (2)2常用连续时间信号的时域波形 (3)2.1单位阶跃信号 (3)2.2单位冲激信号 (3)2.3正弦信号 (4)2.4实指数信号 (5)2.5虚指数信号 (5)2.6复指数信号 (6)3 连续时间信号的时域运算 (7)3.1相加 (7)3.2相乘 (7)3.3数乘 (8)3.4微分 (8)3.5积分 (9)4 连续时间信号的时域变换 (10)4.1反转 (10)4.2时移 (10)4.3展缩 (11)4.4倒相 (11)4.5综合变化 (12)5连续时间信号简单的时域分解 (13)5.1信号的交直流分解 (13)5.2信号的奇偶分解 (14)6连续时间系统的卷积积分的仿真波形 (15)7连续时间系统的冲激响应、阶跃响应的仿真波形 (16)7.1 IMPULSE（）函数 (17)7.2 STEP（）函数 (19)8连续时间系统对正弦信号、实指数信号的零状态响应的仿真波形 (21)8.1正弦信号的零状态响应 (21)8.2实指数信号的零状态响应 (22)9小结即心得体会 (24)致谢 (25)参考文献 (26)附录 (27)摘要MATLAB目前已发展成为由MATLAB 语言、MATLAB 工作环境、MATLAB 图形处理系统、MATLAB 数学函数库和MATLAB 应用程序接口五大部分组成的集数值计算、图形处理、程序开发为一体的功能强大的系统。
实验1 信号的时域描述与运算
实验1 信号的时域描述与运算一、实验目的1、掌握信号的MATLAB 表示及其可视化方法。
2、掌握信号基本时域运算的MATLAB 实现方法。
3、利用MATLAB 分析常用信号，加深对信号时域特性的理解。
二、实验原理1、连续时间信号的MATLAB 表示在MATLAB 中连续时间信号可以用两种方法来表示，即向量表示法和符号对象表示法。
在MATLAB 中连续时间信号是用信号等时间间隔采样后的采样值来近似表示的，当采样间隔足够小时，这些采样值就可以很好的近似表示出连续时间信号，这种表示方法称为向量表示法。
表示一个连续时间信号需要使用两个向量，其中一个向量用于表示信号的时间范围，另一个向量表示连续时间信号在该时间范围内的采样值。
例如一个正弦信号可以表示如下：>>t=0:0.01:10; >>x=sin(t); 利用plot(t,x)命令可以绘制上述信号的时域波形。
如果连续时间信号可以用表达式来描述，则还可以用符号表达式来表示信号。
例如对于上述正弦信号，可以用符号对象表示如下：>>syms t; >>x=sin(t); 利用ezplot(x)命令可以绘制上述信号的时域波形。
2、连续时间信号的时域运算对连续时间信号的运算包括两信号相加、相乘、微分、积分，以及移位、反转、尺度变换等。
（1）相加和相乘对于两个采用向量表示的可以直接使用算术运算的运算符“+”和“*”来计算，此时要求表示两信号的向量时间范围和采样间隔相同。
采用符号对象表示的两个信号，可以直接根据符号对象的运算规则运算。
（2）微分和积分对于向量表示法表示的连续时间信号，可以通过数值计算的方法计算信号的微分和积分。
这里微分是用差分来近似求取的，由时间向量[t 1,t 2,⋯,t N ]和采样值向量[x 1,x 2,⋯,x N ]表示的连续时间信号，其微分可以下式实现1()|,1,2,,1k k k t t x xx t k N t+=-'≈=-∆其中∆t 表示采样间隔。
实验一连续时间信号的MATLAB表示与计算实验报告
信号与系统实验报告实验一连续时间信号得MＡＴLAB表示与计算学院专业班级姓名学号指导教师实验报告评分:＿____＿＿题目:连续时间信号得ＭATＬＡB表示与计算一、实验目得(1)初步学习MATLAB语言，熟悉MATＬAB软件得基本使用。
（2)掌握用ＭＡＴLAB描述连续时间信号方法，能够棉靴MＡLTAB 程序,实现各种信号得时域变换与运算，并且以图形得方式再现各种信号得波形。
二、实验内容与步骤1、编写m文件,将实验原理中所给得单位冲激信号与单位阶跃信号得函数文件在MATＬＡＢ文件编辑器中编写好,并分别以文件名ｃhｏｎgｊｉ与heaｖisｉdｅ存入文件夹中以便于使用。
单位阶跃信号：functioｎ y = hｅaｖiｓｉde(t）y ＝ (ｔ>=０);单位冲激信号:funcｔion chonｇｊi（t1, t2, t0)dｔ =0、０1;ｔ= t1:dt:t2;n = ｌength(ｔ）;x ＝ zeros(1,n);x(1, (t0-t1)/dt +1) = 1／dt;staｉrｓ(t,x);axis([t1, ｔ2，０, １、1/dt]）2、编写m文件,并以Q1＿１为文件名存盘，产生信号f(t)=sｉn(2πｔ)，时间范围在-2～2秒之间，dｔ = 0、01。
若将ｄｔ改为0、2,再执行该程序，保存图形，瞧瞧所得图形得效果如何？代码:t1 = [－2：０、01：2]；f1 = ｓiｎ(2*ｐi*t1）;subplot(1,2，1);ｐｌot(t1,f1)；ｔitｌe（'ｄt = 0、01');ｘlaｂel(＇t１'）；ylabel（'siｎ（2*pi*ｔ1)'）;t2＝［-２:0、２:2]；f2 ＝ sin（2*ｐi*t２)；sｕbplot(１,２,２）；plot(t２,f2）;tｉｔle('dt ＝ 0、２'）;xlａbeｌ（＇t2');ylabel('sin(2*pi＊t2）');3、编写m文件,并以Q1_2为文件名存盘绘出 (u（ｔ)-u(t-3））信号,给图形标出名称与坐标轴名称。
MATLAB 3.连续时间LTI系统的时域分析
1
3.连续时间 LTI 系统的时域分析
作业题
2、某连续系统的输入为 e(t)，输出为 y(t)，系统的微分方程为：y ’’(t)+5y ’(6)+6y(t)=3f ’(t)+2f(t) 绘出该系统的单位冲激响应和单位阶跃响应时域图形 若f t = ������ −2������ ������(������)，绘出系统的零状态响应时域图形
将激励信号和单位冲击响应进行卷积，由下图发现结果一致。
6
3.连续时间 LTI 系统的时域分析
作业题
4、如下图所示的电路中，已知 电压 y(t ) 作为系统输出：
R1 R2 R3 4() ， L1 L2 1( H ) ，如果以电阻 R3 上
请画出系统的冲击响应 h(t) 请画出系统在激励e(������) = 12������(������)（v）作用下的零状态响应，要求写出系统微分方程的求 解步骤。 对所求结果进行验证。
零输入响应： eq1='D2y+4*Dy+4*y=0'; cond='Dy(-0.01)=1,y(-0.01)=1'; answer=dsolve(eq1,cond) answer = 103/(100*exp(1/50)*exp(2*t))+ (3*t)/(exp(1/50)*exp(2*t)) >> 全响应： eq1='D2y+4*Dy+4*y=Df+3*f'; eq2='f=exp(-t)*heaviside(t)'; cond='Dy(-0.01)=1,y(-0.01)=1'; answer=dsolve(eq1,eq2,cond); simplify(answer.y) ans = (3*t+2*exp(t+1/50)*heaviside(t)-2*exp(1/50)*heaviside(t)-t*exp(1/50)* heaviside(t) + 103/100)/exp(2*t + 1/50)
MATLAB信号的表示与时域分析,matlab与信号实验连续时间信号的时域分析
MATLAB信号的表⽰与时域分析,matlab与信号实验连续时间信号的时域分析matlab与信号实验连续时间信号的时域分析上机实验 1 连续时间信号的时域分析 ⼀、 实验⽬的 (1 ) 掌握连续时间信号的时域运算的基本⽅法； (2 ) 掌握相关函数的调⽤格式及作⽤； (3 ) 掌握连续信号的基本运算； (4 ) 掌握利⽤计算机进⾏卷积运算的原理和⽅法； (5 ) 熟悉连续信号卷积运算函数 conv 的应⽤。
⼆、 实验原理 信号的基本运算包括信号的相加(减)和相乘(除) 。
信号的时域变换包括信号的平移、翻 转、倒相、尺度变换等，由以下公式所描述。
(1 ) 加(减)：f(t)=f 1 (t)±f 2 (t) (2 ) 乘：f(t)=f 1 (t)×f 2 (t) (3 ) 延时或平移：f(t) →f(t-t 0 ) t 0 >0 时右移；t 0 0 。
MATLAB 程序： clear all; t=0:0.0001:2; y=sin(2*pi*(t)); y1=sin(2*pi*(t-0.2)); plot(t, y, - ,t, y1, -- ); xlabel( t );ylabel( f(t) );title(信号的移位 ); 信号及其移位结果如下图所⽰。
4)翻转 信号的翻转就是将信号的波形以纵轴为对称轴翻转 180°。
将信号 f(t)中的⾃变量 t 替换成- t 即可得到其翻转信号。
MATLAB 程序： clear all; t=0:0.02:1; t1=-1:0.02:0; g1=3*t; g2=3*(-t1); grid on; plot(t, g1, -- ,t1, g2); xlabel( t );ylabel( g(t) );title( 信号的反折 ); 信号及其反折结果如下图所⽰。
2.程序设计实验 (1)已知信号 f1(t)=(-t+4)[U(t)-U(t-4)]，f2(t)=sin(2πt)，⽤ MATLAB 绘出下列信号的时域波形。
连续时间信号在MATLAB中的运算
课程名称：MATLAB 实验实 验 目 的:学会运用MATLAB 进行连续信号的时移、反折和尺度变换；学会运用MATLAB 进行连续信号的相加、相乘运算；学会运用MATLAB 数值计算方法求连续信号的卷积。
实 验 环 境:硬件：PC 机，Inter （R ）Core （TM ）i5-4210U CPU @1.7GHz ，内存4G ；软件：Matlab 版本2013b实 验 内 容 及 过 程:1、试用MATLAB 命令绘制信号/2()sin(10)sin(9)t t f t e t e t ππ--=+的波形图。
2、已知信号()()(1)(1)[(1)()]f t u t u t t u t u t =--+-+-，画出()f t 、(2)f t +、()f t -、(21)f t -+的波形。
3、求信号1()()(2)f t u t u t =--与2()()(1)(2)(3)f t u t u t u t u t =+-----的卷积结果12()()*()f t f t f t =，并画出12(),()f t f t 和()f t 的波形。
4、求信号1()(0.5)(0.5)f t u t u t =+--与自身的卷积结果11()()*()f t f t f t =，并画出1()f t 和()f t 的波形。
实验结果及分析：卷积就是进行反转，平移，但信号又是连续的，通过时间间隔取足够小的离散时间信号的数值计算方法来实现，间隔越小，卷积的结果图像越接近真实的图像。
实验心得：根据例子，之后在自己写，去了解一些函数，比如subplot（3，2，4），三行两列，第四个位置（第二行第一个）subplot（3，2，[4,6]）第四，五，六个位置合成整体（第二行合成一行）。
多去网上找一些例子来实践，尽量是编成函数文件，并且可移植性强，加上注释之类的。
方便下次的使用。
MATLAB实验一连续时间信号在Matlab中的运算
实验一连续时间信号在Matlab中的运算一、实验目的1、学会运用Matlab进行连续时间信号的时移、反褶和尺度变换；2、学会运用Matlab进行连续时间信号微分、积分运算；3、学会运用Matlab进行连续时间信号相加、相乘运算；4、学会运用Matlab进行连续时间信号卷积运算..二、实验原理及实例分析1、信号的时移、反褶和尺度变换信号的平移、反转和尺度变换是针对自变量时间而言的;其数学表达式和波形变换中存在着一定的变化规律..从数学表达式上来看;信号的上述所有计算都是自变量的替换过程..所以在使用Matlab进行连续时间信号的运算时;只需要进行相应的变量代换即可完成相关工作..2、连续时间信号的微分和积分符号运算工具箱有强大的积分运算和求导功能..连续时间信号的微分运算;可使用diff命令函数来完成;其语句格式为：difffunction; ‘variable’;n其中;function表示需要进行求导运算的函数;或者被赋值的符号表达式；variable为求导运算的独立变量；n为求导阶数;默认值为一阶导数..连续时间信号积分运算可以使用int命令函数来完成;其语句格式为：intfunction; ‘variable’; a; b其中;function表示被积函数;或者被赋值的符号表达式；variable为积分变量；a为积分下限;b为积分上限;a和b默认时则求不定积分..3、信号的相加和相乘运算信号的相加和相乘是信号在同一时刻取值的相加和相乘..因此Matlab对于时间信号的相加和相乘都是基于向量的点运算..4、连续信号的卷积运算卷积积分是信号与系统时域分析的重要方法之一..定义为：⎰+∞∞-=*=-2121d )t (f )(f )t (f )t (f )t (f τττMatlab 进行卷积计算可通过符号运算方法和数值计算方法实现.. 1Matlab 符号运算法求连续信号卷积从卷积定义出发;可以利用Matlab 符号运算法求卷积积分;但要注意积分变量和积分限的选取..例：试用Matlab 符号运算法求卷积yt=ut-ut-1*ut-ut-1..2Matlab 数值计算法求连续信号的卷积例：试用Matlab 数值计算法求信号)2t (u )t (u )t (f 1--=和)t (u e )t (f 3t2-=的卷积..三、实验内容1、已知信号的波形课本P11例题;画出()()()()2332----tftftftf，，，的波形图..2、使用微分命令求xsinxlnxy=关于变量x的一阶导数；使用积分命令计算不定积分dxxaxx⎰⎪⎪⎭⎫⎝⎛+-225;定积分()dxxxe x⎰+121..Hz f 12=Ω=π3、已知()()()t t f t t f Ω=Ω=8sin ,sin 21;使用命令画出两信号和及两信号乘积的波形图..其中;4、。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
用MATLAB实现常用的连续时间信号及其时域运算
信息与通信工程学院通信133班卢承慧
一.引言
1.要求：
1.1用MATLAB语言产生连续时间信号
1.2对连续时间信号进行时域运算
2.任务：
①绘制用于产生以下信号的通用程序，要求对于任意给定的参数都能实现所要求的信号。
调试并运行这些程序，具体产生由指导教师制定的信号并绘制信号波形。
t );
a. f（t）=δ(t-
t )。
b. f (t) = Au(t-
②已知信号波形如图7.6所示，使用MATLAB语言求出下列信号的表达式并绘制出各信号波形。
a.f（-t）；
b.f（t-2）；
c.f（1-2t）。
图7.6任务②中的f（t）
3.思考题
编制一通用程序用于产生信号
)()cos()(0t t u t Ae t f at -=-ω，要求
对于任意给定的参数都能实现所要求的信号。
二.基本原理
1.1连续时间信号
如果在所讨论的时间间隔内，除若干个不连续点之外，对于任意时间值都可以给出确定的函数值，此信号就称为连续信号。
从严格意义上来讲，MATLAB 不能处理连续时间信号。
在MATLAB 中，使用连续时间信号在等时间间隔点的样值来近似表示连续时间信号的。
当取样时间间隔足够小时，这些离散的样值就能较好的近似出连续时间信号。
由于在MATLAB 中，矩阵的元素个数是有限的，因此MATLAB 无法表示无限序列。
MATLAB 的绘图命令有很多种，其中比较常用的绘制连续时间信号的绘图命令有“plot ”，“stairs ”，“ezplot ”等。
“plot ”适用于绘制平滑的曲线，而“stairs ”适合于绘制具有阶跃形式的图形，“ezplot ”只能用于符号函数的绘图。
1.2单位阶跃信号
单位阶跃信号的波形图如图1所示，通常以符号u （t ）表示
⎩
⎨⎧><=)0( 1)0( 0)(t t t u
在跳变点t=0处，函数未定义，或在t=0处规定函数值u （0）=
2
1。
图1
用MATLAB 实现单位阶跃信号
%t1:起始时刻；t2:终止时刻；t0：跳变时刻
function u(t1,t0,t2)
t=t1:0.01:t2; %步长值越小，图形越精确 n=length(t); tt=t0:0.01:t2; n1=length(tt);
x=[zeros(1,n-n1),ones(1,n1)]; %产生单位阶跃信号 stairs(t,x),grid on title('单位阶跃信号')
axis([t1 t2 -0.2 1.1]) %为方便波形顶部避开图 框，改变图框坐标
1.3单位冲激信号
单位冲激信号是持续时间无穷小、瞬间幅度无穷大、涵盖面积恒1的理想信号。
记为)(t δ，其表达式为
⎪⎩
⎪⎨⎧≠==∞
-∞
=⎰)0t ( 0)(1)()(当t dt t t f δδ
严格说来，MATLAB 是不能表示单位冲击信号的，但我们可以用时间宽度为dt 、高度为1/dt 的矩形脉冲来近似地表示冲激信号。
当dt 趋近零时，就较好地近似出冲激信号的实际波形。
用MATLAB 实现单位冲激信号：
clear
t1=input('输入起始时刻t1='); t0=input('输入跳变时刻t0='); t2=input('输入终止时刻t2='); t=t1:0.001:t2; n=length(t); tt=t1:0.001:t0;
n1=length(tt); x=zeros(1,n); x(n1)=1/0.001; plot(t,x),grid on title('单位冲激信号') axis([t1 t2 -0.2 1.1])
1.4连续时间信号的时域运算
如前所述，MATLAB 可以用向量或符号两种方法来表示连续信号，而只能用向量表示离散信号，单用符号运算的方法进行连续信号的时域运算较为方便。
用MATLAB 的符号运算命令来表示连续信号的时域运算，然后用“ezplot ”命令绘制出其结果波形图。
在下面的命令中，f ，f1，f2都是用符号表达式表达的连续信号，s 为运算结果信号的符号表达式，t0和a 均为常数。
1.4.1 相加
)()(21t f t f ±
1.4.2 相乘
)()(21t f t f ⋅
1.4.3 平移
⎩⎨⎧<>+→，右移，左移
00)()(0
00t t t t f t f
若f （t ）表达式的自变量t 更换为（t+0t ）（0t 为正或负实数），则f （t+0t ）相当于f （t ）波形在t 轴上的整体移动，当0t >0时，波形左移，当0t <0时，波形右移。
1.4.4 反褶
)()(t f t f -→
信号反褶表示将f （t ）的自变量t 更换为-t ，此时f （-t ）的波形相当于将f （t ）以t=0为轴反褶过来。
1.4.5 尺度
⎩
⎨⎧<<>→扩展压缩
1,a 0 ,1)()(a at f t f
如果将信号f （t ）的自变量t 乘以正实系数a ，则信号波形f （at ）将是f （t ）波形的压缩(a>1)或扩展（a<1）。
相加：s=f1+f2,ezplot(s) 相乘：s=f1*f2,ezplot(s)
平移：s=subs(f,t,t-t0) ,ezplot(s) 反褶：s=subs(f,t,-t) ,ezplot(s) 尺度：s=subs(f,t,a*t) ,ezplot(s)
三.实现方法
论述具体实现方法及编程思路，以及具体程序的编制。
1.任务①
a. f （t ） =δ(t-0t );
编程思路：建立一个单位冲激信号，对该信号进行平移运算。
代码如下：
clear
t1=input('输入起始时刻t1='); t0=input('输入跳变时刻t0='); t2=input('输入终止时刻t2='); t=t1:0.001:t2; n=length(t); tt=t1:0.001:t0; n1=length(tt); x=zeros(1,n); x(n1)=1/0.001; plot(t,x),grid on title('单位冲激信号') axis([t1 t2 -0.2 1.1])
b ． f （t ） =Au(t-0t );
编程思路：建立一个单位阶跃信号，对该信号进行平移运算并乘以一个幅度
值。
代码如下：
function f=u(t)
f=(t>0);
clear
close all
A=input('输入幅度值A=');
t0=input('输入跳变值t0=');
syms t
f=sym('u(t)');
f1=A*f;
y=subs(f1,t,t-t0)
ezplot(y,[-3,3])
2.任务②：
编程思路：根据图7.6写出f(t)=2u(t)-u(t-1)-u(t-2),用符号运算命令来表示以下连续信号的时域运算，然后用“ezplot”命令绘制出其结果图。
a.f（-t）；
b.f（t-2）；
c.f（1-2t）。
代码如下：
function f=u(t)
f=(t>0);
syms t
f=sym('2*u(t)-u(t-1)-u(t-2)');
subplot(2,2,1),ezplot(f,[-3,3]),grid on ,title('f(t)')
y1=subs(f,t,-t);
subplot(2,2,2),ezplot(y1,[-3,3]),grid on ,title('f(-t)')
y2=subs(f,t,t-2);
subplot(2,2,3),ezplot(y2,[-3,3]),grid on ,title('f(t-2)')
y3=subs(f,t,1-2*t);
subplot(2,2,4),ezplot(y3,[-3,3]),grid on ,title('f(1-2t)')
3.思考题
编程思路：建立一个幅值为A 的指数信号at Ae -，
一个余弦信号)cos(t ω和一个平移值为0t 的单位阶跃信号)(0t t u -，然后对它们进行相乘运算。
代码如下：
四、结果验证
任务①：
a ．f （t ） =δ(t-0t );
运行结果：
图1
图2
t );
b．f（t） =Au(t-
运行结果：
图3
图4 任务②：
a.f（-t）；
b.f（t-2）；
c.f（1-2t）。
图5
图6 3.思考题
图7
图8
参考文献
[1] 郑君里, 应启珩, 杨为理. 信号与系统引论. 北京: 高等教育出版社, 2009.
[2] 李敏,陈兴文.信号分析与处理的软硬件实现.大连:海事大学出版社,2014.。