实验二:连续时间信号的频域分析
连续周期信号的频域分析
三、周期信号的频谱及其特点
3. 频谱的特性
(3) 信号的有效带宽
0~2 / 这段频率范围称为周期矩形脉冲信号的 有效频带宽度,即 2π B
信号的有效带宽与信号时域的持续时间成反比。 即 越大,其B越小;反之, 越小,其B 越大。
三、周期信号的频谱及其特点
3. 频谱的特性
(3) 信号的有效带宽 物理意义:在信号的有效带宽内,集中了信 号绝大部分谐波分量。若信号丢失有效带宽以 外的谐波成分,不会对信号产生明显影响。
n=—4 4
1 T /2 2 P T / 2 f (t )dt 0.2 T 包含在有效带宽(0 ~ 2 / )内的各谐波平均功率为
2 2 C0
2 | Cn | 2 0.1806
n=1
4
P 0.1806 1 90% P 0.200
例3 试求周期矩形脉冲信号在其有效带宽(0~2 /t)内
频谱的特性频谱的特性信号的有效带宽信号的有效带宽这段频率范围称为周期矩形脉冲信号的有效频带宽度有效频带宽度即信号的有效带宽与信号时域的持续时间信号的有效带宽与信号时域的持续时间成反比
连续周期信号的频域分析
周期信号的傅里叶级数展开 傅里叶级数的基本性质 周期信号的频谱及其特点 周期信号的功率谱
三、周期信号的频谱及其特点
三、周期信号的频谱及其特点
4. 相位谱的作用
幅频不变,零相位
幅频为常数,相位不变
四、周期信号的功率谱
帕什瓦尔(Parseval)功率守恒定理
2 1 T P 2T f (t ) dt Cn T 2 n 2
物理意义:任意周期信号的平均功率等于信号所 包含的直流、基波以及各次谐波的平均功率之和。
MATLAB信号与系统实验报告19472[五篇范文]
MATLAB信号与系统实验报告19472[五篇范文]第一篇:MATLAB信号与系统实验报告19472信号与系统实验陈诉(5)MATLAB 综合实验项目二连续系统的频域阐发目的:周期信号输入连续系统的响应可用傅里叶级数阐发。
由于盘算历程啰嗦,最适适用MATLAB 盘算。
通过编程实现对输入信号、输出信号的频谱和时域响应的盘算,认识盘算机在系统阐发中的作用。
任务:线性连续系统的系统函数为11)(+=ωωjj H,输入信号为周期矩形波如图 1 所示,用MATLAB 阐发系统的输入频谱、输出频谱以及系统的时域响应。
-3-2-1 0 1 2 300.511.52Time(sec)图 1要领:1、确定周期信号 f(t)的频谱nF&。
基波频率Ω。
2、确定系统函数 )(Ω jn H。
3、盘算输出信号的频谱n nF jn H Y&&)(Ω=4、系统的时域响应∑∞-∞=Ω=nt jnn eY t y&)(MATLAB 盘算为y=Y_n*exp(j*w0*n“*t);要求(画出 3 幅图):1、在一幅图中画输入信号f(t)和输入信号幅度频谱|F(jω)|。
用两个子图画出。
2、画出系统函数的幅度频谱|H(jω)|。
3、在一幅图中画输出信号y(t)和输出信号幅度频谱|Y(jω)|。
用两个子图画出。
解:(1)阐发盘算:输入信号的频谱为(n)输入信号最小周期为=2,脉冲宽度,基波频率Ω=2π/ =π,所以(n)系统函数为因此输出信号的频谱为系统响应为(2)步伐:t=linspace(-3,3,300);tau_T=1/4;%n0=-20;n1=20;n=n0:n1;%盘算谐波次数20F_n=tau_T*Sa(tau_T*pi*n);f=2*(rectpuls(t+1.75,0.5)+rectpuls(t-0.25,0.5)+rectpuls(t-2.25,0.5));figure(1),subplot(2,1,1),line(t,f,”linewidth“,2);%输入信号的波形 axis([-3,3,-0.1,2.1]);grid onxlabel(”Time(sec)“,”fontsize“,8),title(”输入信号“,”fontweight“,”bold“)%设定字体巨细,文本字符的粗细text(-0.4,0.8,”f(t)“)subplot(2,1,2),stem(n,abs(F_n),”.“);%输入信号的幅度频谱xlabel(”n“,”fontsize“,8),title(”输入信号的幅度频谱“,”fontweight“,”bold“)text(-4.0,0.2,”|Fn|“)H_n=1./(i*n*pi+1);figure(2),stem(n,abs(H_n),”.“);%系统函数的幅度频谱xlabel(”n“,”fontsize“,8),title(”系统函数的幅度频谱“,”fontweight“,”bold“)text(-2.5,0.5,”|Hn|“)Y_n=H_n.*F_n;y=Y_n*exp(i*pi*n”*t);figure(3),subplot(2,1,1),line(t,y,“linewidth”,2);%输出信号的波形 axis([-3,3,0,0.5]);grid onxlabel(“Time(sec)”,“fontsize”,8),title(“输出信号”,“fontweight”,“bold”)text(-0.4,0.3,“y(t)”)subplot(2,1,2),stem(n,abs(Y_n),“.”);%输出信号的幅度频谱xlabel(“n”,“fontsize”,8),title(“输出信号的幅度频谱”,“fontweight”,“bold”)text(-4.0,0.2,“|Yn|”)(3)波形:-3-2-1 0 1 2 300.511.52Time(sec)输入信号f(t)-20-15-10-5 0 5 10 15 2000.10.20.30.4n输入信号的幅度频谱|Fn|-20-15-10-5 0 5 10 15 2000.10.20.30.40.50.60.70.80.91n系统函数的幅度频谱|Hn|-3-2-1 0 1 2 300.10.20.30.4Time(sec)输出信号y(t)-20-15-10-5 0 5 10 15 2000.10.20.30.4n输出信号的幅度频谱|Yn| 项目三连续系统的复频域阐发目的:周期信号输入连续系统的响应也可用拉氏变更阐发。
通信原理实验教程(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
信号与系统实验报告实验一信号与系统的时域分析一、实验目的1、熟悉和掌握常用的用于信号与系统时域仿真分析的MA TLAB函数;2、学会用MA TLAB进行信号基本运算的方法;3、掌握连续时间和离散时间信号的MA TLAB产生,掌握用周期延拓的方法将一个非周期信号进行周期信号延拓形成一个周期信号的MATLAB编程。
二、实验内容Q1-1:修改程序Program1_1,将dt改为0.2,再执行该程序,保存图形,看看所得图形的效果如何?dt = 0.01时的程序clear, % Clear all variablesclose all, % Close all figure windowsdt = 0.01; % Specify the step of time variablet = -2:dt:2; % Specify the interval of timex = sin(2*pi*t); % Generate the signalplot(t,x) % Open a figure window and draw the plot of x(t)title('Sinusoidal signal x(t)')xlabel('Time t (sec)')dt = 0.2时的程序clear, % Clear all variablesclose all, % Close all figure windowsdt = 0.2; % Specify the step of time variablet = -2:dt:2; % Specify the interval of timex = sin(2*pi*t); % Generate the signalplot(t,x) % Open a figure window and draw the plot of x(t)title('Sinusoidal signal x(t)')xlabel('Time t (sec)')dt = 0.01时的信号波形dt = 0.2时的信号波形这两幅图形有什么区别,哪一幅图形看起来与实际信号波形更像?答:dt = 0.01的图形比dt = 0.2的图形光滑,dt = 0.01看起来与实际信号波形更像。
信号与系统实验教程只有答案
信号与系统实验教程(只有答案))(实验报告目录实验一信号与系统的时域分析 (2)三、实验内容及步骤 (2)实验二连续时间信号的频域分析 (14)三、实验内容及步骤 (14)实验三连续时间LTI系统的频域分析 (35)三、实验内容及步骤 (35)实验四通信系统仿真 (42)三、实验内容及步骤 (42)实验五连续时间LTI系统的复频域分析 (51)三、实验内容及步骤 (51)实验一信号与系统的时域分析三、实验内容及步骤实验前,必须首先阅读本实验原理,读懂所给出的全部范例程序。
实验开始时,先在计算机上运行这些范例程序,观察所得到的信号的波形图。
并结合范例程序应该完成的工作,进一步分析程序中各个语句的作用,从而真正理解这些程序。
实验前,一定要针对下面的实验项目做好相应的实验准备工作,包括事先编写好相应的实验程序等事项。
Q1-1:修改程序Program1_1,将dt改为0.2,再执行该程序,保存图形,看看所得图形的效果如何?dt = 0.01时的信号波形dt = 0.2时的信号波形这两幅图形有什么区别,哪一幅图形看起来与实际信号波形更像?答:Q1-2:修改程序Program1_1,并以Q1_2为文件名存盘,产生实指数信号x(t)=e-0.5t。
要求在图形中加上网格线,并使用函数axis()控制图形的时间范围在0~2秒之间。
然后执行该程序,保存所的图形。
修改Program1_1后得到的程序Q1_2如下:信号x(t)=e-0.5t的波形图clear, % Clear all variablesclose all, % Close all figure windowsdt = 0.2; % Specify the step of time variablet = -2:dt:2; % Specify the interval of timex = exp(-0.5*t); % Generate the signalplot(t,x)grid on;axis ([0 2 0 1 ])title('Sinusoidal signal x(t)')xlabel('Time t (sec)')Q1-3:修改程序Program1_1,并以Q1_3为文件名存盘,使之能够仿真从键盘上任意输入的一个连续时间信号,并利用该程序仿真信号x(t)=e-2t。
实验二---连续时间信号的频域分析
实验二---连续时间信号的频域分析实验目的:1. 学习连续时间信号的频域分析方法,掌握傅里叶变换理论。
2. 理解信号的时域与频域之间的转换关系,能够实现信号的频域分析及某些信号处理操作。
3. 了解傅里叶变换的性质和应用,能够应用傅里叶变换对各种周期和非周期信号进行分析。
实验原理:1. 傅里叶变换傅里叶变换是将一个连续时间函数在频域中的频谱与该函数在时域中的波形进行对应的数学变换。
连续时间傅里叶变换(CTFT)是将一个无限长但可积的信号,即绝对可积信号,变换为复频域函数。
如果傅里叶变换是定义在时域上的,那么它的自变量是时间t,而它的函数值是一个关于f的复合函数,即分别为实频谱与虚频谱的函数。
- 傅里叶变换是一个线性变换;- 时域中的卷积在频域中对应为乘积;- 频域中的卷积在时域中对应为乘积;- 时域中的移位在频域中对应为复制效应;- 能量守恒:信号在时域中的总能量等于在频域中的总能量;- Parseval定理:信号在时域和频域中的幅度平方和等于常数。
实验步骤:1)连续时间正弦波$f(t)=A sin(2\pi f_0 t)$其中,$f_0 =1200 Hz$,采样间隔 $\Delta t =5*10^{-6}$ s,数据长度 $N= 150$。
$f(t)=\frac{2A}{T_0} t$($-\frac{T_0}{2}<t<\frac{T_0}{2}$)其中,$T_0$ 为周期,数据长度 $N= 500$。
$f(t) =\frac{A}{2}[sgn(t)+1]$($-1<t<1$)绘制信号的频域幅度谱和相位谱,并分析其特点。
实验结果:正弦波:三角波:方波:实验分析:从时域波形可以看出,正弦信号为一定频率下的振荡信号,具有周期性,幅度相等,相位差为 $\frac{\pi}{2}$ 的两个正弦函数相加而成;三角波和方波均为非周期信号。
从频域幅度谱可以看出,正弦波在频域中只存在一个正弦函数,且其频率与时域信号的频率相同;三角波在频域中存在多个频率成分,且成分包含奇数倍或基波的奇数倍;方波在频域中由越来越多的奇数倍频率成分组成,其频率分量越高,能量越小。
实验二 测试系统的时域响应和频域响应
实验一测试系统的时域响应【实验目的】1.了解MATLAB软件的基本特点和功能,熟悉其界面、菜单和工具条,熟悉MATLAB程序设计结构及M文件的编制;2.掌握线性系统模型的计算机表示方法;3.掌握求线性定常连续系统时域输出响应的方法,求得系统的时域响应曲线;4. 了解Simulink 的使用。
【实验指导】一、模型的建立:在线性系统理论中,一般常用的数学模型形式有:(1)传递函数模型;(2)状态空间模型;(3)零极点增益模型这些模型之间都有着内在的联系,可以相互进行转换.1、传递函数模型若已知系统的传递函数为:对线性定常系统,式中s的系数均为常数,且an不等于零,这时系统在MATLAB中可以方便地由分子和分母系数构成的两个向量唯一地确定出来,这两个向量分别用num和den表示.num=[cm,c,m-1,…,c1,c0]den=[an,an-1,…,a1,a0]注意:它们都是按s的降幂进行排列的.则传递函数模型建立函数为:sys=tf(num,den).2、零极点增益模型(略)3、状态空间模型(略)二、模型的转换在一些场合下需要用到某种模型,而在另外一些场合下可能需要另外的模型,这就需要进行模型的转换.三、模型的连接1、并联:parallel[num,den]=parallel(num1,den1,num2,den2)%将并联连接的传递函数进行相加.2、串联:series[num,den]=series(num1,den1,num2,den2)%将串联连接的传递函数进行相乘.3、反馈:feedback[num,den]=feedback(num1,den1,num2,den2,sign)%可以得到类似的连接,只是子系统和闭环系统均以传递函数的形式表示.当sign=1时采用正反馈;当sign= -1时采用负反馈;sign缺省时,默认为负反馈.4、闭环:cloop(单位反馈)[numc,denc]=cloop(num,den,sign)%表示由传递函数表示的开环系统构成闭环系统,sign意义与上述相同.四、线性连续系统的时域响应1 求取线性连续系统的阶跃响应函数为(step) 基本格式为:step(sys) step(num,den)【实验内容】1. 典型一阶系统的传递函数为 11)(+=s s G τ;τ为时间常数,试绘出当τ=0.5、1、 2、4、6、8、时该系统的单位阶跃响应曲线。
《信号、系统与数字信号处理》第二章 连续时间信号与系统的频域分析
0 21
/4
/2
(b)相位图
图2.1-2例2.1-2的频谱图
二、指数形式的傅里叶级数
利用欧拉公式将三角形式的傅里叶级数,表示为 复指数形式的傅氏级数
其中
f t F n1 e jn1t
n
F n1
1 T
t0 T t0
f t e jn1tdt
F n1 是复常数,通常简写为 Fn 。
21t
5
4
2
sin
1t
1 2
sin
31t
解:将 f t 整理为标准形式
f
(t)
1
2cos 1t来自4cos 21t
5
4
1 2
cos
31t
2
1
2
cos
1t
4
cos
21t
4
1 2
cos
31t
2
振幅谱与相位谱如图2-1所示。
cn
2
1
1
1/2
0 1 21 31
(a) 振幅图
n
/4
31
第二章 连续时间信号与系统的频域分析 ——Fourier变换
2. 1 周期信号的傅里叶级数分析 2. 2 非周期信号的频谱--傅里叶变换 2. 3 傅里叶变换的性质及定理 2. 4 系统的频域分析方法 2. 5 无失真传输系统与滤波
LTI系统分析的一个基本任务,是求解系统对任意 激励信号的响应,基本方法是将信号分解为多个基本信 号元。
一、三角形式傅里叶级数
周期信号: f t f t nT
其中
T
是信号的最小重复时间间隔,f1
1 是信号的基波频率。 T
若 f t 满足狄里赫利条件,则 f t 可以展开为三角形
连续时间信号的频域分析
实验二 连续时间信号的频域分析 实验内容:1、周期信号的傅里叶级数与GIBBS 现象给定如一个周期信号1()x t 如图所示:Q2-1: 分别手工计算x1(t) 的傅里叶级数的系数。
信号x1(t) 在其主周期内的数学表达式为:⎩⎨⎧≤≤-=其他,02.02.0,1)(1t t x计算x1(t) 的傅里叶级数的系数的计算过程如下:解:首先,我们根据前面所给出的公式,计算该信号的傅里叶级数的系数。
⎰-=10011)(1T t jk k dt e t x T a ω⎰--=2.02.0021dt e t jk ω⎰--=-1000)(210t jk d e k j t jk ωωω 02.02.0200ωωωk j e e kj kj --=-ππk k )2.0s i n (=(ππω==T 120) 在MATLAB 命令窗口,依次键入:>> k = -10:10;>> ak = sin(0.2*(k+eps)*pi)./((k+eps)*pi)用MATLAB 帮助你计算出你手工计算的傅里叶级数的系数ak 从-10到10共21个系数。
t 22-12.02.0-从命令窗口上抄写x1(t)的21个系数如下:Columns 1 through 40.0000 -0.0208 -0.0378 -0.0432Columns 5 through 8-0.0312 0.0000 0.0468 0.1009Columns 9 through 120.1514 0.1871 0.2000 0.1871Columns 13 through 160.1514 0.1009 0.0468 0.0000Columns 17 through 20-0.0312 -0.0432 -0.0378 -0.0208Columns210.00000Q2-2:仿照程序Program2_1,编写程序Q2_2,以计算x1(t)的傅里叶级数的系数。
连续时间信号与系统的频域分析实验报告
实验四连续时间信号与系统的频域分析一、实验目的掌握连续时间信号的傅里叶变换及傅里叶逆变换的实现方法,掌握连续时间系统的频域分析方法,熟悉MATLAB 相应函数的调用格式和作用,掌握使用MATLAB 来分析连续时间信号与系统的频域特性及绘制信号频谱图的方法。
二、实验原理(一)连续时间信号与系统的频域分析原理1、连续时间信号的额频域分析 连续时间信号的傅里叶变换为:()()dt e t f j F t j ωω-∞∞-⎰=傅里叶逆变换为:()()ωωπωd e j F t f t j ⎰∞∞-=21()ωj F 称为频谱密度函数,简称频谱。
一般是复函数,可记为:()()()ωϕωωj e j F j F =()ωj F 反映信号各频率分量的幅度随频率ω的变化情况,称为信号幅度频谱。
()ωϕ反映信号各频率分量的相位随频率ω的变化情况,称为信号相位频谱。
2、连续时间系统的频域分析 在n 阶系统情况下,数学模型为:()()()()()()()()t f b dtt df b dt t f d b dt t f d b t y a dtt dy a dt t y d a dt t y d a o m m n m m n o n n n n n n ++++=++++------11111111 令初始条件为零,两端取傅里叶变换,得:()()[]()()()[]()ωωωωωωωωj F b j b j b j b j Y a j a j a j a m n m n n n nn01110111++++=++++----表示为()()()()ωωωωj F j b j Y j a kmk kkn k k∑∑===0则 ()()()()()()()()()∑∑==----=++++++++==nk kk mk kk n n n n m m mm j a j b a j a j a j a b j b j b j b j F j Y j H 0001110111ωωωωωωωωωωω3、系统传递函数 系统传递函数定义为:()()()ωωωj H j Y j H =系统传递函数反映了系统内在的固有的特性,它取决于系统自身的结构及参数,与外部 激励无关,是描述系统特性的一个重要参数。
2010信号与系统实验2
实验二:连续和离散系统的频域分析一:实验目的1:学习傅里叶正变换和逆变换,理解频谱图形的物理含义2:了解连续和离散时间系统的单位脉冲响应3:掌握连续时间系统的频率特性二:实验原理1. 傅里叶正变换和逆变换公式 正变换:()()j t F f t e dt ωω∞--∞=⎰逆变换:1()()2j t f t F e d ωωωπ∞-∞=⎰2. 频域分析t j tj e d d e t e ωωωπωωωπ⎰⎰∞∞-∞∞-E =E =)(21)(21)(将激励信号分解为无穷多个正弦分量的和。
⎰∞∞-H E =ωωωπωd e t r tj zs )()(21)(,R(ω)为)(t r zs 傅里叶变换;πωωd )(E 各频率分量的复数振幅 激励单位冲激响应时的零状态响应→ )(t δ)(t h单位阶跃响应时的零状态响应激励→)(t u )(t g3 各函数说明:(1)impulse 冲激响应函数:[Y ,X,T]=impulse(num,den);)1()2()1()1()2()1()()()(11++++++++==--n a s a s a m b s b s b s A s B s H n n m m num 分子多项式系数; num=[b(1) b(2) … b(n+1)]; den 分母多项式系数; den=[a(1) a(2) … a(n+1)];Y ,X,T 分别表示输出响应,中间状态变量和时间变量; 如:352)(2+++=s s s s H ,等价于)(2)()(3)(5)(t e t e t r t r t r +=++ 定义den=[1 5 3];num=[1 2]; [Y ,X,T]=impulse(num,den);(2)step 阶跃响应函数:[Y,X,T]=step(num,den);num 分子多项式;den 分母多项式 Y ,X,T 分别表示输出响应,中间状态变量和时间变量;如:352)(2+++=s s s s H ,den=[1 5 3];num=[1 2];[Y ,X,T]= step (num,den);(3)impz 数字滤波器的冲激响应 [h,t] = impz(b,a,n) b 分子多项式系数;a 分母多项式系数;n 采样样本h 离散系统冲激响应;t 冲激时间,其中t=[0:n-1]', n=length(t)时间样本数(4)freqs 频域响应 [h,w] = freqs(b,a,f) b,a 定义同上,f 频率点个数 h 频域响应,w 频域变量)1()2()1()1()2()1()()()(11++++++++==--m a s a s a n b s b s b s A s B s H m m n n三.实验内容1 周期信号傅里叶级数 已知连续时间信号()()2/8cos 3/4coscos )(321ππ++++=t A t A wt A t x ,其中321,,A A A 取值如下:(X 为学号的后两位)]10,1[,5.02321∈⎪⎩⎪⎨⎧===X X A X A X A ]20,11[,55321∈⎪⎩⎪⎨⎧+==-=X X A XA X A ⎪⎩⎪⎨⎧=-=-=X A X A X A 32151020,>X 要求画出信号的时域波形和频域波形(幅度谱和相位谱)。
信号_频域分析实验报告(3篇)
第1篇一、实验目的1. 理解信号的频域分析方法及其在信号处理中的应用。
2. 掌握傅里叶变换的基本原理和计算方法。
3. 学习使用MATLAB进行信号的频域分析。
4. 分析不同信号在频域中的特性,理解频域分析在实际问题中的应用。
二、实验原理频域分析是信号处理中一种重要的分析方法,它将信号从时域转换到频域,从而揭示信号的频率结构。
傅里叶变换是频域分析的核心工具,它可以将任何信号分解为不同频率的正弦波和余弦波的线性组合。
三、实验内容及步骤1. 信号生成与傅里叶变换- 使用MATLAB生成一个简单的正弦波信号,频率为50Hz,采样频率为1000Hz。
- 对生成的正弦波信号进行傅里叶变换,得到其频谱图。
2. 频谱分析- 分析正弦波信号的频谱图,观察其频率成分和幅度分布。
- 改变正弦波信号的频率和幅度,观察频谱图的变化,验证傅里叶变换的性质。
3. 信号叠加- 将两个不同频率的正弦波信号叠加,生成一个复合信号。
- 对复合信号进行傅里叶变换,分析其频谱图,验证频谱叠加原理。
4. 窗函数- 使用不同类型的窗函数(如矩形窗、汉宁窗、汉明窗等)对信号进行截取,观察窗函数对频谱的影响。
- 分析不同窗函数的频率分辨率和旁瓣抑制能力。
5. 信号滤波- 设计一个低通滤波器,对信号进行滤波处理,观察滤波器对信号频谱的影响。
- 分析滤波器对信号时域和频域特性的影响。
6. MATLAB工具箱- 使用MATLAB信号处理工具箱中的函数,如`fft`、`ifft`、`filter`等,进行信号的频域分析。
- 学习MATLAB工具箱中的函数调用方法和参数设置。
四、实验结果与分析1. 正弦波信号的频谱分析实验结果显示,正弦波信号的频谱图只有一个峰值,位于50Hz处,说明信号只包含一个频率成分。
2. 信号叠加的频谱分析实验结果显示,复合信号的频谱图包含两个峰值,分别对应两个正弦波信号的频率。
验证了频谱叠加原理。
3. 窗函数对频谱的影响实验结果显示,不同类型的窗函数对频谱的影响不同。
连续时间信号与系统的频域分析
目录
5-12 信号的时域抽样与抽样定理 5-13 调制与解调 5-14 频分复用与时分复用
4
引言
• 用时间作为变量描述信号我们称为信号的时域表示,显 示信号随时间变换的快慢、出现先后、存在时间的长短以 及信号是否按一定的时间间隔重复出现等。 • 用频率作为变量描述信号称为频域描述,揭示了信号各 个频率分量的大小,信号的能量主要集中在哪个频率范 围等特性。 • 信号的时域表示和频域表示是从信号的两个不同方面 对信号进行描述, • 在正交函数的基础上对时域信号的进行分解。最常用的 分解就是傅立叶分解,也称为信号的傅立叶分析。
能使信号 f(t)进行正交分解的基底函数,并且分解后均方 差为零的一组正交基底函数称为完备的正交函数集。
一个信号可用完备的正交函数集表示,正交函数集有许 多,如:
• 正弦函数集 • 指数函数集 • walsh函数集 • ……
正交函数集有许多重要的用途,例如进行频谱分析、信道 编码等。
13
5-2 周期信号的傅立叶级数分析
1
T
t0 T t0
f
(t) cos n1tdt
j1 T
t0 t0
T
sin
n1tdt
(5-22)
Fn
1 T
t0 T f (t)e dt j(n)1t
t0
(n取 ~ 之间的整数)
1
T
t0 T t0
f
(t)[cos n1t
j sin n1t]dt
通过比较可以得到指数形式的傅里叶系数与三角形式的傅 里叶系数有以下关系:
当周期信号 f (t)满足狄里赫利条件时,就可以用复指数函 数集或三角函数集的线性组合来表示,这种线性组合的表 示称为傅立叶级数展开。 狄里赫利条件:
连续时间信号与系统的频域分析实验报告(共9篇)
连续时间信号与系统的频域分析实验报告(共9篇)信号与系统实验五__连续时间信号的频域分析实验名称:连续时间信号的频域分析报告人:姓名班级学号一、实验目的1、熟悉傅里叶变换的性质;2、熟悉常见信号的傅里叶变换;3、了解傅里叶变换的MATLAB实现方法。
二、实验内容及运行结果1、编程实现下列信号的幅度频谱:(1)求出f(t)=u(2t+1)-u(2t-1)的频谱函数F(w);请与f1(t) u(2t+1)-u(2t-1)的频谱函数F1(w)进行比较,说明两者的关系。
%(1)f(t)=u(2t+1)-u(2t-1)与f(t)=u(t+1)-u(t-1) syms t w t1 w1Gt=sym('Heaviside(2*t+1)-Heaviside(2*t-1)');Gt1=sym('Heaviside(t1+1)-Heaviside(t1-1)');Fw=fourier(Gt,t,w);Fw1=fourier(Gt1,t1,w1);FFw=maple('convert',Fw,'piecewise');FFw1=maple('convert',Fw1,'piecewise');FFP=abs(FFw);FFP1=abs(FFw1);subplot(2,1,1);ezplot(FFP,[-10*pi 10*pi]);axis([-10*pi 10*pi 0 1.5]);subplot(2,1,2);ezplot(FFP1,[-10*pi 10*pi]);grid;axis([-10*pi 10*pi 0 2.2]);不同点:F1(w)的图像在扩展,幅值是F(w)的两倍。
(2)三角脉冲f2(t)=1-|t|;|t|=1;ft=sym('(1+t)*Heaviside(t+1)-2*t*Heaviside(t)+(t-1)*Heaviside( t-1)');Fw=fourier(ft);subplot(211)ezplot(abs(Fw)); g2)');ft=ifourier(Fw,w,t)ft =exp(-4*t)*heaviside(t)-exp(4*t)*heaviside(-t)(2)F(w)=((i*w)+5*i*w-8)/((i*w)+6*i*w+5)syms t wFw=sym('((i*w)+5*i*w-8)/((i*w)+6*i*w+5)');ft=ifourier(Fw,w,t)ft =dirac(t)+(2*exp(-5*t)-3*exp(-t))*heaviside(t)三、讨论与总论通过本实验,掌握了信号的傅里叶变换的性质以及方法,对傅里叶变换的性质有进一步的提高。
连续时间系统的频域分析
第三章.连续时间系统的频域分析一、任意信号在完备正交函数系中的表示法(§)信号分解的目的:● 将任意信号分解为单元信号之和,从而考查信号的特性。
●简化电路分析与运算,总响应=单元响应之和。
1.正交函数集任意信号)(t f 可表示为n 维正交函数之和:原函数()()()t g t g t g r Λ21,相互正交:⎩⎨⎧=≠=⋅⎰nm K nm dt t g t g m t t n m ,,0)()(21()t g r 称为完备正交函数集的基底。
一个信号可用完备的正交函数集表示,.正弦函数集有许多方便之处,如易实现等,我们主要讨论如何用正弦函数集表示信号。
2.能量信号和功率和信号(§一)设()t i 为流过电阻R 的电流,瞬时功率为R t i t P )()(2=一般说来,能量总是与某一物理量的平方成正比。
令R = 1Ω,则在整时间域内,实信号()t f 的能量,平均功率为: 讨论上述两个式子,只可能出现两种情况: ✍∞<<W 0(有限值) 0=P✍∞<<P 0(有限值)∞=W满足✍式的称为能量信号,满足✍式称功率信号。
3.帕斯瓦尔定理设{})(t g r 为完备的正交函数集,即信号的能量 基底信号的能量 各分量此式称为帕斯瓦尔定理 P331 式(6-81) (P93, P350) 左边是信号能量,右边是各正交函数的能量。
物理意义:一个信号所含有的能量(功率)恒等于此信号在完备正交函数集中各分量能量(功率)之和。
二、周期信号的频谱分析——傅里叶级数(1) 周期信号傅里叶级数有两种形式三角形式: ()∑∞=++=1110sin cos )(n n nt n b t n aa t f ωω=∑∞=++110)cos(n n nt n cc ϕω指数形式:t jn n e n F t f 1)()(1ωω∑∞-∞==(2) 周期信号的频谱是离散谱,三个性质收敛性()↓↑)(,1ωn F n谐波性:(离散性)谱线只出现在1ωn 处,唯一性:)(t f 的谱线唯一(3)两种频谱图的关系● 三角形式:ω~n c ,ωφ~n 单边频谱● 指数形式:ωω~)(1n F , ωφ~n 双边频谱两者幅度关系 )(1ωn F =()021≠n c n000a c F ==● 指数形式的幅度谱为偶函数 ●指数形式的相位谱为奇函数(4) 引入负频率对于双边频谱,负频率)(1ωn ,只有数学意义,而无物理意义。
信号与系统的实验报告(2)
信号与系统实验报告——连续时间系统的复频域分析班级:05911101学号:**********姓名:***实验五连续时间系统的复频域分析——1120111487 信息工程(实验班)蒋志科一、实验目的①掌握拉普拉斯变换及其反变换的定义,并掌握MA TLAB 实现方法 ②学习和掌握连续时间系统系统函数的定义及其复频域分析方法③掌握系统零极点的定义,加深理解系统零极点分布与系统特性的关系。
二、实验原理与方法 1、拉普拉斯变换连续时间信号x(t)的拉普拉斯变换定义为:X s =x (t )e −st dt +∞−∞拉普拉斯反变换为:x t =12πj X (s )e st ds σ+j ∞σ−j ∞在MA TLAB 中可以采用符号数学工具箱中的laplace 函数和ilaplace 函数进行拉氏变换和拉氏反变换。
L=laplace(F)符号表达式F 的拉氏变换,F 中时间变量为t ,返回变量为s 的结果表达式。
L=laplace(F,t)用t 替换结果中的变量s 。
F=ilaplace(L)以s 为变量的符号表达式L 的拉氏反变换,返回时间变量t 的结果表达式。
F=ilaplace(L,x)用x 替换结果中的变量t 。
2、连续时间系统的系统函数连续时间系统的系统函数是系统单位冲激响应的拉氏变换H s =ℎ(t )e −st dt +∞−∞此外,连续时间系统的系统函数还可以由系统输入和输出信号的拉氏变换之比得到H s =Y(s)/X(s) 单位冲激响应h(t)反映了系统的固有性质,而H(s)从复频域反映了系统的固有性质。
对于H(s)描述的连续时间系统,其系统函数s 的有理函数H s =b M s M +b M−1s M−1+⋯+b 0a n s n +a n −1s M−1+⋯+a 03、连续时间系统的零极点分析系统的零点指使式H s 的分子多项式为零的点,极点指使分母多项式为零的点,零点使系统的值为零,极点使系统函数的值无穷大。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
一、实验目的1、掌握连续时间周期信号的傅里叶级数的物理意义和分析方法;2、观察截短傅里叶级数而产生的“Gibbs 现象”,了解其特点以及产生的原因;3、掌握连续时间傅里叶变换的分析方法及其物理意义;4、学习利用MATLAB 语言编写计算CTFS 和CTFT 的仿真程序。
基本要求:掌握并深刻理傅里叶变换的物理意义,掌握信号的傅里叶变换的计算方法,掌握利用MATLAB 编程完成相关的傅里叶变换的计算。
二、实验原理及方法1、连续时间周期信号的傅里叶级数CTFS 分析任何一个周期为T 1的正弦周期信号,只要满足狄利克利条件,就可以展开成傅里叶级数。
其中三角傅里叶级数为:∑∞=++=1000)]sin()cos([)(k k k t k b t k a a t x ωω 9.1或: ∑∞=++=100)cos()(k k kt k ca t x ϕω 9.2其中102T πω=,称为信号的基本频率,k k b a a ,和,0分别是信号)(t x 的直流分量、余弦分量幅度和正弦分量幅度,k k c ϕ、为合并同频率项之后各正弦谐波分量的幅度和初相位,它们都是频率0ωk 的函数,绘制出它们与0ωk 之间的图像,称为信号的频谱图(简称“频谱”),k c -0ωk 图像为幅度谱,k ϕ-0ωk 图像为相位谱。
指数形式的傅里叶级数为:∑∞-∞==k tjk kea t x 0)(ω 9.3其中,k a 为指数形式的傅里叶级数的系数,按如下公式计算:⎰--=2/2/1110)(1T T tjk k dt e t x T a ω 9.4 假设谐波项数为N ,则上面的和成式为:∑-==NNk tjk kea t x 0)(ω 9.5显然,N 越大,所选项数越多,有限项级数合成的结果越逼近原信号x(t)。
2、连续时间信号傅里叶变换----CTFT傅里叶变换在信号分析中具有非常重要的意义,它主要是用来进行信号的频谱分析的。
傅里叶变换和其逆变换定义如下:⎰∞∞--=dt et x j X tj ωω)()( 9.6⎰∞∞-=ωωπωd e j X t x tj )(21)( 9.7 连续时间傅里叶变换主要用来描述连续时间非周期信号的频谱。
任意非周期信号,如果满足狄里克利条件,那么,它可以被看作是由无穷多个不同频率(这些频率都是非常的接近)的周期复指数信号e j ωt 的线性组合构成的,每个频率所对应的周期复指数信号e j ωt 称为频率分量其相对幅度为对应频率的|X(j ω)|之值,其相位为对应频率的X(j ω)的相位。
给定一个连续时间非周期信号x(t),它的频谱也是连续且非周期的。
对于连续时间周期信号,也可以用傅里变换来表示其频谱,其特点是,连续时间周期信号的傅里叶变换时有冲激序列构成的,是离散的——这是连续时间周期信号的傅里叶变换的基本特征。
3、连续周期信号的傅里叶级数CTFS 的MATLAB 实现 3.1 傅里叶级数的MATLAB 计算设周期信号x(t)的基本周期为T 1,且满足狄里克利条件,则其傅里叶级数的系数可由式9.4计算得到。
式9.4重写如下:⎰--=2/2/1110)(1T T tjk k dt et x T a ω基本频率为: 102T πω=对周期信号进行分析时,我们往往只需对其在一个周期内进行分析即可,通常选择主周期。
假定x 1(t)是x(t)中的主周期,则⎰--=2/2/11110)(1T T tjk k dt e t x T a ω 计算机不能计算无穷多个系数,所以我们假设需要计算的谐波次数为N ,则总的系数个数为2N+1个。
在确定了时间范围和时间变化的步长即T 1和dt 之后,对某一个系数,上述系数的积分公式可以近似为:∑⎰---==ntjk n T T t jk k T dt e t x dt e t x T a 12/2/11/)()(10110ωω 121/],,[)](),(),([02010T dt e e et 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 = in put(‘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;例1:给定一个周期为T 1 = 2s 的连续时间周期方波信号,如图所示,其一个周期内的数学表达式为:⎩⎨⎧<<≤≤=21,010,1)(1t t t x% Program2_1clear, close all T = 2; dt = 0.00001; t = -2:dt:2; x1 = u(t) - u(t-1-dt); x = 0;for m = -1:1x = x + u(t-m*T) - u(t-1-m*T-dt);endw0 = 2*pi/T;N = 10; L = 2*N+1;for k = -N: N;ak(N+1+k) = (1/T)*x1*exp(-j*k*w0*t')*dt;endphi = anglel(ak); % Evaluate the phase of ak 执行程序Program2_1后,就完成了信号的傅里叶级数的系数的计算,在命令窗口键入>> ak图9.1 周期方波信号命令窗口就可以显示傅里叶级数的21个系数:在程序Program2_1中添加相应的计算| ak |和绘图语句,就可以绘制出信号的幅度谱和相位谱的谱线图。
3.2 周期信号的合成以及Gibbs 现象从傅里叶级数的合成式(Synthesis equation )∑∞-∞==k tjk kea t x 0)(ω可以看出,用无穷多个不同频率和不同振幅的周期复指数信号可以合成一个周期信号。
然而,我们无法用计算机实现对无穷多个周期复指数信号的合成。
但是,用有限项来合成却是可行的,在实际应用中,多半也就是这么做的。
然而,这样做的一个必然结果,就是引入了误差。
如果一个周期信号在一个周期有内断点存在,那么,引入的误差将除了产生纹波之外,还将在断点处产生幅度大约为9%的过冲(Overshot ),这种现象被称为吉伯斯现象(Gibbs phenomenon )。
为了能够观察到合成信号与原信号的不同以及Gibbs 现象,我们可以利用前面已经计算出的傅里叶级数的系数,计算出截短的傅里叶级数:∑-==NNk tjk kea t x 0)(ω这个计算可用L = 2N+1次循环来完成:t N r j k e r a x x 0)1()(22ω--⋅+=其中r 作为循环次数,x2在循环之前应先清零。
完成这一计算的MATLAB 程序为:x2 = 0; L = 2*N+1;for r = 1:L;x2 = x2+ak(r)*exp(j*(r-1-N)*w0*t); end;完成了所有的计算之后,就可以用绘图函数:plot()和stem()将计算结果包括x1, x2, abs(ak)和angle(ak)以图形的形式给出,便于我们观察。
例 2:修改程序Program2_1,使之能够用有限项级数合成例 1所给的周期方波信号,并绘制出原始周期信号、合成的周期信号、信号的幅度谱和相位谱。
% Program2_2% This program is used to compute the Fourier series coefficients ak of a periodic square waveclear,close allT = 2; dt = 0.00001; t = -2:dt:2; x1 = u(t)-u(t-1-dt); x = 0; for m = -1:1x = x + u(t-m*T) - u(t-1-m*T-dt); % Periodically extend x1(t) to form a periodic signalendw0 = 2*pi/T;N = input('Type in the number of the harmonic components N = :'); L = 2*N+1; for k = -N:1:N;ak(N+1+k) = (1/T)*x1*exp(-j*k*w0*t')*dt; endphi = angle(ak); y=0;for q = 1:L; % Synthesiz the periodic signal y(t) from the finite Fourier series 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)'), axis([-N,N,-0.1,0.6])subplot(224)stem(k,phi,'r.'), title('The phase phi(k) of x(t)'), axis([-N,N,-2,2]), xlabel('Index k') 4 用MATLAB 实现CTFT 计算MATLAB 进行傅里叶变换有两种方法,一种利用符号运算的方法计算,另一种是数值计算,本实验要求采用数值计算的方法来进行傅里叶变换的计算。
严格来说,用数值计算的方法计算连续时间信号的傅里叶变换需要有个限定条件,即信号是时限信号(Time limited signal ),也就是当时间|t|大于某个给定时间时其值衰减为零或接近于零,这个条件与前面提到的为什么不能用无限多个谐波分量来合成周期信号的道理是一样的。