论文初稿-连续周期信号时域采样频谱分析程序设计

连续周期信号时域采样频谱分析程序设计摘要信号的频谱分析是信号与系统分析的基础，本文分析了用数值计算以及使用MATLAB软件与MATHEMATICA软件的方法实现连续周期信号的时域采样频谱分析。数值计算即使用傅里叶级数展开方法进行连续周期信号的频谱分析，使用软件计算即设计实用的程序进行周期信号的时域采样频谱分析。关键词频谱分析；傅里叶级数；MATLAB软件；MATHEMATICA软件CONTINUOUSPERIODICSIGNALSPECTRUMANALYSISINTIMEDOMAINSAMPLINGPROGRAMDESIGNABSTRACTREPAIRINGSPECTRALANALYSISOFTHESIGNALISTHEBASISOFSIGNALSANDSYSTEMSANALYSIS,THISPAPERANALYZESTHESPECTRALANALYSISOFTIMEDOMAINSAMPLESTOACHIEVECONTINUOUSPERIODICSIGNALSWITHNUMERICALCALCULATIONSANDTHEUSEOFMATLABSOFTWAREANDMATHEMATICASOFTWAREMETHODNUMERICALCALCULATIONSTHATTHEUSEOFFOURIERSERIESEXPANSIONMETHODFORSPECTRALANALYSISOFCONTINUOUSPERIODICSIGNAL,CALCULATEDUSINGTHESOFTWAREPROGRAMTHATISDESIGNEDANDPRACTICALPERIODICSIGNALINTHETIMEDOMAINSAMPLINGSPECTRUMANALYSISKEYWORDSSPECTRUMANALYSISFOURIERSERIESMATLABSOFTWAREMATHEMATICASOFTWARE目录1引言12周期信号频谱分析的傳里叶级数展开方法121傅里叶级数1211周期信号的分解2212指数形式傅里叶级数322连续时间傅里叶级数的基本性质5221线性5222卷积特性5223微分特性523技巧性方法5231频移方法5232求导方法73周期信号的频谱分析方法731已知变化规律的周期连续信号频谱识别方法732未知变化规律的周期连续信号频谱识别方法833未知周期信号的频谱的理论方法和计算技巧84MATHEMATICA软件计算方法105结束语146致谢15参考文献151引言随着计算机的发展，数字信号处理的理论与技术得到飞速发展，20世纪60年代以来，我国形成了一系列的数字信号处理的理论与算法，比如，数字滤波器，快速傅立叶变换（FFT），这些都是数字信号处理的技术基础，随着信息科技的飞速发展，信号处理取得了重大的飞跃。信号的频谱分析是数字信号处理中的一个很重要的研究课题，对信号进行频谱分析，是对其进行傅里叶变换，得到其振幅谱与相位谱。对于信号来说，分模拟信号与数字信号。进行频谱分析时，对于模拟信号来说，首先对其进行抽样，使其离散化，然后利用离散傅里叶变换（DFT）或者快速傅里叶变换（FFT），然后对其幅度（ABS）和相位（ANGLE）的图像进行分析，而对于数字信号来说，则可直接进行离散傅里叶变换或快速傅里叶变换。在现实生活中对于信号进行频谱分析具有重要的意义。通过对信号频谱的分析，可以得到信号的频率结构，了解信号的频率成分或系统的特征。在此基础之上，可实现对信号的跟踪控制，从而实现对系统状态的早期预测，发现潜在的危险并诊断可能发生故障的原因，对系统参数进行识别及校正。因此，频谱分析是揭示信号特征的重要方法，也是处理信号的重要手段。频谱分析工作在许多科学技术领域都会遇到。首先，在无线电技术的许多方面，例如通讯、地震测量、电信、导航、雷达、电子对抗、空间技术等；此外，由于光象波、机械震动、冲击、响声等各种非电量都可以通过所谓换能器转换成电流或电压的变化来方便地分析，所以，频谱分析在各种振动、噪声、电声、发动机、建筑、生物、医学等领域也起重要作用。许多工程设计都要对振动频率进行认真分析。如果外界强迫振动的频率与某物体的固有频率相近会发生共振。如电动机转动时会引起台座的受迫振动，作用力的频率如果接近机器某部分的固有频率，就可能因共振而引起损坏；轮船和码头会经常收到海浪的周期性冲击，冲击频率如果与船体或码头结构的固有频率相近，也会引起船头或码头的损坏。这就需要研究清楚各海域海浪的频率范围，避免共振的发生；在地震过程中，建筑物的固有频率如果与地面振动的固有频率相近，及时地震烈度较低，也可能激发建筑物产生剧烈震动而导致损坏。所以，建筑工程师对于一栋建筑物的动力学特征，如固有频率，振型，阻尼等是需要认真研究清楚的，特别要避开经常出现起主导作用的地震频率；研制发动机的科学家们为了弄清楚发动机活塞在不同运行速度下的性能，他们在活塞边缘事先嵌入一个微型无线电发射机，这个发射机可以在实际运行条件下把活塞所收到的各种应力和振动的详细信息通过无线电信号发射过来，科学家对这些信号进行频谱分析和其他分析就可以对活塞性能进行正确直观的评价；研制飞机的科学家们需经常做一中所谓疲劳实验。因为在飞机飞行过程中，他的机翼一方面收到发动机的振动，一方面受到不均匀气流的不断冲击，引起机翼金属应力不断变化。虽然这成千上万次的小振动幅度和能量均很小，但是累计下来也是很大的能量，为了避免发生意外，科学家们也要对这些振动进行频谱分析，以保证飞机的正常运行。在工程领域中，MATLAB软件和MATHEMATICA软件是一种倍受程序开发人员青睐的语言，对于一些需要做大量数据运算处理的复杂应用以及某些复杂的频谱分析算法MATLAB显得游刃有余1。2周期信号频谱分析的傳里叶级数展开方法21傅里叶级数周期信号是定义在（）区间，每隔一定时间T，按相同规律重复变化的信号，如图21,所示，它可以表示为图21周期信号（21）MTTF式中M为任意整数。时间T被称为该信号的重复周期，简称周期。周期的倒数称为该信号的频率3。由式（22）1JTCTF可知，周期信号FT在区间可以展开成在完备正交信号空间中的无穷级数。如果完备的,0TT正交函数集是三角函数集或指数函数集，那么，周期信号所展开的无穷级数就分别称为“三角型傅里叶级数”或“指数型傅里叶技术”，统称为傅里叶级数。需要指出，只有当周期信号满足狄里赫利条件时，才能展开成傅里叶级数。通常遇到的周期信号都满足该条件。狄里赫利条件是（1）函数在任意有限区间内连续，或只有有限个第一类间断点（当T从左或右趋于这个间断点时，函数有有限的左极限或右极限）；（2）在一周期内，函数有有限个极大值或极小值。211周期信号的分解设有周期信号，他的周期是T，角频率,它可分解为TFTF2COSCOS2COSCOS2110TBTBTAA23110INNN式23中的系数A，B称为傅里叶系数。它可由式242121TIITIIDFKDFC求得。为简便，式22的积分区间取为或。由式23可得傅里叶系数,0T2,T0N0，1，2，252NCOSTDTNTFAN0，1，2，262NITB式中T为函数的周期，为角频率。由式25和式26可见，傅里叶系数和都TFTNAB是N的函数，其中是N的偶函数，即；而是N的奇函数，即有ANABN将式（23）中同频率项合并，可写成如下形式2COSCOS22110TATATF2710NN式中,N1，2，28ARCTNN2N0BA如将式（28）的形式转化为式（211）的形式，他们系数之间的关系为N1，2，29，SIN，CON0ABA由式（28）可见，是N的偶函数，即有；而是N的奇函数，即有。ANAN傅里叶系数的这些重要性质是很有用的。式（29）表明，任何满足狄里赫利条件的周期函数可分解为直流和许多余弦或正弦分量。其中第一项是常数项，他是周期信号中所包含的直流分量；式中第二项称为基波20ACOS11TA或一次谐波，它的角频率与元周期信号相同，是基波振幅，是基波初相角；式中第三项1称为二次谐波，他的频率是基波频率的二倍，是二次谐波的振幅，是其COS2T22初相角。以此类推，还有三次，四次、谐波。一般而言，称为N次谐波，N是N次谐波的振幅，是其初相角。式（27）表明，周期信号可以分解为各次谐波分量。AN212指数形式傅里叶级数三角函数形式的傅里叶级数含义比较明确，但运算常感不便，因而经常采用指数形式的傅里叶级数。由于（210）2COSJXEX所以式（29）可以写为102NTNJTNJEATF（211）TJNNTJNNEA1102将上式第三项中的N用N代换，并考虑到到N是偶函数，即；是N的奇函数，即A，则上式可写为NTJNNTJNNEAEATFNN11022如将上式中的写成（其中0），则上式可以写为0N0（212）TJNNTF令复数量，称其为复傅里叶系数，简称傅里叶系数，其模为，相角NJNJNFEA21NF为，则傅里叶级数的指数形式为N（213）TJNNEFTF21根据式（29），傅里叶系数（214）21SICO21NNNNJNNJBAJAAEAF将式（25）和式（26）代入上式，得DTTFTDTTFTN221SNJNT2SICO1NO,1，2，,2DTEFTJN（215）这就是求指数形式傅里叶级数的复系数的公式。NF式（213）表明，任意周期信号可分解为许多不同频率的虚指数信号（）之和，其TJNE各分量的复数幅度（或相量）为。22连续时间傅里叶级数的基本性质221线性若，则对于任意常数，有11FTF22FTF1A2（216）1FTFA222时移特性若，则（217）00TJEFTF式中，为常数。0T222卷积特性若，则11FF22FTF（218）211FTF图22时域卷积运算图23频域相乘运算223微分特性设是以为周期的周期信号，其傅里叶系数为TF0T（219）NCTF则导数的傅里叶系数是（220）J0若的傅里叶系数是，则的傅里叶系数是ND2210NDC23技巧性方法231频移方法有些周期信号的周期规律不明显或周期变化计算过于复杂，则可通过频移和增加（减小）直流增量来改变信号位置，由此可使信号由非奇非偶信号变成奇信号或者偶信号，然后进行傅里叶级数展开时就会比较快5。某信号在时域中乘以，相当于频域中频谱向右移；TJE00某信号在时域中乘以，相当于频域中频谱向左移；也就是时域频域TXXJE00T0表21频移特性前面两个结论可以这么直观的理解信号由无数不同旋转角速度（）的旋转向量叠加而成。TX图24信号频移当与相乘之后，原来的所有旋转向量的旋转角速度都将增加，达到，对TXTJE000应的频谱也将由原来的位置向右搬移。0图25信号向右频移当与相乘之后，原来的所有旋转向量的旋转角速度都将减少，达到，TXTJE000对应的频谱也将由原来的位置向左搬移。0图26信号向右频移232求导方法有些周期信号的周期规律不明显或周期变化计算过于复杂，则可对其进行求导（求导只是对其一个周期求导，因为信号的所有信息在一个周期均有完全的反应），这样就是将复杂的信号分析改变成一个周期内的函数进行分析，可以减少计算过程，若经过一次求导以后函数依旧复杂，可再次对其进行求导。给出一个周期的表达式，周期信号表达式为TFTTF，102COS2KKTAATF对其求一次导则有12KKT2,2KKTAB转换可得。2KK对其求二次导则有12COSKKTTATFKKTC,2转换可得。KKA23周期信号的频谱分析方法31已知变化规律的周期连续信号频谱识别方法通过上面分析我们已经知道，对于这类信号要获得其频谱，我们只需将此周期连续信号做傅里叶级数展开，即S设FT为一个周期连续信号，T为其周期，F11/T是信号的基波频率，为信号的基波12F角频率，则FT可以表示为，（N0,1,2）。对FT进行傅里叶级数展开NTTF式222110COSKKTA式223式224式225式226通过上面推导即可计算已知变化规律的周期连续信号的频谱和振，从而画出相应信号FT的幅频特性杆状图。对于周期连续信号FT，若要计算其连续频谱，可对FT取傅里叶变换FTFT，其定义为式227其中密度函数FJW是一个复函数，它可写为式228式中和分别是频谱函数的模和相位。发现是由频率位置各不相JFJFJF同的一系列狄拉克K1带权叠加，对于这样的频谱函数是不能直接画幅频特性图的，可在各频率坐标K1处以狄拉克函数前面的权系数为高度画幅频杆状图来图示周期连续信号所含频率成份。32未知变化规律的周期连续信号频谱识别方法未知变化规律的周期连续信号通常是通过对信号采样，然后分析采样信号的频谱来获得连续信号的频谱，这种做法首先要满足采样定理。为了使实际信号在采样后能够不失真的再现，采样频率必须大于信号最高频率的两倍。其次，为获得连续信号的频谱，要对采样数据进行离散傅里叶变换（DFT）来获得离散信号的频谱。我们假设周期连续信号采样得到的序列为XN，设XN为有限长序列，长度为N，但它取自周期信号，我们可以认为在采样参数选择恰当时，对它进行周期延拓所得周期序列的最小周期与被采样信号的周期是一致的。通过采样定理的学习我们得知采样信号的频谱是被采样信号频谱的周期延SI22COS11101KKKKTBAAT21TKDTFA21SINTKTKB2KKBAA/RCTNJFDTEFTFTJXRJE拓，所以我们只需要对一个采样周期进行研究便可得到信号的全部信息。现在我们取定采样序列XN，长度为N（称为一个主值区间）对它进行频谱分析即可。对XN取DFT，计算过程如下具体计算步骤如下式229NTAAXNTXKNNNJJEXEXDFK210式230102NNKJXT式中是待采样得周期连续信号，T是采样周期，是采样后的有限长TXATXA离散序列，与是等同的。由上式可知有限长离散信号频谱的采样是由其连续频谱经过采样得N到的。由此得DTFT是计算离散序列的频谱函数的，而DFT则是计算离散序列频谱的采样的。33未知周期信号的频谱的理论方法和计算技巧例20COS210COS31TTTF3684540试画出的振幅谱和相位谱解的周期信号，题中所给的表达式可视为的傅里叶级数展开式。据TFTFTF10COS2NNA可知，其基波频率（RAD/S），其基本周期T2S，2，3，6分别为二、三、六次谐波频率。且有120A0312240A538606其余0NKNJADFTAAEXKNXTTX2JJAEXKNJEX2图31信号的振幅谱图32信号的相位谱图33信号的双边频谱振幅谱图34信号的双边频谱相位谱4MATHEMATICA软件计算方法例1利用MATHEMATICA软件快速对周期信号进行频谱选择2XT_5COS100T/42COS300T/805COS500TTS0001FS1/TSN1000TSRANGE0,N1TSXNXTSLISTPLOTTRANSPOSETS,XN2042,JOINEDTRUEN20TNTSYNXN120Y2FOURIERYN,FOURIERPARAMETERS1,1Y1/2AYABSY1N/2AMMAXAYTH008FRANGE0,FLOORN/21FS/NDATA1TRANSPOSEF,AYDATA2TRANSPOSEAY,F,ARGY1N/2LISTPLOTDATA1,FILLINGAXISSPSELECTDATA2,1THAMUT_UNITSTEPTT0001Q5TS000001LROUNDT/TSST_4Q/TTUTUTT/8T/2T/16UTT/8UT3T/8T/4UT3T/8UT5T/89T/16T/2UT5T/8UT7T/8TTUT7T/8UTTTSRANGE0,L1TSSNSTSDATATRANSPOSETS,SNP1LISTPLOTDATA,ASPECTRATIO1/3,AXESLABEL“T“,“LADDERT“S2FOURIERSN,FOURIERPARAMETERS1,1S1/2AABSS1L/2/CHOPPRINT“F0,F1,F2,F3“,0,1,2,3/TPRINT“A0,A1,A2,A3“,A14PRINT“0,1,2,3“,14ARGS1L/2LISTPLOTTABLEK/T,AK1,K,0,L/21,ASPECTRATIO1/3,AXESLABEL“F/HZ“,“AF“,FILLINGAXIS,PLOTRANGEALLLISTPLOTTABLEK/T,K1,K,0,L/21,ASPECTRATIO1/3,AXESLABEL“F/HZ“,“F“,FILLINGAXISSCT_SUMAK1COS2K/TTK1,K,0,40P2PLOTSCT,T,3T,5T,AXESLABEL“T“,“SCT“,ASPECTRATIO1/3,PLOTRANGEALL,PLOTSTYLEREDSHOWP2,P1运行程序可得出以下结论信号采样图图44信号采样图频点、振幅、相位分别为振幅频谱图图45振幅频谱图相频图图46相频图信号波形图与采样图叠加图47信号波形图和采样图叠加5结束语本次毕业设计至此已经接近尾声，在这几个月的时间里，我通过利用MATLAB和MATHEMATICA强大的数据运算功能以及图像处理功能对于连续时间的复频域分析进行深入的研究。在整个设计过程中，我首先对于所学的基础信号知识进行温习巩固，比如傅立叶级数、时域采样、信号频谱分析等；其次，整个实现过程是通过MATLAB软件与MATHEMATICA软件完成的，MATLAB的图形功能十分强大，具有良好的设计平台，在此次设计过程中，我熟练了MATLAB的编程方法，掌握了很多函数的表示含义及使用方法；MATHEMATICA软件数学实现方法特别的强，并且内容丰富，操作简便。最后，通过此次毕业设计，我对设计所用到的软件有了更加深刻的认识，MATLAB与MATHEMATICA不仅在数值计算方面的功能十分强大，而且其图形仿真功能能够满足各个领域的需要，因此MATLAB与MATHEMATICA已经成为我们工作学习中不可或缺的软件。由于MATLAB与MATHEMATICA软件是专业性较强的软件，所以刚开始使用是比较困难的，通过老师的指导、翻阅了大量的相关资料，不能说已经通悟了MATLAB与MATHEMATICA应用，但是对一些基本知识还是有所了解的，在学习的过程中学习的过程中我们进一步对MATLAB与MATHEMATICA编程中的常用语句、过程已经初步掌握。这次毕业设计，使我懂得了理论与实际相结合是很重要的，只有理论知识是远远不够的，只有把所学的理论知识与实践相结合起来，从理论中得出结论，才能真正为社会服务，从而提高自己的实际动手能力和独立思考的能力。这次课程设计终于顺利完成，在设计中遇到的运行和调试问题，最后在老师的耐心指导下，终于游逆而解。在以后的学习过程中我要不段的学习，不断丰富自己的知识。6致谢在