傅里叶变换、数字滤波器设计、标准表插值算法

1、傅里叶变换周期函数可表示为：其中：周期函数的周期为。频率，角频率，为正整数。周期函数的直流分量。为各次谐波的频率。周期函数可化为：(三角函数公式：)其中： 即周期函数可表示为不同频率成分的正弦函数的和。其中频率为基波的频率。 根据欧拉公式，有：所以周期函数可表示为： = 而 = = = = 令 为正整数 则当 取整数时，可以合写为一个式子 （n = 0, 1，2，.）所以有 为整数非周期函数，当时，有所以取，当时，。从而亦即令则 = =因此有 (1) (2)称式(1)中函数为函数的傅里叶变换，式(2)中函数为函数的傅里叶逆变换。函数即为函数的频谱。图1 是函数y1和y2的函数图。其中y1=si

2、n(t)。y2=sin(t)+0.5*cos(3*t)+0.2*sin(8*t)+0.35*cos(15*t)。y1是标准的正弦函数，y2中加入了高次谐波分量。图1 谐波分量图图2 是偶次谐波的函数图。图2 偶次谐波图图3 是偶次谐波的频谱图。图3 偶次谐波频谱图图4 是偶次谐波5次谐波含量和20次谐波含量的波形图。图4 偶次谐波5次谐波含量和20次谐波含量的波形图傅里叶分析在电路上的应用函数的傅里叶变换记为 ，函数的傅里叶变换记为，即，。 则有傅里叶变换的线性性质 = 傅里叶变换的微分性质 傅里叶变换的积分性质 电路上的一个例子。有一段rlc电路如图5所示图5 rlc电路求电路的电流 ，列方

3、程有函数的傅里叶变换为，函数的傅里叶变换为，对方程两边做傅里叶变换，有求得求的傅里叶逆变换得代入具体的参数值，即可求得电路的电流。函数的卷积已知函数，则积分称为函数和的卷积，记为按傅里叶变换的定义，有 = = = = = =即两个函数卷积的傅里叶变换等于这两个函数傅里叶变换的乘积。数字低通滤波器的设计模拟二阶低通滤波器的电路如图6所示。图6 模拟二阶低通滤波器电路用傅里叶变换分析电路，可以证明其中，。设则有函数为图6模拟二阶低通滤波器的传递函数。为放大系数，为滤波器的截止角频率，为滤波器的品质因数。取，则，。函数的频谱图如图7所示。图7 函数的频谱图()特别的，取，则, ，函数的频谱图如图8所

4、 示。图8 函数的频谱图()取, 的参数，当， 求得。即为函数的半功率点。函数的零极点图如图9所示。时极点位于左半平面。图9 函数的零极点图()特别的，当，时，。当时，函数的零极点位于右半平面。取，函数的零极点图如图10所示。函数极点位于右半平面。图10 函数的零极点图()函数的相频图如图11所示。图11函数的相频图将函数级联，构成多阶低通滤波器，如图12所示的2阶、4阶、6阶低通滤波器的频谱图。图12 2阶、4阶、6阶低通滤波器的频谱图()由，根据卷积定理得。在频域上对函数采样，并对函数做傅里叶逆变换得 。二阶模拟低通滤波器在时域上的传递函数的图形如图13所示。对函数和函数做卷积运算，求得函

5、数，即通过数字滤波器滤波后的结果。函数和函数的图形如图15所示。图14是函数的基波经过滤波器后产生相位延时的例子。图16是模拟二阶低通滤波器电路运行后的结果。函数和图6中模拟电路给出的结果是一致的。图13时域上的传递函数图14 滤波器的相位延时图15 函数和函数图16 模拟二阶低通滤波器电路运行后的结果标准表的计算公式电压表达式： 电流表达式： 电压有效值 电流有效值瞬时有功功率平均有功功率瞬时无功功率 (注)平均无功功率视在功率功率因数 (注)有功电能无功电能视在电能 脉冲频率离散傅里叶变换(dft)离散傅里叶变换的逆变换(idft)奈奎斯特采样定理标准表的插值算法定频采样的同步问题需要插值

6、算法。对采样到的波形分段插值。一副正弦曲线图用线性分段插值后的图形如下。将线性分段插值的图像局部放大，如下图所示。使用三次样条插值算法，得到的图形如下将三次样条插值的图像局部放大，如下图所示。三次样条插值对于 n+1 个给定点的数据集 xi ，我们可以用 n 段三次多项式在数据点之间构建一个三次样条。如果表示对函数 f 进行插值的样条函数，则样条函数满足以下条件。插值特性：s(xi) = f(xi)样条相互连接：si-1(xi) = si(xi), i=1,.,n-1两次连续可导：si-1(xi) = si(xi) 以及 si-1(xi) = si(xi), i=1,.,n-1。由于每个三次多项式需要四个条件才能确定曲线形状，所以对于组成 s的 n 个三次多项式来说，这就意味着需要 4n 个条件才能确定这些多项式。但是，插值特性只给出了 n + 1 个条件，内部数据点给出 n + 1 2 = n 1 个条件，总计是 4n 2 个条件。我们还需要另外两个条件，根据不同的因素我们可以使用不同的条件。其中一项选择条件可以得到给定 与 的钳位三次样条，s(x0) =s(xn) =另外，可以设s(x0) = s(xn) = 0这样就得到自然三次样条。自然三次样条几乎等同于样条设备生成