(精选)信号与系统实验(MATLAB 西电版)实验10 周期信号的合成与分解课件.ppt

一、实验目的 二、实验原理 三、涉及的MATLAB函数 四、实验内容与方法 五、实验要求 六、思考题;一、 实验目的 　　(1) 在理论学习的基础上，通过本实验熟悉信号的合成、分解原理，加深对傅里叶级数的理解; 　　(2) 了解和认识吉布斯现象 (Gibbs)。二、 实验原理　　任何具有确定性的信号都可以表示为随时间变化的某种物理量，比如电压u(t)和电流i(t)等。信号主要表现在随着时间t的变化，波形幅值的大小、持续时间的长短、变化速率的快慢、波动的速度及重复周期的大小的变化等。信号的这一特性称为信号的时间特性。 ;　　信号还可以分解为一个直流分量和许多不同频率的正弦分量之和。主要表现在各频率正弦分量所占比重的大小不同，主要频率分量所占有的频率范围也不同等，信号的这一特性称为信号的频率特性。　　无论是信号的时间特性，还是信号的频率特性，都包含了信号的全部信息量。　　根据周期信号的傅里叶级数展开式可知，任何非正弦周期信号，只要满足狄里赫利条件，都可以分解为一直流分量和由基波及各次谐波(基波的整数倍)分量的叠加。 ;　　例如一个周期的方波信号f(t)可以分解为如下形式：　　f(t)=如图10.1(a)所示。　　同样，由基波及各次谐波分量也可以叠加出来一个周期方波信号，如图10.1(b)所示。　　 至于叠加出来的信号与原信号的误差，则取决于傅里叶级数的项数。 ;三、 涉及的MATLAB函数略四、 实验内容　　1. 验证性实验　　1) 周期信号的分解MATLAB程序：clf; %周期信号的分解t=0：0.01：2*pi; y=zeros(10，max(size(t))); x=zeros(10，max(size(t))); for k=1：2：9;　 x1=sin(k*t)/k;  x(k，：)=x(k，：)+x1;  y((k+1)/2，：)=x(k，：); endsubplot(2，1，1); plot(t，y(1：9，：)); grid; line(［0，pi+0.5］，［pi/4，pi/4］); text(pi+0.5，pi/4，′pi/4′); halft=ceil(length(t)/2); subplot(2，1，2); mesh(t(1：halft)，［1：10］，y(：，1： halft));周期信号的分解如图10.2所示。;图 10.2 周???信号的分解;　　2) 傅里叶级数逼近MATLAB程序： clf; %宽度为1，高度为1，周期为2的正方波，傅里叶级数逼近 t=-2：0.001：2; %信号的抽样点N=20; c0=0.5; f1=c0*ones(1，length(t)); %计算抽样上的直流分量for n=1：N %偶次谐波为零 f1=f1+cos(pi*n*t)*sinc(n/2); endplot(t，f1); axis(［-2 2 -0.2 0.8］);方波的傅里叶级数逼近如图10.3所示。;图 10.3 方波的傅里叶级数逼近;　　3) 用正弦信号的叠加近似合成一频率为50 Hz，幅值为3的方波MATLAB程序： clear all; fs=10000; t=［0：1/fs：0.1］; f0=50; sum=0; subplot(211)for n=1：2：9; ; plot(t，4/pi*1/n*sin(2*pi*n*f0*t)，′k′);  title(′信号叠加前′);  hold on; endsubplot(212)for n=1：2：9;  sum=sum+4/pi*1/n*sin(2*pi*n*f0*t); endplot(t，sum，′k′); title(′信号叠加后′);正弦信号的叠加如图10.4所示。;图 10.4 正弦信号的叠加;　　4) Gibbs现象　　执行下列程序，令N分别为10，20，30，40，50，观察波形的特点，了解吉布斯现象的特点。MATLAB程序： t=-1.5：0.01：1.5; wo=4，E=1; N=10; xN=0; for n=1：N; an=(E/(n*pi))*(sin(n*pi/2)-sin(n*3*pi/2)) xN=xN+an.*cos(n*wo*t); endsubplot(221); plot(t，xN)xlabel(′time′); 