1、第四次实验报告信号、系统与信号处理实验课程实验报告实验四 :周期信号的傅里叶级数学 院 通信工程 班 级 14083414 学 号 14081401 姓 名 陈卓琳 指导教师 毕美华 2015年 12 月 11日实验四周期信号的傅里叶级数一、实验目的1、分析典型的矩形脉冲信号,了解矩形脉冲信号谐波分量的构成。2、观察矩形脉冲信号通过多个数字滤波器后,分解出各谐波分量的情况。3、掌握用傅里叶级数进行谐波分析的方法。4、观察矩形脉冲信号分解出的各谐波分量可以通过叠加合成出原矩形脉冲信号。二、实验原理1. 信号的时间特性与频率特性信号可以表示为随时间变化的物理量,比如电压和电流等,其特性主要表现为随
2、时间的变化,波形幅值的大小、持续时间的长短、变化速率的快慢、波动的速度及重复周期的大小等变化,信号的这些特性称为时间特性。信号还可以分解为一个直流分量和许多不同频率的正弦分量之和。主要表现在各频率正弦分量所占比重的大小不同;主要频率分量所占的频率范围也不同,信号的这些特性称为信号的频率特性。无论是信号的时间特性还是频率特性都包含了信号的全部信息量。2. 信号的频谱信号的时间特性和频率特性是对信号的两种不同的描述方式。根据傅里叶级数原理,任意一个时域的周期信号,只要满足狄利克莱(Dirichlet)条件,就可以将其展开成三角形式或指数形式的傅里叶级数。例如,对于一个周期为T的时域周期信号,可以用
3、三角形式的傅里叶级数求出它的各次分量,在区间内表示为即将信号分解成直流分量及许多余弦分量和正弦分量,研究其频谱分布情况。图4-1 信号的时域特性和频域特性3. 信号的时间特性与频率特性关系信号的时域特性与频域特性之间有着密切的内在联系,这种联系可以用图4-1来形象地表示。其中图4-1(a)是信号在幅度-时间-频率三维坐标系统中的图形;图4-1(b)是信号在幅度-时间坐标系统中的图形即波形图;把周期信号分解得到的各次谐波分量按频率的高低排列,就可以得到频谱图。反映各频率分量幅度的频谱称为振幅频谱。图4-1(c)是信号在幅度-频率坐标系统中的图形即振幅频谱图。反映各分量相位的频谱称为相位频谱。4.
4、 信号频谱的测量在本实验中只研究信号振幅频谱。周期信号的振幅频谱有三个性质:离散性、谐波性、收敛性。测量时利用了这些性质。从振幅频谱图上,可以直观地看出各频率分量所占的比重。测量方法有同时分析法和顺序分析法。同时分析法的基本工作原理是利用多个滤波器,把它们的中心频率分别调到被测信号的各个频率分量上。当被测信号同时加到所有滤波器上,中心频率与信号所包含的某次谐波分量频率一致的滤波器便有输出。在被测信号发生的实际时间内可以同时测得信号所包含的各频率分量。5. 周期方波信号的傅里叶级数5.1一个周期为的正方波, 如图下图4-2所示,一个周期内从幅值为1,其余为0。图4-2 周期性的正方波信号由傅里叶
5、级数展开式可知,方波信号傅里叶级数系数为:;则该周期信号可以表示为傅里叶级数的形式:因为,当k为偶数时,所以进一步带入得当占空比为0.5时候的方波,即时可以看出方波各频率分量中,直流分量为0.5;偶次谐波分量为零;各奇次谐波分量比值为。5.2方波,如图4-3所示,可以分解为图4-3 周期性的方波信号可以看出方波各频率分量中,直流分量为零;偶次谐波分量为零;各奇次谐波分量比值为。6. 周期信号的合成吉布斯现象(Gibbs)根据傅里叶级数可以将周期信号分解成直流分量、基波以及各次谐波分量,同样,由直流分量、基波和各次谐波分量可以叠加出来一个周期信号。例如前述的方波信号,可以由其基波和各次谐波分量按
6、照比例叠加出来,合成方波信号与原信号的误差取决于傅里叶级数的项数。合成波形所包含的谐波分量越多,它越逼近原方波信号,但是间断点除外。用有限项傅里叶级数表示有间断点的信号时,在间断点附近不可避免的会出现振荡和超量。超量的幅度不会随所取项数的增加而减小。只是随着项数的增多,振荡频率变高,并向间断点处压缩,从而使它所占有的能量减少。这种现象称为吉布斯现象。三、实验内容:1.信号的分解与合成Matlab仿真实验1.1方波信号的分解思考:参考下面程序,将频率为50Hz幅值为3的方波进行分解,给出前5项谐波,并在不同坐标系和同一坐标系下绘制各次谐波波形。参考程序及运行结果如下:t=0:0.01:2*pi;
7、 %时间间隔为0.01y=zeros(10,max(size(t); %10*629(t的长度)的矩阵x=zeros(10,max(size(t);for k=1:2:9 %奇次谐波1,3,5,7,9 x1=sin(k*t)/k; %各次谐波正弦分量 x(k,:)=x(k,:)+x1; %x第k(1,3,5,7,9)行存放k次谐波的629个值 y(k+1)/2,:)=x(k,:); %矩阵非零行向量移至1-5行endsubplot(2,1,1);plot(t,y(1:5,:); %绘制y矩阵中1-5行随时间波形grid;halft=ceil(length(t)/2); %行向量长度减半(由于对
8、称前后段一致)subplot(2,1,2); %绘制三维图形:矩阵y中全部行向量的一半mesh(t(1:halft),1:10,y(:,1:halft);图像:第一题程序:t=0:0.0001:0.02; y=zeros(10,max(size(t); x=zeros(10,max(size(t);for k=1:2:9 x1=3*sin(k*100*pi*t)/k; x(k,:)=x(k,:)+x1; y(k+1)/2,:)=x(k,:);endsubplot(2,1,1);plot(t,y(1:5,:); grid;halft=ceil(length(t)/2); subplot(2,1,
9、2); mesh(t(1:halft),1:10,y(:,1:halft);图像:1.2方波信号的合成思考:(1)参考实验原理内容解释下面程序中“f1=f1+cos(pi*n*t)*sinc(n/2)”;观察N值改变时合成波形的变化。(2)参考下面程序分别对周期为4的方波和周期为2的方波进行合成,注意比较:原方波与合成方波;两个方波合成有何不同;当傅里叶级数项数增加时合成方波的变化。答:(1) “f1=f1+cos(pi*n*t)*sinc(n/2)”的意思是直流分量和各奇次谐波分量的累加,等价于数学公式:;N发生变化:t=-4.5:0.001:5.5;t1=-4.499:0.001:5.5;
10、x=ones(1,1000),zeros(1,1000);x=x,x,x,x,x;subplot(2,2,3);plot(t1,x,linewidth,1.5);axis(-4.5,5.5,-0.5,1.5);N=100;c0=0.5;f1=c0*ones(1,length(t)for n=1:Nf1=f1+cos(pi*n*t)*sinc(n/2);endsubplot(2,2,4);plot(t,f1,r,linewidth,1.5);axis(-4.5,5.5,-0.5,1.5);t=-4.5:0.001:5.5;t1=-4.499:0.001:5.5;x=ones(1,1000),ze
11、ros(1,1000);x=x,x,x,x,x;subplot(2,2,1);plot(t1,x,linewidth,1.5);axis(-4.5,5.5,-0.5,1.5);N=10;c0=0.5;f1=c0*ones(1,length(t)for n=1:Nf1=f1+cos(pi*n*t)*sinc(n/2);endsubplot(2,2,2);plot(t,f1,r,linewidth,1.5);axis(-4.5,5.5,-0.5,1.5);图像:由上图可知,N越大,合成方波波形的上下起伏越小。(2)方波x1(t):t=-4.5:0.001:5.5;t1=-8.999:0.002:1
12、1;x=ones(1,1000),zeros(1,1000);x=x,x,x,x,x;subplot(1,2,1);plot(t1,x,linewidth,1.5);axis(-4.5,5.5,-0.5,1.5);N=10;9ic0=0.5;f1=c0*ones(1,length(t)for n=1:Nf1=f1+cos(0.5*pi*n*t)*sinc(n/2);endsubplot(1,2,2);plot(t,f1,r,linewidth,1.5);axis(-4.5,5.5,-0.5,1.5);图像:方波x2(t):t=-4.5:0.001:5.5;t1=-4.499:0.001:5.5
13、;x=ones(1,1000),zeros(1,1000);x=x,x,x,x,x;subplot(1,2,1);plot(t1,x,linewidth,1.5);axis(-4.5,5.5,-0.5,1.5);N=10;c0=0.5;f1=c0*ones(1,length(t)for n=1:N f1=f1+cos(pi*n*t)*sinc(n/2);endsubplot(1,2,2);plot(t,f1,r,linewidth,1.5);axis(-4.5,5.5,-0.5,1.5);图像:合成方波与原方波的比较如图所示,周期、相位都相同,合成方波在正负幅值处起伏较大。两方波合成后的合成方
14、波只有周期不同。傅里叶级数项数增加相当于N变大,N变大,合成方波上下起伏变小,如(1)中图所示。四、实验小结通过本次实验,我了解了矩形脉冲信号,掌握了用傅里叶级数进行谐波分析的方波,收获了以下几点知识:(1) 矩形脉冲信号谐波分量由直流分量(常数);偶次谐波分量(零);奇次谐波分量(比值为)构成。(2) 方波傅里叶级数的表示:对应程序语言:c0=0.5;f1=c0*ones(1,length(t)for n=1:Nf1=f1+cos(pi*n*t)*sinc(n/2);end(3) 第一题中,频率发生改变则修改t的值。运用公式周期等于频率的倒数进行计算。频率变为50Hz,对应的周期为0.02,所以t的值修改为t=0:0.0001:0.02;
copyright@ 2008-2023 冰点文库 网站版权所有
经营许可证编号:鄂ICP备19020893号-2