非周期信号(方波,锯齿波,三角波)的合成分解以及频谱分析的MATLAB实现
- 格式：wps
- 大小：211.50 KB
- 文档页数：23
MATLAB 在数学类科技应用软件中在数值计算方面首屈一指。MATLAB 可以进行矩 阵运算、绘制函数和数据、实现算法、创建用户界面、连接其他编程语言的程序等， 主要应用于工程计算、控制设计、信号处理与通讯、图像处理、信号检测、金融建模 设计与分析等领域。
1.2 主要功能
1.数值分析 2.数值和符号计算 3.工程与科学绘图 4.控制系统的设计与仿真 5.数字图像处理 6.数字信号处理 7.通讯系统设计与仿真 8.财务与金融工程
2
连续周期信号的傅立叶级数分析及其 MATLAB 实现
２ 连续周期信号的傅立叶级数
频域分析法即傅里叶分析法，它是变换域分析法的基石。其中，傅里叶级数 是变换域分析法的理论基础，傅里叶变换作为频域分析法的重要数学工具，具有 明确的物理意义，在不同的领域得到广泛的应用。
2.1 连续时间周期信号的分解
以高等数学的知识，任何周期为 T 的周期函数，在满足狄里赫利条件时，则 该周期信号可以展开成傅里叶级数。傅里叶级数有三角形式和指数形式两种。
(3-2)
2
2
3
MATLAB 实现程序：
n=7;
6
连续周期信号的傅立叶级数分析及其 MATLAB 实现
T0=2;A=2; T1=2; tn_i=1; for tn=0:0.01:T1*T0
y_t(tn_i)=A* rem (tn,T0)/T0; t_t(tn_i)=tn; tn_i=tn_i+1; end; t=0:0.01:T1*T0; x=A/2; pi=3.1415926; w0=2*pi/T0; for i=1:n fw(i)=i*w0; a(i)=-A/(pi*i); y(i,:)=a(i)*sin(fw(i)*t); x=x+y(i,:); end; subplot(1,3,1); plot(t_t,[y_t;x]); title('锯齿波、锯齿波合成图') subplot(1,3,2); plot(t,[x; y]); title('0-n 次谐波及合成图') subplot(1,3,3); stem(fw,a); title('锯齿波频谱图') 生成图形：
1.１.１三角形式的傅里叶级数
f(t)a 2 0 a 1co t) s ( a 2co t) s ( a 3co t) s ( b 1 si n t) (b 2si n t) (
a a 0 co n s t) ( sin n t)(
b 2 n n 1
n n 1
n 1 ,2 ,3 (2 1 )
连续周期信号的傅立叶级数分析及其 MATLAB 实现
连续周期信号的傅立叶级数分析及其 MATLAB 实现 摘要
现在，MATLAB 已经发展成为适合多学科的大型软件，在世界各高校，MALAB 已经 成为线性代数、数值分析、数理统计、优化方法、自动控制、数字信号处理、动态系 统仿真等高级课程的基本应用工具。本次课程设计则在深入研究连续时间信号傅里叶 级数分析理论知识的基础上，利用 MATLAB 强大的图形处理功能、符号运算功能以及 数值计算功能，通过 MATLAB 编程进行图形功能仿真，从而实现连续时间周期信号频 域分析的仿真波形，包括以下内容：用 MATLAB 实现周期信号的傅里叶级数分解与综 合的波形；用 MATLAB 实现周期信号的单边频谱及双边频谱的波形与分析；用 MATLAB 实现典型周期信号的频谱的波形。
5
连续周期信号的傅立叶级数分析及其 MATLAB 实现
plot(t_t,[y_t;x]); axis([0 4 -3 3]); title('方波、方波合成图') subplot(3,1,2); plot(t,[x; y]); title('0-n 次谐波及合成图') subplot(3,1,3); stem(fw,a); title('方波幅频谱图')
生成图形：
方波、方波合成图
2
0
-2
0
0.5
1
1.5
2
2.5
3
3.5
4
0-n次 谐 波 及 合 成 图 5
0
-5
0
0.5
1
1.5
2
2.5
3
3.5
4
方波幅频谱图 4
2
0
0
5
10
15
20
25
30
35
40
45
图 3-1
3.2 锯齿波
周期锯齿波的傅立叶级数展开式：
u(t) A A (sin wt 1 sin 2wt 1 sin 3wt )
2 2
弦分量合并，则三角形式的傅里叶级数可表示为：
ftA 20n 1Ancosn tn
n1,2,3,
2 3
上式中
A0 An
a0 an2
bn2
,n
1tan
an bn
2 4
3
连续周期信号的傅立叶级数分析及其 MATLAB 实现
a0 A0
an An cosn,n 1,2,
1.1 简介 ..............................................1 1.2 主要功能..........................................1 1.3 特点.............................................2 ２ 连续周期信号的傅立叶级数 .............................3 2.1 连续时间周期信号的分解 ............................3
bn An sinn
可以看出，傅里叶系数 a n 和 b n 都是 n 或 n 的函数， 2 5
（2-5） 其中 a n 和 A n
是n 或 n 的偶函数，即有a n a n ；而b n 和 n 是n 或 n 的奇函数，即有b n bn。
2.１.２指数形式的傅里叶级数
根据欧拉公式：
cos(nt
n )
1 [e j(ntn ) 2
e ] j(ntn )
并考虑和奇偶性可将改写为指数形式的傅里叶级数：
f t Fne jnt , n 0, 1, 2, 3, n
(2 6)
即周期信号可分解为一系列不同频率的虚指数信号之和，式中 Fn 称为傅里
叶复系数，可由下式求得：
Fn
1 T
T
2f
T
2
t e jntdt
2.1.3狄里赫里条件
2 7
1.在一个周期内只有有限个间断点；
2.在一个周期内有有限个极值点；
3.在一个周期内函数绝对可积，即
一t般t00周期T信1号都f满足(这t 些)条.dt
（2-8）
4
连续周期信号的傅立叶级数分析及其 MATLAB 实现
3 连续周期信号的傅里叶级数分解与合成的 MATLAB 实现 3.1 方波
MATLAB 是一种科学计算软件，主要适用于矩阵计算及控制和信息处理领域 的分析设计。它使用方便，输入简捷，运算效率高，内容丰富，并且很容易由用户自 行扩展。Matlab 当前已成为美国和其他发达国家在大学教学和科学研究中最常用而 必不可少的工具
MATLAB 需要较多的高等数学的知识，要随着年级的增加，随着 matlab 在专业技术上应用才能逐渐深入掌握
连续周期信号的傅立叶级数分析及其 MATLAB 实现
1 MATLAB 1.1 简介
MATLAB 是由美国 mathworks 公司发布的主要面对科学计算、可视化以及交互 式程序设计的高科技计算环境。它将数值分析、矩阵计算、科学数据可视化以及非线 性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中，为科学 研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方 案，并在很大程度上摆脱了传统非交互式程序设计语言（如 C、Fortran）的编辑模 式，代表了当今国际科学计算软件的先进水平。
式中系数 a n 、 b n 称为傅里叶系数，可由下式求得
2 T
a 0
1 T
T
2f
T
2
t dt
a n
2 T
T
2
T
f
t cos n t dt
2
f (t) a0/ 2
b n
2 T
T
2
T
f
t sin n t dt
2
其中，为基波频率, n为 n 次谐波频率。如果将 2 1 式中同频率的正弦和余
周期方波的傅里叶级数展开式：
（ 3-1）
由式( 3-1) 可知，周期锯齿波信号，可分解为一个常分量和多个成谐波关系的正 弦波信号，或者说，可用一个常分量和多个成谐波关系的正弦波信号叠加逼近．
MATLAB 实现程序： n=7; T0=2;A=2; T1=2; tn_i=1; for tn=0:0.01:T1*T0
1
连续周期信号的傅立叶级数分析及其 MATLAB 实现
1.3 特点
语言简洁，编程效率高 因为 MATLAB 定义了专门用于矩阵运算的运算符，使得矩阵运算很简单。对于这 些运算，几条语句即可代替数十条甚至上百条 C 语言及 Fortran 语言的功能。 交互性能好，使用方便 在 MATLAB 的命令窗口中，输入一条命令，立即就能看到该命令的执行结果，体 现了良好的交互性。交互方式减少了编程和调试程序的工作量，给使用者带来了极大 的方便。因为不像使用 C 语言和 Fortran 语言那样，首先编写源程序，然后对其进行 编译、连接、待形成可执行文件后，方可运行程序，得出结果。 强大的绘图功能，便于数据可视化 MATLAB 提供了两个层次的图形命令：一种是对图形句柄进行的低级图形命令， 另一种是建立在低级图形命令之上的高级图形命令。利用 MATLAB 的高级图形命令可 以轻而易举地绘制二维、三维乃至四维图形，并可进行图形和坐标的标识、视角和光 照设计、色彩精细控制等等。 学科众多，领域广泛的应用工具箱 其工具箱分为两大类：功能性工具箱和学科性工具箱。功能性工具箱主要用来扩 充其符号计算功能、可视建模仿真功能及文字处理功能等。学科性工具箱专业性比较 强，如控制系统工具箱、信号处理工具箱、神经网络工具箱、最优化工具箱、金融工 具箱等，用户可以直接利用这些工具箱进行相关领域的科学研究。 开放性好，易于扩充 除内部函数外，MATLAB 的其他文件都是公开的、可读可改的源文件，体现了 MATLAB 的开放性特点。用户可修改源文件和加入自己的文件，甚至可以构造自己的 工具箱
1.2 主要功能
1.数值分析 2.数值和符号计算 3.工程与科学绘图 4.控制系统的设计与仿真 5.数字图像处理 6.数字信号处理 7.通讯系统设计与仿真 8.财务与金融工程
2
连续周期信号的傅立叶级数分析及其 MATLAB 实现
２ 连续周期信号的傅立叶级数
频域分析法即傅里叶分析法，它是变换域分析法的基石。其中，傅里叶级数 是变换域分析法的理论基础，傅里叶变换作为频域分析法的重要数学工具，具有 明确的物理意义，在不同的领域得到广泛的应用。
2.1 连续时间周期信号的分解
以高等数学的知识，任何周期为 T 的周期函数，在满足狄里赫利条件时，则 该周期信号可以展开成傅里叶级数。傅里叶级数有三角形式和指数形式两种。
(3-2)
2
2
3
MATLAB 实现程序：
n=7;
6
连续周期信号的傅立叶级数分析及其 MATLAB 实现
T0=2;A=2; T1=2; tn_i=1; for tn=0:0.01:T1*T0
y_t(tn_i)=A* rem (tn,T0)/T0; t_t(tn_i)=tn; tn_i=tn_i+1; end; t=0:0.01:T1*T0; x=A/2; pi=3.1415926; w0=2*pi/T0; for i=1:n fw(i)=i*w0; a(i)=-A/(pi*i); y(i,:)=a(i)*sin(fw(i)*t); x=x+y(i,:); end; subplot(1,3,1); plot(t_t,[y_t;x]); title('锯齿波、锯齿波合成图') subplot(1,3,2); plot(t,[x; y]); title('0-n 次谐波及合成图') subplot(1,3,3); stem(fw,a); title('锯齿波频谱图') 生成图形：
1.１.１三角形式的傅里叶级数
f(t)a 2 0 a 1co t) s ( a 2co t) s ( a 3co t) s ( b 1 si n t) (b 2si n t) (
a a 0 co n s t) ( sin n t)(
b 2 n n 1
n n 1
n 1 ,2 ,3 (2 1 )
连续周期信号的傅立叶级数分析及其 MATLAB 实现
连续周期信号的傅立叶级数分析及其 MATLAB 实现 摘要
现在，MATLAB 已经发展成为适合多学科的大型软件，在世界各高校，MALAB 已经 成为线性代数、数值分析、数理统计、优化方法、自动控制、数字信号处理、动态系 统仿真等高级课程的基本应用工具。本次课程设计则在深入研究连续时间信号傅里叶 级数分析理论知识的基础上，利用 MATLAB 强大的图形处理功能、符号运算功能以及 数值计算功能，通过 MATLAB 编程进行图形功能仿真，从而实现连续时间周期信号频 域分析的仿真波形，包括以下内容：用 MATLAB 实现周期信号的傅里叶级数分解与综 合的波形；用 MATLAB 实现周期信号的单边频谱及双边频谱的波形与分析；用 MATLAB 实现典型周期信号的频谱的波形。
5
连续周期信号的傅立叶级数分析及其 MATLAB 实现
plot(t_t,[y_t;x]); axis([0 4 -3 3]); title('方波、方波合成图') subplot(3,1,2); plot(t,[x; y]); title('0-n 次谐波及合成图') subplot(3,1,3); stem(fw,a); title('方波幅频谱图')
生成图形：
方波、方波合成图
2
0
-2
0
0.5
1
1.5
2
2.5
3
3.5
4
0-n次 谐 波 及 合 成 图 5
0
-5
0
0.5
1
1.5
2
2.5
3
3.5
4
方波幅频谱图 4
2
0
0
5
10
15
20
25
30
35
40
45
图 3-1
3.2 锯齿波
周期锯齿波的傅立叶级数展开式：
u(t) A A (sin wt 1 sin 2wt 1 sin 3wt )
2 2
弦分量合并，则三角形式的傅里叶级数可表示为：
ftA 20n 1Ancosn tn
n1,2,3,
2 3
上式中
A0 An
a0 an2
bn2
,n
1tan
an bn
2 4
3
连续周期信号的傅立叶级数分析及其 MATLAB 实现
a0 A0
an An cosn,n 1,2,
1.1 简介 ..............................................1 1.2 主要功能..........................................1 1.3 特点.............................................2 ２ 连续周期信号的傅立叶级数 .............................3 2.1 连续时间周期信号的分解 ............................3
bn An sinn
可以看出，傅里叶系数 a n 和 b n 都是 n 或 n 的函数， 2 5
（2-5） 其中 a n 和 A n
是n 或 n 的偶函数，即有a n a n ；而b n 和 n 是n 或 n 的奇函数，即有b n bn。
2.１.２指数形式的傅里叶级数
根据欧拉公式：
cos(nt
n )
1 [e j(ntn ) 2
e ] j(ntn )
并考虑和奇偶性可将改写为指数形式的傅里叶级数：
f t Fne jnt , n 0, 1, 2, 3, n
(2 6)
即周期信号可分解为一系列不同频率的虚指数信号之和，式中 Fn 称为傅里
叶复系数，可由下式求得：
Fn
1 T
T
2f
T
2
t e jntdt
2.1.3狄里赫里条件
2 7
1.在一个周期内只有有限个间断点；
2.在一个周期内有有限个极值点；
3.在一个周期内函数绝对可积，即
一t般t00周期T信1号都f满足(这t 些)条.dt
（2-8）
4
连续周期信号的傅立叶级数分析及其 MATLAB 实现
3 连续周期信号的傅里叶级数分解与合成的 MATLAB 实现 3.1 方波
MATLAB 是一种科学计算软件，主要适用于矩阵计算及控制和信息处理领域 的分析设计。它使用方便，输入简捷，运算效率高，内容丰富，并且很容易由用户自 行扩展。Matlab 当前已成为美国和其他发达国家在大学教学和科学研究中最常用而 必不可少的工具
MATLAB 需要较多的高等数学的知识，要随着年级的增加，随着 matlab 在专业技术上应用才能逐渐深入掌握
连续周期信号的傅立叶级数分析及其 MATLAB 实现
1 MATLAB 1.1 简介
MATLAB 是由美国 mathworks 公司发布的主要面对科学计算、可视化以及交互 式程序设计的高科技计算环境。它将数值分析、矩阵计算、科学数据可视化以及非线 性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中，为科学 研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方 案，并在很大程度上摆脱了传统非交互式程序设计语言（如 C、Fortran）的编辑模 式，代表了当今国际科学计算软件的先进水平。
式中系数 a n 、 b n 称为傅里叶系数，可由下式求得
2 T
a 0
1 T
T
2f
T
2
t dt
a n
2 T
T
2
T
f
t cos n t dt
2
f (t) a0/ 2
b n
2 T
T
2
T
f
t sin n t dt
2
其中，为基波频率, n为 n 次谐波频率。如果将 2 1 式中同频率的正弦和余
周期方波的傅里叶级数展开式：
（ 3-1）
由式( 3-1) 可知，周期锯齿波信号，可分解为一个常分量和多个成谐波关系的正 弦波信号，或者说，可用一个常分量和多个成谐波关系的正弦波信号叠加逼近．
MATLAB 实现程序： n=7; T0=2;A=2; T1=2; tn_i=1; for tn=0:0.01:T1*T0
1
连续周期信号的傅立叶级数分析及其 MATLAB 实现
1.3 特点
语言简洁，编程效率高 因为 MATLAB 定义了专门用于矩阵运算的运算符，使得矩阵运算很简单。对于这 些运算，几条语句即可代替数十条甚至上百条 C 语言及 Fortran 语言的功能。 交互性能好，使用方便 在 MATLAB 的命令窗口中，输入一条命令，立即就能看到该命令的执行结果，体 现了良好的交互性。交互方式减少了编程和调试程序的工作量，给使用者带来了极大 的方便。因为不像使用 C 语言和 Fortran 语言那样，首先编写源程序，然后对其进行 编译、连接、待形成可执行文件后，方可运行程序，得出结果。 强大的绘图功能，便于数据可视化 MATLAB 提供了两个层次的图形命令：一种是对图形句柄进行的低级图形命令， 另一种是建立在低级图形命令之上的高级图形命令。利用 MATLAB 的高级图形命令可 以轻而易举地绘制二维、三维乃至四维图形，并可进行图形和坐标的标识、视角和光 照设计、色彩精细控制等等。 学科众多，领域广泛的应用工具箱 其工具箱分为两大类：功能性工具箱和学科性工具箱。功能性工具箱主要用来扩 充其符号计算功能、可视建模仿真功能及文字处理功能等。学科性工具箱专业性比较 强，如控制系统工具箱、信号处理工具箱、神经网络工具箱、最优化工具箱、金融工 具箱等，用户可以直接利用这些工具箱进行相关领域的科学研究。 开放性好，易于扩充 除内部函数外，MATLAB 的其他文件都是公开的、可读可改的源文件，体现了 MATLAB 的开放性特点。用户可修改源文件和加入自己的文件，甚至可以构造自己的 工具箱
