实验三-周期信号的频谱分析.doc
《实验三-周期信号的频谱分析.doc》由会员分享,可在线阅读,更多相关《实验三-周期信号的频谱分析.doc(10页珍藏版)》请在咨信网上搜索。
1、实验三_周期信号的频谱分析 作者: 日期:10 个人收集整理 勿做商业用途实验三 信号的频谱分析一、实验目的1、掌握连续时间周期信号的傅里叶级数的物理意义和分析方法;2、观察截短傅里叶级数而产生的“Gibbs现象”,了解其特点以及产生的原因;3、掌握各种典型的连续时间非周期信号的频谱特征二、原理说明:1、连续时间周期信号的傅里叶级数分析任何一个周期为T1的正弦周期信号,只要满足狄利克利条件,就可以展开成傅里叶级数。其中三角傅里叶级数为: 2.1或: 2.2其中,称为信号的基本频率(Fundamental frequency),分别是信号的直流分量、余弦分量幅度和正弦分量幅度,为合并同频率项之后
2、各正弦谐波分量的幅度和初相位,它们都是频率的函数,绘制出它们与之间的图像,称为信号的频谱图(简称“频谱),图像为幅度谱,图像为相位谱.三角形式傅里叶级数表明,如果一个周期信号x(t),满足狄里克利条件,那么,它就可以被看作是由很多不同频率的互为谐波关系(harmonically related)的正弦信号所组成,其中每一个不同频率的正弦信号称为正弦谐波分量 (Sinusoid component),其幅度(amplitude)为。也可以反过来理解三角傅里叶级数:用无限多个正弦谐波分量可以合成一个任意的非正弦周期信号.指数形式的傅里叶级数为: 2.3其中,为指数形式的傅里叶级数的系数,按如下公式
3、计算: 2。4指数形式的傅里叶级数告诉我们,如果一个周期信号x(t),满足狄里克利条件,那么,它就可以被看作是由很多不同频率的互为谐波关系(harmonically related)的周期复指数信号所组成,其中每一个不同频率的周期复指数信号称为基本频率分量,其复幅度(complex amplitude)为。这里“复幅度(complex amplitude)”指的是通常是复数。上面的傅里叶级数的合成式说明,我们可以用无穷多个不同频率的周期复指数信号来合成任意一个周期信号。然而,用计算机(或任何其它设备)合成一个周期信号,显然不可能做到用无限多个谐波来合成,只能取这些有限个谐波分量来近似合成。假设
4、谐波项数为N,则上面的和成式为: 2。5显然,N越大,所选项数越多,有限项级数合成的结果越逼近原信号x(t)。本实验可以比较直观地了解傅里叶级数的物理意义,并观察到级数中各频率分量对波形的影响包括“Gibbs”现象:即信号在不连续点附近存在一个幅度大约为9%的过冲,且所选谐波次数越多,过冲点越向不连续点靠近.这一现象在观察周期矩形波信号和周期锯齿波信号时可以看得很清楚。2、傅里叶级数的MATLAB计算设周期信号x(t)的基本周期为T1,且满足狄里克利条件,则其傅里叶级数的系数可由式2.4计算得到。式2.4重写如下:基本频率为: 对周期信号进行分析时,我们往往只需对其在一个周期内进行分析即可,通
5、常选择主周期(Principle period)。假定x1(t)是x(t)中的主周期,则计算机不能计算无穷多个系数,所以我们假设需要计算的谐波次数为N,则总的系数个数为2N+1个.在确定了时间范围和时间变化的步长即T1和dt之后,对某一个系数,上述系数的积分公式可以近似为: 对于全部需要的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)=
6、);N = input(Type in the number N=); k = -N:N; L = 2N+1; ak = x1exp(-jk*w0t)dt/T;需要强调的是,时间变量的变化步长dt的大小对傅里叶级数系数的计算精度的影响非常大,dt越小,精度越高,但是,计算机计算所花的时间越长。例题1:给定一个周期为T1 = 2s的连续时间周期方波信号,如图所示,其一个周期内的数学表达式为:112-1-2x(t)t0图1 周期方波信号 解:首先,我们根据前面所给出的公式,计算该信号的傅里叶级数的系数。 因为:w0 = 2/T1 = ,代入上式得到:在MATLAB命令窗口,依次键入: k = -1
7、0:10; ak = (-j).k). (sin(k+eps)pi/2)。/(k+eps)*pi) % The expression of akeps 叫做机器的浮点运算误差限!PC机上eps的默认值为2.2204*10(-16)。若某个量的绝对值小于eps,就认为这个量为0ak = Columns 1 through 4 0。0000 0 + 0。0354i 0.0000 0 + 0。0455i Columns 5 through 8 -0。0000 0 + 0。0637i -0。0000 0 + 0.1061i Columns 9 through 12 -0.0000 0 + 0.3183
8、i 0。5000 0 - 0。3183i Columns 13 through 16 -0.0000 0 - 0.1061i 0.0000 0 - 0。0637i Columns 17 through 20 -0。0000 0 0。0455i 0。0000 0 - 0.0354i Column 21 -0.0000 从MATLAB命令窗口,我们得到了该周期信号从到共21个系数.紧接着再键入以下命令: subplot(221) stem(k,abs(ak),k.) title(The Fourier series coefficients) xlabel(Frequency index k)就得
9、到一幅如右图所示的描述与k之间的关系的图形。以上是我们通过手工计算得到的这个周期信号的傅里叶级数表达式及其频谱图,下面给出完成傅里叶级数系数计算的相应MATLAB范例程序。 Program3_1% This program is used to evaluate the Fourier series coefficients ak of a periodic square waveclear, close allT = 2; dt = 0。0001; t = -2:dt:2;x1 = ut(t) ut(t-1-dt); x = 0;for m = -1:1 % Periodically ext
10、end x1(t) to form a periodic signal x = x + ut(tm*T) - ut(t1m*T-dt);endw0 = 2*pi/T;N = 10; 谐波次数L = 2*N+1;for k = N: N; Evaluate the Fourier series coefficients ak ak(N+1+k) = (1/T)*x1*exp(j*kw0t)*dt;endphi = angle(ak); % ak相位执行程序Program3_1后,就完成了信号的傅里叶级数的系数的计算,在命令窗口键入 ak命令窗口就可以显示傅里叶级数的21个系数:ak = Colu
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 实验 周期 信号 频谱 分析
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【天****】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【天****】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。