DSP谐波分析的理论基础ccs仿真结果及分析.docx

基于DSP的方波信号谐波分析一切周期非正弦信号都是由某一频率的正弦波（基波）和若干数倍于此频率的正弦波（高次谐波）叠加而成。谐波分析的工作就是了解一个信号的基波和谐波的组成及其占比。该设计使用基于ICETEK—VC5509—A的DSP实验箱，软件环境为CCS3.3。一·理论部分1信号采集首先是对信号的采集。主要是对输入信号进行AD转换，采样率越高，则分析精度越高，但硬件消耗同时加大。模拟信号经过 A/D 变换转换为数字信号的过程称之为采样，信号采样后其频谱产生了周期延拓，每隔一个采样频率ωs，重复出现一次。（1）.频混现象频混现象又称为频谱混叠效应，它是由于采样信号频谱发生变化，而出现高、低频成分发生混淆的一种现象，如图3.1.1所示。信号x(t)的傅里叶变换为X(ω)，其频带范围为-ωm~+ωm；采样信号x(t)的傅里叶变换是一个周期谱图，其周期为ωs，并且：ωs=2π／TsTs为时域采样周期．当采样周期Ts较小时，ωs＞2ωm，周期谱图相互分离如图1中（b）所示；当Ts较大时，ωs＜2ωm，周期谱图相互重叠，即谱图之间高频与低频部分发生重叠，如图3.1中(c)所示，此即为频混现象，这将使信号复原时丢失原始信号中的高频信息。图1采样信号的频混现象下面从时域信号波形来看这种情况。图2(a)是频率正确的情况，以及其复原信号；(b)是采样频率过低的情况，复原的是一个虚假的低频信号。图2发生频混现象的时域信号波形当采样信号的频率低于被采样信号的最高频率时，采样所得的信号中混入了虚假的低频分量，这种现象叫做频率混叠。（2）. 采样定理上述情况表明，如果ωs＞2ωm，就不发生频混现象，因此对采样脉冲序列的间隔Ts须加以限制，即采样频率ωs（2π／Ts）或 fs(1／Ts)必须大于或等于信号x(t)中的最高频率ωm的两倍，即ωs＞2ωm，或 fs＞2fm。为了保证采样后的信号能真实地保留原始模拟信号的信息，采样信号的频率必须至少为原信号中最高频率成分的2倍。这是采样的基本法则，称为采样定理。需要注意的是，在对信号进行采样时，满足了采样定理，只能保证不发生频率混叠，对信号的频谱作逆傅立叶变换时，可以完全变换为原时域采样信号，而不能保证此时的采样信号能真实地反映原信号。工程实际中采样频率通常大于信号中最高频率成分的3到5倍。2 FFT变换与谐波分析（1）FFT变换采样到信号后进行调整，主要是化为以弧度为单位的数值（比例关系）。之后进行FFT变换即可得到信号频谱图。FFT变换即快速傅立叶变换的算法是蝶型算法，即交叉排序后进行乘加运算。傅立叶分析是联系时域和频域的桥梁,是频域数字化分析的理论基础。一个波形的傅立叶变换的实质是:把这个波形分解成许多不同频率的正弦波之和。如果这些正弦波加起来成为原来的波形,那么,就确定了这个波形的傅立叶变换。傅立叶变换可以辨别或区分开组成任意波形的一些不同频率的正弦波和它们各自的振幅,在数学上,这种关系被表示为:式中,χ( t)是被分解为正弦函数之后的波形; F (ω)是χ( t)的傅立叶变换。傅立叶变换可以看作是时间函数在频率域上的表示,由傅立叶变换给出的频率域包含的信息和原函数时间域内所包含的完全相同,不同的仅是信息的表示方法。自从1965年FFT算法问世以来,已深刻地改变了作数字谱分析的途径。它代替了过去计算功率谱密度时先计算相应相关函数的离散估计,然后求出相关函数的DFT的谱而得到功率谱的计算途径。如果给定的周期函数满足狄里赫利条件,则该周期函数定可展开为一个收敛的正弦函数级数。这样:其中,两式中的各个系数的计算公式为式中, T为周期;ω为角频率,；A0 为直流分量; A1 cos (ωt +ψ1 ) 称为一次谐波; Ak cos ( kωt +ψk )成为k次谐波。因此,可以看出,一个周期振动可看成是基频是的谐振动和其他高次谐振动的叠加。频谱是由等距离的频谱线组成的,谐波频率具有简单的倍数关系。它们的振幅却随值k变大而减小。整个频率范围就是信号的频带,频谱是离散频谱。通过上述傅里叶变换,方波可表示为:送入r（n）、MDIT-FFT程序框图根据DIT-FFT原理和过程，DIT-FFT的完整程序框图包括以下几部分：(1)倒序:输入自然顺序序列x(n)，根据倒序规律，进行倒序处理；(2)循环层1:确定运算的级数，L=1M (N=2M)；确定一蝶形两输入数据距离B=2L-1(3)循环层2:确定L级的(B=)2L-1个旋转因子；旋转因子指数p=2M-LJ，J=0B-1；(4)循环层3:对于同一旋转因子，用于同一级2M-L个蝶形运算中：k的取值从J到N-1，步长为2L (使用同一旋转因子的蝶形相距的距离) (5)完成一个蝶形运算。3 频谱分析信号的频谱分析是信号与系统分析的基础。用数值计算的方法实现确知信号的频谱分析，