实验一方波分解为多次正弦波

1、实验一方波分解为多次正弦波之和的设计一、实验目的1、了解信号分解与正交函数集；2、了解三角函数的正交性；3、了解基波与谐波的概念与关系；二、实验原理1、信号分解与正交函数集信号通常以时间函数表示，所以信号的分量及其分解指的就是函数的分量及 其分解。可利用与矢量分解相类比的方法来来研究如何将一信号分解为其分量。与矢量用另一矢量上的分量表示原矢量类似。在一定的时间区间tl,t2内，若用函数f2（t）中的G2f2（t）来近似的表示原函数fl（t），将存在一误差函数 （t）， 且有：（t） fl G2f2（t）（ 1）在矢量近似中最佳系数选择的依据，是使得误差矢量的长度的平方最小；而5的选择，则是要求

2、使误差函数的方均值最小。误差函数的方均值为：荀 仆 t2）F 2（t）dt（ 2）此值最小时有：C12 t2 fi（t）f2（t）dt t2f；（t）dt（ 3）系数C（2是在最小方均误差的意义上代表二函数存、f2（t）的相关联的程度的度量。当C12 0时，由式（3）可知，此时有：t2 f1（t）f2（t）dt 0（ 4）如果满足这一条件，则称的与f2（t）在区间魚龙）内正交。此时f1（t）与f2（t） 就构成一个正交函数集。f,（t）与f2（t）两函数正交时，f1（t）在f2（t）中的分量一个函数可以在另一个函数中具有分量，则和矢量的情况类似，可以将一代表信号的函数表示为该函数在一正交函数集

3、中的分量的加权和。在区间(tit)内互相正交的n个函数g1(t)> g2(t)、g3(t)、gn(t)组成一个n维的正交信号空间。此函数集中的函数之间，在区间(tit)内具有如下关系：2gm(t)dt km( 5)：gi(t)gm(t)dt 0 l m( 6)其中km为一常数。若km 1，则上述函数集就称为是归一化正交的。任意- 个代表信号的的函数f(t)，在区间(tit)内，可以用组成信号空间的n个正交函 数的线性组合来近似的表示为:f (t) qgi(t) C2g2(t) . Crgr(t)qgn(t)( 7)若要使近似值的方均误差最小，则f (t)在函数gr (t)中的分量系数为：

4、Cr$f (t)gr(t)dt： $g2(t)dt 1 kr ti2 f (t)gr(t)dt( 8)若用一正交函数集中的分量和各次谐波分量之和，那么该矢量集必须是一完备的正交矢量集。与此相似，用一正交函数集中的分量去代表任意一个函数，该函数集也必须是一完备的正交函数集。完备的正交函数集往往都是由无穷多的函 数组成。任意一信号表示为正交函数集中的分量之和时，所取分量函数的项数越多精确度越高，即方均误差最小。当所取项数无限增大时，方均误差趋于零，这是的正交函数集也成为完备的。对于一个在区间(t1,t2)内的完备正交函数集中的所有函数，不可能找到另外一个异于零的函数能在同一区间内和它们相正交。 反

5、 之，若存在正交函数集以外的函数， 与正交函数集中的所有函数正交， 则该正交 函数集必不是完备的。信号在正交函数集中的分解是多样。在矢量分解中，坐标轴经过变换，可以 有不同的选取方法；同样，表示信号的正交函数集也可以经过变换而有不同选取 方法。如同坐标变换不影响矢量本身一样，正交函数集的变换也不影响所表示的 函数本身。故可以用一个正交函数集变换到另一个正交函数集去表示一个函数。在各种正交函数集中，傅里叶级数是既方便又很有用的。 除傅里叶级数外，其他 如沃尔什函数、勒让德函数、切比雪夫函数等。2、三角函数的正交性正弦函数与余弦函数满足如下关系：t1 T cosm tsinn tdt 0 m,n

6、为任意整数(9)t T cosm t cosn tdt 1 Tsinm tsinn tdt 0 m n(10)£ T cos2 n tdtTsin2n tdt T 2( 11)这表示上述各余弦函数与正弦函数，在区间(ti,ti T)内，均相互正交。即1、 cos t、 cos2 t、cos3 t、 cosn t、 sin t、 sin2 t、 sin3 t、sinn t、(其中T 2 ；'为上述三角函数的公共周期) 构成一组完备的正交函数集。故任意一个周期信号都可以展开成为正弦函数及余 弦函数的无穷级数。3、连续时间周期信号分解为三角函数之和由前一部分可知，对于任意一个周期为

7、T的周期信号f(t),都可以求出它在 上述三角函数集中各函数中的分量，从而可将 f(t)在区间(t1, t1 T)内表示为上 述三角函数集中各函数的加权和。即：f(t) a 2ancosn t bnsinn t(12)n 1其中a 2实际就是f(t)在区间内的平均值，意即直流分量。设定其为'2是 为了使得后面出计算ao的公式与其他下标不等于零的an 一致。由公式(8)、( 9)、(10)、( 11)可知，各正弦函数余弦函数的分量系数为：ant1T f (t)cos n tdt /T cos2 n tdt 2/T , T f (t )cos n tdt(13)bnt1T f(t)sinn

8、 tdt/1Tsin2n tdt2/TT f(t)sinn tdt(14)但是，要将周期信号分解为谐波分量，代表该周期信号的函数f(t)应当满足 狄利克雷条件。即：1) 在一周期内，函数式绝对可积的，即：2 T|f(t)|dt应为有限值；2）在一周期内，函数的极值数目有限；3）在一周期内，函数f （t）或者为连续的，或者具有有限个这样的间断点，即当t从较大的时间值和较小的时间值分别趋向于间断点时，函数具有两个不同的有限的极限值。实际工程中的周期信号，大多都满足狄利克雷条件。周期性方波也满足上述 条件，即方波可以展开为三角函数的加权和。4、方波分解为多次正弦波之和的原理由前一部分可知。代表周期性

9、方波信号的函数f（t）满足狄利克雷条件，即方 波可以表示为多次正弦波之和。如图1所示方波信号，其周期为2且正半周期负 半周期是形状全同的矩形，在区间（0，2）内可用函数表示为：方液1-煙0-1J100511.S图1周期为2的方波信号若将f（t）展开为三角傅里叶级数，即将f（t）分解为多次正弦波之和，则有式（13）、式（14）可知，在区间（0, 2）内，如图1所示的周期为2的方波信号 的a，an，bn的值分别为：ao = 22 o f (t)dt 0an = 2 2 o f (t)cosn tdt 0bn = 2 2 o f (t)sinn tdt 0 当 n为偶数bn22 0 f (t)si

10、nn tdt 4/n 当 n 为奇数则在区间(0, 2)内f(t)可表示为：f(t) 4 (sin t 'sin3 t sin5 t Jsin7 t )(15)即周期为2s的方波信号中含有大量的正弦波，其频率分别为1/2,3/2,5/2 , 7/2 其中频率为1/2的正弦波称为基波，其他频率的正弦波称为谐波。即 一周期性方波，可表示为基波与无穷多谐波之和。实用中进行信号分析时，不可能取无穷多次谐波之和，而只能用有限项来近 似表示。这样就无法避免有一误差 (t)，如果将基波加到n次谐波之和后的函数 表示为 fn(t)，则有 f(t) = fn(t)+ (t)，即(t)=f(t)- (t)

11、。5、正弦波合成并与原始方波进行比较模型的建立由原理部分可知，方波中含有无穷多次的正弦波，即方波可以分解为无穷多 次正弦波之和。反过来可用无限多次正弦波相加可以合成方波，从而完成设计。 但是在实际中不可能取无穷中多次的正弦波，取有限次正弦波合成，相较于原始方波信号，必然存在一定的误差。若用误差函数表示误差。若方波信号函数表示为f(t)，多次正弦波合成后函数表示为fn(t)，则误差函数(t)=f(t)- fn(t)。 然后利用相关绘图函数画出f(t)、fn(t)、(t)即可。故设计建立的模型可表示如下：图2正弦波合成程序设计模型图示描述其中产生方波可以调用square函数。对于多次正弦波的合成可

12、以使用一循环语句，通过使用in put输入循环的最大值N,来控制合成的谐波的数量。然后 通过subplot函数建立三个子窗口，运用 plot函数在三个子窗口中分别画出原 始方波、合成波以及误差函数的波形图即可。同时，利用subplot函数建立三个子窗口，将原始方波，合成波及误差函数画在同一窗口下，可以更加直观的比较合成波与原始方波及观察误差函数。三、实验内容1、正弦波合成并与原始方波比较的源程序代码正弦波合成并与原始方波进行比较的源程序代码如下：%多次正弦波合成并与原始方波比较的设计t=0:0.000111:6;n=1;fn=0;y=square(pi气50);%周期为2s的方波N=input

13、( 'N=' ); %输入 N值for num=1:Nfn=fn+4/(n*pi)*sin(n*t*pi);%加至U n次谐波正弦波之和n=n+2;endfm=y-fn; %误差函数subplot(3,1,1)%在第一个子窗口中画出原始方波图像plot(t,y, 'k' ),grid onaxis(0 6 -2 2)xlabel('时间 s' );ylabel('振幅');title('方波')subplot(3,1,2)%在第二个子窗口中画出加到n次谐波的正弦波合成图像plot(t,fn, 'k'

14、 ),grid onxlabel('时间 s ' );ylabel('振幅');title('正弦波合成')subplot(3,1,3)%在第三个子窗口中画出误差函数图像plot(t,fm, 'k'),axis(0 6 -2 2)xlabel（'时间 s' ）;ylabel（'振幅'）;title（'误差函数'）源程序中，通过输入的N值控制for循环的次数，从而控制合成的谐波次数 n以及合成的谐波数量。通过改变输入的N可以改变合成波的谐波的最高次数即 相加的谐波次数，从而改变合成波的

15、波形。通过比较加到不同次数的合成波的波 形及其相应的误差函数波形，可以知道正弦波合成随级数增大的趋势情况。2、程序运行结果1）N=3时，原始方波，合成波及误差函数图像如下：图3原始方波，加到五次谐波的合成波及相应误差函数图像2）N=10时，原始方波，合成波及误差函数图像如下：图4原始方波，加到19次谐波的合成波及相应误差函数图像3) N=100时，原始方波，合成波及误差函数图像如下：图5原始方波，加到199次谐波的合成波及相应误差函数图像由图上可以看出，当输入N=100时，合成波波形已经非常接近方波波形。 除 在方波跃变点附近合成波的波形与原始方波波形之间有较大误差之外， 其他各处 波形误差为

16、零。4) N=10000时，原始方波，合成波及误差函数图像如下：图6原始方波，加到19999次谐波的合成波及相应误差函数图像分析随着N值得增大，正弦波合成方波的效果如何变化？分析问什么会出现 这样的显现3、正弦波合成趋势图源程序代码1)正弦波合成趋势图二维图源程序代码如下：%皆波合成趋势二维图t=0:0.0111:6;y=zeros(1000,max(size(t);%建立全零矩阵x=zeros(size(t);n=1;N=input( 'N=' ); %俞入 N值for k=1:Nx=x+(4/( n*pi)*si n(n *pi*t);y(n, : )=x;n=n+2;en

17、dplot(t,y(1:101,: ),gridonaxis(0 6 -2 2)xlabel('时间s' );ylabel('振幅');title('正弦波合成趋势二维图')该程序中，首先通过zeros函数建立相关全零矩阵并赋给x、y。然后利用 for循环语句求出加到不同次数谐波的合成波。最后利用所建立的矩阵及plot函数画出合成波随级数增加合成波改变的二维趋势图。通过输入N值控制级数的项数。2)正弦波合成趋势图三维图源程序代码如下：%皆波合成趋势三维图t=0:0.0111:2;y=zeros(1000,max(size(t);x=zeros(s

18、ize(t);for k=1:2:200x=x+sin(k*pi*t)/k; y(k+1)/2, : )=x;endhalft=ceil(le ngth(t)/2);mesh(t(1:halft),1:1000,y(: ,1:halft)xlabel('时间 s' );ylabel('谐波次数');zlabel('振幅');title('谐波合 成趋势三维图)基本原理与上述二维图类似，通过建立全零矩阵，然后利用所建立的矩阵以 及mesh函数画出三维图。为使图形更加美观，只画出了合成波形的半个周期， 即区间(0,1)内的合成波波形。该程序中也通过输入N值控制级数的项数。并且在图形中利用y轴表示级数项数的增加，这样可以更直观看到合成波随级数项数 的变化情况。4