信号分析与处理实验报告(基于MATLAB)资料

  • 格式:doc
  • 大小:712.82 KB
  • 文档页数:30

下载文档原格式

  / 30
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

武汉工程大学电气信息学院

title('幅度')

xlabel('t')

axis([-0.5,20.5,0,1.1])

grid on

subplot(2,2,4) %将当前窗口分成2行2列个子窗口,并在第4个子窗口绘图plot(t,fg)

title('相位')

xlabel('t')

axis([-0.5,20.5,-3.5,3.5])

grid on

三、实验数据与结果分析

1、

2、

四、思考:

1、为什么图二中t=0处曲线是间断的,如何使其成为连续的曲线?

因为axis函数对纵坐标的的上边界限定过小,使图形在边界处不能完整的显示。

xlabel('t')

axis([-7,7,-1,2])

grid on

三、实验数据与结果分析

1.

2.

3.

四、思考:

1、代数运算符号*和.*的区别是?

*是矩阵相乘,是矩阵A行元素与B的列元素相乘的和

.*是数组相乘,表示数组A和数组B中的对应元素相乘

实 验 内 容

实验三 连续时间信号的卷积

一、实验内容

1、已知两连续时间信号如下图所示,绘制信号f 1(t)、f 2(t)及卷积结果f(t)的波形;设时间变化步长dt 分别取为0.5、0.1、0.01,当dt 取多少时,程序的计算结果就是连续时间卷积的较好近似?

2、、计算信号

()()()11==-a t u e t f at 和()()t tu t f sin 2=的卷积f(t),f 1(t)、f 2(t)的时间

范围取为0~10,步长值取为0.1。绘制三个信号的波形。

二、实验方法与步骤

1、绘制信号f 1(t)、f 2(t)及卷积结果f(t)的波形,当dt 取0.01时程序的计算结果就是连续时间卷积的较好近似 程序代码如下:

clear all close all clc dt=0.01 t1=0:dt:2; t2=-1:dt:1;

三、实验数据与结果分析

1.

2.

三、实验数据与结果分析1.

2.

实 验 内 容

实验五 连续时间信号的频域分析

一、实验内容

1、如图5.4所示的奇谐周期方波信号,周期为T1=1,幅度为A=1,将该方波信号展开成三角形式Fourier 级数并分别采用频域矩形窗和Hanning 窗加权,绘制两种窗函数加权后的方波合成图像。时间范围取为-2~2,步长值取为0.01。

2、将图5.5中的锯齿波展开为三角形式Fourier 级数,按(2)式求出Fourier 级数的系数,并在频域分别采用矩形窗、Hanning 窗和三角窗加权,观察其Gibbs 效应及其消除情况。时间范围取为-2~2,步长值取为0.01。

3、选做:编程计算连续时间周期信号的三角形式傅里叶级数展开的系数

二、实验方法与步骤

1、将方波信号展开成三角形式Fourier 级数并分别采用频域矩形窗和Hanning 窗加权 方波展开的三角式傅立叶级数为:

()()t k k t x L

k 1,5,3,1sin 4ωπ

⋅∑

=∞

=

采用频域矩形窗加权,则展开式变为:

()()()[]t k k t x K

k 10

12sin 124

ωπ

+⋅+∑==

a0=2/T*int(f,t,0,T); %求函数f对t从0到T的定积分

a0=simplify(a0) %得出结果

syms k

fa=t*cos(k*w*t);

fb=t*sin(k*w*t);

ak=2/T*int(fa,t,0,T); %求函数fa对t从0到T的定积分bk=2/T*int(fb,t,0,T); %求函数fb对t从0到T的定积分ak=simplify(ak)

bk=simplify(bk)

三、实验数据与结果分析

1.

2.

3.

根据绘制的幅频特性曲线,系统具有低通滤波特性2.

根据绘制的幅频特性曲线,系统具有带通滤波特性