第5章_信号与系统_MATLAB
信号与系统第五章习题答案
５．２ 教材习题同步解析
5.1 设信号 f (t ) 为包含 0~ ω m 的频带有限信号，试确定 f (3t ) 的抽样频率。 【知识点窍】主要考察奈奎斯特频率的概念。 【逻辑推理】时域的信号的压缩，在频域中将会扩展。时域中压缩多少培，频域中将扩展多少 培。另外，抽样频率等于奈奎斯特频率与 2π 之比。 解：因为信号 f (t ) 为包含 0~ ω m 的频带有限信号，则信号 f (3t ) 为包含 0~3ω m 的频带有限信 号。
g [n ] = 0 。
g [n ] = 0.8 g [n − 1] + ε [n]
可依次迭代得
g [0] = 0.8 g[− 1] + ε [0] = 1
g [1] = 0.8 g [0] + ε [1] = 0.8 + 1 L
g [2 ] = 0.8 g [ 1] + ε [2] = 0.8 2 + 0.8 + 1 g [n ] = 0.8 g [n − 1] + ε [n ] = 0.8 n + 0.8 n−1 + L + 0.82 + 0.8 + 1 1 − 0.8 n+1 1 − 0 .8 = 5 1 − 0.8 n+1 =
273
②将序列 f 2 [− i ] 沿正 n 轴平移 n 个单位，成为 f 2 [n − i ] ； ③求乘积 f 1 [i ] f 2 [n − i ] ； ④按式 f 1 [n] ∗ f 2 [n ] = 2）阵列表法 3）解析法：利用卷积和定义求解。 解： f [n] ∗ h[n] = 上式是公比为
第5章_信号与系统_MATLAB
1 0.8
相加用算术运算符“+”实现 相乘用数组运算符“.*”实现
0.6 0.4 0.2 0
例:画信号Aeatcos(w0t+f)的波形
-0.2 -0.4 -0.6 -0.8 0
t=0:0.001:8; A=1; a=-0.4; w0=2*pi;phi=0; ft1=A*exp(a*t).*sin(w0*t+phi); plot(t,ft1)
分析LTI系统可采用时域方法或变换域方法，如 傅立叶变换、拉普Leabharlann Baidu斯变换和Z变换。 LTI系统可分为连续时间系统和离散时间系统 描述连续LTI系统方法：用常系数微分方程、系 统的传递函数或状态方程。
一.利用MATLAB进行系统的时域分析
连续时间系统零状态响应的求解
连续时间系统冲激响应和阶跃响应的求解
y(t)
1
2
Time(sec)
3
4
5
例2 求系统 y" (t)+2y ' (t)+100y(t)=10x(t) 的零状 态响应，已知x(t) =d (t) 。
%连续时间系统的冲激响应 ts=0;te=5;dt=0.01; 1 sys=tf([10],[1 2 100]); 0.5 t=ts:dt:te; y=impulse(sys,t); 0 plot(t,y); -0.5 xlabel('Time(sec)') ylabel('h(t)') -1
信号与系统matlab实验及答案
产生离散衰减正弦序列()π0.8sin 4n x n n ⎛⎫
= ⎪⎝⎭
, 010n ≤≤，并画出其波形图。
n=0:10;
x=sin(pi/4*n).*0.8.^n;
stem(n,x);xlabel( 'n' );ylabel( 'x(n)' );
用MATLAB 生成信号()0sinc at t -, a 和0t 都是实数,410t -<<，画波形图。观察并分析a 和0t 的变化对波形的影响。
t=linspace(-4,7); a=1;
t0=2;
y=sinc(a*t-t0); plot(t,y);
t=linspace(-4,7); a=2;
t0=2;
y=sinc(a*t-t0); plot(t,y);
t=linspace(-4,7); a=1;
t0=2;
y=sinc(a*t-t0); plot(t,y);
三组对比可得a 越大最大值越小，t0越大图像对称轴越往右移
某频率为f 的正弦波可表示为()()cos 2πa x t ft =，对其进行等间隔抽样，得到的离散样值序列可表示为()()a t nT x n x t ==，其中T 称为抽样间隔，代表相邻样值间的时间间隔，1
s f T
=
表示抽样频率，即单位时间内抽取样值的个数。抽样频率取40 Hz s f =，信号频率f 分别取5Hz, 10Hz, 20Hz 和30Hz 。请在同一张图中同时画出连续信号()a x t t 和序列()x n nT 的波形图，并观察和对比分析样值序列的变化。可能用到的函数为plot, stem, hold on 。
matlab信号与系统课程设计
matlab信号与系统课程设计
一、引言
信号与系统是电子信息类专业中的一门重要课程，Matlab作为一种常用的计算工具，可以帮助学生更好地理解和应用信号与系统的相关知识。本文将介绍如何使用Matlab进行信号与系统课程设计。
二、实验目的
本次实验旨在通过使用Matlab软件，帮助学生深入理解信号与系统的相关知识，并掌握Matlab软件在信号与系统中的应用。
三、实验内容
本次实验分为两个部分：信号处理和系统分析。
1. 信号处理
1.1 生成离散时间序列信号
使用Matlab中的discrete函数生成一个离散时间序列信号。要求该信号包含10个采样点，采样频率为100Hz，幅度随机取值。
1.2 时域分析
对生成的离散时间序列信号进行时域分析。计算出该信号的均值、方差、标准差和自相关函数，并画出该信号及其自相关函数的图像。1.3 频域分析
对生成的离散时间序列信号进行频域分析。计算出该信号的功率谱密度，并画出该功率谱密度函数图像。
2. 系统分析
2.1 系统建模
使用Matlab中的tf函数建立一个一阶低通滤波器系统模型。该系统
的传递函数为H(s)=1/(s+1)。
2.2 系统分析
对建立的一阶低通滤波器系统进行分析。计算出该系统的单位脉冲响应、单位阶跃响应和零极点图，并画出相应的图像。
四、实验步骤
4.1 生成离散时间序列信号
首先，打开Matlab软件，新建一个m文件，命名为
“signal_processing.m”。
在m文件中输入以下代码：
t = 0:0.01:0.09;
x = rand(1,10);
stem(t,x);
信号与系统及MATLAB实现
《信号与系统及MATLAB实现》实验指导书
前言
长期以来，《信号与系统》课程一直采用单一理论教学方式，同学们依靠做习题来巩固和理解教学内容，虽然手工演算训练了计算能力和思维方法，但是由于本课程数学公式推导较多，概念抽象，常需画各种波形，作题时难免花费很多时间，现在，我们给同学们介绍一种国际上公认的优秀科技应用软件MATLAB，借助它我们可以在电脑上轻松地完成许多习题的演算和波形的绘制。
MATLAB的功能非常强大，我们此处仅用到它的一部分，在后续课程中我们还会用到它，在未来地科学研究和工程设计中有可能继续用它，所以有兴趣的同学，可以对MATLAB 再多了解一些。
MATLAB究竟有那些特点呢？
1．高效的数值计算和符号计算功能，使我们从繁杂的数学运算分析中解脱出来；
2．完备的图形处理功能，实现计算结果和编程的可视化；
3．友好的用户界面及接近数学表达式的自然化语言，易于学习和掌握；
4．功能丰富的应用工具箱，为我们提供了大量方便实用的处理工具；
MATLAB的这些特点，深受大家欢迎，由于个人电脑地普及，目前许多学校已将它做为本科生必须掌握的一种软件。正是基于这些背景，我们编写了这本《信号与系统及MATLAB实现》指导书，内容包括信号的MA TLAB表示、基本运算、系统的时域分析、频域分析、S域分析、状态变量分析等。通过这些练习，同学们在学习《信号与系统》的同时，掌握MATLAB的基本应用，学会应用MATLAB的数值计算和符号计算功能，摆脱烦琐的数学运算，从而更注重于信号与系统的基本分析方法和应用的理解与思考，将课程的重点、难点及部分习题用MATLAB进行形象、直观的可视化计算机模拟与仿真实现，加深对信号与系统的基本原理、方法及应用的理解，为学习后续课程打好基础。另外同学们在进行实验时，最好事先预习一些MATLAB的有关知识，以便更好地完成实验，同时实验中也可利用MATLAB的help命令了解具体语句以及指令的使用方法。
信号与系统MATL实验及代码
实验一、MATLAB编程基础及典型实例
一、实验目的
（1）熟悉MATLAB软件平台的使用；
（2）熟悉MATLAB编程方法及常用语句；
（3）掌握MATLAB的可视化绘图技术；
（4）结合《信号与系统》的特点，编程实现常用信号及其运算。
示例一：在两个信号进行加、减、相乘运算时，参于运算的两个向量要有相同的维数，并且它们的时间变量范围要相同，即要对齐。编制一个函数型m文件，实现这个功能。
function [f1_new,f2_new,n]=duiqi(f1,n1,f2,n2)
a=min(min(n1),min(n2));
b=max(max(n1),max(n2));
n=a:b;
f1_new=zeros(1,length(n));
f2_new=zeros(1,length(n));
tem1=find((n>=min(n1))&(n<=max(n1))==1);
f1_new(tem1)=f1;
tem2=find((n>=min(n2))&(n<=max(n2))==1);
f2_new(tem2)=f2;
四、实验内容与步骤
− 2 t （2）绘制信号x(t)= e
sin( t=0:0.1:30; 2
3
t ) 的曲线，t的范围在0~30s，取样时间间隔为0.1s。
y=exp(-sqrt(2)*t).*sin(2*t/3); plot(t,y);
（3）在n=[-10:10]范围产生离散序列：x
( n)
=
⎧2 n ,
− 3
⎨
≤ n ≤ 3
,并绘图。
⎩0, Other
n=-10:1:10;
信号与系统MATLAB实验课后习题答案
读书破万卷下笔如有神
3-1
a=[1,1,1];
b=[1,1];
sys=tf(b,a);
t=[0:0.01:10];
figure;
subplot(2,2,1);
step(sys);
subplot(2,2,2);
x_step=zeros(size(t));
x_step(t>0)=1;
x_step(t==0)=1/2;
lsim(sys,x_step,t);
subplot(2,2,3);
impulse(sys,t);
title('Impulse Response');
xlabel('Time(sec)');
ylabel('Amplitude');
subplot(2,2,4);
x_delta=zeros(size(t));
x_delta(t==0)=100;
[y1,t]=lsim(sys,x_delta,t);
y2=y1;
plot(t,y2);
title('Impulse Response');
xlabel('Time(sec)');
ylabel('Amplitude');
Step ResponseLinear Simulation Results1.51.5
11eedduuttiillppmm0.50.5AA0001005510Time (sec)Time (sec)Impulse Response1Impulse Response
1
0.5ed0.5u ed ti u lp til m p0m A0A-0.5-0.510500510Time(sec)Time(sec) (sec)
读书破万卷下笔如有神
(完整版)信号与系统Matlab实验作业
实验一 典型连续时间信号和离散时间信号
一、实验目的
掌握利用Matlab 画图函数和符号函数显示典型连续时间信号波形、典型时间离散信号、连续时间信号在时域中的自变量变换。
二、实验内容
1、典型连续信号的波形表示（单边指数信号、复指数信号、抽样信号、单位阶跃信号、单位冲击信号）
1）画出教材P28习题1-1(3) ()[(63)(63)]t f t e u t u t =----的波形图。 function y=u(t) y=t>=0;
t=-3:0.01:3;
f='exp(t)*(u(6-3*t)-u(-6-3*t))';
ezplot(f,t);
grid on;
2）画出复指数信号()()j t f t e σω+=当0.4, 8σω==(0
波形图。 t=0:0.01:10;
f1='exp(0.4*t)*cos(8*t)';
f2='exp(0.4*t)*sin(8*t)';
figure(1)
ezplot(f1,t);
grid on;
figure(2)
ezplot(f2,t);
grid on;
t=-10:0.01:10; f='sin(t)/t'; ezplot(f,t); grid on;
t=0:0.01:10;
f='(sign(t-3)+1)/2'; ezplot(f,t);
grid on;
5）单位冲击信号可看作是宽度为∆，幅度为1/∆的矩形脉冲，即t=t 1处的冲击信号为
11111 ()()0 t t t x t t t other
δ∆⎧<<+∆⎪=-=∆⎨⎪⎩
画出0.2∆=, t 1=1的单位冲击信号。 t=0:0.01:2;
信号与系统实验__matlab
信号与系统部分
01.分别用MATLAB 表示并绘出下列连续时间信号的波形：
2()(2)()t f t e u t -=-
02.分别用MATLAB 表示并绘出下列连续时间信号的波形：
[]()cos(
)()(4)2
t
f t u t u t π=--
03.分别用MATLAB 表示并绘出下列离散时间信号的波形：
()12()()k
f k u k =-
04.分别用MATLAB 表示并绘出下列离散时间信号的波形：
[]()()(8)f t k u k u k =--
05.已知信号f (t)的波形如下图所示，试用MATLAB 绘出满足下列要求的信号波形。
()f t -
06.已知信号f (t)的波形如下图所示，试用MATLAB 绘出满足下列要求的信号波形。
()f at
a ＝0.5
07.已知信号f (t)的波形如下图所示，试用MATLAB 绘出满足下列要求的信号波形。
()f at ，
a ＝2
08.已知信号f (t)的波形如下图所示，试用MATLAB 绘出满足下列要求的
信号波形。
(0.51)f t +
09.已知两信号1()(1)()f t u t u t =+-，2()()(1)f t u t u t =--，求卷积积分12()()()g t f t f t =*。
10.已知两信号1()()f t tu t =，20()()0
t t
t te u t f t t e
-≥⎧=⎨<⎩，求卷积积分12()()()g t f t f t =*。
11.已知{}{}12()1,1,1,2,()1,2,3,4,5f k f k ==，求两序列的卷积和。
《信号与系统》MATLAB仿真实验讲义
《信号与系统》MATLAB仿真实验讲义
（第二版）
肖尚辉编写
宜宾学院电信系电子信息教研室《信号与系统》课程
2004年3月 宜宾
使用对象：电子专业02级3/4班（本科）
实验一 产生信号波形的仿真实验
一、实验目的：熟悉MATLAB软件的使用，并学会信号的表示和以及用MATLAB
来产生信号并实现信号的可视化。
二、实验时数：3学时+3学时（即两次实验内容）
三、实验内容：
信号按照自变量的取值是否连续可分为连续时间信号和离散时间信号。对信号进行时域分析，首先需要将信号随时间变化的规律用二维曲线表示出来。对于简单信号可以通过手工绘制其波形，但对于复杂的信号，手工绘制信号波形显得十分困难，且难以绘制精确的曲线。
在MATLAB中通常用三种方法来产生并表示信号，即（1）用MATLAB软件的funtool符合计算方法（图示化函数计算器）来产生并表示信号；（2）用MATLAB软件的信号处理工具箱（Signal Processing Toolbox）来产生并表示信号；（3）用MATLAB软件的仿真工具箱Simulink中的信号源模块。
（一） 用MATLAB软件的funtool符合计算方法（图示化函数计算器）来产生并表示信号
在MATLAB环境下输入指令funtool，则回产生三个视窗。即
figure No.1：可轮流激活，显示figure No.3的计算结果。
figure No.2：可轮流激活，显示figure No.3的计算结果。
figure No.3：函数运算器，其功能有：f，g可输入函数表达式；x是自变量，在缺省时在[－2pi，2pi]的范围内；自由参数是a；在分别输入完毕后，按下面四排的任一运算操作键，则可在figure No.1或figure No.2产生相应的波形。
基于MATLAB的信号与系统仿真实验毕业设计(含源文件)
第 5 章 系统主要功能............................................................................................................................ 14 5.1 连续时间信号时域分析.......................................................................................................... 14 5.2 连续系统时域分析.................................................................................................................. 20 5.3 离散系统时域分析.................................................................................................................. 22 5.4 傅里叶变换.............................................................................................................................. 24
信号与系统MATLAB仿真——信号及其运算
信号与系统MATLAB仿真——信号及其运算1. 知识回顾
（1）信号的分类：
确定信号与随机信号；
周期信号与⾮周期信号；
周期信号在时间上必须是⽆始⽆终的
f(t)=f(t+T)
f[k]=f[k+N]
连续时间信号和离散时间信号；
连续信号是指在信号的定义域内，除若⼲个第⼀类间断点外，对于任意时刻都由确定的函数值的信号
离散信号是指在信号的定义域内，只在某些不连续规定的时刻给出函数值，⽽在其他时刻没有给出函数的信号能量信号、功率信号与⾮功率⾮能量信号；
时限与频限信号；
物理可实现信号。
（2）信号能量：
E=lim
T→∞∫T−T f2(t)dt
P=lim
T→∞
1
2T∫T
−T
f2(t)dt
E=
lim
N→∞
N
∑
k=−N|f[k]|2
P=
lim
N→∞
1
2N+1
N
∑
k=−N|f[k]|2
能量信号：0<E<∞，P=0；
功率信号：0<P<∞，E=∞。
（3）冲激函数的性质
加权特性（筛选特性）：
f(t)δ(t−t0)=f(t0)δ(t−t0)取样特性：
∫+∞−∞f(t)δ(t−t0)=f(t0)偶函数：
f(t)=f(−t)展缩特性：
δ(at)=1
|a|δ(t)
δ(at−t0)=1
|a|δ(t−
t0
a)
导数及其特性。
（4）正弦
两个频率相同的正弦信号相加，即使其振幅和相位各不相同，但相加后结果仍是原频率的正弦信号；
若⼀个正弦信号的频率是另⼀个正弦信号频率的整数倍时，则合成信号是⼀个⾮正弦周期信号，其周期等于基波的周期。正弦型序列：
f[k]=A sin(Ω0k+φ)
2π/Ω0是正整数：周期序列，周期为N；
信号与系统分析图文 (5)
第5章 离散信号与系统的时域分析 图 5-5 anU(n)的波形
第5章 离散信号与系统的时域分析
5.1.3 典型离散信号的MATLAB
在MATLAB中，离散信号用一个行向量或一个列向 量表示。在MATLAB中向量是从1开始编导的，即x(1)是 x向量的第1个元素。在表示信号或信号运算时，如果这 些编号与所需要的信号标号不能对应，可以创建另外一 个标号向量， 使信号的标号与实际情况一致。MATLAB
第5章 离散信号与系统的时域分析 3. 图形形式 图形形式即信号的波形。例如上面f1(n)、f3(n)分别 如图5-1(a)、(b)
图 5-1 离散信号的波形
第5章 离散信号与系统的时域分析
5.1.2
1. 单位样值(Unit Sample)信号δ(n) (5-1)
δ(n)的波形如图5-2(a)所示。
第5章 离散信号与系统的时域分析 若令相继时刻tn与tn+1之间的间隔为T, 则离散信号只在 均匀离散时刻t=…, -2T, -T, 0, T, 2T, …时有定义，它可以表 示为f(nT)。为了方便，不妨把f(nT)简记为f(n), 这样的离散
一个离散时间信号f(n)可以用三种方法来描述。
第5章 离散信号与系统的时域分析
第5章 离散信号与系统的时域分析 图 5-9 斜变序列
第5章 离散信号与系统的时域分析
5.2 离散信号的基本运算及MATLAB实现
电子教案-信号与系统第四版(含习题解答)-信号与系统的MATLAB仿真
信号与系统的MATLAB 仿真
一、信号生成与运算的实现
1.1 实现)3
(sin )()(
π±==
=t t
t
t S t f a )(sin )sin()sin(sin )()(t c t t t t t t t S t f a '='
'==
==πππ
π
ππ m11.m
t=-3*pi:0.01*pi:3*pi; % 定义时间范围向量t f=sinc(t/pi); % 计算Sa(t)函数 plot(t,f); % 绘制Sa(t)的波形 运行结果：
1.2 实现)10()
sin()(sin )(±==
=t t
t t c t f ππ m12.m
t=-10:0.01:10; % 定义时间范围向量t f=sinc(t); % 计算sinc(t)函数
plot(t,f); % 绘制sinc(t)的波形 运行结果：
1.3 信号相加：t t t f ππ20cos 18cos )(+=
m13.m
syms t; % 定义符号变量t
f=cos(18*pi*t)+cos(20*pi*t); % 计算符号函数f(t)=cos(18*pi*t)+cos(20*pi*t) ezplot(f,[0 pi]); % 绘制f(t)的波形 运行结果：
1.4 信号的调制：t t t f ππ50cos )4sin 22()(+=
m14.m
syms t; % 定义符号变量t
f=(2+2*sin(4*pi*t))*cos(50*pi*t) % 计算符号函数f(t)=(2+2*sin(4*pi*t))*cos(50*pi*t) ezplot(f,[0 pi]); % 绘制f(t)的波形 运行结果：
Matlab 信号与系统
%integration t= -3:0.1:3; for x=1:length(t)
y2(x)=quad('f2_2', -3,t(x)); end plot(t,y2)
35
三角波f(t)微分与积分的波形
1
0.9
0.8
0.7
0.6
0.5
0.4
0.3
0.2
0.1
0
-3
-2
-1
0
1
2
3
0.
dx(t)/dt
信号基本运算的MATLAB实现
尺度变换、翻转、时移、 相加、相乘、 差分与求和、微分与积分
24
一、基本信号的MATLAB表示
➢ 指数信号Aeat ➢ 指数序列ak ➢ 正弦型信号 ➢ 抽样函数Sa(t) ➢ 矩形脉冲信号 ➢ 三角波脉冲信号
y = A*exp(a*t); 幂运算a.^k实现 内部函数cos( ) 和sin( ) sinc(t) y = rectpuls(t,width) y = tripuls(t, width,skew)
1
0.9
0.8
0.7
0.6
0.5
0.4
0.3
0.2
0.1
0
-50 -40 -30 -20 -10
0
10 20 30 40 50
信号与系统——Matlab部分习题【范本模板】
课程设计(论文）
课程名称:信号与系统
课程论文题目：信号与系统——MATLAB课程设计姓名:
系：
专业：物理学（电子信息工程方向）年级：大学二年级
学号：
指导教师：
职称：
2017年 5 月22日
目录
摘要---------------------------------------------------------------------- 1关键词-------------------------------------------------------------------- 1
一、M2—7 -------------------------------------------------------------- 2
1、问题重述--------------------------------------------------------------- 2
2、问题分析--------------------------------------------------------------- 2
3、仿真程序与仿真结果----------------------------------------------------- 2
（1）仿真程序---------------------------------------------------------------------------------------------------------- 2频率为262Hz的正弦信号声音波形 ------------------------------------------------------------------------------- 2频率为294Hz的正弦信号声音波形 ------------------------------------------------------------------------------- 3频率为330Hz的正弦信号声音波形 ------------------------------------------------------------------------------- 3频率为349Hz的正弦信号声音波形 ------------------------------------------------------------------------------- 4频率为392Hz的正弦信号声音波形 ------------------------------------------------------------------------------- 4频率为440Hz的正弦信号声音波形 ------------------------------------------------------------------------------- 5频率为494Hz的正弦信号声音波形 ------------------------------------------------------------------------------- 5频率为524Hz的正弦信号声音波形 ------------------------------------------------------------------------------- 6结果分析：--------------------------------------------------------------------------------------------------------------- 6 (2）仿真程序：--------------------------------------------------------------------------------------------------------- 6仿真结果：--------------------------------------------------------------------------------------------------------------- 7结果分析：--------------------------------------------------------------------------------------------------------------- 7 (3）仿真程序：--------------------------------------------------------------------------------------------------------- 7仿真结果：--------------------------------------------------------------------------------------------------------------- 7结果分析：--------------------------------------------------------------------------------------------------------------- 8
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
t=-4:0.001:4; ft=tripuls(t,4,0.5); % 原始信号 subplot(2,2,1); plot(t,ft); title('x(t)') ft1=tripuls(2*t,4,0.5); % 压缩 subplot(2,2,2); plot(t,ft1); title('x(2t)') ft2=tripuls(-t,4,0.5); % 翻转 subplot(2,2,3); plot(t,ft2); title('x(-t)') ft3=tripuls(t+1,4,0.5); % 平移 subplot(2,2,4); plot(t,ft3); title('x(t+1)')
信号的MATLAB表示
基本信号的MATLAB表示
指数信号Aeat 、指数序列ak 、抽样函数Sa(t)、
正弦型信号、矩形脉冲信号、三角脉冲信号
信号基本运算的MATLAB实现
尺度变换、翻转、时移、 相加、相乘、 差分与求和、微分与积分
精品课件
一、基本信号的MATLAB表示
➢ 指数信号Aeat
y = A*exp(a*t);
0
0.5
1
1.5
2
2.5
3
3.5
4
精品课件
一、基本信号的MATLAB表示
% tripuls t=-3:0.001:3; ft=tripuls(t,4,0.5); plot(t,ft)
ft=tripuls(t,4,1);
精品课件
1
0.9
0.8
0.7
0.6
0.5
0.4
0.3
0.2
0.1
0
-3
-2
-1
0
%decaying exponential 1 0.9
t=0:001:10;
0.8
0.7
A=1;
0.6
0.5
a=-0.4;
0.4
0.3
ft=A*exp(a*t);
0.2
0.1
plot(t,ft)
0
0
1
2
3
4
5
6
7
8
9
10
t=0:0.1:10; A=1; a=-0.4; ft=A*exp(a*t); stem(t,f精t品)课件
-50 -40 -30 -20 -10
0
10 20 30 40 50
%产生 f[k]=u(k-k0);k1<=k<=k2
k=[k1:k2];f=[(k-k0)>=0];
k0=0;k1=-50;k2=50;
[f,k]=stepseq(k0,k1,k2);
stem(k,f) 精品课件
二、信号基本运算的MATLAB实现
精品课件
x(t) 1
x(2t) 1
0.8
0.8
0.6
0.6
0.4
0.4
0.2
0.2
0
-4
-2
0
2
4
0
-4
-2
0
2
4
x(-t) 1
x(t+1) 1
0.8
0.8
0.6
0.6
0.4
0.4
0.2
0.2
0
0
-4
-2
0
2
精4品课件
-4
-2
0
2
4
二、信号基本运算的MATLAB实现
2. 信号的相加与相乘
1
0.8
stem(k,f) 精品课件
一、基本信号的MATLAB表示
% unit step sequence
1
0.9
k=-50:50;
0.8
0.7
uk=[zeros(1,50), ones(1,51)]; 0.6
stem(k,uk)
0.5 0.4
0.3
0.2
0.1
function [f,k]=stepseq(k0,k1,k2) 0
✓ 相加用算术运算符“+”实现 0.6
0.4
✓ 相乘用数组运算符“.*”实现 0.2
0
例:画信号Aeatcos(w0t+f)的波形
-0.2 -0.4
t=0:0.001:8;
-0.6
A=1; a=-0.4;
-0.8
0
1
2
3
4
5
6
7
8
w0=2*pi;phi=0;
ft1=A*exp(a*t).*sin(w0*t+phi);
0.3
0.2
0.1
function [f,k]=impseq(k0,k1,k2) 0
-50 -40 -30 -20 -10
0
10 20 30 40 50
%产生 f[k]=delta(k-k0)；k1<=k<=k2
k=[k1:k2];f=[(k-k0)==0];
k0=0;k1=-50;k2=50;
[f,k]=impseq(k0,k1,k2);
1
2
3
1
0.9
0.8
0.7
0.6
0.5
0.4
0.3
0.2
0.1
0
-3
-2
-1
0
1
2
3
一、基本信号的MATLAB表示
% unit impuls sequence
1 0.9
k=-50:50;
0.8 0.7
delta=[zeros(1,50),1,zeros(1,50)00]..56;
0.4wenku.baidu.com
stem(k,delta)
for x=1:length(t) % 积分
y2(x)=integral(@(t)
tripuls(t,4,0.5),t(1),t(x));
plot(t,ft1) 精品课件
二、信号基本运算的MATLAB实现
3. 离散序列的差分与求和
➢ 差分
y=diff(f );
y = [f(2)-f(1) f(3)-f(2) ... f(m)-f(m-1)]
➢ 求和
y=sum(f(k1:k2));
y = f(k1)+f(k1+1)+..+f(k2)
连续信号的微分与积分
1
0.9
0.8
0.7
0.6
0.5
0.4
0.3
0.2
0.1
0
0
1
2
3
4
5
6
7
8
9
10
一、基本信号的MATLAB表示
% rectpuls
t=0:0.001:4; T=1; ft=rectpuls(t-2*T,T); plot(t,ft) axis([0,4,-0.5,1.5])
1.5
1
0.5
0
-0.5
➢ 指数序列ak
幂运算a.^k实现
➢ 正弦型信号
内部函数cos( ) 和sin( )
➢ 抽样函数Sa(t)
sinc(t)
➢ 矩形脉冲信号 y = rectpuls(t,width)
➢ 三角波脉冲信号 y = tripuls(t,width,skew) %skew:斜
度
精品课件
一、基本信号的MATLAB表示
➢ 微分 y=diff(f)/h; h为数值计算所取时间间隔 ➢ 定积分 integral(function_handle,a,b); function_handle被积函数句柄，a和b定积分区间。
精品课件
二、信号基本运算的MATLAB实现
例:已知三角波x(t)，画出其微分与积分的波
形
% 原始信号 h=0.01; t=-4:h:4; ft=tripuls(t,4,0.5); subplot(2,1,1);plot(t,ft); title('x(t)');grid on