实验二_连续时间信号的频域分析
信号与系统实验总结1
实验总结班级：10电子班学号：1039035 姓名：田金龙这学期的实验都有：信号的时域分析、线性时不变系统的时域分析、连续时间信号系统的频域分析、连续时间在连续时间信号的频域LTI系统的复频域分析、连续时间LTI系统的频域分析。
在这学期的学习中学习了解到很多关于信号方面的处理方法加上硬件动手的实践能力，让我对课堂上所学到的知识有了更深层次的理解也加深了所学知识的印象。
下面则是对每次实验的分析和总结：实验一：信号的时域分析在第一次试验中进行信号的时域分析还有的就是学会使用MATLAB软件来利用它实现一些相关的运算并且绘制出相关的信号图。
在时域分析中掌握连续时间信号和离散时间信号的描述方法，并能够实现各种信号的时域变化和运算。
了解单位阶跃信号和单位冲激信号的拓展函数，以便于熟悉这两种函数在之后的程序中的应用。
在能够对简单信号的描述的前提下，通过一些简单的程序，实现信号的分析，时域反相，时域尺度变换和周期信号的描述。
clear,close alldt=0.01;t=-2:dt:2;x=u(t);plot(t,x)title('u signal u(t)')grid on连续时间信号的时域分析后，则是离散时间信号的仿真。
通过对连续时间信号的描述和对离散时间信号的描述，发现它们的不同之处在于对时间的定义和对函数的图形描述。
在离散时间信号的图形窗口描述时，使用的是stem（n，x）函数。
在硬件实验中，使用一些信号运算单元，加法器，减法器，倍乘器，反相器，积分器和微分器。
输入相应的简单信号，观察通过不同运算单元输出的信号。
实验二：线性时不变系统的时域分析在线性时不变系统的时域分析中主要研究的就是信号的卷积运算，学会进行信号的卷积运算和MATLAB对卷积运算的实现。
而系统则通常是由若干部件或单元组成的一个整体，根据系统所处理的信号不同，系统又有多种不同的分类。
而在学习总最常研究的则是线性时不变系统，而线性时不变系统则是形同同时满足齐次性和叠加性。
连续时间信号的频域分析及Matlab实现
function CTF3()
1/2 exp(-2 t) heaviside(t)
syms t v w x;
F = fourier(x);
0.4
x = 1/2*exp(-2*t)*sym('Heaviside(t)');0.2
fliplr例子
6 4
>> n = 0:4; >> a = [5 4 3 2 1]; >> subplot(2,1,1),stem(n,a);
2
>> b = fliplr(a);
>> k = -4:4; >> c = [b,a(2:end)];
0 -5
-4
-3
-2
-1
0
1
2
3
4
5
6
4
>> subplot(2,1,2),stem(k,c);
0 0 0.5 1 t 1/2/abs(2+i w) 0.25 0.2 0.15 0.1 -6 -4 -2 0 w 2 4 6 1.5 2
subplot(2,1,1);
ezplot(x); subplot(2,1,2);
ezplot(abs(F));
f(t) = u(t+1) - u(t-1) 1
function [A_sym,B_sym] = CTF2()
syms t n k x T = 5; tao = T/5; a = 0; Nf = 16; Nn = 32; x1 = sym('Heaviside(t+0.5)')*h; x = x1 - sym('Heaviside(t-0.5)')*h; A0 = 2*int(x,t,-a,T-a)/T;%求出三角函数展开系数A0
实验二 信号的频域分析
实验二信号的频域分析一．实验目的1.通过MATLAB编程观察周期信号的合成过程，进一步了解周期信号的傅里叶级数分解特性。
2.学习用MATLAB绘制周期信号的频谱的方法，观察周期信号的离散性、谐波性和收敛性。
3.用MATLAB研究矩形波的频谱，观察周期变化和脉冲变化对频谱的影响。
4.学习用MATLAB的符号运算功能计算傅里叶变换和反变换的方法。
5.掌握用MATLAB绘制非周期信号频谱的数值方法。
6.掌握信号的采集方法与过程以及信号恢复的原理与方法。
二．实验内容1. 题目：2-1已知信号周期如图所示。
其中T=4s，A=1，试：1）画出两个周期信号的频谱图。
2）求傅里叶级数最高谐波次数20的部分和波形，并比较。
源程序：T=4N=20t=linspace(-T/2,T/2,1000)w0=2*pi/Tf=sawtooth(w0*t,1)for k=0:Na(k+1)=2/T*trapz(t,f.*cos(k*w0*t))b(k+1)=2/T*trapz(t,f.*sin(k*w0*t))endAn=a-j*bAn(1)=a(1)tt=[t-T,t,t+T]ff=[f,f,f]subplot(2,3,1);plot(tt,ff,'LineWidth',1.5);ylabel('f(t)')title('周期信号波形','Fontsize',8)subplot(2,3,2);h=stem(0:N,abs(An),'.');hold onh=stem(0:-1:-N,abs(An),'.')title('幅度频谱','Fontsize',8)xlabel('\omega');ylabel('F(j\omega)')subplot(2,3,3);h=stem(0:N,angle(An),'.');hold onh=stem(0:-1:-N,-angle(An),'.')title('相位频谱','Fontsize',8)xlabel('\omega');ylabel('\phi(j\omega)')f=sawtooth(w0*t,0.5)for k=0:Na(k+1)=2/T*trapz(t,f.*cos(k*w0*t)) b(k+1)=2/T*trapz(t,f.*sin(k*w0*t)) end An=a-j*b An(1)=a(1) tt=[t-T,t,t+T] ff=[f,f,f]subplot(2,3,4);plot(tt,ff,'LineWidth',1.5);ylabel('f(t)') title('周期信号波形','Fontsize',8)subplot(2,3,5);h=stem(0:N,abs(An),'.');hold on h=stem(0:-1:-N,abs(An),'.')title('幅度频谱','Fontsize',8)xlabel('\omega');ylabel('F(j\omega)')subplot(2,3,6);h=stem(0:N,angle(An),'.');hold on h=stem(0:-1:-N,-angle(An),'.')title('相位频谱','Fontsize',8)xlabel('\omega');ylabel('\phi(j\omega)')程序运行结果（截图）：-1010-1-0.500.51f (t )周期信号波形幅度频谱ωF (j ω)-2020相位频谱ωφ(j ω)-1010-1-0.500.51f (t )周期信号波形幅度频谱ωF (j ω)-20020相位频谱ωφ(j ω)2. 题目：2-2求下列信号的傅里叶变换。
通信原理实验教程(MATLAB)
实验教程目录实验一：连续时间信号与系统的时域分析-------------------------------------------------6一、实验目的及要求---------------------------------------------------------------------------6二、实验原理-----------------------------------------------------------------------------------61、信号的时域表示方法------------------------------------------------------------------62、用MATLAB仿真连续时间信号和离散时间信号----------------------------------73、LTI系统的时域描述-----------------------------------------------------------------11三、实验步骤及内容--------------------------------------------------------------------------15四、实验报告要求-----------------------------------------------------------------------------26 实验二：连续时间信号的频域分析---------------------------------------------------------27一、实验目的及要求--------------------------------------------------------------------------27二、实验原理----------------------------------------------------------------------------------271、连续时间周期信号的傅里叶级数CTFS---------------------------------------------272、连续时间信号的傅里叶变换CTFT--------------------------------------------------283、离散时间信号的傅里叶变换DTFT -------------------------------------------------284、连续时间周期信号的傅里叶级数CTFS的MATLAB实现------------------------295、用MATLAB实现CTFT及其逆变换的计算---------------------------------------33三、实验步骤及内容----------------------------------------------------------------------34四、实验报告要求-------------------------------------------------------------------------48 实验三：连续时间LTI系统的频域分析---------------------------------------------------49一、实验目的及要求--------------------------------------------------------------------------49二、实验原理----------------------------------------------------------------------------------491、连续时间LTI系统的频率响应-------------------------------------------------------492、LTI系统的群延时---------------------------------------------------------------------503、用MATLAB计算系统的频率响应--------------------------------------------------50三、实验步骤及内容----------------------------------------------------------------------51四、实验报告要求-------------------------------------------------------------------------58 实验四：调制与解调以及抽样与重建------------------------------------------------------59一、实验目的及要求--------------------------------------------------------------------------59二、实验原理----------------------------------------------------------------------------------591、信号的抽样及抽样定理---------------------------------------------------------------592、信号抽样过程中的频谱混叠----------------------------------------------------------623、信号重建--------------------- ----------------------------------------------------------624、调制与解调----------------------------------------------------------------------------------645、通信系统中的调制与解调仿真---------------------------------------------------------66三、实验步骤及内容------------------------------------------------------------------------66四、实验报告要求---------------------------------------------------------------------------75 实验五：连续时间LTI系统的复频域分析----------------------------------------------76一、实验目的及要求------------------------------------------------------------------------76二、实验原理--------------------------------------------------------------------------------761、连续时间LTI系统的复频域描述--------------------------------------------------762、系统函数的零极点分布图-----------------------------------------------------------------773、拉普拉斯变换与傅里叶变换之间的关系-----------------------------------------------784、系统函数的零极点分布与系统稳定性和因果性之间的关系------------------------795、系统函数的零极点分布与系统的滤波特性-------------------------------------------806、拉普拉斯逆变换的计算-------------------------------------------------------------81三、实验步骤及内容------------------------------------------------------------------------82四、实验报告要求---------------------------------------------------------------------------87 附录：授课方式和考核办法-----------------------------------------------------------------88实验一信号与系统的时域分析一、实验目的1、熟悉和掌握常用的用于信号与系统时域仿真分析的MA TLAB函数；2、掌握连续时间和离散时间信号的MATLAB产生，掌握用周期延拓的方法将一个非周期信号进行周期信号延拓形成一个周期信号的MATLAB编程；3、牢固掌握系统的单位冲激响应的概念，掌握LTI系统的卷积表达式及其物理意义，掌握卷积的计算方法、卷积的基本性质；4、掌握利用MA TLAB计算卷积的编程方法，并利用所编写的MA TLAB程序验证卷积的常用基本性质；掌握MATLAB描述LTI系统的常用方法及有关函数，并学会利用MATLAB求解LTI系统响应，绘制相应曲线。
连续时间信号与系统的频域分析报告
连续时间信号与系统的频域分析报告1. 引言连续时间信号与系统的频域分析是信号与系统理论中的重要分支，通过将信号和系统转换到频域，可以更好地理解和分析信号的频谱特性。
本报告将对连续时间信号与系统的频域分析进行详细介绍，并通过实例进行说明。
2. 连续时间信号的频域表示连续时间信号可以通过傅里叶变换将其转换到频域。
傅里叶变换将信号分解成一系列不同频率的正弦和余弦波的和。
具体来说，对于连续时间信号x(t)，其傅里叶变换表示为X(ω)，其中ω表示频率。
3. 连续时间系统的频域表示连续时间系统可以通过频域中的频率响应来描述。
频率响应是系统对不同频率输入信号的响应情况。
通过系统函数H(ω)可以计算系统的频率响应。
系统函数是频域中系统输出与输入之比的函数，也可以通过傅里叶变换来表示。
4. 连续时间信号的频域分析频域分析可以帮助我们更好地理解信号的频谱特性。
通过频域分析，我们可以获取信号的频率成分、频谱特性以及信号与系统之间的关系。
常用的频域分析方法包括功率谱密度估计、谱线估计等。
5. 连续时间系统的频域分析频域分析也可以用于系统的性能评估和系统设计。
通过分析系统的频响特性，我们可以了解系统在不同频率下的增益和相位变化情况，进而可以对系统进行优化和设计。
6. 实例分析以音频信号的频域分析为例，我们可以通过对音频信号进行傅里叶变换，将其转换到频域。
通过频域分析，我们可以获取音频信号的频谱图，从而了解音频信号的频率成分和频率能量分布情况。
进一步，我们可以对音频信号进行系统设计和处理，比如对音乐进行均衡、滤波等操作。
7. 结论连续时间信号与系统的频域分析是信号与系统理论中重要的内容，通过对信号和系统进行频域分析，可以更好地理解和分析信号的频谱特性。
频域分析也可以用于系统的性能评估和系统设计，对于音频信号的处理和优化具有重要意义。
总结：通过本报告，我们了解了连续时间信号与系统的频域分析的基本原理和方法。
频域分析可以帮助我们更好地理解信号的频谱特性和系统的频响特性，对系统设计和信号处理具有重要意义。
实验二--连续时间信号的频域分析
实验二连续时间信号的频域分析专业班级通信1601 姓名宁硕学号 20 评分：实验日期： 2017 年 12 月 13日指导教师: 张鏖峰一、实验目的1、掌握连续时间周期信号的傅里叶级数的物理意义和分析方法；2、观察截短傅里叶级数而产生的“Gibbs现象”，了解其特点以及产生的原因；3、掌握连续时间傅里叶变换的分析方法及其物理意义；4、掌握各种典型的连续时间非周期信号的频谱特征以及傅里叶变换的主要性质；5、学习掌握利用MATLAB语言编写计算CTFS、CTFT和DTFT的仿真程序，并能利用这些程序对一些典型信号进行频谱分析，验证CTFT、DTFT的若干重要性质。
基本要求：掌握并深刻理傅里叶变换的物理意义，掌握信号的傅里叶变换的计算方法，掌握利用MATLAB编程完成相关的傅里叶变换的计算。
以看得很清楚。
二、实验原理及方法任何一个周期为T1的正弦周期信号，只要满足狄利克利条件，就可以展开成傅里叶级数。
其中三角傅里叶级数为：2.1或：2.2指数形式的傅里叶级数为：2.3其中，为指数形式的傅里叶级数的系数，按如下公式计算：2.4傅里叶变换在信号分析中具有非常重要的意义，它主要是用来进行信号的频谱分析的。
傅里叶变换和其逆变换定义如下：2.52.6连续时间傅里叶变换主要用来描述连续时间非周期信号的频谱。
按照教材中的说法，任意非周期信号，如果满足狄里克利条件，那么，它可以被看作是由无穷多个不同频率（这些频率都是非常的接近）的周期复指数信号ejt的线性组合构成的，每个频率所对应的周期复指数信号ejt称为频率分量（frequency component），其相对幅度为对应频率的|X(j)|之值，其相位为对应频率的X(j)的相位三、实验内容和要求Q2-1 编写程序Q2_1，绘制下面的信号的波形图：其中，0 = 0.5π，要求将一个图形窗口分割成四个子图，分别绘制cos(0t)、cos(30t)、cos(50t) 和x(t) 的波形图，给图形加title，网格线和x 坐标标签，并且程序能够接受从键盘输入的和式中的项数。
实验二 连续时间信号的频域分析
实验二连续时间信号的频域分析专业班级通信1601 姓名宁硕学号20 评分：实验日期： 2017 年 12 月 13日指导教师:张鏖峰一、实验目的1、掌握连续时间周期信号的傅里叶级数的物理意义和分析方法；2、观察截短傅里叶级数而产生的“Gibbs现象”，了解其特点以及产生的原因；3、掌握连续时间傅里叶变换的分析方法及其物理意义；4、掌握各种典型的连续时间非周期信号的频谱特征以及傅里叶变换的主要性质；5、学习掌握利用MATLAB语言编写计算CTFS、CTFT和DTFT的仿真程序，并能利用这些程序对一些典型信号进行频谱分析，验证CTFT 、DTFT 的若干重要性质。
基本要求：掌握并深刻理傅里叶变换的物理意义，掌握信号的傅里叶变换的计算方法，掌握利用MATLAB 编程完成相关的傅里叶变换的计算。
以看得很清楚。
二、实验原理及方法任何一个周期为T 1的正弦周期信号，只要满足狄利克利条件，就可以展开成傅里叶级数。
其中三角傅里叶级数为：∑∞=++=1000)]sin()cos([)(k k k t k b t k a a t x ωω或： ∑∞=++=100)cos()(kk k t k A A t x ϕω指数形式的傅里叶级数为：∑∞-∞==kt jk k e F t x 0)(ω 其中，k F 为指数形式的傅里叶级数的系数，按如下公式计算：⎰--=2/2/111)(1T Tt jk k dt e t x T F ω傅里叶变换在信号分析中具有非常重要的意义，它主要是用来进行信号的频谱分析的。
傅里叶变换和其逆变换定义如下：⎰∞∞--=dt e t x j X t j ωω)()(⎰∞∞-=ωωπωd e j X t x tj )(21)( 连续时间傅里叶变换主要用来描述连续时间非周期信号的频谱。
按照教材中的说法，任意非周期信号，如果满足狄里克利条件，那么，它可以被看作是由无穷多个不同频率（这些频率都是非常的接近）的周期复指数信号e j?t 的线性组合构成的，每个频率所对应的周期复指数信号e j?t 称为频率分量（frequency component ），其相对幅度为对应频率的|X(j?)|之值，其相位为对应频率的X(j?)的相位三、实验内容和要求Q2-1 编写程序Q2_1，绘制下面的信号的波形图：Λ-+-=)5cos(51)3cos(31)cos()(000t t t t x ωωω∑∞==10)cos()2sin(1n t n n nωπ其中，?0 = π，要求将一个图形窗口分割成四个子图，分别绘制cos(?0t)、cos(3?t)、cos(5?t) 和x(t) 的波形图，给图形加title，网格线和x坐标标签，并且程序能够接受从键盘输入的和式中的项数。
实验二 连续信号与系统的频域分析
(2)绘出f(t)的时域波形及频谱图。
f(t) 1
(2)电路的系统函数为 H(jω)
1 j 1 j
(b)用MATLAB求系统的单位冲激响应。
(c) 当输入为 f (t) sint sin3t t
求系统的稳态响应。
程序清单: syms w t; Hw=(1-j*w)/(1+j*w); ht=ifourier(Hw,t); ft=sin(t)+sin(3*t); Fw=fourier(ft); Yw=Fw*Hw; yt=ifourier(Yw,t);
函数fourier()——傅立叶正变换 函数ifourier()——傅立叶逆变换
3、连续时间系统的响应
已知某电路的系统传递函数为 H(jω)=1/(0.08(jω)2+0.4jω+1)
用MATLAB绘制系统的幅频特性曲线和相频特性曲线，并分
析该系统的频率特性。 系统频率特性
H(j)=|H(j)|ej() |H(j)|——系统的幅频特性
for i=1:9
a(i)=subs(an,n,i); % 计算系数a1~a9,存于数组a中 end
a0=double(a0);a=double(a); %转换成数值型
stem(0,a0,i,a); %绘f(t)的频谱图
2、非周期信号的分析 (1)已知某一连续时间信号为
f t e2 t
试绘出它的时域波形及相应的频谱图。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
实验二 连续时间信号的频域分析一、实验目的1、掌握连续时间周期信号的傅里叶级数的物理意义和分析方法；2、观察截短傅里叶级数而产生的“Gibbs 现象”，了解其特点以及产生的原因；3、掌握连续时间傅里叶变换的分析方法及其物理意义；4、掌握各种典型的连续时间非周期信号的频谱特征以及傅里叶变换的主要性质；5、学习掌握利用Matlab 语言编写计算CTFS 、CTFT 和DTFT 的仿真程序，并能利用这些程序对一些典型信号进行频谱分析，验证CTFT 、DTFT 的若干重要性质。
基本要求：掌握并深刻理傅里叶变换的物理意义，掌握信号的傅里叶变换的计算方法，掌握利用Matlab 编程完成相关的傅里叶变换的计算。
二、原理说明1、连续时间周期信号的傅里叶级数CTFS 分析任何一个周期为T 1的正弦周期信号，只要满足狄利克利条件，就可以展开成傅里叶级数。
三角傅里叶级数为：∑∞=++=1000)]sin()cos([)(k k k t k b t k a a t x ωω 2.1或： ∑∞=++=100)cos()(k k kt k ca t x ϕω 2.2其中102T πω=，称为信号的基本频率（Fundamental frequency ），k k b a a ，和，0分别是信号)(t x 的直流分量、余弦分量幅度和正弦分量幅度，k k c ϕ、为合并同频率项之后各正弦谐波分量的幅度和初相位，它们都是频率0ωk 的函数，绘制出它们与0ωk 之间的图像，称为信号的频谱图（简称“频谱”），k c －0ωk 图像为幅度谱，k ϕ－0ωk 图像为相位谱。
三角形式傅里叶级数表明，如果一个周期信号x(t)，满足狄里克利条件，就可以被看作是由很多不同频率的互为谐波关系（harmonically related ）的正弦信号所组成，其中每一个不同频率的正弦信号称为正弦谐波分量 (Sinusoid component)，其幅度（amplitude ）为k c 。
也可以反过来理解三角傅里叶级数：用无限多个正弦谐波分量可以合成一个任意的非正弦周期信号。
指数形式的傅里叶级数为：∑∞-∞==k tjk kea t x 0)(ω 2.3其中，k a 为指数形式的傅里叶级数的系数，按如下公式计算：⎰--=2/2/1110)(1T T tjk k dt e t x T a ω 2.4 指数形式的傅里叶级数告诉我们，如果一个周期信号x(t)，满足狄里克利条件，那么，它就可以被看作是由很多不同频率的互为谐波关系（harmonically related ）的周期复指数信号所组成，其中每一个不同频率的周期复指数信号称为基本频率分量，其复幅度（complex amplitude ）为k a 。
这里“复幅度（complex amplitude ）”指的是k a 通常是复数。
上面的傅里叶级数的合成式说明，我们可以用无穷多个不同频率的周期复指数信号来合成任意一个周期信号。
然而，用计算机（或任何其它设备）合成一个周期信号，显然不可能做到用无限多个谐波来合成，只能取这些有限个谐波分量来近似合成。
假设谐波项数为N ，则上面的和成式为：∑-==NNk tjk kea t x 0)(ω 2.5显然，N 越大，所选项数越多，有限项级数合成的结果越逼近原信号x(t)。
本实验可以比较直观地了解傅里叶级数的物理意义，并观察到级数中各频率分量对波形的影响包括“Gibbs ”现象：即信号在不连续点附近存在一个幅度大约为9%的过冲，且所选谐波次数越多，过冲点越向不连续点靠近。
这一现象在观察周期矩形波信号和周期锯齿波信号时可以看得很清楚。
2、连续时间信号傅里叶变换----CTFT傅里叶变换在信号分析中具有非常重要的意义，它主要是用来进行信号的频谱分析的。
傅里叶变换和其逆变换定义如下：⎰∞∞--=dt et x j X tj ωω)()( 2.6⎰∞∞-=ωωπωd e j X t x tj )(21)( 2.7 连续时间傅里叶变换主要用来描述连续时间非周期信号的频谱。
按照教材中的说法，任意非周期信号，如果满足狄里克利条件，那么，它可以被看作是由无穷多个不同频率（这些频率都是非常的接近）的周期复指数信号e j t 的线性组合构成的，每个频率所对应的周期复指数信号e j t 称为频率分量（frequency component ），每个频率的幅度为对应频率的|X(j ω)|之值，其相位为对应频率的X(j ω)的相位。
X(j ω)通常为复函数，可以按照复数的极坐标表示方法表示为：X(j ω)=| X(j ω)|e j X(j ω)其中，| X(j ω)|称为x(t)的幅度谱，而X(j ω)则称为x(t)的相位谱。
给定一个连续时间非周期信号x(t)，它的频谱也是连续且非周期的。
对于连续时间周期信号，也可以用傅里叶变换来表示其频谱，其特点是，连续时间周期信号的傅里叶变换由冲激序列构成的，是离散的——这是连续时间周期信号的傅里叶变换的基本特征。
4.1 傅里叶级数的Matlab 计算设周期信号x(t)的基本周期为T 1，且满足狄里克利条件，则其傅里叶级数的系数可由式2.4计算得到。
式2.4重写如下：⎰--=2/2/1110)(1T T tjk k dt et x T a ω基本频率为： 102T πω=对周期信号进行分析时，我们往往只需对其在一个周期进行分析即可，通常选择主周期（Principle period ）。
假定x 1(t)是x(t)中的主周期，则⎰--=2/2/11110)(1T T tjk k dt e t x T a ω 计算机不能计算无穷多个系数，所以我们假设需要计算的谐波次数为N ，则总的系数个数为2N+1个。
在确定了时间围和时间变化的步长即T 1和dt 之后，对某一个系数，上述系数的积分公式可以近似为： ∑⎰---==n tjk nT T t jk k T dt e t x dt e t x T a 12/2/11/)()(10110ωω121/],,[)](),(),([02010T dt e e e t x t x t x M t jk t jk t jk M ⋅⋅=---ωωωΛΛ对于全部需要的2N+1个系数，上面的计算可以按照矩阵运算实现。
Matlab 实现系数计算的程序如下：dt = 0.01;T = 2; t = -T/2:dt:T/2; w0 = 2*pi/T;x1 = input(‘Type in the periodic signal x(t) over one period x1(t)=’); N = input(‘Type in the number N=’); k = -N:N; L = 2*N+1;ak = x1*exp(-j*k*w0*t’)*dt/T;需要强调的是，时间变量的变化步长dt 的大小对傅里叶级数系数的计算精度的影响非常大，dt 越小，精度越高，但是，计算机计算所花的时间越长。
5、用Matlab 实现CTFT 及其逆变换的计算5.1 用Matlab 实现CTFT 的计算Matlab 进行傅里叶变换有两种方法，一种利用符号运算的方法计算，另一种是数值计算，本实验要求采用数值计算的方法来进行傅里叶变换的计算。
严格来说，用数值计算的方法计算连续时间信号的傅里叶变换需要有个限定条件，即信号是时限信号（Time limited signal ），也就是当时间|t|大于某个给定时间时其值衰减为零或接近于零，这个条件与前面提到的为什么不能用无限多个谐波分量来合成周期信号的道理是一样的。
计算机只能处理有限大小和有限数量的数。
采用数值计算算法的理论依据是：⎰∞∞-Ω-=dt et x j X tj )()(ω∑∞-∞=-→=k Tjk T T ekT x ω)(lim若信号为时限信号，当时间间隔T 取得足够小时，上式可演变为：∑-=-=NNk Tjk ekT x Tj X ωω)()(T e e e t x t x t x N t j t j t j N ],,,[)](,),(),([12211221+---+⋅=ωωωΛΛ上式用Matlab 表示为：X=x*exp(j*t’*w)*T其中X 为信号x(t)的傅里叶变换，w 为频率Ω，T 为时间步长。
相应的Matlab 程序：T = 0.01; dw = 0.1; %时间和频率变化的步长 t = -10:T:10;w = -4*pi:dw:4*pi;X(j ω)可以按照下面的矩阵运算来进行：X=x*exp(-j*t’*ω)*T; %傅里叶变换 X1=abs(X); %计算幅度谱 phai=angle(X); %计算相位谱为了使计算结果能够直观地表现出来，还需要用绘图函数将时间信号x(t)，信号的幅度谱|X(j ω)|和相位谱∠ X(j ω)分别以图形的方式表现出来，并对图形加以适当的标注。
四、实验步骤及容6、编写程序，绘制下面的信号的波形图：Λ-+-=)5cos(51)3cos(31)cos()(000t t t t x ωωω∑∞==10)cos()2sin(1n t n n nωπ其中，0 = 0.5π，要求将一个图形窗口分割成四个子图，分别绘制cos(0t)、cos(30t)、cos(50t) 和x(t) 的波形图，给图形加title，网格线和x坐标标签。
dt=0.002;%定义时间间隔t=-2:dt:2;%设置时间起始和截止w0=0.5*pi;%题目要求，定义wox1=cos(w0*t);%输入信号x1subplot(221);%将信号1的图形放在两行两列中的第一个plot(t,x1),grid on;%输出信号1的图形，并加网格x2=cos(3*w0*t);%输入信号x2subplot(222);%将信号2的图形放在两行两列中的第二个plot(t,x2),grid on;%输出信号2的图形，并加网格x3=cos(5*w0*t);%输入信号3subplot(223);%将信号3的图形放在两行两列中的第三个plot(t,x3),grid on;%输出信号3的图形，并加网格ak=0;%设傅里叶系数的初始值为0for n=1:9%循环计算9次傅里叶系数ak=ak+(-1).^(n+1).*(1/(2*n-1)).*cos((2*n-1).*w0*t);%计算傅里叶系数endsubplot(224);%将该图形放在两行两列中的第四个plot(t,ak),grid on;%输入该图形，以t为横坐标，ak为纵坐标，并加网格7、给定如下两个周期信号：编写程序，计算x 1(t)和x 2(t)的傅里叶级数的系数，分别写出ak 从-10到10共21个系数；仿照程序Program2_2，计算并绘制出原始信号x 1(t)和x 2(t)的波形图，画出用有限项级数合成的y 1(t)和y 2(t)的波形图，以及x 1(t) 和x 2(t)的幅度频谱和相位频谱的谱线图。