周期三角波激励下的动态电路响应解法研究

周期三角波激励下的动态电路响应解法研究摘要由于直接对周期三角波激励下的动态电路求解有一定的难度，而通过分析，正弦信号经过任何复杂的系统，系统输出和输入波形一致，只是幅度被放大，被延时，波形不变，周期三角波是一系列不同频率，不同相位，不同振幅正弦信号的叠加，所以把三角波通过傅里叶级数展开，化为正弦信号求解方法求解。但是传统的手工计算较为困难，因此我们设计一个程序，用计算机快速的解决这类问题。关键词周期三角波动态电路傅里叶级数展开CYCLEOFTRIANGULARWAVEEXCITATION,DYNAMICCIRCUITRESPONSEMETHODRESEARCHABSTRACTASADIRECTDYNAMICCYCLEOFTHETRIANGULARWAVEEXCITATIONCIRCUITFORSOLVINGACERTAINDEGREEOFDIFFICULTY,ANDTHROUGHANALYSIS,SINUSOIDALSIGNALTHROUGHANYCOMPLEXSYSTEM,THESYSTEMOUTPUTANDINPUTTHESAME,BUTTHEMAGNITUDEISAMPLIFIEDBYTHEDELAY,WAVEFORMCONSTANT,PERIODICTRIANGULARWAVEISARANGEOFDIFFERENTFREQUENCIES,DIFFERENTPHASES,DIFFERENTAMPLITUDESINUSOIDALSIGNALSUPERIMPOSED,SOTHETRIANGULARWAVEBYFOURIERSERIESEXPANSIONINTOSOLVINGMETHODFORSOLVINGASINUSOIDALSIGNALHOWEVER,THETRADITIONALMANUALCALCULATIONSMOREDIFFICULT,SOWEDESIGNEDAPROGRAMTOSOLVETHISPROBLEMQUICKLYBYCOMPUTERKEYWORDSPERIODICTRIANGULARWAVEDYNAMICCIRCUITFOURIERSERIESEXPANSION目录1引言12周期信号的频谱分析方法121傅里叶级数展开法122周期信号离散化频谱分析法13动态电路频域系统函数计算方法231系统函数HS2311由微分方程求解系统函数2312在动态电路系统中求解系统函数232系统的频域系统函数HJ计算方法44正弦信号激励下的动态电路零状态响应计算方法441单频率正弦信号激励下系统的零状态响应求解442多频率正弦信号激励下系统的零状态响应求解55周期三角波激励下的动态电路响应551周期三角脉冲信号的傅立叶级数求解552周期三角脉激励下动态电路响应频域系统求解66周期三角波激励下的动态电路响应解法程序设计761程序设计思路762程序清单77程序应用实例8结束语10参考文献101引言众所周知，我们若想直接对周期三角波激励下的动态电路求解，将会有一定的难度，而通过我们所学知识分析了解到单一频率的正弦信号经过任何复杂的系统，系统输出和输入的波形一致，只是幅度被放大，被延时，波形不变，但多频率信号因电路对不同频率成分放大延时的量不同，因而输出信号的波形不同于输入信号。周期三角波是一系列不同频率，不同相位，不同振幅正弦信号的叠加，所以把三角波通过傅里叶级数展开，化为正弦信号叠加求解。我们以三角波激励下的动态响应为研究对象，探讨周期三角波信号这种特殊的多频率正弦信号激励下的动态电路响应求解方法及输出信号的特点。研究发现我们讨论的时域解法手工计算很繁琐，如果用MATHEMATICA编程求解则效率显著提高。2周期信号的频谱分析方法21傅里叶级数展开法任何一个满足狄里赫利DIRICHLET条件的非正弦周期信号函数都可以分解为一个恒定分量与无穷多个频率为非正弦周期信号频率的整数倍、不同幅值的正弦分量的和1。2111110SINCO2NTBTATF傅里叶系数212TDTF0213NTNA1COS2214TTFB0I频率相同的余弦项与正弦项合并为一个正弦函数215COSSINCO111NNNTATTA傅里叶级数又可写成下列形式2216110SNNTTF21720AA218NNB219IRG式中为直流（常数）分量，为基波FUNDAMENTALWAVE或一次谐波FIRST0ASI11THARMONICWAVE，为N次谐波NTHHARMONICWAVE。SIN1TN22周期信号离散化频谱分析法对于复杂周期信号的傅里叶级数展开中，、的积分运算过于复杂，手工难以算出其FANB结果，所以我们需要运用周期信号离散化频谱分析方法，对周期信号进行合理采样，利用采样数据计算周期频谱。对于已知周期信号，其周期，最高幅度，我们对其离散化分析步骤为3TF01T5E1合理选取采样间隔对于周期函数，采样区间我们通常选取周期函数的一个周期来进行研S究。理论上，周期信号有无限多个频率成分，但频率越高，振幅越小，所以技术层面上，可以忽略高频成分，而视为周期信号含有限个频率成分，若我们需要保留个频率成分，则采样间隔的表M达式为221TS42写出给在一个周期的函数关系式。TF3时间离散化采样点数，。STN1,0NKTSS4对周期信号作数学采样。STF5对信号的采样数据作快速傅里叶变换，实数形式傅里叶级数展开式的频谱称为实频谱。复数形式傅里叶级数展开式的频谱称为复频谱。实频谱中角频率的取值范围从0到是单边频谱，而复频谱中的取值范围从一到,是双边离散频谱。实频谱的幅值是相应复频谱的两倍，初相位与复频谱的初相位互为余角。6利用采样数据快速傅里叶变换算出周期信号的前个频率成分的振幅、频率、初相位。若M以或为纵坐标，为横坐标就可以分别画出幅值或相位随频率变化的关系曲线。称FAFF为幅频或相频特性，二者合称为周期信号的频谱。TF3动态电路频域系统函数计算方法31系统函数HS系统的响应一方面与激励有关，也与系统自身结构有关4。系统函数就是试图描述系统的本身特性，而与输入无关。若为激励，为零状态响应，且令系统的单位冲激响应为，则TFTYTH有311TFHT对上式等号两端同时求拉氏变换，并设，YLTSYTFLTSFTLTSH则有TFLTSF312H故有313S称为复频域系统函数，简称为系统函数。可见系统函数就是系统零状态响应的象SHSTY函数与激励的象函数之比，也是系统单位冲击响应的拉式变换。YTFSFTH311由微分方程求解系统函数（314）0111TYBDTDTYBTYDATTATNMNMNN系统为零状态，对方程两边取单边拉氏变换，可得5（315）011SFBSBYASASNMN所以微分方程描述的LTI系统的系统函数为（316）011ASSASFYHNN312在动态电路系统中求解系统函数在动态电路S域中，我们很难通过电路图写出系统的微分方程，因而我们只能根据电路的特性，由电路KCL、KVL关系，通过运算电路法，映射在S域中，写出系统的系统函数。研究电路问题的基本依据是描述互连各支路（或元件）电流、电压相互关系的基尔霍夫定律（KCL和KVL）和电路元件端电压与流经该元件电流的电压电流关系VCR6。现讨论它们在S域的形式。KCL方程描述了在任意时刻流入（或流出）任一结点（或割集）各电流关系的方程，0TI它是各电流的一次函数线性函数），若各电流的象函数为称其为象电流，则由线性TIJIJ性质有3170SI上式表明，对任一结点（或割集）,流入（或流出）该结点的象电流的代数和恒为零。在此仍称式317为KCL。同理,KVL方程也是回路中各支路电压的一次函数，若各支路电压的象函数为0TUTUJ称其为象电压,则由线性性质有SUJ3180SU上式表明，对任一回路，各支路象电压的代数和恒等于零。在此仍称式318为KVL。对于线性时不变二端元件R、L、C，若规定其端电压与电流为关联参考方向，其相应TUTI的象函数分别为和,那么由拉普拉斯变换的线性性质及微分性质、积分性质可得到他们SI的S域模型。零状态条件时元件的S域模型7R,L,C元件的时域电压电流关系为319DTUCTITILVRTCC零状态条件时元件的S域网络模型31101SICSVLIR由上式可见（1）R、LS,1/CS是复频域阻抗，S域的电压电流关系满足复频域（广义）的欧姆定律。（2）S域的模型将元器件时域的微、积分运算关系变为代数运算关系。举例说明用S域等效模型，可以得到网络的系统函数。例311图311（A）所示电路系统，输入为，输出为，试求系统函数HS1TV2TV解1先画出S域等效电路图，如图311（B）因为求系统函数，所以系统响应为零状态响应，因而元件的S域模型为无初始条件时的S域模型2根据S域模型，直接表达出输入和输出，由系统函数的概念得到系统函数621/42/1SSVH3111ARC时域电路图BRCS域电路图图311电路图由以上讨论可见，经过拉普拉斯变换，可以将时域中用微分、积分形式描述的元件端电压UT与电流IT的关系，变换为S域中用代数方程描述US与IS的关系，而且在S域中KCL、KVL也成立。这样，在分析电路的各种问题时，将原电路中已知电压源、电流源都变换为相应象函数；未知电压、电流也用象函数表示；各电路元件都用其S域模型替代（初始状态变换为相应的内部象电源），则可画出原电路的S域电路模型。对该S域电路而言，用以分析计算正弦稳态电路的各种方法（如无源支路的串、并联。电压源与电流源的等效变换，等效电源定理以及回路法，终点法等等）都适用。这样。可按S域的电路模型解出所需位置响应的象函数，取其逆变换就得到所需的时域响应。需注意的是，在做电路的S域模型时，应画出其所有的内部象电源，并特别注意其参考方向。32系统的频域系统函数HJ计算方法频率响应是指系统在正弦信号激励下稳态响应随信号频率的变化情况，系统未定时，令,则JS得稳定系统的频率响应特性函数321根据具体输入周期信号的傅里叶级数展开式，我们依次可以得到频域系统函数所对应的角频率，模值，主幅角。我们代入列321中，得3224正弦信号激励下的动态电路零状态响应计算方法41单频率正弦信号激励下系统的零状态响应求解设正弦信号激励下的系统为稳定系统且为零状态，今设系统的激励为单频率正弦信号，TF且在时刻作用于系统8，即0T411COS1NNTATF且令该系统的频域系统函数为。1JH则可得频域系统函数模值为4121JNH频域系统函数主幅角为4131JARGN23612JJHJSS于是得系统零状态响应函数为COS1NNNHTHATY414S11JNARGTJ可见系统的正弦零状态响应仍为与激励同频率的正弦函数，但其振幅却增大了TNF倍，相位增加了。1JNH1JRG42多频率正弦信号激励下系统的零状态响应求解由上述我们知道单频率正弦信号激励下系统的零状态我们很容易可以求解，而实际中，我们需要解决诸多由一系列单频率正弦信号叠加起来的多频率正弦信号激励系统的零状态响应求解9，我们今令一个系统的激励为多频率正弦信号，且令其表达式为TS42111COSNNTATS其该系统的频域系统函数为1JNH模值及主幅角为4221JNHA423RGN由式414，我们可以得出其系统零状态响应函数为11COSNNNTY4245周期三角波激励下的动态电路响应51周期三角脉冲信号的傅立叶级数求解FTTT/2T图51三角脉冲波波形图如图51所示，有一周期三角波脉冲信号，我们先对其进行傅里叶级数展开。TTTTETF2,1JHARGTJHA由式5111110SINCOS2NTBTATF512TDTFA0513NT1COS514TNTFB0I2我们对作偶拓展，知可展开成周期为的余弦级数，且由式512和513可得TFTF515NTDTNTFA01COS2COS4421201TDTNETTE5162为奇数N517EDTFTA004518422为奇数NABANN519为奇数）（所以的傅里叶级数展开为TF为奇数5110COS14212TTNETFNN52周期三角脉激励下动态电路响应频域系统求解根据对三角波激励信号的傅里叶级数展开，且我们令一个动态电路的频域系统函数为521所以频域系统函数的模值为522211364NJHAN频域系统函数的主幅角为5236ARCT1JGN所以三角脉冲激励下的响应函数为236JJ5246ARCTNCOS36414COS112121NEHTHATYNNNN6周期三角波激励下的动态电路响应解法程序设计61程序设计思路程序设计是用精确的语言告诉计算机告诉计算机该做什么，首先我们因明确我们要解决的问题是设计一个周期三角波激励下的动态电路响应求解，以周期三角波激励下的动态电路响应为研究对象，探讨周期三角波信号激励下的动态电路响应求解方法及输出信号的特点。且我们使用的编程软件是MATHEMATICA，根据我们所学知识，以及对编程语言的应用，在MATHEMATIA中编程，最后验证我们所编写的程序有无错误，能否正确执行并是否可以解决我们所要求得的问题。具体程序我们将在下面具体展示。62程序清单程序如下所示10CLEAR“F“,“T“,U,S,T,TS,LUT_UNITSTEPTT0001Q5TS000001LROUNDT/TSST_2Q/TTUTUTT/2TTUTT/2UTTTSRANGE0,L1TSSNSTSDATATRANSPOSETS,SNP1LISTPLOTDATA,ASPECTRATIO1/3,AXESLABEL“T“,“LADDERT“S2FOURIERSN,FOURIERPARAMETERS1,1S1/2AABSS1L/2/CHOPFRANGE0,L/21/TARGS1L/2FE0,1,3,5,7,9/TAEA1,2,4,6,8,10E1,2,4,6,8,10PRINT“F0,F1,F3,F5,F7,F9“,FEPRINT“A0,A1,A3,A5,A7,A9“,AEPRINT“0,1,3,5,7,9“,ELISTPLOTTABLEK/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_SUMAEKCOS2FEKTEK,K,1,6P2PLOTSCT,T,3T,5T,AXESLABEL“T“,“SCT“,ASPECTRATIO1/3,PLOTRANGEALL,PLOTSTYLEREDSHOWP2,P1YT02YT5YT2FT5FTHS_25S/S202S5HFF_HI2FHARGHFFEAHABSHFFEYT_SUMAEKAHKCOS2FEKTEKHK,K,1,6PLOTYT,T,3T,5T,AXESLABEL“T“,“YT“,ASPECTRATIO1/3,PLOTRANGEALL,PLOTSTYLERED7实例例证RLC三阶电路具有电阻电感电容的无源二端网络如图71所示，其中R1，L1H,C2F。现已电压为输入，电压为输出，分析RLC三阶系统。1UT2UT图71时域电路图我们先利用拉普拉斯变换进行分析，将时域映射到复频域（S域）电阻不变，电感，电容，时域电路变成S域电路，LS1CS即如图72所示图72S域电路图由复频域电路图建立复频域代数方程,代入R,L,C的值。得出输出象函数7112132USUSS由输出象函数得系统函数7121/23SSH根据我们所求的系统函数，将函数中的各项参数带入程序，可得振幅频谱图图73幅频图相频图图74相频图信号波形图与采样图叠加000200020004T12345SCT图75叠加图输出信号波形图000200020004T125125125125YT图76输出波形图结语这次毕业设计让我受益匪浅。通过这次设计我对自己在大学四年时间里所学的知识得到了全面的回顾，并充分发挥对所学知识的理解和对毕业设计的思考及书面表达能力，最终完成了这份论文。通过这次毕业设计我发现，只有理论水平提高了，才能够将课本知识与实践相整合，理论知识服务于教学实践，以增强自己的动手能力。这个设计十分有意义，我获得很深刻的经验。通过这次毕业设计，我们知道了理论和实际的距离，也知道了理论和实际想结合的重要性，也从中得知了很多书本上无法得知的知识。我们的学习不但要立足于书本，以解决理论和实际教学中的实际问题为目的，还要以实践相结合，理论问题即实践课题，解决问题即课程研究，学生自己就是一个专家，通过自己的手来解决问题比用脑子解决问题更加深刻。学习就应该采取理论与实践结合的方式，理论的问题，也就是实践性的课题。这种做法既有助于完成理论知识的巩固，又有助于带动实践，解决实际问题，加强我们的动手能力