文档视界 最新最全的文档下载
当前位置:文档视界 › 基于Matlab的信号与系统实验指导-2012新版

基于Matlab的信号与系统实验指导-2012新版

基于Matlab 的信号和系统实验指导

实验一 连续时间信号在Matlab 中的表示

一、实验目的

1、学会运用Matlab 表示常用连续时间信号的方法

2、观察并熟悉这些信号的波形和特性 二、实验原理及实例分析

1、信号的定义和分类

2、如何表示连续信号?

连续信号的表示方法有两种;符号推理法和数值法。

从严格意义上讲,Matlab 数值计算的方法不能处理连续时间信号。然而,可利用连续信号在等时间间隔点的取样值来近似表示连续信号,即当取样时间间隔足够小时,这些离散样值能被Matlab 处理,并且能较好地近似表示连续信号。

3、Matlab 提供了大量生成基本信号的函数。如: (1)指数信号:K*exp(a*t)

(2)正弦信号:K*sin(w*t+phi)和K*cos(w*t+phi) (3)复指数信号:K*exp((a+i*b)*t) (4)抽样信号:sinc(t*pi)

注意:在Matlab 中用和Sa(t)类似的sinc(t)函数表示,定义为:sinc(t)sin(t)/(t)ππ= (5)矩形脉冲信号:rectpuls(t,width)

(6)周期矩形脉冲信号:square(t,DUTY),其中DUTY 参数表示信号的占空比DUTY%,即在一个周期脉冲宽度(正值部分)和脉冲周期的比值。占空比默认为0.5。

(7)三角波脉冲信号:tripuls(t, width, skew),其中skew 取值范围在-1~+1之间。 (8)周期三角波信号:sawtooth(t, width)

(9)单位阶跃信号:y=(t>=0) 三、实验内容

1、验证实验内容 直流及上述9个信号

2、程序设计实验内容

(1)利用Matlab 命令画出下列连续信号的波形图。 (a ))4/3t (2cos π+

(b )

)t (u )e 2(t

-- (c ))]2()(u )][t (cos 1[--+t u t π

(2)利用Matlab 命令画出复信号j(t /4)f (t)2e π+=的实部、虚部、模和辐角。 四、实验报告要求

1、格式:实验名称、实验目的、实验原理、实验环境、实验内容、实验思考等

2、实验内容:程序设计实验部分源代码及运行结果图示。

实验二 连续时间信号在Matlab 中的运算

一、实验目的

1、学会运用Matlab 进行连续时间信号的时移、反褶和尺度变换;

2、学会运用Matlab 进行连续时间信号微分、积分运算;

3、学会运用Matlab 进行连续时间信号相加、相乘运算;

4、学会运用Matlab 进行连续时间信号卷积运算。

二、实验原理及实例分析

1、信号的时移、反褶和尺度变换

信号的平移、反转和尺度变换是针对自变量时间而言的,其数学表达式和波形变换中存在着一定的变化规律。从数学表达式上来看,信号的上述所有计算都是自变量的替换过程。所以在使用Matlab 进行连续时间信号的运算时,只需要进行相应的变量代换即可完成相关工作。

2、连续时间信号的微分和积分

符号运算工具箱有强大的积分运算和求导功能。

连续时间信号的微分运算,可使用diff 命令函数来完成,其语句格式为: diff(function, ‘variable ’,n)

其中,function 表示需要进行求导运算的函数,或者被赋值的符号表达式;variable 为求导运算的独立变量;n 为求导阶数,默认值为一阶导数。

连续时间信号积分运算可以使用int 命令函数来完成,其语句格式为: int(function, ‘variable ’, a, b)

其中,function 表示被积函数,或者被赋值的符号表达式;variable 为积分变量;a 为积分下限,b 为积分上限,a 和b 默认时则求不定积分。

3、信号的相加和相乘运算

信号的相加和相乘是信号在同一时刻取值的相加和相乘。因此Matlab 对于时间信号的相加和相乘都是基于向量的点运算。

4、连续信号的卷积运算

卷积积分是信号和系统时域分析的重要方法之一。定义为:

?+∞

-=*=-2121d )t (f )(f )t (f )t (f )t (f τττ

Matlab 进行卷积计算可通过符号运算方法和数值计算方法实现。 (1)Matlab 符号运算法求连续信号卷积

从卷积定义出发,可以利用Matlab 符号运算法求卷积积分,但要注意积分变量和积分限的选取。

例:试用Matlab 符号运算法求卷积y(t)=[u(t)-u(t-1)]*[u(t)-u(t-1)]。

(2)Matlab 数值计算法求连续信号的卷积

例:试用Matlab 数值计算法求信号)2t (u )t (u )t (f 1--=和)t (u e )t (f 3t 2-=的卷积。

三、实验内容

1、已知信号的波形(课本P13例题),画出()

()()()2332----t f t f t f t f ,,,的波形图。

2、使用微分命令求y sin ln x x x =

**关于变量x 的一阶导数;使用积分命令计

算不定积分 dx x ax x ????

? ??+-22

5,定积分()dx x xe x ?+1

021。

3、已知()()()t t f t t f Ω=Ω=8sin ,sin 21,使用命令画出两信号和及两信号乘积的波形图。其中,Hz f 12=Ω

=

π

4

四、实验报告要求

1、格式:实验名称、实验目的、实验原理、实验环境、实验内容(上述4部分代码及结果图形)、实验思考等。

实验三 连续时间LTI 系统的时域分析

一、实验目的

1、学会使用符号法求解连续系统的零输入响应和零状态响应

2、学会使用数值法求解连续系统的零状态响应

3、学会求解连续系统的冲激响应和阶跃响应 二、实验原理及实例分析

1、连续时间系统零输入响应和零状态响应的符号求解

连续时间系统可以使用常系数微分方程来描述,其完全响应由零输入响应和零状态响应组成。MATLAB 符号工具箱提供了dsolve 函数,可以实现对常系数微分方程的符号求解,其调用格式为:

dsolve(‘eq1,eq2…’,’cond1,cond2,…’,’v ’)

其中参数eq 表示各个微分方程,它和MATLAB 符号表达式的输入基本相同,微分和导数的输入是使用Dy ,D2y ,D3y 来表示y 的一价导数,二阶导数,三阶导数;参数cond 表示初始条件或者起始条件;参数v 表示自变量,默认是变量t 。通过使用dsolve 函数可以求出系统微分方程的零输入响应和零状态响应,进而求出完全响应。 [实例1]试用Matlab 命令求齐次微分方程0)()(2)(='+''+'''t y t y t y 的零输入响应,已知起始条件为2)0(,1)0(,1)0(=''='=---y y y 。

注意,程序中绘图的时间区间一定要t>0,本程序中取[0, 8],程序运行后结果如下。

2、连续时间系统零状态响应的数值求解

在实际工程中使用较多的是数值求解微分方程。对于零输入响应来说,其数值解可以通过函数initial来实现,而该函数中的参量必须是状态变量所描述的系统模型,由于现在还没有学习状态变量相关内容,所以此处不做说明。对于零状态响应,MATLAB控制系统工具箱提供了对LTI系统的零状态响应进行数值仿真的函数lsim,利用该函数可以求解零初始条件下的微分方程的数值解。其调用格式为:

y=lsim(sys,f,t),其中t表示系统响应的时间抽样点向量,f是系统的输入向量;sys表示LTI系统模型,用来表示微分方程、差分方程或状态方程。在求解微分方程时,sys是有tf函数根据微分方程系数生成的系统函数对象,其语句格式为:sys=tf(a,b)。其中,a和b分别为微分方程右端和左端的系数向量。例如,对于微分方程

)()()()()()()()(0'1''2'''30'1''2'''3t f b t f b t f b f f b t y a t y a t y a t y a +++=+++

可以使用),(];,,,[];,,,[01230123a b tf sys b b b b b a a a a a ===获得其LTI 模型。注意,如果微分方程的左端或者右端表达式有缺项,则其向量a 或者b 中对应元素应该为零,不能省略不写。

3、连续时间系统冲激响应和阶跃响应的求解

在连续时间LTI系统中,冲激响应和阶跃响应是系统特性的描述。在MATLAB中,对于冲激响应和阶跃响应的数值求解,可以使用控制工具箱中提供的函数impulse和step来求解。

)

, (

) ,

(

t

sys

step

y t

sys

impulse

y

=

=

其中t表示系统响应的时间抽样点向量,sys表示LTI系统模型。

三、实验内容

1、已知系统的微分方程和激励信号,使用MATLAB命令画出系统的零状态响应和零输入响应(零状态响应分别使用符号法和数值法求解,零输入响应只使用符号法求解)。要求题目2必做,题目1选做。

2、已知系统的微分方程,使用MATLAB命令画出系统的冲激响应和阶跃响应(数值法)。要求题目2必做,题目1选做。

四、实验报告要求

1、格式:实验名称、实验目的、实验原理、实验环境、实验内容(上述几部分代码及结果图形)、实验思考等。

实验四傅里叶变换(FT)及其性质

一、实验目的

1、学会运用Matlab求连续时间信号的傅里叶

2、学会运用Matlab求连续时间信号的频谱图

3、学会运用Matlab分析连续时间信号的傅里叶变换的性质

二、实验原理及实例分析

(一)傅里叶变换的实现

例1:用Matlab 符号运算求解法求单边指数信号)()(2t u e t f t -=的FT 。

例2:用Matlab 符号运算求解法求2

11

)(ω

ω+=

j F 的IFT 。

例3:用Matlab 命令绘出例1中单边指数数信号的频谱图。

例4:用Matlab命令求图示三角脉冲的FT,并画出其幅度谱。

例5:用Matlab数值计算法求例3的三角脉冲幅度频谱图。

(二)FT 的性质

1、尺度变换

例6:设矩形信号)5.0()5.0()(--+=t u t u t f ,利用Matlab 命令绘出该信号及其频谱图。同时绘出)2()2/(t f t f 和的频谱图,并加以比较。

下面利用Matlab将常规矩形脉冲信号的频谱和其调制信号(课本例3-4信号)频谱进行比较。Matlab源程序如下:

傅里叶变换的其它性质可用类似的方法验证,希望大家课下练习。

三、实验内容

[注意:(1)写代码时j i]

3、分别利用Matlab 符号运算求解法和数值计算法求下图所示信号的FT ,并画出

其频谱图。

4、已知门函数自身卷积为三角波信号,试用Matlab 命令验证FT 的时域卷积定理。 四、实验报告要求

实验名称、实验目的、实验原理、实验环境、实验内容(上述几部分代码及结果图形)、实验思考等。

实验五 信号抽样及抽样定理

一、实验目的

1. 学会运用MATLAB 完成信号抽样以及对抽样信号的频谱进行分析

2. 学会运用MATLAB 改变抽样时间间隔,观察抽样后信号的频谱变化

3. 学会运用MATLAB 对抽样后的信号进行重建

二、实验原理及实例分析 (一)信号抽样

信号抽样是利用抽样脉冲序列

)(t p 从连续信号)(t f 中抽取一系列的离散值,通过

抽样过程得到的离散值信号称为抽样信号,记为)(t f s 。从数学上讲,抽样过程就是信

号相乘的过程,即)()()(t p t f t f s ?=

因此,可以使用傅里叶变换的频域卷积性质来求抽样信号)(t f s 的频谱。常用的抽样脉冲序列有周期矩形脉冲序列和周期冲激脉冲序列。

上式表明,信号在时域被抽样后,它的频谱是原连续信号频谱以抽样角频率为间隔周期的延拓,即信号在时域抽样或离散化,相当于频域周期化。在频谱的周期重复过程中,其频谱幅度受抽样脉冲序列的傅里叶系数加权,即被n P 加权。

可以看出,)(ωs F 是以s ω为周期等幅地重复。

信号与系统实验题目及答案

第一个信号实验的题目 1实现下列常用信号 (1)(5)u t +;(2)(1)t δ-;(3)cos(3)sin(2)t t +;(4)()[(1)(2)]f t t u t t u t t =?---; (5)0.5()4cos(),010t f t e t t π-=?= 2连续信号的基本运算与波形变换 已知信号2 2,2 1 ()33 t t f t ? -+-≤≤?=???,试画出下列各函数对时间t 的波形: (1)()f t -(2)(2)f t -+(3)(2)f t (4)1 (1)2 d f t dt +(5)(2)t f d ττ-∞-? 3连续信号的卷积运算 实现12()()f t f t *,其中1()f t 、2()f t 从第2个题目中任选3对组合。 4连续系统的时域分析 (1) 描述某连续系统的微分方程为()2()()()2()y t y t y t f t f t ''''++=+,求当输入信号为 2()2()t f t e u t -=时,该系统的零状态响应()y t 。 (2) 已知描述某连续系统的微分方程为2()()3()()y t y t y t f t '''+-=,试用MATLAB 绘出 该系统的冲激响应和阶跃响应的波形。 实验一答案: (1)(5)u t +在MATLAB 软件的输入程序及显示波形如下:

(2)(1)t δ-在MATLAB 软件的输入程序及显示波形如下: (3)cos(3)sin(2)t t +在MATLAB 软件的输入程序及显示波形如下: (4)()[(1)(2)]f t t u t t u t t =?---在MATLAB 软件的输入程序及显示波形如下: (5)0.5()4cos(),010t f t e t t π-=?=在MATLAB 软件的输入程序及显示波形如下:

信号与系统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 。 fs = 40; t = 0 : 1/fs : 1 ; % ?μ?ê·?±e?a5Hz,10Hz,20Hz,30Hz f1=5; xa = cos(2*pi*f1*t) ; subplot(1, 2, 1) ;

信号与系统实验指导书

信号与系统实验指导书 通信教研室 辽宁工业大学 2009年8月

目录 实验一时域卷积积分- 1 - 实验二连续系统时域分析- 10 - 实验三离散系统时域分析- 13 - 实验四二阶低通滤波器的响应- 16 -

1 MATLAB 在信号与系统中的应用 1.1.1 实验目的 (1) 练习连续信号的产生。 (2) 练习傅里叶变换的分析。 (3) 练习连续系统分析。 1.1.2 实验仪器 计算机、MATLAB 软件环境。 1.1.3 实验内容 在下面的实验操作中,认真保存、记录每项操作的作用和目的。 (一) 练习连续信号的产生 已知连续信号()()sin()j t x t e t αω+=+。要求编写程序文件siggen.m ,完成以下功能: (1) 在0≤t ≤5之间,产生该信号。其中0.6,5αω=-=。 (2) 在3个子图上分别画出该信号、信号的实部和虚部,并对图形进行标注。 此外, (3) 将(2)中产生的图形文件以bmp 格式保存到桌面。 (二) 练习连续傅里叶变换的分析 已知信号12()sin(2)2cos(2)s t f t f t ππ=+,其中f 1=47Hz ,f 2=88Hz 。要求: (1) 在0≤t ≤5范围内,步长增量为0.001,求出该信号的傅里叶变换; (2) 在2个子图上,分别绘制该信号的波形和幅频、相频响应曲线图。 (三) 练习连续系统分析 某LTI 系统输入信号为信号1110()0 t u t ≤≤?=? ?其它 ,系统的冲激响应为0.2()t h t e -=,

长度为15。要求: (1) 在时间间隔为0.5前提下,完成系统的卷积计算; (2) 在2个子图上,绘出输入信号和输出信号曲线,并进行标注。 1.1.4 实验预习 (1) 读懂各例题实验程序,了解MA TLAB基本操作方法。 (2) 根据实验内容预先编写实验程序。 1.1.5 实验报告 (1) 列写实验内容和已调试通过的实验程序,并按实验记录完成实验报告,打印实验程序产生的曲线图形。 (2) 思考题 ①若通信信号由有用信号和信道噪声组成,该信号该如何产生? ②连续系统分析的方法

matlab信号与系统实验报告

实验一 基本信号的产生与运算 一、 实验目的 学习使用MATLAB 产生基本信号、绘制信号波形、实现信号的基本运算。 二、 实验原理 MATLAB 提供了许多函数用于产生常用的基本信号:如阶跃信号、脉冲信号、指数信号、正弦信号和周期方波等等。这些信号是信号处理的基础。 1、 利用MATLAB 产生下列连续信号并作图。 (1)51),1(2)(<<---=t t u t x (2)300),3 2sin()(3.0<<=-t t e t x t (3)1.01.0,3000cos 100cos )(<<-+=t t t t x (4)2000),8.0cos()1.0cos()(<<=t t t t x ππ 答:(1)、 >> t=-1:0.02:5; >> x=(t>1); >> plot(t,-2*x); >> axis([-1,5,-3,1]); >> title('杨婕婕 朱艺星'); >> xlabel('x(t)=-2u(t-1)'); (2)、 >> t=0:0.02:30; >> x=exp(-0.3*t).*sin(2/3*t); >> plot(t,x);

>> title('杨婕婕朱艺星'); >> xlabel('x(t)=exp(-0.3*t).*sin(2/3*t)'); 因为原函数在t=15后x(t)取值接近于零,所以将横坐标改成0到15,看得更清晰 axis([0,15,-0.2,0.6]); (3)>> t=-0.1:0.01:0.1;x=cos(100*t)+cos(3000*t);plot(t,x); >> title('杨婕婕朱艺星'); >>xlabel('x=cos(100*t)+cos(3000*t)'); 因为t的间隔取太大,以至于函数不够准确,缩小t的间隔: t=-0.1:0.002:0.2;x=cos(100*t)+cos(3000*t); plot(t,x);title('杨婕婕') >> t=-0.1:0.0001:0.1; x=cos(100*t)+cos(3000*t); >> plot(t,x);title('杨婕婕朱艺星'); >> xlabel('x=cos(100*t)+cos(3000*t)');

信号与系统实验报告1

学生实验报告 (理工类) 课程名称:信号与线性系统专业班级:M11通信工程 学生学号:1121413017 学生姓名:王金龙 所属院部:龙蟠学院指导教师:杨娟

20 11 ——20 12 学年第 1 学期 金陵科技学院教务处制 实验报告书写要求 实验报告原则上要求学生手写,要求书写工整。若因课程特点需打印的,要遵照以下字体、字号、间距等的具体要求。纸张一律采用A4的纸张。 实验报告书写说明 实验报告中一至四项内容为必填项,包括实验目的和要求;实验仪器和设备;实验内容与过程;实验结果与分析。各院部可根据学科特点和实验具体要求增加项目。 填写注意事项 (1)细致观察,及时、准确、如实记录。 (2)准确说明,层次清晰。 (3)尽量采用专用术语来说明事物。 (4)外文、符号、公式要准确,应使用统一规定的名词和符号。 (5)应独立完成实验报告的书写,严禁抄袭、复印,一经发现,以零分论处。 实验报告批改说明 实验报告的批改要及时、认真、仔细,一律用红色笔批改。实验报告的批改成绩采用百分制,具体评分标准由各院部自行制定。 实验报告装订要求

实验批改完毕后,任课老师将每门课程的每个实验项目的实验报告以自然班为单位、按学号升序排列,装订成册,并附上一份该门课程的实验大纲。

实验项目名称:常用连续信号的表示 实验学时: 2学时 同组学生姓名: 无 实验地点: A207 实验日期: 11.12.6 实验成绩: 批改教师: 杨娟 批改时间: 一、实验目的和要求 熟悉MATLAB 软件;利用MATLAB 软件,绘制出常用的连续时间信号。 二、实验仪器和设备 586以上计算机,装有MATLAB7.0软件 三、实验过程 1. 绘制正弦信号)t Asin t (f 0?ω+=(),其中A=1,πω2=,6/π?=; 2. 绘制指数信号at Ae t (f =),其中A=1,0.4a -=; 3. 绘制矩形脉冲信号,脉冲宽度为2; 4. 绘制三角波脉冲信号,脉冲宽度为4;斜度为0.5; 5. 对上题三角波脉冲信号进行尺度变换,分别得出)2t (f ,)2t 2(f -; 6. 绘制抽样函数Sa (t ),t 取值在-3π到+3π之间; 7. 绘制周期矩形脉冲信号,参数自定; 8. 绘制周期三角脉冲信号,参数自定。 四、实验结果与分析 1.制正弦信号)t Asin t (f 0?ω+=(),其中A=1,πω2=,6/π?= 实验代码: A=1;

信号(MATLAB)实验指导书

《信号与系统》实验指导书 张建奇骆崇编写 浙江工业大学之江学院信息工程分院 2012年2月

目录 实验一MATLAB的基本使用 (1) 一、实验目的 (1) 二、实验原理 (1) 三、实验内容与要求 (8) 四、实验报告 (9) 实验二时域波形的MATLAB实现 (10) 一、实验目的 (10) 二、预习要求 (10) 三、实验原理 (10) 四、实验内容与要求 (18) 五、实验报告 (19) 实验三用MATLAB对系统时域分析 (20) 一、实验目的 (20) 二、预习要求 (20) 三、实验原理 (20) 四、实验内容与要求 (29)

实验一MATLAB的基本使用 一、实验目的 1、了解和掌握MATLAB的基本操作 2、了解MATLAB的库函数 3、会用MATLAB进行简单的操作。 二、实验原理 1、界面操作 MATLAB是“MATrix LABoratory”的缩写(矩阵实验室),它是由美国Mathworks公司于1984年正式推出的一种科学计算软件,由于其强大的功能,在欧美的一些大学里MATLAB已经成为许多诸如数字信号处理、自动控制理论等高级教程的主要工具软件,同时也成为理工科学生,必须掌握的一项基本技能。 当需要运行程序时,只需选择桌面上(或开始)中的MATLAB6.5应用程序图标即可 通常情况下,MATLAB的工作环境主要由一下几个窗口组成: 命令窗口(Command Window)

工作区间浏览器(Workspace) 历史命令窗口(Command History) 图形窗口(Figure) 文本编辑窗口(Editor) 当前路径窗口(Current Directory) MATLAB的命令窗与命令操作 当用户使用命令窗口进行工作时,在命令窗口中可以直接输入相应的命令,系统将自动显示信息。 例如在命令输入提示符“>>”后输入指令: >>t=[1,2,3;4,5,6;7,8,9]; 按回车键(Enter)后,系统即可完成对变量t的赋值。 MATALB提供了非常方便的在线帮助命令(help),它可提供各个函数的用法指南,包括格式、参数说明、注意事项及相关函数等内容。 2、图形窗 MATLAB图形窗(Figure)主要用于显示用户所绘制的图形。 通常,只要执行了任意一种绘图命令,图形窗就会自动产生。

信号与系统MATLAB实验报告

《信号与系统》MATLAB实验报告 院系:专业: 年级:班号: 姓名:学号: 实验时间: 实验地点:

实验一 连续时间信号的表示及可视化 实验题目: )()(t t f δ=;)()(t t f ε=;at e t f =)((分别取00<>a a 及); )()(t R t f =;)()(t Sa t f ω=;)2()(ft Sin t f π=(分别画出不同周期个数 的波形)。 解题分析: 以上各类连续函数,先运用t = t1: p:t2的命令定义时间范围向量,然后调用对应的函数,建立f 与t 的关系,最后调用plot ()函数绘制图像,并用axis ()函数限制其坐标范围。 实验程序: (1) )()(t t f δ= t=-1:0.01:3 %设定时间变量t 的范围及步长 f=dirac(t) %调用冲激函数dirac () plot(t,f) %用plot 函数绘制连续函数 axis([-1,3,-0.5,1.5]) %用axis 函数规定横纵坐标的范围 (2) )()(t t f ε= t=-1:0.01:3 %设定时间变量t 的范围及步长 f=heaviside(t) %调用阶跃函数heaviside () plot(t,f) %用plot 函数绘制连续函数 title('f(t)=heaviside(t)') %用title 函数设置图形的名称 axis([-1,3,-0.5,1.5]) %用axis 函数规定横纵坐标的范围 (3) at e t f =)( a=1时: t=-5:0.01:5 %设定时间变量t 的范围及步长 f=exp(t) %调用指数函数exp ()

实验1 用MATLAB进行信号频谱分析(推荐文档)

实验1 用MATLAB 进行信号频谱分析 一、实验目的 ㈠ 初步掌握MATLAB 产生常用离散时间信号的编程方法。 ㈡ 学习编写简单的FFT 算法程序,对离散信号进行幅频谱分析。 ㈢ 观察离散时间信号频谱的特点。 二、实验原理 ㈠ 常用的离散时间信号 在 MATLAB 语言主要是研究离散信号的。常用的离散信号有: 1.单位取样序列 ???≠==000 1)(n n n δ 2.单位阶跃序列 ?? ?<≥=0 01 )(n n n u 3.实指数序列 R a n a n x n ∈?=;)( 4.复指数序列 n e n x n j ?=+)(0)(ωσ 5.正(余)弦序列 )c o s ()(0θω+=n n x n ? 6.周期序列 n N n x n x ?+=)()( ㈡ 离散信号的产生 离散信号的图形显示使用stem 指令。 在 MATLAB 中的信号处理工具箱中,主要提供的信号是离散信号。

由于MATLAB 对下标的约定为从1开始递增,例如x=[5,4,3,2,1,0],表示x(1)=5,x(2)=4,X(3)=3… 因此要表示一个下标不由1开始的数组x(n),一般应采用两个矢量,如 n=[-3,-2,-1,0,l ,2,3,4,5]; x=[1,-l ,3,2,0,4,5,2,1]; 这表示了一个含9个采样点的矢量:X(n)={x(-3),x(-2),x(-1),x(0),x(1),x(2),x(3),x(4),x(5)}。 1.单位取样序列 ?? ?≠==δ0 001)(n n n 这一函数实现的方法有二: 方法一:可利用MATLAB 的zeros 函数。 x=zeros(1,N); %建立一个一行N 列的全零数组 x(1)=1; %对X (1)赋1 方法二:可借助于关系操作符实现 n=1:N; x=[n==1]; %n 等于1时逻辑关系式结果为真,x=1;n 不等于1时为假,x=0 如要产生 ?????≤<<=≤≤=-δ2 0210 100)(10)(n n n n n n n n n n n n 则可采用MATLAB 实现: n=n1:n2; x=[(n-n0)==0];%n=n0时逻辑关系式结果为真,x=1;n ≠n0时为

信号与系统实验报告_1(常用信号的分类与观察)

实验一:信号的时域分析 一、实验目的 1.观察常用信号的波形特点及产生方法 2.学会使用示波器对常用波形参数的测量 二、实验仪器 1.信号与系统试验箱一台(型号ZH5004) 2.40MHz双踪示波器一台 3.DDS信号源一台 三、实验原理 对于一个系统特性的研究,其中重要的一个方面是研究它的输入输出关系,即在一特定的输入信号下,系统对应的输出响应信号。因而对信号的研究是对系统研究的出发点,是对系统特性观察的基本手段与方法。在本实验中,将对常用信号和特性进行分析、研究。 信号可以表示为一个或多个变量的函数,在这里仅对一维信号进行研究,自变量为时间。常用信号有:指数信号、正弦信号、指数衰减正弦信号、复指数信号、Sa(t)信号、钟形信号、脉冲信号等。 1、信号:指数信号可表示为f(t)=Ke at。对于不同的a取值,其波形表现为不同的形式,如下图所示: 图1―1 指数信号 2、信号:其表达式为f(t)=Ksin(ωt+θ),其信号的参数:振幅K、角频率ω、与初始相位θ。其波形如下图所示:

图1-2 正弦信号 3、指数衰减正弦信号:其表达式为其波形如下图: 图1-3 指数衰减正弦信号 4、Sa(t)信号:其表达式为:。Sa(t)是一个偶函数,t= ±π,±2π,…,±nπ时,函数值为零。该函数在很多应用场合具有独特的运用。其信号如下图所示:

图1-4 Sa(t)信号 5、钟形信号(高斯函数):其表达式为:其信号如下图所示: 图1-5 钟形信号 6、脉冲信号:其表达式为f(t)=u(t)-u(t-T),其中u(t)为单位阶跃函数。其信号如下图所示: 7、方波信号:信号为周期为T,前T/2期间信号为正电平信号,后T/2期间信号为负电平信号,其信号如下图所示 U(t)

(完整word版)信号与系统matlab实验

习题三 绘制典型信号及其频谱图 1.更改参数,调试程序,绘制单边指数信号的波形图和频谱图。观察参数a对信号波形 及其频谱的影响。 程序代码: close all; E=1;a=1; t=0:0.01:4; w=-30:0.01:30; f=E*exp(-a*t); F=1./(a+j*w); plot(t,f);xlabel('t');ylabel('f(t)'); figure; plot(w,abs(F));xlabel('\omega');ylabel('|F(\omega)|';

E=1,a=1,波形图频谱图更改参数E=2,a=1;

更改参数a,对信号波形及其频谱的影响。(保持E=2)上图为a=1图像 a=2时

a=4时 随着a的增大,f(t)曲线变得越来越陡,更快的逼近0,而对于频谱图,随着a增大,图像渐渐向两边张开,峰值减小,陡度减小,图像整体变得更加平缓。 2.矩形脉冲信号 程序代码: close all; E=1;tao=1; t=-4:0.1:4; w=-30:0.1:30;

f=E*(t>-tao/2&tao/2)+0*(t<=-tao/2&t>=tao/2); F=(2*E./w).*sin(w*tao/2); plot(t,f);xlabel('t');ylabel('f(t)'); figure; plot(w,abs(F));xlabel('\omega');ylabel('|F(\omega)|') ; figure; plot(w,20*log10(abs(F))); xlabel('\omega');ylabel('|F(\omega)| in dB'); figure; plot(w,angle(F));xlabel('\omega');ylabel('\phi(\omega )');

matlab信号与系统实验报告

实验一 基本信号的产生与运算 一、 实验目的 学习使用MATLAB 产生基本信号、绘制信号波形、实现信号的基本运算。 二、 实验原理 MATLAB 提供了许多函数用于产生常用的基本信号:如阶跃信号、脉冲信号、指数信号、正弦信号和周期方波等等。这些信号是信号处理的基础。 1、 利用MATLAB 产生下列连续信号并作图。 (1)51),1(2)(<<---=t t u t x (2)300),3 2 sin()(3.0<<=-t t e t x t (3)1.01.0,3000cos 100cos )(<<-+=t t t t x (4)2000),8.0cos()1.0cos()(<<=t t t t x ππ 答:(1)、 >> t=-1:0.02:5; >> x=(t>1); >> plot(t,-2*x); >> axis([-1,5,-3,1]); >> title('杨婕婕 朱艺星'); >> xlabel('x(t)=-2u(t-1)');

(2)、 >> t=0:0.02:30; >> x=exp(-0.3*t).*sin(2/3*t); >> plot(t,x); >> title('杨婕婕朱艺星'); >> xlabel('x(t)=exp(-0.3*t).*sin(2/3*t)');

因为原函数在t=15后x(t)取值接近于零,所以将横坐标改成0到15,看得更清晰 axis([0,15,-0.2,0.6]);

(3)>> t=-0.1:0.01:0.1;x=cos(100*t)+cos(3000*t);plot(t,x); >> title('杨婕婕朱艺星'); >>xlabel('x=cos(100*t)+cos(3000*t)'); 因为t的间隔取太大,以至于函数不够准确,缩小t的间隔:t=-0.1:0.002:0.2;x=cos(100*t)+cos(3000*t); plot(t,x);title('杨婕婕')

北京理工大学信号与系统实验实验报告

实验1 信号的时域描述与运算 一、实验目的 1. 掌握信号的MATLAB表示及其可视化方法。 2. 掌握信号基本时域运算的MA TLAB实现方法。 3. 利用MA TLAB分析常用信号,加深对信号时域特性的理解。 二、实验原理与方法 1. 连续时间信号的MATLAB表示 连续时间信号指的是在连续时间范围内有定义的信号,即除了若干个不连续点外,在任何时刻信号都有定义。在MATLAB中连续时间信号可以用两种方法来表示,即向量表示法和符号对象表示法。 从严格意义上来说,MATLAB并不能处理连续时间信号,在MATLAB中连续时间信号是用等时间间隔采样后的采样值来近似表示的,当采样间隔足够小时,这些采样值就可以很好地近似表示出连续时间信号,这种表示方法称为向量表示法。表示一个连续时间信号需要使用两个向量,其中一个向量用于表示信号的时间范围,另一个向量表示连续时间信号在该时间范围内的采样值。例如一个正弦信号可以表示如下: >> t=0:0.01:10; >> x=sin(t); 利用plot(t,x)命令可以绘制上述信号的时域波形,如图1所示。 如果连续时间信号可以用表达式来描述,则还可以采用符号表达式來表示信号。例如对于上述正弦信号,可以用符号对象表示如下: >> x=sin(t); >> ezplot(X); 利用ezplot(x)命令可以绘制上述信号的时域波形 Time(seconds) 图1 利用向量表示连续时间信号

t 图 2 利用符号对象表示连续时间信号 sin(t) 2.连续时间信号的时域运算 对连续时间信号的运算包括两信号相加、相乘、微分、积分,以及位移、反转、尺度变换(尺度伸缩)等。 1)相加和相乘 信号相加和相乘指两信号对应时刻的值相加和相乘,对于两个采用向量表示的可以直接使用算术运算的运算符“+”和“*”来计算,此时要求表示两信号的向量时间范围和采样间隔相同。采用符号对象表示的两个信号,可以直接根据符号对象的运算规则运算。 2)微分和积分 对于向量表示法表示的连续时间信号,可以通过数值计算的方法计算信号的微分和积分。这里微分使用差分来近似求取的,由时间向量[N t t t ,,,21?]和采样值向量[N x x x ,,,21?]表示的连续时间信号,其微分可以通过下式求得 1,,2,1,|)('1-?=?-≈ +=N k t x x t x k k t t k 其中t ?表示采样间隔。MA TLAB 中用diff 函数来计算差分 k k x x -+1。 连续时间信号的定积分可以由MATLAB 的qud 函数实现,调用格式为 quad ('function_name',a,b) 其中,function_name 为被积函数名,a 、b 为积分区间。

信号与系统MATLAB实验

2016-2017学年第一学期 信号与系统实验报告 班级: 姓名: 学号: 成绩: 指导教师:

实验一常见信号的MATLAB 表示及运算 一.实验目的 1.熟悉常见信号的意义、特性及波形 2.学会使用MATLAB 表示信号的方法并绘制信号波形 3. 掌握使用MATLAB 进行信号基本运算的指令 4. 熟悉用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 所定义的时间点上的样值。 说明:plot 是常用的绘制连续信号波形的函数。 严格说来,MATLAB 不能表示连续信号,所以,在用plot()命令绘制波形时,要对自变量t 进行取值,MATLAB 会分别计算对应点上的函数值,然后将各个数据点通过折线连接起来绘制图形,从而形成连续的曲线。因此,绘制的只是近似波形,而且,其精度取决于t 的取样间隔。t 的取样间隔越小,即点与点之间的距离越小,则近似程度越好,曲线越光滑。例如:图1-1是在取样间隔为p=0.5时绘制的波形,而图1-2是在取样间隔p=0.1时绘制的波形,两相对照,可以看出图1-2要比图1-1光滑得多。

信号与系统课程设计应用MATLAB实现连续信号的采样与重构仿真

设计题目应用MATLAB实现连续信号的采样与重构仿真 1、设计目的 信号与系统课程设计是学习《信号与系统》课程必要的教学环节。由于该课程是专业基础课,需要通过实践了巩固基础知识,为使学生取得最现代化的设计技能和研究方法,课程设计训练也就成为了一个重要教学环节。通过对信号与系统一书的重新认识,我们将学习如何利用MATLAB软件进行仿真与重构并加深对滤波器的理解,这样的课程设计出了对我们的学习起着只关重要的作用,还可以很好的培养我们自己的动手能力。本次课程设计,我们会引入一个模拟的信号,通过MATLAB软件的防真技术来实现对它的分析、理解与学习。 MATLAB软件是今年来比较长用的一种数学软件,它有很强大的功能,主要侧重于某些理论知识的灵活运用。本次课程设计的目的是:增加对仿真软件MATLAB的感性认识,熟悉MATLAB软件平台的使用和MATLAB编程方法及常用语句;、初步掌握MATLAB的编程方法和特点;加深理解采样与重构的概念,应用MATLAB编程实现对信号的采样与重构;分别计算在临界采样、过采样、欠采样三种不同条件下重构信号的误差,并由此总结采样频率对信号重构误差的影响;学生需要自拟题目,根据自己手中的资料独立思考与分析,明确实习内容,制定实习步骤与方案,独立完成作业。 2、原理说明 2.1.1MATLAB MATLAB是美国Math Works公司产品,MATLAB现已被广泛于数学、通信、信号处理、自动控制、神经网络、图形处理等许多不同学科的研究中。并越来越多的应用到我们的学习生活中来,是目前通信工程上最广泛应用的软件之一。最初的MATLAB 只是一个数学计算工具。但现在的MATLAB已经远不仅仅是一个“矩阵实验室”,它已经成为一个集概念设计、算法开发、建模仿真,实时实现于一体的集成环境,它拥有许多衍生子集工具。 沈阳大学

《信号与系统》课程教学大纲

《信号与系统》课程教学大纲 课程编码:A0303051 总学时:64 理论学时:64 实验学时:0 学分:4 适用专业:通信工程 先修课程:电路,高等数学,复变函数与积分变换,线性代数 一、课程的性质与任务 《信号与系统》是电类专业的一门重要的专业课程。它的任务是研究信号和线性非时变系统的基本理论和基本分析方法,要求掌握最基本的信号变换理论,并掌握线性非时变系统的分析方法,为学习后续课程,以及从事相关领域的工程技术和科学研究工作奠定坚实的理论基础。通过本课程的学习,学生将理解信号的函数表示与系统分析方法,掌握连续时间系统和离散时间系统的时域分析和频域分析,连续时间系统的S域分析和离散时间系统的Z域分析,以及状态方程与状态变量分析法等相关内容。通过实验,使学生掌握利用计算机进行信号与系统分析的基本方法,加深对信号与线性非时变系统的基本理论的理解,训练学生的实验技能和科学实验方法,提高分析和解决实际问题的能力。

二、课程学时分配 教学章节理论实践 第一章:信号与系统导论6 第二章:连续系统的时域分析8 第三章:信号与系统的频域分析18 第四章:连续系统的复频域分析10 第五章:系统函数的零、极点分析8 第六章:离散系统的时域分析6 第七章:离散系统的Z域分析8 总计64 三、课程的基本教学内容及要求 第一章信号与系统导论(6学时) 1.教学内容 (1)历史的回顾,应用领域,信号的概念 (2)系统的概念,常用的基本信号 (3)信号的简单处理,单位冲激函数 2.重点及难点 教学重点:信号的描述、阶跃信号与冲激信号;信号的运算;线性时不变系统判据;系统定义 教学难点:信号及其分类,信号分析与处理,系统分析 3.课程教学要求

答案-信号与系统实验报告.

大连理工大学 本科实验报告 课程名称:___信号与系统实验学院:信息与通信工程学院专业:电子信息工程 班级: 学号: 学生姓名: 2012年12月11日

信号与系统实验 项目列表 信号的频谱图 Signals Frequency Spectrum 连续时间系统分析 Analysis for Continuous-time System 信号抽样 Signal Sampling 离散时间LTI系统分析 Analysis for Discrete-time LTI System 语音信号的调制解调 Modulation and Demodulation for Audio Signals Simulink?模拟信号的调制解调 Modulation and Demodulation for Analog Signals in Simulink ?

实验1信号的频谱图 一、实验目的 1. 掌握周期信号的傅里叶级数展开; 2. 掌握周期信号的有限项傅里叶级数逼近; 3. 掌握周期信号的频谱分析; 4. 掌握连续非周期信号的傅立叶变换; 5. 掌握傅立叶变换的性质。 二、实战演练(5道题) 1.已知周期三角信号如下图1-5所示,试求出该信号的傅里叶级数,利用MA TLAB编程 实现其各次谐波的叠加,并验证其收敛性。 解: 调试程序如下: clc clear t=-2:0.001:2; omega=pi; y=-(sawtooth(pi*t,0.5)/2+0.5)+1; plot(t,y),grid on; xlabel('t'),ylabel('周期三角波信号'); axis([-2 2 -0.5 1.5]) n_max=[1 3 5 11 47]; N=length(n_max); for k=1:N n=1:2: n_max(k); c=n.^2; b=4./(pi*pi*c); x=b*cos(omega*n'*t)+0.5; figure; plot(t,y,'b'); hold on; plot(t,x,'r'); hold off; xlabel('t'),ylabel('部分和的波形'); axis([-2 2 -0.5 1.5]);grid on; title(['最大谐波数=',num2str(n_max(k))]) end 运行结果如下:

信号与系统MATLAB实验总汇

实验一、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)绘制信号x(t)=)3 2sin(2t e t ?的曲线,t 的范围在0~30s ,取样时间间隔为0.1s 。t=0:0.1:30; y=exp(-sqrt(2)*t).*sin(2*t/3); plot(t,y);

(3)在n=[-10:10]范围产生离散序列:?? ?≤≤?=Other n n n x ,033,2)(,并绘图。n=-10:1:10; z1=((n+3)>=0); z2=((n-3)>=0); x=2*n.*(z1-z2); stem(n,x);(4)编程实现如下图所示的波形。 t=-2:0.001:3; f1=((t>=-1)&(t<=1)); f2=((t>=-1)&(t<=2)); f=f1+f2; plot(t,f); axis([-2,3,0,3]);

MATLAB实验教程

目录实验一、MATLAB 基本应用 实验二、实验二信号的时域表示 实验三、连续信号卷积 实验四、典型周期信号的频谱表示 实验五、傅立叶变换性质研究 实验六、系统的零极点分析 实验七离散信号分析

实验一 MATLAB 基本应用 一、实验目的:学习MATLAB的基本用法,了解 MATLAB 的目录结构和基本功能以及MATLAB在信号与系统中的应用。 二、实验内容: 例一已知x的取值范围,画出y=sin(x)的图型。 参考程序:x=0:0.05:4*pi; y=sin(x); plot(y) 例二计算y=sin(π/5)+4cos(π/4) 例三已知z 取值范围,x=sin(z);y=cos(z);画三维图形。 z=0:pi/50:10*pi; x=sin(z); y=cos(z); plot3(x,y,z)

xlabel('x') ylabel('y') zlabel('z') 例四已知x的取值范围,用subplot函数绘图。 参考程序:x=0:0.05:7; y1=sin(x); y2=1.5*cos(x); y3=sin(2*x); y4=5*cos(2*x); subplot(2,2,1),plot(x,y1),title('sin(x)') subplot(2,2,2),plot(x,y2),title('1.5*cos(x)') subplot(2,2,3),plot(x,y3),title('sin(2*x)') subplot(2,2,4),plot(x,y4),title('5*cos(2*x)')

单位阶跃信号 u(t+3)-2u(t) 其中 u(t)=1/2+(1/2)sign(t) Sigh(t)是符号函数t>0时为1,t<0时为-1; clear all %利用maple中的heaviside函数来实现u(t+3)-2u(t) syms t f=maple('Heaviside(t+3)-2*Heaviside(t)'); t1=-3*pi:0.01:3*pi;%注意精度问题 ff=subs(f,t,t1); figure(1); plot(t1,ff); axis([-5,5,-1.2,1.2]); title('u(t+3)-2u(t)'); %利用自己编写的heaviside函数来实现u(t+3)-2u(t) t=-5:0.01:5; f=Heaviside(t+3)-2*Heaviside(t); figure(2); plot(t,f); axis([-5,5,-1.2,1.2]); title('u(t+3)-2u(t)') %利用符号函数来实现u(t) t=-5:0.05:5; f=sign(t); ff=1/2+1/2*f; figure(3); plot(t,ff); axis([-5 5 -0.1 1.1]); title('u(t)') %利用符号函数来实现u(t+3)-2u(t) t=-5:0.01:5; f=(1/2+1/2*sign(t+3))-2*(1/2+1/2*sign(t)); figure(4) plot(t,f),axis([-5,5,-1.2,1.2]),title('u(t+3)-2u(t)')

基带信号眼图实验——matlab仿真

基带信号眼图实验——matlab 仿真

————————————————————————————————作者:————————————————————————————————日期: ?

数字基带信号的眼图实验——matla b仿真 一、实验目的 1、掌握无码间干扰传输的基本条件和原理,掌握基带升余弦滚降系统的实现方法; 2、通过观察眼图来分析码间干扰对系统性能的影响,并观察在输入相同码率的NRZ 基带信号下,不同滤波器带宽对输出信号码间干扰大小的影响程度; 3、熟悉MATL AB 语言编程。 二、实验预习要求 1、复习《数字通信原理》第七章7.1节——奈奎斯特第一准则内容; 2、复习《数字通信原理》第七章7.2节——数字基带信号码型内容; 3、认真阅读本实验内容,熟悉实验步骤。 三、实验原理和电路说明 1、基带传输特性 基带系统的分析模型如图3-1所示,要获得良好的基带传输系统,就应该 () n s n a t nT δ-∑() H ω() n s n a h t nT -∑基带传输抽样判决 图3-1?基带系统的分析模型 抑制码间干扰。设输入的基带信号为()n s n a t nT δ-∑,s T 为基带信号的码元周期,则经过基 带传输系统后的输出码元为 ()n s n a h t nT -∑。其中 1 ()()2j t h t H e d ωωωπ +∞ -∞ = ? ?(3-1) 理论上要达到无码间干扰,依照奈奎斯特第一准则,基带传输系统在时域应满足: 10()0,s k h kT k =?=? ? , 为其他整数 ?? ?(3-2) 频域应满足:

信号与系统实验报告

信号与系统实验教程(实验报告) 班级: 姓名: 学号: XXXX大学信息科学与技术学院二〇一六年五月十五日

实验一:连续时间信号的频域分析 实验目的: 1、掌握连续时间周期信号的傅里叶级数的物理意义和分析方法; 2、观察截短傅里叶级数产生的Gibbs现象,了解其特点及产生的原因; 3、掌握连续时间傅里叶变换的分析方法及其物理意义; 4、掌握各种典型的连续时间非周期信号的频谱特征以及傅里叶变换的主要性质; 5、学习掌握利用MATLAB语言编写计算CTFS、CTFT的程序,并能利用这些程序对一些典型信号进行频谱分析,验证CTFT的若干重要性质。 实验要求: 掌握并深刻理傅里叶变换的物理意义,掌握信号的傅里叶变换的计算方法,掌握利用MATLAB编程完成相关的傅里叶变换的计算。 实验原理: 1.傅里叶级数:任何周期函数都可以用正弦函数和余弦函数构成的无穷级数来表示。 2.吉布斯现象:将具有不连续点的周期函数(如矩形脉冲)进行傅立叶级数展开后,选取有限项进行合成。当选取的项数越多,在所合成的波形中出现的峰起越靠近原信号的不连续点。当选取的项数很大时,该峰起值趋于一个常数,大约等于总跳变值的9%。 实验内容: 1.参照例2-1程序,上机验证周期方波信号的傅里叶级数ck,并画出幅度谱|ck|。 1.1 程序代码 (1)准备:定义单位阶跃函数和delta函数 % filename u.m function y = u(t) y = (t>=0); % filename delta.m function y = delta(t) dt = 0.001; y = (u(t)-u(t-dt))/dt; 将u.m和delta.m分别保存到work文件夹中,或者保存在自己建立的文件夹中并将此文件夹设为工作路径(file->set path...) (2)验证方波信号的傅里叶级数ck并画频谱图: clear, close all T = 2; dt = 0.00001; t = -2:dt:2; x1 = u(t) - u(t-1-dt); x = 0;

相关文档
相关文档 最新文档