武波形合成与分解.docx

武汉大学教学实验报告 电子信息学院电子信息科学与技术专业 2010年04月25日 实验名称波形合成与分解 实验名称波形合成与分解 指导教师 卜方玲 姓名戴永珊 年级2008级学号2008301200097成绩 预习部分 实验目的 实验基本原理 主要仪器设备（含必要的元器件、工具） 实验目的 在理论学习的基础上，通过本实验熟悉信号的合成、分解原理，了解信号频谱的含义， 加深对傅里叶变换性质和作用的理解。 实验基本原理 根据傅里叶分析的原理，任何周期信号都可以用一组三角函数 {Sin（ ，0t）； C0S（ n，0t）}的 组合表示，即： x（t） = a。*1 cos（ ? pt） - 0 sin（ ? ^t） - a? cos（ 2 ^t） - b? sin（ 2 ^t） 即可以用一组正弦波和余弦波来合成任意形状的周期信号。 主要仪器设备（含必要的元器件、工具） MATLAB^ 件 实验内容： （1）方波的合成 图示方波是一个奇谐信号， 由傅里叶级数可知，它是由无穷个奇次谐波分量合成的， 本实 验用图形的方式来表示它的合成。方波信号可以分解为 ： 1? ,n =1,3,5,7,9,*nx(t) =2A sin( 2 1 ? ,n =1,3,5,7,9,* n 二 n 土 所以可以用前5项谐波近似合成50Hz,幅值为3的方波 （2）设计谐波合成三角波的实验 三角波是一个奇谐信号，只含有正弦分量 4 E处 n兀 1 x(t) ： 一2 a si n( )si n(2 二 n f°t)二，n = 1,3, 5, 7,9,… H n A 2 n 设计分析方波、三角波频谱的分析实验，写出实验步骤，并完成实验(并比较二者频 谱的特点)。 二、实验操作部分 实验数据、表格及数据处理 实验操作过程(可用图表示) 实验结论 实验数据、表格及数据处理 (1)用前5项谐波近似合成50Hz,幅值为3的方波的代码如下：(代码一) clear all close all f0=5; %基波频率 n=10; %最高谐波次数 f=n *f0; %最大频率 T=1/(20*f); %采样间隔 t=0:T:1; %观察0-1秒的范围 A=4.7; %采样幅度=2A/pi=3 Y1=2*A*si n(2*pi*f0*t)/pi; plot(t,Y1); hold on; Y2=2*A*si n(2*pi*f0*t)/pi+2*A*si n(2*3*pi*f0*t)/(3*pi); plot(t,Y2); hold on; Y3=2*A*si n(2*pi*f0*t)/pi+2*A*si n(2*3*pi*f0*t)/(3*pi)+2*A*si n(5*2*pi*f 0*t)/(5*pi); plot(t,Y3); hold on; Y4=2*A*si n(2*pi*f0*t)/pi+2*A*si n(2*3*pi*f0*t)/(3*pi)+2*A*si n(5*2*pi*f 0*t)/(5*pi)+2*A*si n(7*2*pi*f0*t)/(7*pi); plot(t,Y4); hold on; Y5=2*A*si n(2*pi*f0*t)/pi+2*A*si n(2*3*pi*f0*t)/(3*pi)+2*A*si n(5*2*pi*f 0*t)/(5*pi)+2*A*si n(7*2*pi*f0*t)/(7*pi)+2*A*si n(9*2*pi*f0*t)/(9*pi); plot(t,Y5); hold on; (2)用前5项谐波近似合成50Hz,幅值为1.5的三角波的代码如下：(代码二) clear all close all f=5; %基波频率 n=10; %最高谐波次数 fm=n*f; %最大频率 T=1/(20*f); %采样间隔 t=0:T:1; %观察0-1秒的范围 E=3; %控制采样幅度 Y1=4*E*si n(1*pi/2)*si n(1*2*pi*f*t)/(1*1*pi*pi); plot(t,Y1); hold on; Y2=4*E*si n(1*pi/2)*si n(1*2*pi*f*t)/(1*1*pi*pi)+4*E*si n( 3*pi/2)*si n(3* 2*pi*f*t)/(3*3*pi*pi); plot(t,Y2); hold on; Y3=4*E*si n(1*pi/2)*si n(1*2*pi*f*t)/(1*1*pi*pi)+4*E*si n( 3*pi/2)*si n(3* 2*pi*f*t)/(3*3*pi*pi)+4*E*si n(5*pi/2)*si n(5*2*pi*f*t)/(5*5*pi*pi); plot(t,Y3); hold on; Y4=4*E*si n(1*pi/2)*si