实验二连续时间信号卷积及MATLAB实现

实验二连续时间信号卷积及MATLAB实现

一、实验目的：熟悉使用MATLAB软件来分析连续时间信号的卷积积分运算并用图形可视化相关结果.

二、实验内容：

1. 卷积积分

卷积积分在信号与线形系统分析中具有非常重要的意义，市信号与系统分析的基本方法

之一.

用MATLAB实现连续信号f i(t)与f2(t)卷积的过程如下：

(1) 将连续时间信号f i(t)与f2(t)以时间间隔△进行取样，得到离散序列f i(k △)和f2(k △)；

(2) 构造与f i(k △)和f 2(k △)相对应的时间向量k i和k;

(3) 调用conv()函数计算卷积积分f(t)的近似向量f(n △) 2

(4) 构造f(n △)对应的时间向量k.

2. 已知两连续信号如下图所示，用MATLAB f(t)=f1(t)*f2(t), 并绘出f(t)的时域波

形图。(设定取样时间间隔为p)

当p=0.5时，其参考程序为：

p=0.5;

k1=0:p:2;

f1=0.5*k1;

k2=k1;

f2=f1;

[f,k]=sconv(f1,f2,k1,k2,p)

所以其图形如下：

当P=0.1时，其参考程序为：

p=0.1；

k1=0:p:2;

f1=0.5*k1;

k2=k1;

f2=f1;

[f,k]=sconv(f1,f2,k1,k2,p)

所以其图形如下：

同理两信号如下所示，用MATLAB求f(t)=f1(t)*f2(t). 当P=0.1时，参考程序为：

p=0.1；

k1=0:p:2;

f1=2*ones(1,length(k1));

k2=0:p:4;

f2=ones(1,length(k2));

[f,k]=sconv(f1,f2,k1,k2,p)

当P=0.5时，参考程序为:

p=0.5;

k1=0:p:2;

f1=2*ones(1,length(k1));

k2=0:p:4;

f2=ones(1,length(k2));

[f,k]=sconv(f1,f2,k1,k2,p)

信号与系统matlab实验及答案

产生离散衰减正弦序列()π0.8sin 4n x n n ?? = ??? , 010n ≤≤，并画出其波形图。 n=0:10; x=sin(pi/4*n).*0.8.^n; stem(n,x);xlabel( 'n' );ylabel( 'x(n)' ); 用MATLAB 生成信号()0sinc at t -, a 和0t 都是实数,410t -<<，画波形图。观察并分析a 和0t 的变化对波形的影响。 t=linspace(-4,7); a=1;

t0=2; y=sinc(a*t-t0); plot(t,y); t=linspace(-4,7); a=2; t0=2; y=sinc(a*t-t0); plot(t,y);

t=linspace(-4,7); a=1; t0=2; y=sinc(a*t-t0); plot(t,y);

三组对比可得a 越大最大值越小，t0越大图像对称轴越往右移 某频率为f 的正弦波可表示为()()cos 2πa x t ft =，对其进行等间隔抽样，得到的离散样值序列可表示为()()a t nT x n x t ==，其中T 称为抽样间隔，代表相邻样值间的时间间隔，1 s f T = 表示抽样频率，即单位时间内抽取样值的个数。抽样频率取40 Hz s f =，信号频率f 分别取5Hz, 10Hz, 20Hz 和30Hz 。请在同一张图中同时画出连续信号()a x t t 和序列()x n nT 的波形图，并观察和对比分析样值序列的变化。可能用到的函数为plot, stem, hold on 。 fs = 40; t = 0 : 1/fs : 1 ; % ?μ?ê·?±e?a5Hz,10Hz,20Hz,30Hz f1=5; xa = cos(2*pi*f1*t) ; subplot(1, 2, 1) ;

移动通信实验线性分组码卷积码实验

实验二抗衰落技术实验（4学时） 1.线性分组码实验 2.卷积码实验 姓名： 学号： 班级： 日期： 成绩：

1、线性分组码实验 一、实验目的 了解线性分组码在通信系统中的意义。 掌握汉明码编译码及其检错纠错原理，理解编码码距的意义。二、实验模块 主控单元模块 2号数据终端模块 4号信道编码模块 5号信道译码模块 示波器 三、实验原理

汉明码编译码实验框图 2、实验框图说明 汉明码编码过程：数字终端的信号经过串并变换后，数据进行了分组，分组后的数据再经过汉明码编码，数据由4bit变为7bit。 注：为方便对编码前后的数据进行对比观测，本实验中加入了帧头指示信号。帧头指示信号仅用于线性分组码编码时将输入信号的比特流进行分组，其上跳沿指示了分组的起始位置。 四、实验步骤 （注：实验过程中，凡是涉及到测试连线改变或者模块及仪器仪表的更换时，都需先停止运行仿真，待连线调整完后，再开启仿真进行后续调节测试。） 任务一汉明码编码规则验证 概述：本项目通过改变输入数字信号的码型，观测延时输出，编码输出及译码输出，验证汉明码编译码规则。 1、登录e-Labsim仿真系统，创建实验文件，选择实验所需模块和示波器。 2、按表格所示进行连线。 3、调用示波器观测2号模块的DoutMUX和4号模块的编码输出TH4编码数据，

6、此时系统初始状态为：2号模块提供32K编码输入数据，4号模块进行汉明码编码，无差错插入模式，5号模块进行汉明码译码。 7、实验操作及波形观测。 0000 0001 0010

0100 0101

0111 1000

matlab信号与系统实验报告

实验一 基本信号的产生与运算 一、 实验目的 学习使用MATLAB 产生基本信号、绘制信号波形、实现信号的基本运算。 二、 实验原理 MATLAB 提供了许多函数用于产生常用的基本信号：如阶跃信号、脉冲信号、指数信号、正弦信号和周期方波等等。这些信号是信号处理的基础。 1、 利用MATLAB 产生下列连续信号并作图。 （1）51),1(2)(<<---=t t u t x （2）300),3 2sin()(3.0<<=-t t e t x t （3）1.01.0,3000cos 100cos )(<<-+=t t t t x （4）2000),8.0cos()1.0cos()(<<=t t t t x ππ 答：（1）、 >> t=-1:0.02:5; >> x=(t>1); >> plot(t,-2*x); >> axis([-1,5,-3,1]); >> title('杨婕婕 朱艺星'); >> xlabel('x(t)=-2u(t-1)'); （2）、 >> t=0:0.02:30; >> x=exp(-0.3*t).*sin(2/3*t); >> plot(t,x);

>> title('杨婕婕朱艺星'); >> xlabel('x(t)=exp(-0.3*t).*sin(2/3*t)'); 因为原函数在t=15后x(t)取值接近于零，所以将横坐标改成0到15，看得更清晰 axis([0,15,-0.2,0.6]); （3）>> t=-0.1:0.01:0.1;x=cos(100*t)+cos(3000*t);plot(t,x); >> title('杨婕婕朱艺星'); >>xlabel('x=cos(100*t)+cos(3000*t)'); 因为t的间隔取太大，以至于函数不够准确，缩小t的间隔： t=-0.1:0.002:0.2;x=cos(100*t)+cos(3000*t); plot(t,x);title('杨婕婕') >> t=-0.1:0.0001:0.1; x=cos(100*t)+cos(3000*t); >> plot(t,x);title('杨婕婕朱艺星'); >> xlabel('x=cos(100*t)+cos(3000*t)');

西电数字信号处理上机实验报告

数字信号处理上机实验报告 14020710021 张吉凯 第一次上机 实验一： 设给定模拟信号()1000t a x t e -=，t 的单位是ms 。 (1) 利用MATLAB 绘制出其时域波形和频谱图（傅里叶变换），估计其等效带宽（忽略谱分量降低到峰值的3%以下的频谱）。 (2) 用两个不同的采样频率对给定的()a x t 进行采样。 ○1()()15000s a f x t x n =以样本秒采样得到。 ()()11j x n X e ω画出及其频谱。 ○2()()11000s a f x t x n =以样本秒采样得到。 ()() 11j x n X e ω画出及其频谱。 比较两种采样率下的信号频谱，并解释。 （1）MATLAB 程序： N=10; Fs=5; T s=1/Fs; n=[-N:T s:N]; xn=exp(-abs(n)); w=-4*pi:0.01:4*pi; X=xn*exp(-j*(n'*w)); subplot(211) plot(n,xn); title('x_a(t)时域波形'); xlabel('t/ms');ylabel('x_a(t)'); axis([-10, 10, 0, 1]); subplot(212); plot(w/pi,abs(X)); title('x_a(t)频谱图'); xlabel('\omega/\pi');ylabel('X_a(e^(j\omega))');

ind = find(X >=0.03*max(X))*0.01; eband = (max(ind) -min(ind)); fprintf('等效带宽为%fKHZ\n',eband); 运行结果： 等效带宽为12.110000KHZ

卷积信号实验报告

信号与系统上机实验报告一连续时间系统卷积的数值计算 140224 班张鑫学号 14071002 一、实验原理 计算两个函数的卷积 卷积积分的数值运算实际上可以用信号的分段求和来实现，即： 如果我们只求当 t = n? t1 是r ( t )的值，则由上式可以得到： ?t足够小时，r(t2)就是e(t)和f(t)卷积积分的数值近似值由上面的公式可 当1 以得到卷积数值计算的方法如下： (1)将信号取值离散化，即以为周期，对信号取值，得到一系列宽度间隔为 的矩形脉冲原信号的离散取值点，用所得离散取值点矩形脉冲来表示原来的连续时间信号； (2)将进行卷积的两个信号序列之一反转，与另一信号相乘，并求积分，所得为t=0时的卷积积分的值。以为单位左右移动反转的信号，与另一信号相乘求积 分，求的t<0和t>0时卷积积分的值； (3)将所得卷积积分值与对应的t标在图上，连成一条光滑的曲线，即为所求卷积积分的曲线。 1

信号与系统上机实验报告一二、处理流程图 三、C程序代码 #include"stdafx.h" #include"stdio.h" //#include "stdilb.h" float u(float t) { while (t>= 0) return(1); while (t<0) return(0); } float f1(float t) { return(u(t+2)-u(t-2)); } float f2(float t) { return(t*(u(t)-u(t-2))+(4-t)*(u(t-2)-u(t-4))); } int_tmain(int argc, _TCHAR* argv[]) {

信号(MATLAB)实验指导书

《信号与系统》实验指导书 张建奇骆崇编写 浙江工业大学之江学院信息工程分院 2012年2月

目录 实验一MATLAB的基本使用 (1) 一、实验目的 (1) 二、实验原理 (1) 三、实验内容与要求 (8) 四、实验报告 (9) 实验二时域波形的MATLAB实现 (10) 一、实验目的 (10) 二、预习要求 (10) 三、实验原理 (10) 四、实验内容与要求 (18) 五、实验报告 (19) 实验三用MATLAB对系统时域分析 (20) 一、实验目的 (20) 二、预习要求 (20) 三、实验原理 (20) 四、实验内容与要求 (29)

实验一MATLAB的基本使用 一、实验目的 1、了解和掌握MATLAB的基本操作 2、了解MATLAB的库函数 3、会用MATLAB进行简单的操作。 二、实验原理 1、界面操作 MATLAB是“MATrix LABoratory”的缩写（矩阵实验室），它是由美国Mathworks公司于1984年正式推出的一种科学计算软件，由于其强大的功能，在欧美的一些大学里MATLAB已经成为许多诸如数字信号处理、自动控制理论等高级教程的主要工具软件，同时也成为理工科学生，必须掌握的一项基本技能。 当需要运行程序时，只需选择桌面上（或开始）中的MATLAB6.5应用程序图标即可 通常情况下，MATLAB的工作环境主要由一下几个窗口组成： 命令窗口（Command Window）

工作区间浏览器（Workspace） 历史命令窗口（Command History） 图形窗口(Figure) 文本编辑窗口(Editor) 当前路径窗口（Current Directory） MATLAB的命令窗与命令操作 当用户使用命令窗口进行工作时，在命令窗口中可以直接输入相应的命令，系统将自动显示信息。 例如在命令输入提示符“>>”后输入指令： >>t=[1,2,3;4,5,6;7,8,9]; 按回车键(Enter)后，系统即可完成对变量t的赋值。 MATALB提供了非常方便的在线帮助命令(help)，它可提供各个函数的用法指南，包括格式、参数说明、注意事项及相关函数等内容。 2、图形窗 MATLAB图形窗(Figure)主要用于显示用户所绘制的图形。 通常，只要执行了任意一种绘图命令，图形窗就会自动产生。

信号与系统MATLAB实验报告

《信号与系统》MATLAB实验报告 院系：专业： 年级：班号： 姓名：学号： 实验时间： 实验地点：

实验一 连续时间信号的表示及可视化 实验题目： )()(t t f δ=；)()(t t f ε=；at e t f =)(（分别取00<>a a 及）； )()(t R t f =；)()(t Sa t f ω=；)2()(ft Sin t f π=（分别画出不同周期个数 的波形）。 解题分析： 以上各类连续函数，先运用t = t1: p:t2的命令定义时间范围向量，然后调用对应的函数，建立f 与t 的关系，最后调用plot （）函数绘制图像，并用axis （）函数限制其坐标范围。 实验程序： （1） )()(t t f δ= t=-1:0.01:3 %设定时间变量t 的范围及步长 f=dirac(t) %调用冲激函数dirac （） plot(t,f) %用plot 函数绘制连续函数 axis([-1,3,-0.5,1.5]) %用axis 函数规定横纵坐标的范围 （2） )()(t t f ε= t=-1:0.01:3 %设定时间变量t 的范围及步长 f=heaviside(t) %调用阶跃函数heaviside （） plot(t,f) %用plot 函数绘制连续函数 title('f(t)=heaviside(t)') %用title 函数设置图形的名称 axis([-1,3,-0.5,1.5]) %用axis 函数规定横纵坐标的范围 （3） at e t f =)( a=1时： t=-5:0.01:5 %设定时间变量t 的范围及步长 f=exp(t) %调用指数函数exp （）

卷积编码实验报告

实验名称：___ 卷积编码_______ 1、使用MATLAB进行卷积编码的代码编写、运行、仿真等操作； 2、熟练掌握MATLAB软件语句； 3、理解并掌握卷积编码的原理知识。 二、实验原理 卷积码是由Elias于1955 年提出的，是一种非分组码，通常它更适用于前向纠错法，因为其性能对于许多实际情况常优于分组码，而且设备较简单。 卷积码的结构与分组码的结构有很大的不同。具体地说，卷积码并不是将信息序列分成不同的分组后进行编码，而是将连续的信息比特序列映射为连续的编码器输出符号。卷积码在编码过程中，将一个码组中r 个监督码与信息码元的相关性从本码组扩展到以前若干段时刻的码组，在译码时不仅从此时刻收到的码组中提取译码信息，而且还可从与监督码相关的各码组中提取有用的译码信息。这种映射是高度结构化的，使得卷积码的译码方法与分组译码所采用的方法完全不同。可以验证的是在同样复杂度情况下，卷积码的编码增益要大于分组码的编码增益。对于某个

特定的应用，采用分组码还是卷积码哪一种更好则取决于这一应用的具体情况和进行比较时可用的技术。 （一）卷积编码的图形表示 卷积码的编码器是由一个有k 个输人位，n 个输出位，且有m 个移位寄存器构成的有限状态的有记忆系统，其原理如图1所示。 图1 卷积码编码器的原理图 描述这类时序网络的方法很多，它大致可分为两大类型：解析表示法与图形表示法。在解析法中又可分为离散卷积法、生成矩阵法、码多项式法等；在图形表示法中也可分为状态图法、树图法和网络图法等。 图2给出的是一个生成编码速率为1／2 卷积码的移位寄存器电路。输人比特在时钟触发下从左边移人到电路中，每输入一位，分别去两个模2加法器的输出值并复用就得到编码器的输出。对这一编码，每输入一比特就产生两个输出符号，故编码效率为

实验1 用MATLAB进行信号频谱分析(推荐文档)

实验1 用MATLAB 进行信号频谱分析 一、实验目的 ㈠ 初步掌握MATLAB 产生常用离散时间信号的编程方法。 ㈡ 学习编写简单的FFT 算法程序，对离散信号进行幅频谱分析。 ㈢ 观察离散时间信号频谱的特点。 二、实验原理 ㈠ 常用的离散时间信号 在 MATLAB 语言主要是研究离散信号的。常用的离散信号有： 1．单位取样序列 ???≠==000 1)(n n n δ 2．单位阶跃序列 ?? ?<≥=0 01 )(n n n u 3．实指数序列 R a n a n x n ∈?=;)( 4．复指数序列 n e n x n j ?=+)(0)(ωσ 5．正(余)弦序列 )c o s ()(0θω+=n n x n ? 6．周期序列 n N n x n x ?+=)()( ㈡ 离散信号的产生 离散信号的图形显示使用stem 指令。 在 MATLAB 中的信号处理工具箱中,主要提供的信号是离散信号。

由于MATLAB 对下标的约定为从1开始递增，例如x=[5,4,3,2,1,0],表示x(1)=5,x(2)=4,X(3)=3… 因此要表示一个下标不由1开始的数组x(n)，一般应采用两个矢量，如 n=[-3，-2，-1，0，l ，2，3，4，5]； x=[1，-l ，3，2，0，4，5，2，1]； 这表示了一个含9个采样点的矢量：X(n)={x(-3)，x(-2)，x(-1)，x(0)，x(1)，x(2)，x(3)，x(4)，x(5)}。 1．单位取样序列 ?? ?≠==δ0 001)(n n n 这一函数实现的方法有二： 方法一：可利用MATLAB 的zeros 函数。 x=zeros(1，N)； %建立一个一行N 列的全零数组 x(1)=1； %对X （1）赋1 方法二：可借助于关系操作符实现 n=1:N; x=[n==1]; %n 等于1时逻辑关系式结果为真，x=1；n 不等于1时为假，x=0 如要产生 ?????≤<<=≤≤=-δ2 0210 100)(10)(n n n n n n n n n n n n 则可采用MATLAB 实现： n=n1:n2； x=[(n-n0)==0]；%n=n0时逻辑关系式结果为真，x=1；n ≠n0时为

(完整word版)信号与系统matlab实验

习题三 绘制典型信号及其频谱图 1.更改参数，调试程序，绘制单边指数信号的波形图和频谱图。观察参数a对信号波形 及其频谱的影响。 程序代码： close all; E=1;a=1; t=0:0.01:4; w=-30:0.01:30; f=E*exp(-a*t); F=1./(a+j*w); plot(t,f);xlabel('t');ylabel('f(t)'); figure; plot(w,abs(F));xlabel('\omega');ylabel('|F(\omega)|';

E=1，a=1,波形图频谱图更改参数E=2，a=1;

更改参数a，对信号波形及其频谱的影响。（保持E=2）上图为a=1图像 a=2时

a=4时 随着a的增大，f(t)曲线变得越来越陡，更快的逼近0，而对于频谱图，随着a增大，图像渐渐向两边张开，峰值减小，陡度减小，图像整体变得更加平缓。 2.矩形脉冲信号 程序代码： close all; E=1;tao=1; t=-4:0.1:4; w=-30:0.1:30;

f=E*(t>-tao/2&tao/2)+0*(t<=-tao/2&t>=tao/2); F=(2*E./w).*sin(w*tao/2); plot(t,f);xlabel('t');ylabel('f(t)'); figure; plot(w,abs(F));xlabel('\omega');ylabel('|F(\omega)|') ; figure; plot(w,20*log10(abs(F))); xlabel('\omega');ylabel('|F(\omega)| in dB'); figure; plot(w,angle(F));xlabel('\omega');ylabel('\phi(\omega )');

卷积码实验报告

苏州科技大学天平学院电子与信息工程学院 信道编码课程设计报告 课设名称卷积码编译及译码仿真 学生姓名圣鑫 学号1430119232 同组人周妍智 专业班级通信1422 指导教师潘欣欲 一、实验名称 基于MAATLAB的卷积码编码及译码仿真 二、实验目的 卷积码就是一种性能优越的信道编码。它的编码器与译码器都比较容易实现,同时它具有较强的纠错能力。随着纠错编码理论研究的不断深入,卷积码的实际应用越来越广泛。本实验简明地介绍了卷积码的编码原理与Viterbi译码原理。并在SIMULINK模块设计中,完成了对卷积码的编码与译码以及误比特统计整个过程的模块仿真。最后,通过在仿真过程中分别改变卷积码的重要参数来加深理解卷积码的这些参数对卷积码的误码性能的影响。经过仿真与实测,并对测试结果作了分析。 三、实验原理

1、卷积码编码原理 卷积码就是一种性能优越的信道编码,它的编码器与解码器都比较易于实现,同时还具有较强的纠错能力,这使得它的使用越来越广泛。卷积码一般表示为(n,k,K)的形式,即将 k个信息比特编码为 n 个比特的码组,K 为编码约束长度,说明编码过程中相互约束的码段个数。卷积码编码后的 n 各码元不经与当前组的 k 个信息比特有关,还与前 K-1 个输入组的信息比特有关。编码过程中相互关联的码元有 K*n 个。R=k/n 就是编码效率。编码效率与约束长度就是衡量卷积码的两个重要参数。典型的卷积码一般选 n,k 较小,K 值可取较大(>10),但以获得简单而高性能的卷积码。 卷积码的编码描述方式有很多种:冲激响应描述法、生成矩阵描述法、多项式乘积描述法、状态图描述,树图描述,网格图描述等。 2、卷积码Viterbi译码原理 卷积码概率译码的基本思路就是:以接收码流为基础,逐个计算它与其她所 有可能出现的、连续的网格图路径的距离,选出其中可能性最大的一条作为译码估值输出。概率最大在大多数场合可解释为距离最小,这种最小距离译码体现的正就是最大似然的准则。卷积码的最大似然译码与分组码的最大似然译码在原理上就是一样的,但实现方法上略有不同。主要区别在于:分组码就是孤立地求解单个码组的相似度,而卷积码就是求码字序列之间的相似度。基于网格图搜索的译码就是实现最大似然判决的重要方法与途径。用格图描述时,由于路径的汇聚消除了树状图中的多余度,译码过程中只需考虑整个路径集合中那些使似然函数最大的路径。如果在某一点上发现某条路径已不可能获得最大对数似然函数,就放弃这条路径,然后在剩下的“幸存”路径中重新选择路径。这样一直进行到最后第 L 级(L 为发送序列的长度)。由于这种方法较早地丢弃了那些不可能的路径,从而减轻了译码的工作量,Viterbi 译码正就是基于这种想法。对于(n, k, K )卷积码,其网格图中共 2kL 种状态。由网格图的前 K-1 条连续支路构成的路径互不相交,即最初 2k_1 条路径各不相同,当接收到第 K 条支路时,每条路径都有 2 条支路延伸到第 K 级上,而第 K 级上的每两条支路又都汇聚在一个节点上。在Viterbi译码算法中,把汇聚在每个节点上的两条路径的对数似然函数累加

matlab信号与系统实验报告

实验一 基本信号的产生与运算 一、 实验目的 学习使用MATLAB 产生基本信号、绘制信号波形、实现信号的基本运算。 二、 实验原理 MATLAB 提供了许多函数用于产生常用的基本信号：如阶跃信号、脉冲信号、指数信号、正弦信号和周期方波等等。这些信号是信号处理的基础。 1、 利用MATLAB 产生下列连续信号并作图。 （1）51),1(2)(<<---=t t u t x （2）300),3 2 sin()(3.0<<=-t t e t x t （3）1.01.0,3000cos 100cos )(<<-+=t t t t x （4）2000),8.0cos()1.0cos()(<<=t t t t x ππ 答：（1）、 >> t=-1:0.02:5; >> x=(t>1); >> plot(t,-2*x); >> axis([-1,5,-3,1]); >> title('杨婕婕 朱艺星'); >> xlabel('x(t)=-2u(t-1)');

（2）、 >> t=0:0.02:30; >> x=exp(-0.3*t).*sin(2/3*t); >> plot(t,x); >> title('杨婕婕朱艺星'); >> xlabel('x(t)=exp(-0.3*t).*sin(2/3*t)');

因为原函数在t=15后x(t)取值接近于零，所以将横坐标改成0到15，看得更清晰 axis([0,15,-0.2,0.6]);

（3）>> t=-0.1:0.01:0.1;x=cos(100*t)+cos(3000*t);plot(t,x); >> title('杨婕婕朱艺星'); >>xlabel('x=cos(100*t)+cos(3000*t)'); 因为t的间隔取太大，以至于函数不够准确，缩小t的间隔：t=-0.1:0.002:0.2;x=cos(100*t)+cos(3000*t); plot(t,x);title('杨婕婕')

信 卷积实验报告

信号与系统实验报告学院：电子信息与电气工程学院 班级: 13级电信<1>班 学号: 20131060104 姓名:李重阳

实验三 信号卷积实验 一、实验目的 1、理解卷积的概念及物理意义； 2、通过实验的方法加深对卷积运算的图解方法及结果的理解。 二、实验原理说明 卷积积分的物理意义是将信号分解为冲激信号之和，借助系统的冲激响应，求解系统对任意激励信号的零状态响应。设系统的激励信号为x （t ），冲激响应为h （t ），则系统的零状态响应为()()()*y t x t h t ==()()x t h t d ττ∞-∞-?。 1、两个矩形脉冲信号的卷积过程 两信号x （t ）与h （t ）都为矩形脉冲信号，如图3-1所示。下面由图解的方法（图3-1）给出两个信号的卷积过程和结果，以便与实验结果进行比较。 图3-1 两矩形脉冲的卷积积分的运算过程与结果 2、矩形脉冲信号与锯齿波信号的卷积 信号f1（t ）为矩形脉冲信号， f2（t ）为锯齿波信号，如图3-2所示。根据卷积积分的运算方法得到f1（t ）和f2（t ）的卷积积分结果f （t ），如图3-2（c ）所示。 图3-2 矩形脉冲信号与锯齿脉冲信号的卷积积分的结果 3、本实验进行的卷积运算的实现方法 在本实验装置中采用了DSP 数字信号处理芯片，因此在处理模拟信号的卷积积分运算时，是先通过A/D 转换器把模拟信号转换为数字信号，利用所编写的相应程序控制DSP 芯片实现数字信号的卷积运算，再把运算结果通过D/A 转换为模拟信号输出。结果与模拟信号的直接运算结果是一致的。数字信号处理系统逐步和完全取代模拟信号处理系统是科学技术发展的必然趋势。图3-3为信号卷积的流程图。 图3-3 信号卷积的流程图 三、实验内容 1、检测矩形脉冲信号的自卷积结果。 用双踪示波器同时观察输入信号和卷积后的输出信号，把输入信号的幅度峰峰值调节为4V ，再调节输入信号的频率或占空比使输入信号的时间宽度满足表中的要求，观察输出信号有何变化，判断卷积的结果是否正确，并记录表3-1。 实验步骤如下： ①将函数发生器的SW702置于“方波”上。 ②连接函数发生器H701与数字滤波器的PB01，在TPB01上可观察到输入波形。将示波器接在TPB01上观测输入波形，并调节函数发生器模块上的频率旋钮与幅度旋钮，使信号频率为1KHz ，幅度为4V 。（注意：输入波形的频率幅度要在H701与PB01连接后，在TPB01上测试。） ③将红色拨动开关SWB01调整为“0001”。 ④按下复位键S1。 ⑤将示波器的CH1接于TP901；CH2接于TP903。可分别观察到输入信号的波形与卷积后的输出信号的波形。 表3-1 输入信号卷积后的输出信号

信号与系统MATLAB实验

2016-2017学年第一学期 信号与系统实验报告 班级： 姓名： 学号： 成绩： 指导教师：

实验一常见信号的MATLAB 表示及运算 一．实验目的 1．熟悉常见信号的意义、特性及波形 2．学会使用MATLAB 表示信号的方法并绘制信号波形 3. 掌握使用MATLAB 进行信号基本运算的指令 4. 熟悉用MATLAB 实现卷积积分的方法 二．实验原理 信号一般是随时间而变化的某些物理量。按照自变量的取值是否连续，信号分为连续时间信号和离散时间信号，一般用()f t 和()f k 来表示。若对信号进行时域分析，就需要绘制其波形，如果信号比较复杂，则手工绘制波形就变得很困难，且难以精确。MATLAB 强大的图形处理功能及符号运算功能，为实现信号的可视化及其时域分析提供了强有力的工具。 根据MATLAB 的数值计算功能和符号运算功能，在MATLAB 中，信号有两种表示方法，一种是用向量来表示，另一种则是用符号运算的方法。在采用适当的MATLAB 语句表示出信号后，就可以利用MATLAB 中的绘图命令绘制出直观的信号波形了。下面分别介绍连续时间信号和离散时间信号的MATLAB 表示及其波形绘制方法。 1.连续时间信号 所谓连续时间信号，是指其自变量的取值是连续的，并且除了若干不连续的点外，对于一切自变量的取值，信号都有确定的值与之对应。从严格意义上讲，MATLAB 并不能处理连续信号。在MATLAB 中，是用连续信号在等时间间隔点上的样值来近似表示的，当取样时间间隔足够小时，这些离散的样值就能较好地近似出连续信号。在MATLAB 中连续信号可用向量或符号运算功能来表示。 ⑴ 向量表示法 对于连续时间信号()f t ，可以用两个行向量f 和t 来表示，其中向量t 是用形如12::t t p t 的命令定义的时间范围向量，其中，1t 为信号起始时间，2t 为终止时间，p 为时间间隔。向量f 为连续信号()f t 在向量t 所定义的时间点上的样值。 说明：plot 是常用的绘制连续信号波形的函数。 严格说来，MATLAB 不能表示连续信号，所以，在用plot()命令绘制波形时，要对自变量t 进行取值，MATLAB 会分别计算对应点上的函数值，然后将各个数据点通过折线连接起来绘制图形，从而形成连续的曲线。因此，绘制的只是近似波形，而且，其精度取决于t 的取样间隔。t 的取样间隔越小，即点与点之间的距离越小，则近似程度越好，曲线越光滑。例如：图1-1是在取样间隔为p=0.5时绘制的波形，而图1-2是在取样间隔p=0.1时绘制的波形，两相对照，可以看出图1-2要比图1-1光滑得多。

实验十二 信号卷积实验报告有数据)

实验十二信号卷积实验 一、实验目的 1、理解卷积的概念及物理意义。 2、通过实验的方法加深对卷积运算的图解方法及结果的理解。 二、实验仪器 1、双踪示波器 1台 2、信号源及频率计模块S2 1块 3、数字信号处理模块S4 1块 三、实验原理 卷积积分的物理意义是将信号分解为冲激信号之和，借助系统的冲激响应，求解系统对任意激励信号的零状态响应。设系统的激励信号为)t(x，冲激响应为)t(h，则系统的零状态响应为： 对于任意两个信号)t(f 1和)t( f 2 ，两者做卷积运算定义为： ?∞∞-τ τ - =d) t( f)t(f )t(f 2 1=)t(f 1 *)t( f 2 =)t( f 2 *)t(f 1 表12-1 常用信号卷积表 （一）两个矩形脉冲信号的卷积过程 两信号)t(x与)t(h都为矩形脉冲信号，如图12-1所示。下面由图解的方法（图12-1）给出两个信号的卷积过程和结果，以便与实验结果进行比较。 图解法的一般步骤为：

图12-1两矩形脉冲的卷积积分的运算过程与结果（二）矩形脉冲信号与锯齿波信号的卷积 信号)(t f 1为锯齿波信号，)t( f 2 为矩形脉冲信号，如图12-2所示。根据卷积积分的运算 方法得到)t(f 1和)t( f 2 的卷积积分结果)(t y，如图12-2(i)所示。 图12-2矩形脉冲信号与锯齿脉冲信号的卷积积分的结果 （三）本实验进行的卷积运算的实现方法 在本实验装置中采用了DSP数字信号处理芯片，因此在处理模拟信号的卷积积分运算时，是先通过A/D转换器把模拟信号转换为数字信号，利用所编写的相应程序控制DSP芯片实现数字信号的卷积运算，再把运算结果通过D/A转换为模拟信号输出。结果与模拟信

信号与系统MATLAB实验总汇

实验一、MATLAB 编程基础及典型实例 一、实验目的 （1）熟悉MATLAB 软件平台的使用； （2）熟悉MATLAB 编程方法及常用语句； （3）掌握MATLAB 的可视化绘图技术； （4）结合《信号与系统》的特点，编程实现常用信号及其运算。 示例一：在两个信号进行加、减、相乘运算时，参于运算的两个向量要有相同的维数，并且它们的时间变量范围要相同，即要对齐。编制一个函数型m 文件，实现这个功能。function [f1_new,f2_new,n]=duiqi(f1,n1,f2,n2) a=min(min(n1),min(n2)); b=max(max(n1),max(n2)); n=a:b; f1_new=zeros(1,length(n)); f2_new=zeros(1,length(n)); tem1=find((n>=min(n1))&(n<=max(n1))==1); f1_new(tem1)=f1; tem2=find((n>=min(n2))&(n<=max(n2))==1); f2_new(tem2)=f2; 四、实验内容与步骤 （2）绘制信号x(t)=)3 2sin(2t e t ?的曲线，t 的范围在0~30s ，取样时间间隔为0.1s 。t=0:0.1:30; y=exp(-sqrt(2)*t).*sin(2*t/3); plot(t,y);

（3）在n=[-10:10]范围产生离散序列：?? ?≤≤?=Other n n n x ,033,2)(,并绘图。n=-10:1:10; z1=((n+3)>=0); z2=((n-3)>=0); x=2*n.*(z1-z2); stem(n,x);（4）编程实现如下图所示的波形。 t=-2:0.001:3; f1=((t>=-1)&(t<=1)); f2=((t>=-1)&(t<=2)); f=f1+f2; plot(t,f); axis([-2,3,0,3]);

MATLAB实验报告卷积

实验报告 学院：机电班级：姓名：学号： 实验名称：连续时间信号卷积运算的MATLAB实现 1.实验目的：掌握卷积的概念及计算方法 2.熟悉通过调用conv()函数求解连续时间信号卷积的数值分析 法 实验环境：MATLAB 6.5.1软件 实验要求： 1、已知信号f1（t）=t/2*[ε(t)- ε(t-2)], f2（t）= [ε (t)- ε(t-1)]，通过调用conv（）函数编程实现卷积计算y(t)= f1（t）* f2（t）,画出波形。 2、已知信号f（t）=e-t *ε(t), h（t）= t2 *e-2t *ε(t),y(t)= f（t）* h（t） （1）用符号分析法编程实现计算y（t）的理论解； （2）过调用conv（）函数编程实现卷积计算y(t)的数值解并画图 实验程序及结果： 第一题： M文件 (1) function f=uCT(t) f=(t>=0); 主程序：

k1=0:p:2; k2=0:p:1; f1=k1/2.*[uCT(k1)-uCT(k1-2)]; f2=uCT(k2)-uCT(k2-1); y=conv(f1,f2)*p; k0=k1(1)+k2(1); k3=length(f1)+length(f2)-2; k=k0:p:k3*p+k0; subplot(311) plot(k1,f1); xlabel('t') ylabel('f1(t)') axis([-0.5 2.5 -0.5 1.5]) grid on subplot(312); plot(k2,f2) grid on axis([-0.5 2.5 -0.5 1.5]) xlabel('t') ylabel('f2(t)') subplot(313)

基带信号眼图实验——matlab仿真

基带信号眼图实验——matlａb 仿真

————————————————————————————————作者：————————————————————————————————日期： ?

数字基带信号的眼图实验——matla ｂ仿真 一、实验目的 1、掌握无码间干扰传输的基本条件和原理，掌握基带升余弦滚降系统的实现方法； 2、通过观察眼图来分析码间干扰对系统性能的影响，并观察在输入相同码率的ＮRZ 基带信号下，不同滤波器带宽对输出信号码间干扰大小的影响程度； 3、熟悉MATL ＡB 语言编程。 二、实验预习要求 1、复习《数字通信原理》第七章７.1节——奈奎斯特第一准则内容； 2、复习《数字通信原理》第七章7．2节——数字基带信号码型内容； 3、认真阅读本实验内容,熟悉实验步骤。 三、实验原理和电路说明 １、基带传输特性 基带系统的分析模型如图3-1所示,要获得良好的基带传输系统，就应该 () n s n a t nT δ-∑() H ω() n s n a h t nT -∑基带传输抽样判决 图3-１?基带系统的分析模型 抑制码间干扰。设输入的基带信号为()n s n a t nT δ-∑，s T 为基带信号的码元周期,则经过基 带传输系统后的输出码元为 ()n s n a h t nT -∑。其中 1 ()()2j t h t H e d ωωωπ +∞ -∞ = ? ?（3-1） 理论上要达到无码间干扰,依照奈奎斯特第一准则,基带传输系统在时域应满足： 10()0,s k h kT k =?=? ? ， 为其他整数 ?? ?(３－２) 频域应满足：

离散时间信号的产生及信号的卷积和运算实验报告2

离散时间信号的产生及信号的卷积和运算 实验报告 班级：___________ 姓名：__________ 学号:____________ 一、实验目的和原理 实验原理： （一）DTFT 和DFT 的定义及其相互关系： 序列x[n] 的DTFT 定义：∑=∞ -∞ =-n jn ωj ω x[n]e )X(e 它是关于自变量ω的复函数，且是以π2为周期的连续函数。)X(e j ω 可以表示为： )(e jX )(e X )X(e j ωim j ωre j ω+= 其中，)(e X j ω re 和)(e X j ωim 分别是)X(e j ω的实部和虚部；还可以表示为： )(ωj j ωj ωe )X(e )X(e θ= 其中，)X(e j ω 和}arg{)()X(e j ω=ωθ分别是)X(e j ω的幅度函数和相位函数； 它们都是ω的实函数，也是以π2为周期的周期函数。 序列x[n]的N 点DFT 定义： ∑∑-=-=-===10 1 22][][)(][N n kn N N n kn N j k N j W n x e n x e X k X ππ ][k X 是周期为N 的序列。 )X(e j ω与][k X 的关系：][k X 是对)X(e j ω在一个周期中的谱的等间隔N 点采样，即： k N j ω)X(e k X πω2| ][= =， 而)X(e j ω 可以通过对][k X 内插获得，即：

]2/)1)][(/2([1 ) 22sin() 22sin( ] [1----=?--= ∑N N k j N k j ω e N k N k N k X N )X(e πωπωπω (二) 线性时不变离散时间系统的变换域表示： LTI 离散时间系统的时域差分方程为： ∑∑==-=-M k k N k k k n x p k n y d )()( (1) 传递函数： 对上面的差分方程两边求z 变换，得： ∑∑∑∑=-=-=-=-=? =N k k k M k k k M k k k N k k k z d z p z X z Y z p z X z d z Y 0 00 ) () ()()( 我们定义LTI 离散时间系统的输出的Z 变换Y(z)与输入的Z 变换X(z)的比值为该系统的传递函数，即) () ()(z X z Y z H = 为系统的传递函数。 N N M M z d z d d z p z p p z D z p z H ----++++++= =......)()()(110110 分解因式 ∏-∏-=∑∑= =-=-=-=-N i i M i i N i i k M i i k z z K z d z p z H 11 11 0)1()1()(λξ ，其中i ξ和i λ称为零、极点。 利用系统的传递函数)(z H ，我们可以分析系统的零极点，稳定性及实现结构等特点。 (2) 频率响应： 因为大多数离散时间信号都可以分解为n j e ω的线性组合，所以研究输入n j e ω-的响应具有极大的意义，即当输入为n j e n x ω=][时，输出为： )()()(][) (ωωωωωj n j m m j n j m n j m e H e e m h e e m h n y === ∑∑∞ -∞ =--∞ -∞ = 这里，∑∞-∞ =-= n n j j e n h e H ωω )()(是h(n)的DTFT ，称为LTI 离散时间系统的频率

根据Matlab的信号与系统实验指导材料

基于Matlab 的信号与系统实验指导 实验一 连续时间信号在Matlab 中的表示 一、实验目的 1、学会运用Matlab 表示常用连续时间信号的方法 2、观察并熟悉这些信号的波形和特性 二、实验原理及实例分析 1、信号的定义与分类 2、如何表示连续信号？ 连续信号的表示方法有两种；符号推理法和数值法。 从严格意义上讲，Matlab 数值计算的方法不能处理连续时间信号。然而，可利用连续信号在等时间间隔点的取样值来近似表示连续信号，即当取样时间间隔足够小时，这些离散样值能被Matlab 处理，并且能较好地近似表示连续信号。 3、Matlab 提供了大量生成基本信号的函数。如： （1）指数信号：K*exp(a*t) （2）正弦信号：K*sin(w*t+phi)和K*cos(w*t+phi) （3）复指数信号：K*exp((a+i*b)*t) （4）抽样信号：sin(t*pi) 注意：在Matlab 中用与Sa(t)类似的sinc(t)函数表示，定义为：)t /()t (sin )t (sinc ππ= （5）矩形脉冲信号：rectpuls(t,width) （6）周期矩形脉冲信号：square(t,DUTY)，其中DUTY 参数表示信号的占空比

DUTY%，即在一个周期脉冲宽度（正值部分）与脉冲周期的比值。占空比默认为0.5。 （7）三角波脉冲信号：tripuls(t, width, skew)，其中skew 取值范围在-1~+1之间。 （8）周期三角波信号：sawtooth(t, width) （9）单位阶跃信号：y=(t>=0) 三、实验内容 1、验证实验内容 直流及上述9个信号 2、程序设计实验内容 （1）利用Matlab 命令画出下列连续信号的波形图。 （a ）)4/3t (2cos π+ （b ） )t (u )e 2(t -- （c ）)]2()(u )][t (cos 1[--+t u t π （2）利用Matlab 命令画出复信号) 4/t (j 2e )t (f π+=的实部、虚部、模和辐角。 四、实验报告要求 1、格式：实验名称、实验目的、实验原理、实验环境、实验内容、实验思考等 2、实验内容：程序设计实验部分源代码及运行结果图示。