武汉大学信号与系统仿真实验报告
- 格式:doc
- 大小:578.59 KB
- 文档页数:37
实验一:连续时间信号的表示及可视化
一、实验名称:连续时间信号的表示及可视化
()(),()(),()at f t t f t t f t e δε===(分别取a>0和a<0)
; ()(),()(),()(2)f t R t f t Sa t f t Sin ft ωπ===(分别画出不同周期个数的波
形)
二、实验目的:
(1)掌握应用matlab 绘制连续时间信号图的基本方法 (2)复习《信号与系统》课程中有关连续时间信号的相关知识
(3)通过观察实验结果,对几个典型的连续时间信号图形有直观的理解 三、解题分析:
连续时间函数与离散时间函数在编程中的区别主要体现在如下两个方面:第一,自变量的取值范围不同,离散时间函数的自变量是整数,而连续时间函数的自变量为一定范围内的实数;第二,绘图所用的函数不同,连续函数图形的绘制不止一个,下面将以fplot 函数为主进行编程和绘图。 四、实验程序:
1
)()(t t f δ=
t=-1000:0.1:1000;
y=dirac(t); plot(y,t);
2
)()(t t f ε=
t=-1:0.1:10; y=heaviside(t); plot(t,y)
3
at e t f =)(
t=-1:0.1:2; a=1;
y=exp(a*t); plot(t,y);
4
)()(t R t f =
t=0:0.1:10;
y=heaviside(t)-heaviside(t-5); plot(t,y)
5
)()(t Sa t f ω=
t=-100:0.1:100; y=sin(t)./t; plot(t,y)
6
)2()(ft Sin t f π=
t=-1:0.01:1; f=3;
y=sin(2*pi*f*t); plot(t,y)
五、结果分析及实验体会
结果分析:
连续函数的自变量t 为一定范围内的连续值,函数波形图为连续不间断的。
实验二:离散时间信号的表示及可视化
一、实验名称:离散时间信号的表示及可视化
()(),()(),()an f n n f n n f n e δε===(分别取a<0和a>0)
()()N f n R n =(分别取不同的N 值),()(),()()f n Sa n f n Sin n ωω==(分别取不同的ω值) 二、实验目的:
(1)掌握应用matlab 绘制离散时间信号图的基本方法 (2)复习《信号与系统》课程中有关离散时间信号的相关知识 (3)通过观察实验结果,对几个典型的离散时间信号图形有直观的理解 三、解题分析
本实验中要求绘制离散时间信号图,可以应用matlab 中的stem 函数来实现。用matlab 表示一离散序列x[k]时,可用两个向量来表示。其中一个向量表示自变量k 的取值范围,另一个向量表示序列x[k]的值。之后可用stem(k,f)画出序列波形。当序列是从k=0开始时,可以只用一个向量x 来表示序列。由于计算机内寸的限制,matlab 无法表示一个无穷长的序列。对于典型的离散时间信号,可用逻辑表达式来实现不同自变量时的取值。 四、实验程序
1
)()(n n f δ=
n=-10:10; y=[n==0]; stem(n,y);
2
)()(n n f ε=
n=-10:10; y=[n>=0]; stem(y);
3
an
e
n
f=
)
((分别取0
0<
>a
a及)
n=0:8;
y=exp(0.5*n); stem(y);
grid on;
n=0:8;
y=exp(-0.5*n); stem(y);
grid on;
4
)()(n R n f N (分别取不同的N 值)
N=5
n=0:10; y=[n<=5]; stem(y); grid on ;
N=3
n=0:10; y=[n<=3]; stem(y); grid on ;
5
)()(ωn Sa n f =
ω=50
n=0:30;
y=sin(n*50)./(n*50); stem(y); grid on ;
ω=100
n=0:30;
y=sin(n*100)./(n*100); stem(y); grid on ;
6
)()(ωn Sin n f =
ω=100
n=0:20;
y=sin(n*100); stem(y); grid on ;
ω=50
n=0:20;
y=sin(n*50); stem(y); grid on ;
五、结果分析及实验体会:
结果分析:
通过与上次实验进行对比,可明显地看出离散函数和连续函数的区别。离散时间信号只在n取整数的时候才取值,表现在图上就是一系列不连续的点。由(5)和(6)能明显地看出函数值的轮廓,可见他们是连续函数的抽样取值。对于()()
=,当ω取不同值是,
f n Sin nω
函数在图形上表现的周期就不同。同样的n取值范围内,在ω取0.64时,其波形比ω取0.32时多出一个波形。
实验三:系统的时域求解
一、实验名称:系统的时域求解
(1)设()(0.9)(),()()(10)n
hn un xn un un ===
-
,求:()()*()y n x n h n =,
并画出()x n 、()h n 、()y n 的波形。
(2)求因果线性移不变系统()0.81(2)()(2)y n y n x n x n =-+--的单位抽样响应()h n ,并绘出()j H e ω的幅频及相频特性曲线。 二、实验目的:
(1)对《信号与系统》课程中有关时域求解的相关知识进行回忆巩固。
(2)学习用matlab 对时域系统求解的方法。 (3)加深对卷积求解、幅频特性和相频特性的理解。 三、解题分析:
离散卷积是数字信号处理中的一个基本运算,MTLAB 提供的计算两个离散序列卷积的函数是conv ,其调用方式为 y=conv(x,h) 。其中调用参数x,h 为卷积运算所需的两个序列,返回值y 是卷积结果。
matlab 函数conv 的返回值y 中只有卷积的结果,没有y 的取值范围。由离散序列卷积的性质可知,当序列x 和h 的起始点都为k=0时,y 的取值范围为k=0至length(x)+length(h)-2。
许多离散LTI 都可用如下的线性常系数的差分方程描述
[][]N N
n
n
n n a y k n b x k n ==-=-∑∑
其中x[k]、y[k]分别系统的输入和输出。在已知差分方程的N 个初始状态y[k],和输入x[k],就可由下式迭代计算出系统的输出