实验二 信号分解与合成
--谢格斯 110701336 聂楚飞110701324
一、实验目的
1、观察电信号的分解。
2、掌握带通滤波器的有关特性测试方法。
3、观测基波和其谐波的合成。
二、实验内容
1、观察信号分解的过程及信号中所包含的各次谐波。
2、观察由各次谐波合成的信号。
三、预备知识
1、了解李沙育图相关知识。
2、课前务必认真阅读教材中周期信号傅里叶级数的分解以及如何将各次谐波进行叠加等相关内容。
四、实验仪器
1、信号与系统实验箱一台(主板)。
2、电信号分解与合成模块一块。
3、20M 双踪示波器一台。
五、实验原理
任何电信号都是由各种不同频率、幅度和初相的正弦波迭加而成的。对周期信号由它的傅里叶级数展开式可知,各次谐波为基波频率的整数倍。而非周期信号包含了从零到无穷大的所有频率成份,每一频率成份的幅度均趋向无限小,但其相对大小是不同的。
通过一个选频网络可以将电信号中所包含的某一频率成份提取出来。本实验采用性能较佳的有源带通滤波器作为选频网络,因此对周期信号波形分解的实验方案如图2-3-1所示。
将被测方波信号加到分别调谐于其基波和各次奇谐波频率的一系列有源带通滤波器电路上。从每一有源带通滤波器的输出端可以用示波器观察到相应频率的正弦波。本实验所用的被测信号是Hz 531=ω左右的周期信号,而用作选频网络的五种有源带通滤波器的输出
频率分别是543215432ωωωωω、、、、
,因而能从各有源带通滤波器的两端观察到基波和各次谐波。其中,在理想情况下,如方波的偶次谐波应该无输出信号,始终为零电平,而奇次谐波则具有很好的幅度收敛性,理想情况下奇次谐波中一、三、五、七、九次谐波的幅度比应为1:(1/3):(1/5):(1/7):(1/9)。但实际上因输入方波的占空比较难控制在50%,且方波可能有少量失真以及滤波器本身滤波特性的有限性都会使得偶次谐波分量不能达到理想零的情况。
六、实验步骤
1、把系统时域与频域分析模块插在主板上,用导线接通此模块“电源接入”和主板上的电源(看清标识,防止接错,带保护电路),并打开此模块的电源开关。
2、调节函数信号发生器,使其输出Hz 53左右(其中在Hz Hz 56~50之间进行选择,
使其合成的效果更好)的方波(要求方波占空比为50%,这个要求较为严格),峰峰值为5V 左右。将其接至该实验模块的各带通滤波器的“输入”端,用示波器观察各带通滤波器的输出。(注:观察频率时,可打开实验箱上的频率计实验模块。即按下该模块电源开关S2。) 3、用示波器的两个探头,直接观察基波与三次谐波的相位关系,或者采用李沙育图的方法,同时考察其幅度关系,看其相位差是否为零,幅度之比是否为3:1(可以用相应带通滤波器中的调幅和调相电位器进行相关的调节,保证了相位和幅度满足实验的要求,以下的步骤中均可用到调相和调幅,使我们认识到调相和调幅在信号分解和合成的重要性)。
4、将方波分解所得基波和三次谐波,用导线与其对应的插孔相连,观测加法器的输出“合成”波形,并记录所得的波形。
5、同时考察基波、三次谐波、五次谐波的相位和幅度的关系,还是用李沙育图观察其相位关系,用观察法使其幅度关系为5:3:1,。
6、验证各高次谐波与基波之间的相位差是否为零。可用李沙育图形法进行测量,其方法如下:
用导线将函数发生器的方波输出端与带通滤波器输入端连接起来,即把方波信号分先后送入各带通滤波器,如图2-3-1所示。
图2-3-1 信号分解的过程
具体方法一:基波与标准同频同相信号相位比较(李沙育相位测量法) 把函数信号发生器模块产生的正弦波电压调至5V (峰峰值),使其送入示波器的X 轴,再把BPF-1ω的基波送入Y 轴,示波器采用X-Y 方式显示,观察李沙育图形。(注:当滤波器的增益不为1时,即X 轴和Y 轴信号幅度不一致时,在0
90=Φ时其李沙育图形并不为圆,而是椭圆,但其是垂直椭圆,与0
900<Φ<时的椭圆并不相同。)
当两信号相位差为00时,波形为一条直线;当两信号相位差为0
90时,波形为一个圆;当两信号相位差为0
900<Φ<时,波形为椭圆,如图3-2所示。
00900<Φ<时:??
?
??=ΦB A arcsin
B
A 00=Φ 090=Φ 900<Φ<
图2-3-2 李沙育图形
具体方法二:基波与各高次谐波相位比较(李沙育频率测试法)
把BPF-1ω处的基波送入示波器的X 轴,再分别把BPF-13ω、BPF-15ω处的高次谐波送入Y 轴,示波器采用X-Y 方式显示,观察李沙育图形。
当基波与三次谐波相位差为0
0(即过零点重合)、0
90、1800时,波形分别如图3-3所示。
00=Φ 090=Φ 0180=Φ
图2-3-3 基波与三次谐波相位的观察 以上是三次谐波与基波产生的典型的Lissajous 图,通过图形上下端及两旁的波峰个数,确定频率比,即3:1,实际上可用同样的方法观察五次谐波与基波的相移和频比,其应为5:1。
7、方波波形合成
(1)将函数发生器输出的Hz 53左右(其中在Hz Hz 56~50之间进行选择,使其输出的效果更好)方波信号送入各带通滤波器输入端。
(2)在五个带通滤波器输出端逐个测量各谐波输出幅度,
(3)用示波器观察并记录加法器输出端基波与各奇次谐波的叠加波形,如图2-3-4所示。
图2-3-4 基波与三次和五次谐波叠加后的波形七、实验报告
比较基波与三次谐波不同相位差的波形和频率幅度
相位差为φ=0o相位差为φ=90o
相位差Φ=180o
基波、三次谐波、五次谐波间的合成。
基波与三次谐波的合成基波与五次谐波的合成
基波、三次谐波、五次谐波的合
三、
八、分析相位、幅值在波形合成中的作用
相位对单个波形来说影响波的位置,即改变初相位能使波在时域坐标轴上左右移动。在波的叠加过程中,波的相位则会使得叠加波形的形状发生很大改变
而幅度的话影响合成波的幅度大小
实验一周期信号的分解与合成 一、实验目的 1.用同时分析法观测50Hz 非正弦周期信号的频谱。 2.观测基波和其谐波的合成。 二、实验原理 1.一个非正弦周期函数可以用一系列频率成整数倍的正弦函数来表示,其中与非正弦具有相同频率的成分称为基波或一次谐波,其它成分则根据其频率为基波频率的2、3、4、...、n 等倍数分别称二次、三次、四次、...、n 次谐波,其幅度将随谐波次数的增加而减小,直至无穷小。 2.不同频率的谐波可以合成一个非正弦周期波,反过来,一个非正弦周期波也可以分解为无限个不同频率的谐波成分。 3.一个非正弦周期函数可用傅里叶级数来表示,级数各项系数之间的关系可用一各个频谱来表示,不同的非正弦周期函数具有不同的频谱图,各种不同波形及其傅氏级数表达式见表1-1 表1-1 各种不同波形的傅里叶级数表达式(下) 1.方波
2.三角波 3.半波 4.全波 5.矩形波 三、预习要求 在做实验前必须认真复习教材中关于周期性信号傅利叶级数分解的有关内容。 四、实验内容 1. 50HZ方波信号的频谱。 2. 周期矩形脉冲的频谱;脉冲宽度为1;周期为4;则基波角频率为0.5pi 3. 使用不同频率的谐波合成方波信号;注意观察随着谐波数的增加合成的波形发生的变化。 4. 使用不同频率的谐波合成矩形脉冲信号;注意观察随着谐波数的增加合成的波形。 五、思考题 1.什么样的周期性函数没有直流分量和余弦项?
附: 1. 50HZ方波信号的频谱。 >> w1= ; %基波角频率 >> n=0:1:30; >>bn= ; %三角级数中系数bn,参考书p122 >> stem(n*w1,bn),grid on >> xlabel('\omega(rad/s)'),ylabel('bn') >> title('方波信号频谱分析图') 2. 周期矩形脉冲的频谱;脉冲宽度为1;周期为4;则基波角频率为0.5pi tao= ; w1= ; n=-15:1:15; fn= ; %矩形脉冲级数系数fn,参考书p130,用matlab自带函数sinc stem(n,fn),grid on xlabel('n'); ylabel('Fn'); title('周期矩形脉冲的频谱图'); 3. %使用不同频率的谐波合成方波信号;注意观察随着谐波数的增加合成的波形 %发生的变化。 t=-1:0.001:1; omega=2*pi; y=square(2*pi*t,50); plot(t,y);grid on xlabel('t'); ylabel('周期方波信号'); axis([-1 1 -1.5 1.5]); n_max=[1 3 5 11 47]; N=length(n_max); for k=1:N n=1:2:n_max(k); b=4./(pi*n); x=b*sin(omega*n'*t); figure; plot(t,y) hold on; plot(t,x); hold off; xlabel('t'); ylabel('部分和的波形');
实验二信号分解与合成 --谢格斯110701336 聂楚飞110701324 一、实验目的 1、观察电信号的分解。 2、掌握带通滤波器的有关特性测试方法。 3、观测基波和其谐波的合成。 二、实验内容 1、观察信号分解的过程及信号中所包含的各次谐波。 2、观察由各次谐波合成的信号。 三、预备知识 1、了解李沙育图相关知识。 2、课前务必认真阅读教材中周期信号傅里叶级数的分解以及如何将各次谐波进行叠加等相关内容。 四、实验仪器 1、信号与系统实验箱一台(主板)。 2、电信号分解与合成模块一块。 3、20M双踪示波器一台。 五、实验原理 任何电信号都是由各种不同频率、幅度和初相的正弦波迭加而成的。对周期信号由它的 傅里叶级数展开式可知,各次谐波为基波频率的整数倍。而非周期信号包含了从零到无穷大的所有频率成份,每一频率成份的幅度均趋向无限小,但其相对大小是不同的。 通过一个选频网络可以将电信号中所包含的某一频率成份提取出来。本实验采用性能较 佳的有源带通滤波器作为选频网络,因此对周期信号波形分解的实验方案如图2-3-1所示。 将被测方波信号加到分别调谐于其基波和各次奇谐波频率的一系列有源带通滤波器电路上。从每一有源带通滤波器的输出端可以用示波器观察到相应频率的正弦波。本实验所用 的被测信号是 1 53Hz左右的周期信号,而用作选频网络的五种有源带通滤波器的输出 频率分别是「2 2、3 3、4 4、5 5,因而能从各有源带通滤波器的两端观察到基波和各 次谐波。其中,在理想情况下,如方波的偶次谐波应该无输出信号,始终为零电平,而奇次谐波则具有很好的幅度收敛性,理想情况下奇次谐波中一、三、五、七、九次谐波的幅度比应为1: (1/3):(1/5):(1/7):(1/9)。但实际上因输入方波的占空比较难控制在50%,且方 波可能有少量失真以及滤波器本身滤波特性的有限性都会使得偶次谐波分量不能达到理想零的情况。 六、实验步骤 1、把系统时域与频域分析模块插在主板上,用导线接通此模块“电源接入”和主板上 的电源(看清标识,防止接错,带保护电路),并打开此模块的电源开关。 2、调节函数信号发生器,使其输出53Hz左右(其中在50Hz ~ 56Hz之间进行选择,
数字信号处理作业提交日期:2016年7月15日
实验一 维纳滤波器的设计 第一部分 设计一维纳滤波器。 (1)产生三组观测数据,首先根据()(1)()s n as n w n =-+产生信号()s n ,将其加噪(信噪比分别为20,10,6dB dB dB ),得到观测数据123(),(),()x n x n x n 。 (2)估计()i x n ,1,2,3i =的AR 模型参数。假设信号长度为L ,AR 模型阶数为N ,分析实验结果,并讨论改变L ,N 对实验结果的影响。 1 实验原理 滤波技术是信号分析、处理技术的重要分支,无论是信号的获取、传输,还是信号的处理和交换都离不开滤波技术,它对信号安全可靠和有效灵活地传递是至关重要的。信号分析检测与处理的一个十分重要的内容就是从噪声中提取信号,实现这种功能的有效手段之一是设计一种具有最佳线性过滤特性的滤波器,当伴有噪声的信号通过这种滤波器的时候,它可以将信号尽可能精确地重现或对信号做出尽可能精确的估计,而对所伴随噪声进行最大限度地抑制。维纳滤波器就是这种滤波器的典型代表之一。 维纳(Wiener )是用来解决从噪声中提取信号的一种过滤(或滤波)方法。这种线性滤波问题,可以看做是一种估计问题或一种线性估计问题。 设一线性系统的单位样本响应为()h n ,当输入以随机信号()x n ,且 ()() () x n s n v n =+,其中()s n 表示原始信号,即期望信号。()v n 表示噪声,则输出()y n 为()=()()m y n h m x n m -∑,我们希望信号()x n 经过线性系统()h n 后得到的()y n 尽可能接近 于()s n ,因此称()y n 为估计值,用?()s n 表示。 则维纳滤波器的输入-输出关系可用下面表示。 设误差信号为()e n ,则?()()()e n s n s n =-,显然)(n e 可能是正值,也可能是负值,并且它是一个随机变量。因此,用它的均方误差来表达误差是合理的,所谓均方误差最小即 它的平方的统计期望最小:222?[|()|][|()()|][|()()|]E e n E s n s n E s n y n =-=-=min 。而要使均方误差最小,则需要满足2[|()|]j E e n h ?=0. 进一步导出维纳-霍夫方程为:()()()()*(),0,1,2...xs xx xx i R m h i R m i R m h m m =-==∑ 写成矩阵形式为:xs xx R R h =,可知:1xs xx h R R -=。表明已知期望信号与观测数据的互相关函数以及观测信号的自相关函数时,可以通过矩阵求逆运算,得到维纳滤波器的
电子工程系 信号与系统课程实验报告 2011-----2012学年第一学期 专业: 电子信息工程技术班级: 学号 : 姓名: 指导教师: 实常用连续时间信号的实现
一、实验目的 (1)了解连续时间信号的特点; (2)掌握连续时间信号表示的向量法和符号法; (3)熟悉MATLAB Plot函数等的应用。 二、实验原理 1、信号的定义 信号是随时间变化的物理量。信号的本质是时间的函数。 2、信号的描述 1)时域法 时域法是将信号表示成时间的函数f(t)来对信号进行描述的方法。信号的时间特性指的是信号的波形出现的先后,持续时间的长短,随时间变化的快慢和大小,周期的长短等。 2)频域(变换域)法 频域法是通过正交变换,将信号表示成其他变量的函数来对信号进行描述的方法。一般常用的是傅立叶变换。信号的频域特性包括频带的宽窄、频谱的分布等。 信号的频域特性与时域特性之间有着密切的关系。 3、信号的分类 按照特性的不同,信号有着不同的分类方法。 (1)确定性信号:可以用一个确定的时间函数来表示的信号。 随机信号:不可以用一个确定的时间函数来表示,只能用统计特性加以描述的信号。 (2)连续信号:除若干不连续的时间点外,每个时间点在t上都有对应的数值信号。离散信号:只在某些不连续的点上有数值,其他时间点上信号没有定义的信号。 (3)周期信号:存在T,使得等式f(t+T)=f(t)对于任意时间t都成立的信号。非周期信号:不存在使得等式f(t+T)=f(t)对于任意时间t都成立的信号。 绝对的周期信号是不存在的,一般只要在很长时间内慢走周期性就可以了。 (4)能量信号:总能量有限的信号。 功率信号:平均功率有限切非零的信号。 (5)奇信号:满足等式f(t)=--f(--t)的信号。偶信号:满足等式f(t)=f(--t)的信号。 三、涉及的MATLAB函数 1、plot函数 功能:在X轴和Y轴方向都按线性比例绘制二维图形。 调用格式: Plot(x,y):绘出相x对y的函数线性图。 Plot(x1,y1,x2,y2,…..):会出多组x对y的线性曲线图。 2、ezplot函数 功能:绘制符号函数在一定范围内的二维图形。简易绘制函数曲线。 调用格式: Ezplot (fun):在[-2π,2π]区间内绘制函数。 Ezplot (fun,[min,max]):在[min,max]区间内绘函数。 Ezplot (funx,funy):定义同一曲面的函数,默认的区间是[0, 2π]。】 3、sym函数 功能:定义信号为符号的变量。 调用格式:sym(fun):fun为所要定义的表达式。 4、subplot函数
实验指导书 实验项目名称:典型信号的合成和分解 实验项目性质:普 通 所属课程名称:工程测试技术 实验计划学时:2 一.实验目的 通过本实验熟悉信号的合成、分解原理,了解信号频谱的含义和特点。 二.实验内容和要求 1.周期信号的合成和分解 在有限区间内,凡满足狄里赫利条件的周期信号x(t)都可以展开傅里叶三角函数级数。 001001 ()(cos sin )2 cos()(1,2,3,)2n n n n n n n a x t a n t b n t a A n t n ωωω?∞=∞==++=+-=∑∑ 式中 0a ——常值分量 00/20/202()T T a x t dt T -=? n a ——余弦分量的幅值
00/20/202()cos T n T a x t n tdt T ω-=? n b ——正弦分量的幅值 00/20/202()sin T n T b x t n tdt T ω-=? n A ——n 次谐波的振幅,是n 的偶函数 n A = n ?——n 次谐波的相角,是n 的奇函数 arctan n n n a b ?= 可见,周期信号是由周期信号是由一个或几个、乃至无穷多个不同频率的谐波叠加而成的。也就是说,复杂周期信号是由几个乃至无穷多个简单的周期信号组成的,这些组成的周期信号的频率具有公约数,周期具有公共的周期。 因此,周期信号可以分解成多个乃至无穷多个谐波信号。反过来说,我们可以用一组谐波信号来合 成任意形状的周期信号。 例如对于如右图所示的方 波,其时域描述表达式为 000()()02()02x t x t nT T A t x t T A t =+????<?=???--<?? 其傅立叶三角函数展开式为
实验一MATLAB语言的基本使用方法 实验类别:基础性实验 实验目的: (1)了解MATLAB程序设计语言的基本方法,熟悉MATLAB软件运行环境。 (2)掌握创建、保存、打开m文件的方法,掌握设置文件路径的方法。 (3)掌握变量、函数等有关概念,具备初步的将一般数学问题转化为对应计算机模型并进行处理的能力。 (4)掌握二维平面图形的绘制方法,能够使用这些方法进行常用的数据可视化处理。 实验内容和步骤: 1、打开MATLAB,熟悉MATLAB环境。 2、在命令窗口中分别产生3*3全零矩阵,单位矩阵,全1矩阵。 3、学习m文件的建立、保存、打开、运行方法。 4、设有一模拟信号f(t)=1.5sin60πt,取?t=0.001,n=0,1,2,…,N-1进行抽样,得到 序列f(n),编写一个m文件sy1_1.m,分别用stem,plot,subplot等命令绘制32 点序列f(n)(N=32)的图形,给图形加入标注,图注,图例。 5、学习如何利用MATLAB帮助信息。 实验结果及分析: 1)全零矩阵 >> A=zeros(3,3) A = 0 0 0 0 0 0 0 0 0 2)单位矩阵 >> B=eye(3) B = 1 0 0 0 1 0 0 0 1 3)全1矩阵 >> C=ones(3) C = 1 1 1 1 1 1 1 1 1 4)sy1_1.m N=32; n=0:N-1; dt=0.001; t=n*dt; y=1.5*sin(60*pi*t); subplot(2,1,1), plot(t,y); xlabel('t'); ylabel('y=1.5*sin(60*pi*t)'); legend('正弦函数'); title('二维图形'); subplot(2,1,2), stem(t,y) xlabel('t'); ylabel('y=1.5*sin(60*pi*t)'); legend('序列函数'); title('条状图形'); 00.0050.010.0150.020.0250.030.035 t y = 1 . 5 * s i n ( 6 * p i * t ) 二维图形 00.0050.010.0150.020.0250.030.035 t y = 1 . 5 * s i n ( 6 * p i * t ) 条状图形
实验一:信号的时域分析 一、实验目的 1.观察常用信号的波形特点及产生方法 2.学会使用示波器对常用波形参数的测量 二、实验仪器 1.信号与系统试验箱一台(型号ZH5004) 2.40MHz双踪示波器一台 3.DDS信号源一台 三、实验原理 对于一个系统特性的研究,其中重要的一个方面是研究它的输入输出关系,即在一特定的输入信号下,系统对应的输出响应信号。因而对信号的研究是对系统研究的出发点,是对系统特性观察的基本手段与方法。在本实验中,将对常用信号和特性进行分析、研究。 信号可以表示为一个或多个变量的函数,在这里仅对一维信号进行研究,自变量为时间。常用信号有:指数信号、正弦信号、指数衰减正弦信号、复指数信号、Sa(t)信号、钟形信号、脉冲信号等。 1、信号:指数信号可表示为f(t)=Ke at。对于不同的a取值,其波形表现为不同的形式,如下图所示: 图1―1 指数信号 2、信号:其表达式为f(t)=Ksin(ωt+θ),其信号的参数:振幅K、角频率ω、与初始相位θ。其波形如下图所示:
图1-2 正弦信号 3、指数衰减正弦信号:其表达式为其波形如下图: 图1-3 指数衰减正弦信号 4、Sa(t)信号:其表达式为:。Sa(t)是一个偶函数,t= ±π,±2π,…,±nπ时,函数值为零。该函数在很多应用场合具有独特的运用。其信号如下图所示:
图1-4 Sa(t)信号 5、钟形信号(高斯函数):其表达式为:其信号如下图所示: 图1-5 钟形信号 6、脉冲信号:其表达式为f(t)=u(t)-u(t-T),其中u(t)为单位阶跃函数。其信号如下图所示: 7、方波信号:信号为周期为T,前T/2期间信号为正电平信号,后T/2期间信号为负电平信号,其信号如下图所示 U(t)
实验十三 信号分解及合成 一、 实验目的 1、 了解和熟悉波形分解与合成原理。 2、 了解和掌握用傅里叶级数进行谐波分析的方法。 二、 实验仪器 1、 双踪示波器 2、 数字万用表 3、 信号源及频率计模块S2 4、 数字信号处理模块S4 三、 实验原理 (一)信号的频谱与测量 信号的时域特性和频域特性是对信号的两种不同的描述方式。对于一个时域的周期信号 ()f t ,只要满足狄利克菜(Dirichlet)条件,就可以将其展开成三角形式或指数形式的傅里 叶级数。 例如,对于一个周期为T 的时域周期信号()f t ,可以用三角形式的傅里叶级数求出它的 各次分量,在区间11(,)t t T +内表示为 () 01 ()cos sin 41,3,5,7,n n n f t a a n t b n t A k Tk ω ∞ ==+Ω+Ω=??? ∑ ()01 ()cos sin n n n f t a a n t b n t ∞ ==+Ω+Ω∑ 即将信号分解成直流分量及许多余弦分量和正弦分量,研究其频谱分布情况。 图1 c a
信号的时域特性和频域特性 信号的时域特性与频域特性之间有着密切的内在联系,这种联系可以用图13-1来形象地表示。其中图(a)是信号在幅度—时间—频率三维坐标系统中的图形;图(b)是信号在幅度一时间坐标系统中的图形即波形图:把周期信号分解得到的各次谐波分量按频率的高低排列,就可以得到频谱图。反映各频率分量幅度的频谱称为振幅频谱。图(c)是信号在幅度—频率坐标系统中的图形即振幅频谱图。反映各分量相位的频谱称为相位频谱。在本实验中只研究信号振幅频谱。周期信号的振幅频谱有三个性质:离散性、谐波性、收敛性。测量时利用了这些性质。从振幅频谱图上,可以直观地看出各频率分量所占的比重。测量方法有同时分析法和顺序分析法。 同时分析法的基本工作原理是利用多个滤波器,把它们的中心频率分别调到被测信号的各个频率分量上。当被测信号同时加到所有滤波器上,中心频率与信号所包含的某次谐波分景频率-致的滤波器便有输出。在被测信号发生的实际时间内可以同时测得信号所包含的各频率分量。在本实验中采用同时分析法进行频谱分析,如图132所示。 (二)方波的分解 我们以下图的方波为例:占空比为50% 方波在一个周期内的解析式为:0()2 A t T f t T A t T <≤?? =? -<≤?? 故有 () 01 ()cos sin 41,3,5,7,n n n f t a a n t b n t A k Tk ω ∞ ==+Ω+Ω=??? ∑ 于是,所求级数 b
武汉工程大学 数字信号处理实验报告 姓名:周权 学号:1204140228 班级:通信工程02
一、实验设备 计算机,MATLAB语言环境。 二、实验基础理论 1.序列的相关概念 2.常见序列 3.序列的基本运算 4.离散傅里叶变换的相关概念 5.Z变换的相关概念 三、实验内容与步骤 1.离散时间信号(序列)的产生 利用MATLAB语言编程产生和绘制单位样值信号、单位阶跃序列、指数序列、正弦序列及随机离散信号的波形表示。 四实验目的 认识常用的各种信号,理解其数字表达式和波形表示,掌握在计算机中生成及绘制数字信号波形的方法,掌握序列的简单运算及计算机实现与作用,理解离散时间傅里叶变换,Z变换及它们的性质和信号的频域分
实验一离散时间信号(序列)的产生 代码一 单位样值 x=2; y=1; stem(x,y); title('单位样值 ') 单位阶跃序列 n0=0; n1=-10; n2=10; n=[n1:n2]; x=[(n-n0)>=0]; stem(n,x); xlabel('n'); ylabel('x{n}'); title('单位阶跃序列');
实指数序列 n=[0:10]; x=(0.5).^n; stem(n,x); xlabel('n'); ylabel('x{n}'); title('实指数序列');
正弦序列 n=[-100:100]; x=2*sin(0.05*pi*n); stem(n,x); xlabel('n'); ylabel('x{n}'); title('正弦序列');
随机序列 n=[1:10]; x=rand(1,10); subplot(221); stem(n,x); xlabel('n'); ylabel('x{n}'); title('随机序列');
深圳大学实验报告课程名称:信号与系统 实验项目名称:信号的分解与合成实验 学院:信息工程工程学院 专业: 电子信息工程 指导教师: 报告人:学号:班级: 实验时间: 实验报告提交时间: 教务处制
电位器W01、W02、W03可以将基波,三次谐波,五次谐波,七次谐波的幅度调节成1:1/3 : 1/5 : 1/7,通过导线将其连接至信号的合成的输入插座IN01、IN02、IN03、IN04J ,通过测试勾可以观察到合成后的波形。 2、验证三次谐波与基波之间的相位差是否为180,五次谐波与基波之间的相位差是否为0.可用李沙育图形法进行测量,其测量方法如下:用导线将函数发生器的方便输出端与带通滤波器输入端连接起来,即把方波信号分先后送入各带通滤波器,如图(1)所示. 具体方法:基波与各高次谐波相位比较(李沙育频率测试法) 把BFP-1ω处的基波送入示波器的X 轴,再分别把BFP-31ω、BFP-51ω处的高次谐波送入Y 轴,示波器采用X —Y 方式显示,观察李沙育图。 当基波与三次谐波相位差为0、90、180时,波形分别如图所示. 以上是三次谐波与基波产生的典型的李沙育图,通过图形上下端及两旁的波峰个数,确定频率比.
五、实验步骤与相应实验结果: 1、把电信号分解与合成模块插在主板上,用导线接通此模块“电源插入”和主板上的电源,并打开此模块的电源开关. 2、调节函数信号发生器,使其输出10KHz左右的方波,占空比为50%,峰峰值为6V左右,如图(2)所示。将其接至该实验模块的“输入端",用示波器观察各次谐波的输出即各次谐波,分别如图(3)、图(4)、图(5)、图(6)所示. 图(2)输出方波信号 图(3)基次谐波图(4)三次谐波 图(5)五次谐波图(6)七次谐波
东南大学电工电子实验中心 实验报告 课程名称:电子线路实践 第七次实验 实验名称:信号的产生、分解与合成 院(系):电子科学与工程学院专业: 姓名:姜勖学号:06A11315 实验室:104实验组别:27 同组人员:徐媛媛实验时间:年月日 评定成绩:审阅教师: 实验四信号的产生、分解与合成 一、实验内容及要求 设计并安装一个电路使之能够产生方波,并从方波中分离出主要谐波,再将这些谐波合成为原始信号或其他周期信号。 1.基本要求 (1)设计一个方波发生器,要求其频率为1kHz,幅度为5V; (2)设计合适的滤波器,从方波中提取出基波和3次谐波; (3)设计一个加法器电路,将基波和3次谐波信号按一定规律相加,将合成后的信号与原始信号比较,分析它们的区别及原因。 2.提高要求 设计5次谐波滤波器或设计移相电路,调整各次谐波的幅度和相位,将合成后的信号与原始信号比较,并与基本要求部分作对比,分析它们的区别及原因。 3.创新要求 用类似方式合成其他周期信号,如三角波、锯齿波等。 分析项目的功能与性能指标: 功能:通过振荡电路产生一个方波,并将其通过滤波得到1、3、5次谐波,最后通过加法电路合成新的波形。 性能指标: (1)方波:频率1KHz,幅度5V。 (2)滤波器:基础要求从方波中提取基波和三次谐波,提高要求提取五次谐波。 (3)移相电路:通过移相电路调节滤出来的1、3、5次谐波相位,使得其与原方波相位差近似为0。
(4)加法器电路:将基波和3次谐波和5次谐波信号按一定规律相加。 1、信号的产生 通过震荡电路产生1kHz ,幅度为5V 的方波信号。 2、滤波器的设计 根据方波的傅里叶展开式: 可知原信号分解只包含奇次谐波分量。因此设计不同中心频率的带通滤波器,可将各次谐波滤出。 3、相位校正电路 由于滤波器用到了对不同频率有不同响应的储能元件,对于滤除的波形会产生附加相位。若要让各次谐波叠加出原有信号,必须调节其相位使之同相。用全通滤波器可在不影响相对幅度的前提下改变相位。 4、加法电路 将滤除的基波、3次谐波、5次谐波相加,得到近似的方波信号。对于滤波器对不同频率分量不成比例的衰减,可在加法电路中选择合适的比例给予响应的补偿。 二、电路设计(预习要求) (1) 电路设计思想(请将基本要求、提高要求、创新要求分别表述): 1、信号发生电路: 利用运放和RC 回路构成振荡电路,通过分别调节正反向RC 回路的时间常数和运放同相输入端的参考电压来调节震荡电路的频率以及占空比。用一对稳压二极管限制输出电压幅度,并对稳压管导通压降进行一定的补偿。 2、有源带通滤波器: 根据实验要求,设计有源带通滤波器,将所需频率的信号以尽量小的衰减输出,同时对其它频率有非常大的衰减。因此需要增加滤波器的阶数。初步选择采用二阶有源带通滤波器,通过理论计算,调节其中一个电阻来改变中心频率。根据实际搭出的电路效果,可尝试使用四阶有源带通滤波器,以求获得更好的滤波效果。 3、相移电路: 由于滤波器难免对滤出的谐波分量产生附加相位,需要在选频电路之后加一全通网络校正相位,抵消相位差。移向电路有两种,分为正向移向和反向移向。 4、加法电路 将所得到的各次谐波分量叠加,得到近似的方波。同时,加法电路可对滤波对原信号分量的衰减进行补偿。 (2) 电路结构框图(请将基本要求、提高要求、创新要求分别画出): 基础要求:因基础要求与提高要求相比,除缺少5次滤波与移相电路外,其余部分均相同,其结构框图已包含在提高要求的框图中,故不单独列出。 提高要求: (3)电路原理图(各单元电路结构、工作原理、参数计算和元器件选择说明): 分工:徐媛媛(滤波电路的设计、搭建和调试);姜勖(方波产生、相移及加法电路设计搭建和调试) 方波振荡及鉴幅电路: 采用迟滞比较及RC 反馈回路以及比较器鉴幅电路,总电路图如下: 设从输出端的对输入端的负反馈电阻分别为1f R 和2f R ,则前部分方波的振荡周期为111222 ln(12)ln(12)f f R R T R C R C R R =+++,通过电位器分别调节1f R 和2f R 的阻值使方波的频率为1kHz ,占空比为50%。
实验一 MATLAB 仿真软件的基本操作命令和使用方法 实验容 1、帮助命令 使用 help 命令,查找 sqrt (开方)函数的使用方法; 2、MATLAB 命令窗口 (1)在MATLAB 命令窗口直接输入命令行计算3 1)5.0sin(21+=πy 的值; (2)求多项式 p(x) = x3 + 2x+ 4的根; 3、矩阵运算 (1)矩阵的乘法 已知 A=[1 2;3 4], B=[5 5;7 8],求 A^2*B
(2)矩阵的行列式 已知A=[1 2 3;4 5 6;7 8 9],求A (3)矩阵的转置及共轭转置 已知A=[1 2 3;4 5 6;7 8 9],求A' 已知B=[5+i,2-i,1;6*i,4,9-i], 求B.' , B' (4)特征值、特征向量、特征多项式 已知A=[1.2 3 5 0.9;5 1.7 5 6;3 9 0 1;1 2 3 4] ,求矩阵A的特征值、特征向量、特征多项式;
(5)使用冒号选出指定元素 已知:A=[1 2 3;4 5 6;7 8 9];求A 中第3 列前2 个元素;A 中所有列第2,3 行的元素; 4、Matlab 基本编程方法 (1)编写命令文件:计算1+2+…+n<2000 时的最大n 值;
(2)编写函数文件:分别用for 和while 循环结构编写程序,求 2 的0 到15 次幂的和。
5、MATLAB基本绘图命令 (1)绘制余弦曲线 y=cos(t),t∈[0,2π]
(2)在同一坐标系中绘制余弦曲线 y=cos(t-0.25)和正弦曲线 y=sin(t-0.5), t∈[0,2π] (3)绘制[0,4π]区间上的 x1=10sint 曲线,并要求: (a)线形为点划线、颜色为红色、数据点标记为加号; (b)坐标轴控制:显示围、刻度线、比例、网络线 (c)标注控制:坐标轴名称、标题、相应文本; >> clear;
实验二 信号分解与合成 --谢格斯 110701336 聂楚飞110701324 一、实验目的 1、观察电信号的分解。 2、掌握带通滤波器的有关特性测试方法。 3、观测基波和其谐波的合成. 二、实验内容 1、观察信号分解的过程及信号中所包含的各次谐波。 2、观察由各次谐波合成的信号。 三、预备知识 1、了解李沙育图相关知识. 2、课前务必认真阅读教材中周期信号傅里叶级数的分解以及如何将各次谐波进行叠加等相关内容. 四、实验仪器 1、信号与系统实验箱一台(主板)。 2、电信号分解与合成模块一块。 3、20M双踪示波器一台. 五、实验原理 任何电信号都是由各种不同频率、幅度和初相的正弦波迭加而成的。对周期信号由它的傅里叶级数展开式可知,各次谐波为基波频率的整数倍。而非周期信号包含了从零到无穷大的所有频率成份,每一频率成份的幅度均趋向无限小,但其相对大小是不同的. 通过一个选频网络可以将电信号中所包含的某一频率成份提取出来。本实验采用性能较佳的有源带通滤波器作为选频网络,因此对周期信号波形分解的实验方案如图2-3—1所示。 将被测方波信号加到分别调谐于其基波和各次奇谐波频率的一系列有源带通滤波器电路上。从每一有源带通滤波器的输出端可以用示波器观察到相应频率的正弦波。本实验所用的被测信号是Hz 531=ω左右的周期信号,而用作选频网络的五种有源带通滤波器的输出频 率分别是543215432ωωωωω、、、、 ,因而能从各有源带通滤波器的两端观察到基波和各次谐波.其中,在理想情况下,如方波的偶次谐波应该无输出信号,始终为零电平,而奇次谐波则具有很好的幅度收敛性,理想情况下奇次谐波中一、三、五、七、九次谐波的幅度比应为1:(1/3):(1/5):(1/7):(1/9)。但实际上因输入方波的占空比较难控制在50%,且方波可能有少量失真以及滤波器本身滤波特性的有限性都会使得偶次谐波分量不能达到理想零的情况。 六、实验步骤 1、把系统时域与频域分析模块插在主板上,用导线接通此模块“电源接入"和主板上的电源(看清标识,防止接错,带保护电路),并打开此模块的电源开关. 2、调节函数信号发生器,使其输出Hz 53左右(其中在Hz Hz 56~50之间进行选择,
学生实验报告 (理工类) 课程名称:信号与系统实验专业班级:电子信息(1)班学生学号:1005101058 学生姓名:严生生 所属院部:信息技术学院指导教师:杨婧 20 11 ——20 12 学年第 1 学期 金陵科技学院教务处制
实验报告书写要求 实验报告原则上要求学生手写,要求书写工整。若因课程特点需打印的,要遵照以下字体、字号、间距等的具体要求。纸张一律采用A4的纸张。 实验报告书写说明 实验报告中一至四项内容为必填项,包括实验目的和要求;实验仪器和设备;实验内容与过程;实验结果与分析。各院部可根据学科特点和实验具体要求增加项目。 填写注意事项 (1)细致观察,及时、准确、如实记录。 (2)准确说明,层次清晰。 (3)尽量采用专用术语来说明事物。 (4)外文、符号、公式要准确,应使用统一规定的名词和符号。 (5)应独立完成实验报告的书写,严禁抄袭、复印,一经发现,以零分论处。 实验报告批改说明 实验报告的批改要及时、认真、仔细,一律用红色笔批改。实验报告的批改成绩采用百分制,具体评分标准由各院部自行制定。 实验报告装订要求 实验批改完毕后,任课老师将每门课程的每个实验项目的实验报告以自然班为单位、按学号升序排列,装订成册,并附上一份该门课程的实验大纲。
实验项目名称:常用连续信号的表示实验学时: 1 同组学生姓名:实验地点: B402 实验日期:实验成绩: 批改教师:杨婧批改时间: 一、实验目的和要求 熟悉MATLAB软件,利用MATLAB软件,绘制出常用的连续时间信号。 二、实验仪器和设备 586以上计算机,装有MATLAB7.0软件。 三、实验过程 1,绘制正弦信号f(t)=Asin(ωt+ψ),其中A=1,ω=2π, ψ=π/6; 2,绘制指数信号f(t)=Ae^at,其中A=1,a=-0.4; 3,绘制矩形脉冲信号,脉冲宽度为2; 4,绘制三角波脉冲信号,脉冲宽度为4;斜度为0.5; 5,对上题三角波脉冲信号进行尺度变换,分别得出f(2t),f(2-2t); 6,绘制抽样函数Sa(t),t取值在-3π到+3π之间; 7,绘制周期矩形脉冲信号,参数自定; 8,绘制周期三角脉冲信号,参数自定; 1,打开MATLAB界面,建立新文件。 2,根据实验要求,编写程序。
信息科学与技术学院本科三年级 数字信号处理实验报告 2011 年12 月21日
实验一 序列的傅立叶变换 实验目的 进一步加深理解DFS,DFT 算法的原理;研究补零问题;快速傅立叶变换 (FFT )的应用。 实验步骤 1. 复习DFS 和DFT 的定义,性质和应用; 2. 熟悉MATLAB 语言的命令窗口、编程窗口和图形窗口的使用;利用提供的 程序例子编写实验用程序;按实验内容上机实验,并进行实验结果分析;写出完整的实验报告,并将程序附在后面。 实验内容 1. 周期方波序列的频谱试画出下面四种情况下的的幅度频谱,并分析补零后,对信号频谱的影响。 实验结果: 60 ,7)4(;60,5)3(; 40,5)2(;20,5)1()] (~[)(~,2,1,01 )1(,01,1)(~=========±±=???-+≤≤+-+≤≤=N L N L N L N L n x DFS k X m N m n L m N L m N n m N n x ) 52.0cos()48.0cos()(n n n x ππ+=
2. 有限长序列x(n)的DFT (1) 取x(n)(n=0:10)时,画出x(n)的频谱X(k) 的幅度; (2) 将(1)中的x(n)以补零的方式,使x(n)加长到(n:0~100)时,画出 x(n)的频谱X(k) 的幅度; (3) 取x(n)(n:0~100)时,画出x(n)的频谱X(k) 的幅度。利用FFT 进行谱分析 已知:模拟信号 以t=0.01n(n=0:N-1)进行采样,求N 点DFT 的幅值谱。 请分别画出N=45; N=50;N=55;N=60时的幅值曲线。 实验结果: ) 8cos(5)4sin(2)(t t t x ππ+=
实验二 波形合成与分解 1.实验目的 在理论学习的基础上,通过本实验熟悉信号的合成、分解原理,了解信号频谱的含义,加深对傅里叶变换性质和作用的理解。 2.实验原理 根据傅里叶分析的原理,任何周期信号都可以用一组三角函数)}cos();{sin(00t n t n ωω的组合表示,即: )2sin()2cos()sin()cos()(020201010t b t a t b t a a t x ωωωω++++= 即可以用一组正弦波和余弦波来合成任意形状的周期信号。 3.实验内容 (1) 方波的合成 图示方波是一个奇谐信号,由傅里叶级数可知,它是由无穷个奇次谐波分量 合成的,本实验用图形的方式来表示它的合成。方波信号可以分解为: ,9,7,5,3,1,1)2sin(2)(10=?=∑∞ =n n t nf A t x n ππ 用前5项谐波近似合成50Hz,幅值为3的方波,写出实验步骤。 a.只考察从 0=t s 到10=t s 这段时间内的信号。 b.画出基波分量)sin()(t t y =。 c.将三次谐波加到基波之上,并画出结果,并显示。 3/)*3sin()sin()(t t t y += d.再将一次、三次、五次、七次和九次谐波加在一起。 9/)*9sin(7/)*7sin(5/)*5sin(3/)*3sin()sin()(t t t t t t y ++++= e.合并从基波到十九次谐波的各奇次谐波分量。 f.将上述波形分别画在一幅图中,可以看出它们逼近方波的过程。注意“吉布斯现象”。周期信号傅里叶级数在信号的连续点收于该信号,在不连续点收敛于信号左右极限的平均值。如果我们用周期信号傅里叶级数的部分和来近似周期信号,在不连续点附近将会出现起伏和超量。在实际中,如果应用这种近似,就应该选择足够大的N ,以保证这些起伏拥有的能量可以忽略。 (2) 设计谐波合成三角波的实验,写出实验步骤,并完成实验。
实验三 常见信号的MATLAB 表示及运算 一、实验目的 1.熟悉常见信号的意义、特性及波形 2.学会使用MATLAB 表示信号的方法并绘制信号波形 3. 掌握使用MATLAB 进行信号基本运算的指令 4. 熟悉用MATLAB 实现卷积积分的方法 二、实验原理 根据MATLAB 的数值计算功能和符号运算功能,在MA TLAB 中,信号有两种表示方法,一种是用向量来表示,另一种则是用符号运算的方法。在采用适当的MA TLAB 语句表示出信号后,就可以利用MA TLAB 中的绘图命令绘制出直观的信号波形了。 1.连续时间信号 从严格意义上讲,MATLAB 并不能处理连续信号。在MATLAB 中,是用连续信号在等时间间隔点上的样值来近似表示的,当取样时间间隔足够小时,这些离散的样值就能较好地近似出连续信号。在MATLAB 中连续信号可用向量或符号运算功能来表示。 ⑴ 向量表示法 对于连续时间信号()f t ,可以用两个行向量f 和t 来表示,其中向量t 是用形如12::t t p t =的命令定义的时间范围向量,其中,1t 为信号起始时间,2t 为终止时间,p 为时间间隔。向量f 为连续信号()f t 在向量t 所定义的时间点上的样值。 ⑵ 符号运算表示法 如果一个信号或函数可以用符号表达式来表示,那么我们就可以用前面介绍的符号函数专用绘图命令ezplot()等函数来绘出信号的波形。 ⑶ 常见信号的MATLAB 表示 单位阶跃信号 单位阶跃信号的定义为:10 ()0 t u t t >?=? 方法一: 调用Heaviside(t)函数 首先定义函数Heaviside(t) 的m 函数文件,该文件名应与函数名同名即Heaviside.m 。 %定义函数文件,函数名为Heaviside,输入变量为x,输出变量为y function y= Heaviside(t) y=(t>0); %定义函数体,即函数所执行指令
3.(1)用双线性变换法设计一个Chebyshev型高通滤波器程序如下 Rp=1.2;Rs=20;T=0.001;fp=300;fs=200; wp=2*pi*fp*T;ws=2*pi*fs*T; wp1=(2/T)*tan(wp/2);ws1=(2/T)*tan(ws/2); [n,wn]=cheb1ord(wp1,ws1,Rp,Rs,'s'); [b,a]=cheby1(n,Rp,wn,'high','s'); [bz,az]=bilinear(b,a,1/T); [db,mag,pha,grd,w]=freqz_m(bz,az);plot(w/pi,db); axis([0,1,-30,2]); 3.(2) a用双线性变换法设计一个Butterworth型数字低通滤波器程序如下Rp=1;Rs=25;T=0.001;fp=300;fs=200; wp=2*pi*fp*T;ws=2*pi*fs*T; wp1=(2/T)*tan(wp/2);ws1=(2/T)*tan(ws/2); [n,wn]=buttord(wp1,ws1,Rp,Rs,'s'); [b,a]=butter(n,wn,'low','s'); [bz,az]=bilinear(b,a,1/T); [db,mag,pha,grd,w]=freqz_m(bz,az);plot(w/pi,db); axis([0,1,-30,2]); b用脉冲响应不变法设计一个Butterworth数字低通滤波器的程序如下:wp=400*pi;ws=600*pi;Rp=1;Rs=25; [n,wn]=buttord(wp,ws,Rp,Rs,'s') [b,a]=butter(n,wn,'s') [db,mag,pha,w]=freqs_m(b,a,500*2*pi);
数字信号处理实验 利用FFT对信号进行频谱分析
一 实验目的 学习用FFT FFT 。 二 实验原理 用FFT 对信号作频分析是学习数字信号处理的重要内容,经常需要进行分析的信号是模拟信号的时域离散信号。对信号进行谱分析的重要问题是频谱分辨率D 和分析误差。频谱分辨率直接和FFT 的变换区间N 有关,因为FFT 能够实现的频率分辨率是2π/N ,因此要求2π/N 小于等于D 。可以根据此式选择FFT 的变换区间N 。误差主要来自于用FFT 作频谱分析时,得到的是离散谱,而信号(周期信号除外)是连续谱,只有当N 较大时,离散谱的包络才能逼近连续谱因此N 要适当选择大一些。 三 实验内容 1.模拟信号)8cos(5)4sin(2)(t t t x ππ+=,以)1:0(01.0-==N n n t 进行采样,求: (1)N =40点FFT 的幅度频谱,从图中能否观察出信号的2个频谱分量? (2)提高采样点数,如N =128,256,512,再求该信号的幅度频谱,此时幅度频谱发生了什么变化?信号的2个模拟频率和数字频率各为多少?FFT 频谱分析结果与理论上是否一致? 实验代码: clc;clear all; N=40; % N=128;%%%%%%对N 的值进行改变 % N=256; % N=512; n=0:N-1; t=0.01*n; x=2*sin(4*pi*t)+5*cos(8*pi*t);
x1=x(1:N);X1=fft(x1,2048); figure, subplot(211),plot(0:N-1,x1);xlabel('n');ylabel('x(n)');title('时域波形');grid; subplot(212),plot(abs(X1));xlabel('k');ylabel('|X(k)|');title('幅频特性');grid; set(gcf,'color','w'); N=40 N=128