信号与系统实验二

信 号 与 系 统

实 实验名称：连续时间信号的频域分析

指导老师：

班 级：

学 号：姓 名：

验 报 告

实验二

苏永新

09通信工程1班

2009963924 王维

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

一、实验目的

1、掌握连续时间周期信号的傅里叶级数的物理意义和分析方法； 2、观察截短傅里叶级数而产生的“Gibbs现象”，了解其特点以及产生的原因； 3、掌握连续时间傅里叶变换的分析方法及其物理意义；

4、掌握各种典型的连续时间非周期信号的频谱特征以及傅里叶变换的主要性质； 5、学习掌握利用MATLAB语言编写计算CTFS、CTFT和DTFT的仿真程序，并能利用这些程序对一些典型信号进行频谱分析，验证CTFT、DTFT的若干重要性质。

基本要求：掌握并深刻理傅里叶变换的物理意义，掌握信号的傅里叶变换的计算方法，掌握利用MATLAB编程完成相关的傅里叶变换的计算。

二、实验原理及方法

1、连续时间周期信号的傅里叶级数CTFS分析

任何一个周期为T1的正弦周期信号，只要满足狄利克利条件，就可以展开成傅里叶级数。

其中三角傅里叶级数为：

x(t)?a0??[akcos(k?0t)?bksin(k?0t)] 2.1

k?1??或： x(t)?a0??ck?1kcos(k?0t??k) 2.2

其中?0?2?，称为信号的基本频率（Fundamental frequency），a0，ak，和bk分别是信T1号x(t)的直流分量、余弦分量幅度和正弦分量幅度，ck、?k为合并同频率项之后各正弦谐波分量的幅度和初相位，它们都是频率k?0的函数，绘制出它们与k?0之间的图像，称为信号的频谱图（简称“频谱”），ck－k?0图像为幅度谱，?k－k?0图像为相位谱。

三角形式傅里叶级数表明，如果一个周期信号x(t)，满足狄里克利条件，那么，它就可以被看作是由很多不同频率的互为谐波关系（harmonically related）的正弦信号所组成，其中每一个不同频率的正弦信号称为正弦谐波分量 (Sinusoid component)，其幅度（amplitude）为ck。也可以反过来理解三角傅里叶级数：用无限多个正弦谐波分量可以合成一个任意的非正弦周期信号。

指数形式的傅里叶级数为：

x(t)?k????aek?jk?0t 2.3

其中，ak为指数形式的傅里叶级数的系数，按如下公式计算：

1 ak?T1T1/2?T1/2?x(t)e?jk?0tdt 2.4

指数形式的傅里叶级数告诉我们，如果一个周期信号x(t)，满足狄里克利条件，那么，它就可以被看作是由很多不同频率的互为谐波关系（harmonically related）的周期复指数信号所组成，其中每一个不同频率的周期复指数信号称为基本频率分量，其复幅度（complex amplitude）为ak。这里“复幅度（complex amplitude）”指的是ak通常是复数。

上面的傅里叶级数的合成式说明，我们可以用无穷多个不同频率的周期复指数信号来合成任意一个周期信号。然而，用计算机（或任何其它设备）合成一个周期信号，显然不可能做到用无限多个谐波来合成，只能取这些有限个谐波分量来近似合成。

假设谐波项数为N，则上面的和成式为：

x(t)?k??N?aekNjk?0t 2.5

显然，N越大，所选项数越多，有限项级数合成的结果越逼近原信号x(t)。本实验可以比较直观地了解傅里叶级数的物理意义，并观察到级数中各频率分量对波形的影响包括“Gibbs”现象：即信号在不连续点附近存在一个幅度大约为9%的过冲，且所选谐波次数越多，过冲点越向不连续点靠近。这一现象在观察周期矩形波信号和周期锯齿波信号时可以看得很清楚。

三、实验内容和要求

实验前，必须首先阅读本实验原理，读懂所给出的全部范例程序。实验开始时，先在计算机上运行这些范例程序，观察所得到的信号的波形图。并结合范例程序应该完成的工作，进一步分析程序中各个语句的作用，从而真正理解这些程序。 实验前，一定要针对下面的实验项目做好相应的实验准备工作，包括事先编写好相应的实验程序等事项。

Q2-1 编写程序Q2_1，绘制下面的信号的波形图：

?111n?(0t)?cos3(?0t)?cos5(?0t)????sin()cosn( x(t)?cos??0t)

35n2n?1其中，?0 = 0.5π，要求将一个图形窗口分割成四个子图，分别绘制cos(?0t)、cos(3?0t)、cos(5?0t)

和x(t) 的波形图，给图形加title，网格线和x坐标标签，并且程序能够接受从键盘输入的和式中的项数。

抄写程序Q2_1如下：

clear, close all

dt = 0.00001; t = -2:dt:4; w0=0.5*pi; x1=cos(w0.*t); x2=cos(3*w0.*t); x3=cos(5*w0.*t);

N=input ( 'Type in the number of the harmonic components N ='); x=0;

for q=1:N;

x=x+(sin(q*(pi/2)).*cos(q*w0*t))/q; end

subplot(221), plot(t,x1)

axis([-2 4 -2 2]); grid on,

title('signal cos(w0.*t)'), subplot(222), plot(t,x2)

axis([-2 4 -2 2]); grid on

title('signal cos(3*w0.*t))'), subplot(223) plot(t,x3)

axis([-2 4 -2 2]); grid on

title('signal cos(5*w0.*t))'), xlabel ('time t (sec)') subplot(224) plot(t,x)

axis([-2 4 -2 2]); grid on

title('signal x(t)'), xlabel ('time t (sec)')

执行程序Q2_1所得到的图形如下：

signal cos(w0.*t)210-1-2-2024210-1-2-2signal cos(3*w0.*t))0signal x(t)24signal cos(5*w0.*t))210-1-2-202time t (sec)4210-1-2-202time t (sec)4

Q2-2 给程序Program2_1增加适当的语句，并以Q2_2存盘，使之能够计算例题2-1中的

周期方波信号的傅里叶级数的系数，并绘制出信号的幅度谱和相位谱的谱线图。

通过增加适当的语句修改Program2_1而成的程序Q2_2抄写如下：

clear,close all

T = 2; dt = 0.00001; t = -2:dt:2; x1 = u(t+0.2)-u(t-0.2-1-dt); x = 0; for m = -1:1

x = x + u(t+0.2-m*T) - u(t-0.2-1-m*T-dt); end w0 = 2*pi/T; N = 10; L = 2*N+1; for k = -N:1:N;

ak(N+1+k) = (1/T)*x1*exp(-j*k*w0*t')*dt; end

phi = angle(ak); y=0; for q = 1:L;

y = y+ak(q)*exp(j*(-(L-1)/2+q-1)*2*pi*t/T); end;

subplot(221),

plot(t,x), title('The original signal x(t)'), axis([-2,2,-0.2,1.2]), subplot(223),

plot(t,y), title('The synthesis signal y(t)'), axis([-2,2,-0.2,1.2]), xlabel('Time t'), subplot(222)

k=-N:N; stem(k,abs(ak),'k.'), title('The amplitude |ak| of x(t)'),

搜索“diyifanwen.net”或“第一范文网”即可找到本站免费阅读全部范文。收藏本站方便下次阅读，第一范文网，提供最新IT计算机信号与系统实验二 全文阅读和word下载服务。