(完整版)信号与系统Matlab实验作业
实验一 典型连续时间信号和离散时间信号
一、实验目的
掌握利用Matlab 画图函数和符号函数显示典型连续时间信号波形、典型时间离散信号、连续时间信号在时域中的自变量变换。
二、实验内容
1、典型连续信号的波形表示（单边指数信号、复指数信号、抽样信号、单位阶
跃信号、单位冲击信号） 1）画出教材P28习题1-1(3) ()[(63)(63)]t f t e u t u t =----的波形图。 function y=u(t) y=t>=0; t=-3:0.01:3;
f='exp(t)*(u(6-3*t)-u(-6-3*t))'; ezplot(f,t); grid on;
2）画出复指数信号()()j t f t e σω+=当0.4, 8σω==(0 波形图。 t=0:0.01:10; f1='exp(0.4*t)*cos(8*t)'; f2='exp(0.4*t)*sin(8*t)'; figure(1) ezplot(f1,t); grid on; figure(2) ezplot(f2,t); grid on; t=-10:0.01:10; f='sin(t)/t'; ezplot(f,t); grid on; t=0:0.01:10; f='(sign(t-3)+1)/2'; ezplot(f,t); grid on; 5）单位冲击信号可看作是宽度为?，幅度为1/?的矩形脉冲，即t=t 1处的冲击信号为 11111 ()()0 t t t x t t t other δ??<<+? ?=-=???? 画出0.2?=, t 1=1的单位冲击信号。 t=0:0.01:2; f='5*(u(t-1)-u(t-1.2))'; ezplot(f,t); grid on; axis([0 2 -1 6]); 2、典型离散信号的表示（单位样值序列、单位阶跃序列、实指数序列、正弦序 列、复指数序列） 编写函数产生下列序列： 1）单位脉冲序列，起点n0，终点n f，在n s处有一单位脉冲。 2）单位阶跃序列，起点n0，终点n f，在n s前序列值为0，在n s后序列值为1。 对于1）、2）小题，最后以参数n0= -10，n f=10，n s= -3为例，画出各自 波形。 、 n0=-10;nf=10;ns=-3;n=n0:nf; x1=[zeros(1,ns-n0),1,zeros(1,nf-ns)]; figure(1); stem(n,x1); title('单位脉冲序列'); x2=[zeros(1,ns-n0),1,ones(1,nf-ns)]; figure(2); stem(n,x2); title('单位阶跃序列'); 3）画出教材P21图1-26，即[][]n x n a u n =当a =1.2, 0.6, -1.5, -0.8的单边指数序 列（-2≤n ≤5）。 n=-2:5; subplot(2,2,1) x1=1.2.^n.*u(n);stem(n,x1); title('1.2^n*u(n)'); subplot(2,2,2) x2=0.6.^n.*u(n);stem(n,x2); title('0.6^n*u(n)'); subplot(2,2,3) x3=(-1.5).^n.*u(n);stem(n,x3); title('(-1.5)^n*u(n)'); subplot(2,2,4) x4=(-0.8).^n.*u(n);stem(n,x4); title('(-0.8)^n*u(n)'); 4）画出教材P21图1-27，即00[]sin(), 7 x n n π =ΩΩ= 的正弦序列（-7≤n ≤14）。 n=-7:14; x=sin(pi/7*n); stem(n,x); title('x[n]=sin(\Omega_0n) 正弦序列'); 5）画出复指数序列/6[]j n x n e π=和3[]j n x n e =的实部和虚部（-50≤n ≤50）。 n=-50:50; figure(1) x1=cos(pi/6*n);stem(n,x1); title('cos(n\pi/6) 实部'); figure(2) x2=sin(pi/6*n);stem(n,x2); title('sin(n\pi/6) 虚部'); figure(3) x3=cos(3*n);stem(n,x3); title('cos(3*n) 实部'); figure(4) x4=sin(3*n);stem(n,x4); title('sin(3*n) 虚部'); 3、信号的自变量变换 1）编写程序（函数），画出教材P10图1-13(a)即f(t)的波形(-6 2）利用1）中建立的函数，通过自变量替换方式依次画出图1-13(b)、(c)、 (d)即f(t+5)、f(-t+5)、f(-2t+5)的波形(-6 syms t; f='u(t)-u(t-2)'+(1+t)*'u(t+1)-u(t)'; subplot(2,2,1);ezplot(f,[-2,3]); axis([-2 3 -0.2 1.2]);title('f(t)');grid on; f1=subs(f,t,t+5); subplot(2,2,2);ezplot(f1,[-7,-2]); axis([-7 -2 -0.2 1.2]);title('f(t+5)');grid on; f2=subs(f,t,-t+5); subplot(2,2,3);ezplot(f2,[2,7]); axis([2 7 -0.2 1.2]);title('f(-t+5)');grid on; f3=subs(f,t,-2*t+5); subplot(2,2,4);ezplot(f3,[-1,4]); axis([-1 4 -0.2 1.2]);title('f(-2t+5)');grid on; 实验二连续和离散时间LTI系统的响应及卷积一、实验目的 掌握利用Matlab工具箱求解连续时间系统的冲激响应、阶跃响应，离散时间系统的单位样值响应，理解卷积概念。 二、实验内容 1、连续时间系统的冲击响应、阶跃响应 a. 利用impulse函数画出教材P44例2-15: LTI系统 () 3()2() dy t y t x t dt +=的冲击 响应的波形。 a=[0 1 3]; b=[0 2]; impulse(b,a); b. 利用step函数画出教材P45例2-17: LTI系统 1 ''()3'()2()'()2()2 y t y t y t x t x t ++=+的阶跃响应的波形。 a=[1 3 2]; b=[0.5 2]; step(b,a); 2、离散时间系统的单位样值响应 利用impz 函数画出教材P48例2-21: []3[1]3[2][3][]y n y n y n y n x n --+---=的单位样值响应的图形。 a=[1 -3 3 -1]; b=[0 1]; impz(b,a); 3、连续时间信号卷积 画出函数f1(t)＝(1+t)[u(t)-u(t-1)]和f2(t)=u(t-1)-u(t-2)的图形，并利用附在后面的sconv.m函数画出卷积积分f1(t)* f2(t)图形。 function sconv(f1,f2,k1,k2) f3=conv(f1,f2); ks=k1(1)+k2(1); ke=k1(end)+k2(end); k=length(k1)+length(k2)-1; k3=linspace(ks,ke,k); subplot(2,2,1) plot(k1,f1) title('f1(t)') xlabel('t') ylabel('f1(t)') subplot(2,2,2) plot(k2,f2) title('f2(t)') xlabel('t') ylabel('f2(t)') subplot(2,2,3) plot(k3,f3); h=get(gca,'position'); h(3)=2.5*h(3); set(gca,'position',h) title('f(t)=f1(t)*f2(t)') xlabel('t') ylabel('f(t)') t=-1:0.01:3; f1=(1+t).*(0.5*sign(t)-0.5*sign(t-1)); f2=(0.5*sign(t-1)-0.5*sign(t-2)); sconv(f1,f2,t,t); 4、画出教材P60例2-28中h[n]、x[n]的图形（图2-14(a)(b)），并利用conv函数求出 卷积x[n]*h[n]并画出图形（图2-14(f)）。 function dconv(x1,x2,k1,k2) x3=conv(x1,x2); 实验三连续时间周期信号的傅里叶级数 一、实验目的 掌握连续时间周期信号的傅里叶级数的展开和合成，理解吉布斯现象，掌握周期矩形脉冲信号的频谱及脉冲宽度、周期对周期信号频谱的影响。 二、实验内容 1、周期信号的傅里叶级数的展开和合成 画出如下图对称方波（取E=1、T=1），并采用有限项傅里叶级数对原函数进行逼近，画出对称方波的1、3、5、7、9、11次谐波的傅里叶级数合成波形，观察吉布斯现象。 function F_series(m) sum=0; t=-3:0.01:3; E=1;T=1; ta=T/2;w=2*pi/T; for n=1:2*m-1 fn=(2*E*ta/T)*sin(w*ta*n/2)/(w*ta*n/2); f=(E*ta/T)+cos(n*w*t)*fn-E/2; sum=sum+f; end figure(m) plot(t,sum);grid on; title([num2str(2*m-1) '次谐波的傅里叶级数合成波形']); for i=1:6 F_series(i); end 第一个信号实验的题目 1实现下列常用信号 （1）(5)u t +；（2）(1)t δ-；（3）cos(3)sin(2)t t +；（4）()[(1)(2)]f t t u t t u t t =?---； （5）0.5()4cos(),010t f t e t t π-=?= 2连续信号的基本运算与波形变换 已知信号2 2,2 1 ()33 t t f t ? -+-≤≤?=???，试画出下列各函数对时间t 的波形： （1）()f t -（2）(2)f t -+（3）(2)f t （4）1 (1)2 d f t dt +（5）(2)t f d ττ-∞-? 3连续信号的卷积运算 实现12()()f t f t *，其中1()f t 、2()f t 从第2个题目中任选3对组合。 4连续系统的时域分析 （1） 描述某连续系统的微分方程为()2()()()2()y t y t y t f t f t ''''++=+，求当输入信号为 2()2()t f t e u t -=时，该系统的零状态响应()y t 。 （2） 已知描述某连续系统的微分方程为2()()3()()y t y t y t f t '''+-=，试用MATLAB 绘出 该系统的冲激响应和阶跃响应的波形。 实验一答案： （1）(5)u t +在MATLAB 软件的输入程序及显示波形如下： （2）(1)t δ-在MATLAB 软件的输入程序及显示波形如下： （3）cos(3)sin(2)t t +在MATLAB 软件的输入程序及显示波形如下： （4）()[(1)(2)]f t t u t t u t t =?---在MATLAB 软件的输入程序及显示波形如下： （5）0.5()4cos(),010t f t e t t π-=?=在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) ; 《信号与系统及MATLAB实现》实验指导书 前言 长期以来，《信号与系统》课程一直采用单一理论教学方式，同学们依靠做习题来巩固和理解教学内容，虽然手工演算训练了计算能力和思维方法，但是由于本课程数学公式推导较多，概念抽象，常需画各种波形，作题时难免花费很多时间，现在，我们给同学们介绍一种国际上公认的优秀科技应用软件MA TLAB，借助它我们可以在电脑上轻松地完成许多习题的演算和波形的绘制。 MA TLAB的功能非常强大，我们此处仅用到它的一部分，在后续课程中我们还会用到它，在未来地科学研究和工程设计中有可能继续用它，所以有兴趣的同学，可以对MA TLAB 再多了解一些。 MA TLAB究竟有那些特点呢？ 1．高效的数值计算和符号计算功能，使我们从繁杂的数学运算分析中解脱出来； 2．完备的图形处理功能，实现计算结果和编程的可视化； 3．友好的用户界面及接近数学表达式的自然化语言，易于学习和掌握； 4．功能丰富的应用工具箱，为我们提供了大量方便实用的处理工具； MA TLAB的这些特点，深受大家欢迎，由于个人电脑地普及，目前许多学校已将它做为本科生必须掌握的一种软件。正是基于这些背景，我们编写了这本《信号与系统及MA TLAB实现》指导书，内容包括信号的MA TLAB表示、基本运算、系统的时域分析、频域分析、S域分析、状态变量分析等。通过这些练习，同学们在学习《信号与系统》的同时，掌握MA TLAB的基本应用，学会应用MA TLAB的数值计算和符号计算功能，摆脱烦琐的数学运算，从而更注重于信号与系统的基本分析方法和应用的理解与思考，将课程的重点、难点及部分习题用MA TLAB进行形象、直观的可视化计算机模拟与仿真实现，加深对信号与系统的基本原理、方法及应用的理解，为学习后续课程打好基础。另外同学们在进行实验时，最好事先预习一些MA TLAB的有关知识，以便更好地完成实验，同时实验中也可利用MA TLAB的help命令了解具体语句以及指令的使用方法。 学生实验报告 （理工类） 课程名称：信号与线性系统专业班级：M11通信工程 学生学号：1121413017 学生姓名：王金龙 所属院部：龙蟠学院指导教师：杨娟 20 11 ——20 12 学年第 1 学期 金陵科技学院教务处制 实验报告书写要求 实验报告原则上要求学生手写，要求书写工整。若因课程特点需打印的，要遵照以下字体、字号、间距等的具体要求。纸张一律采用A4的纸张。 实验报告书写说明 实验报告中一至四项内容为必填项，包括实验目的和要求；实验仪器和设备；实验内容与过程；实验结果与分析。各院部可根据学科特点和实验具体要求增加项目。 填写注意事项 （1）细致观察，及时、准确、如实记录。 （2）准确说明，层次清晰。 （3）尽量采用专用术语来说明事物。 （4）外文、符号、公式要准确，应使用统一规定的名词和符号。 （5）应独立完成实验报告的书写，严禁抄袭、复印，一经发现，以零分论处。 实验报告批改说明 实验报告的批改要及时、认真、仔细，一律用红色笔批改。实验报告的批改成绩采用百分制，具体评分标准由各院部自行制定。 实验报告装订要求 实验批改完毕后，任课老师将每门课程的每个实验项目的实验报告以自然班为单位、按学号升序排列，装订成册，并附上一份该门课程的实验大纲。 实验项目名称：常用连续信号的表示 实验学时： 2学时 同组学生姓名： 无 实验地点： A207 实验日期： 11.12.6 实验成绩： 批改教师： 杨娟 批改时间： 一、实验目的和要求 熟悉MATLAB 软件；利用MATLAB 软件，绘制出常用的连续时间信号。 二、实验仪器和设备 586以上计算机，装有MATLAB7.0软件 三、实验过程 1. 绘制正弦信号)t Asin t (f 0?ω+=（），其中A=1，πω2=，6/π?=； 2. 绘制指数信号at Ae t (f =），其中A=1，0.4a -=； 3. 绘制矩形脉冲信号，脉冲宽度为2； 4. 绘制三角波脉冲信号，脉冲宽度为4；斜度为0.5； 5. 对上题三角波脉冲信号进行尺度变换，分别得出)2t (f ，)2t 2(f -； 6. 绘制抽样函数Sa （t ）,t 取值在-3π到+3π之间； 7. 绘制周期矩形脉冲信号，参数自定； 8. 绘制周期三角脉冲信号，参数自定。 四、实验结果与分析 1．制正弦信号)t Asin t (f 0?ω+=（），其中A=1，πω2=，6/π?= 实验代码： A=1; 信号与系统综合实验项目doc 信号与系统综合实验项目 (竞 实 验 指 导 项目一 用MATLAB 验证时域抽样定理 目的： 通过MATLAB 编程实现对时域抽样定理的验证，加深抽样定理的明白得。同时训练应用运算机分析咨询题的能力。 任务： 连续信号f(t)=cos(8*pi*t)+2*sin(40*pi*t)+cos(24*pi*t)，通过理想抽样后得到抽样信号fs(t)，通过理想低通滤波器后重构信号f(t)。 方法： 1、确定f(t)的最高频率fm 。关于无限带宽信号，确定最高频率fm 的方法：设其频谱的模降到10-5左右时的频率为fm 。 2、确定Nyquist 抽样间隔T N 。选定两个抽样时刻：T S 《信号与系统》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 （） plot(t,f) %用plot 函数绘制连续函数 title('f=exp(t)') %用title 函数设置图形的名称 axis([-5,5,-1,100]) %用axis 函数规定横纵坐标的围 a=2时： t=-5:0.01:5 f=exp(2*t) % 调用指数函数exp （） plot(t,f) title('f=exp(2*t)') axis([-5,5,-1,100]) a=-2时： t=-5:0.01:5 f=exp(-2*t) plot(t,f) title('f=exp(-2*t)') axis([-5,5,-1,100]) （4）)()(t R t f = t=-5:0.01:5 f=rectpuls(t,2) % 用rectpuls(t,a)表示门函数，默认以零点为中心，宽度为a plot(t,f) title('f=R(t)') axis([-5 5 -0.5 1.5]) （5）)()(t Sa t f ω= ω=1时： t=-20:0.01:20 f=sin(t)./t % 调用正弦函数sin （）,并用sin （t ）./t 实现抽 样函数 plot(t,f) 实验一：信号的时域分析 一、实验目的 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.更改参数，调试程序，绘制单边指数信号的波形图和频谱图。观察参数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 )'); 广州大学 综合设计性实验 报告册 实验项目选频网络的设计及应用研究 学院物电学院年级专业班电子131 姓名朱大神学号成绩 实验地点电子楼316 指导老师 《综合设计性实验》预习报告 实验项目：选频网络的设计及应用研究 一 引言： 选频网络在信号分解、振荡电路及其收音机等方面有诸多应用。比如，利用选频网络可以挑选出一个周期信号中的基波和高次谐波。选频网络的类型和结构有很多，本实验将通过设计有源带通滤波器实现选频。 二 实验目的： （1）熟悉选频网络特性、结构及其应用，掌握选频网络的特点及其设计方法。 （2）学会使用交流毫伏表和示波器测定选频网络的幅频特性和相频特性。 （3）学会使用Multisim 进行电路仿真。 三 实验原理： 带通滤波器： 这种滤波器的作用是只允许在某一个通频带范围内的信号通过，而比通频带下限频率低和比上限频率高的信号均加以衰减和抑制。 典型的带通滤波器可以从二阶低通滤波器中将其中一级改成高通而成，如图1所示。 电路性能参数可由下面各式求出。 通带增益：CB R R R R A f vp 144+= 其中B 为通频带宽。 中心频率：)1 1(121 3 12 20R R C R f += π 通带宽度：)2 1(14 321R R R R R C B f -+= 品质因数：B f Q 0 = 此电路的优点是，改变f R 和4R 的比值，就可以改变通带宽度B 而不会影响中心频率0f 。 四 实验内容： 设计一个中心频率Hz f 20000=，品质因数5>Q 的带通滤波器。 五 重点问题： （1）确定带通滤波器的中心频率、上限频率及下限频率。 （2）验证滤波器是否能筛选出方波的三次谐波。 六 参考文献： [1]熊伟等.Multisim 7 电路设计及仿真应用.北京：清华大学出版社，2005. [2]吴正光，郑颜.电子技术实验仿真与实践.北京：科学出版社，2008. [4]童诗白等.模拟电子技术基础（第三版）.北京：高等教育出版社， 2001. 图1 二阶带通滤波器 实验一 名称：连续时间信号分析 姓名：王嘉琦 学号：201300800636 班级：通信二班 一、实验目的 （一）掌握使用Matlab 表示连续时间信号 1、学会运用Matlab 表示常用连续时间信号的方法 2、观察并熟悉常用信号的波形和特性 （二）掌握使用Matlab 进行连续时间信号的相关运算 1、学会运用Matlab 进行连续时间信号的时移、反褶和尺度变换 2、学会运用Matlab 进行连续时间信号微分、积分运算 3、学会运用Matlab 进行连续时间信号相加、相乘运算 4、学会运用Matlab 进行连续时间信号卷积运算 二、实验条件 Matlab 三、实验内容 1、利用Matlab 命令画出下列连续信号的波形图。 （1）)4/3t (2cos π+ 代码： k=2;w=3;phi=pi/4; t=0:0.01:3; ft=k*cos(w*t+phi); plot(t,ft),grid on; axis([0,3,-2.2,2.2]) title('余弦信号') （2） )t (u )e 2(t -- 代码： k=-1;a=-1; t=0:0.01:3; ft=2-k*exp(a*t); plot(t,ft),grid on axis([0,3,2,3]) title('指数信号') （3）)]2()(u )][t (cos 1[--+t u t π 代码： k=1;w=pi;phi=0; t=0:0.01:2; ft=1+k*cos(w*t+phi); plot(t,ft),grid on; axis([0,3,0,2]) title('余弦信号') 实验1 信号的时域描述与运算 一、实验目的 1. 掌握信号的MATLAB表示及其可视化方法。 2. 掌握信号基本时域运算的MA TLAB实现方法。 3. 利用MA TLAB分析常用信号，加深对信号时域特性的理解。 二、实验原理与方法 1. 连续时间信号的MATLAB表示 连续时间信号指的是在连续时间范围内有定义的信号，即除了若干个不连续点外，在任何时刻信号都有定义。在MATLAB中连续时间信号可以用两种方法来表示，即向量表示法和符号对象表示法。 从严格意义上来说，MATLAB并不能处理连续时间信号，在MATLAB中连续时间信号是用等时间间隔采样后的采样值来近似表示的，当采样间隔足够小时，这些采样值就可以很好地近似表示出连续时间信号，这种表示方法称为向量表示法。表示一个连续时间信号需要使用两个向量，其中一个向量用于表示信号的时间范围，另一个向量表示连续时间信号在该时间范围内的采样值。例如一个正弦信号可以表示如下： >> t=0:0.01:10; >> x=sin(t); 利用plot(t,x)命令可以绘制上述信号的时域波形，如图1所示。 如果连续时间信号可以用表达式来描述，则还可以采用符号表达式來表示信号。例如对于上述正弦信号，可以用符号对象表示如下： >> x=sin(t); >> ezplot(X); 利用ezplot(x)命令可以绘制上述信号的时域波形 Time(seconds) 图1 利用向量表示连续时间信号 t 图 2 利用符号对象表示连续时间信号 sin(t) 2.连续时间信号的时域运算 对连续时间信号的运算包括两信号相加、相乘、微分、积分，以及位移、反转、尺度变换（尺度伸缩）等。 1）相加和相乘 信号相加和相乘指两信号对应时刻的值相加和相乘，对于两个采用向量表示的可以直接使用算术运算的运算符“+”和“*”来计算，此时要求表示两信号的向量时间范围和采样间隔相同。采用符号对象表示的两个信号，可以直接根据符号对象的运算规则运算。 2）微分和积分 对于向量表示法表示的连续时间信号，可以通过数值计算的方法计算信号的微分和积分。这里微分使用差分来近似求取的，由时间向量[N t t t ,,,21?]和采样值向量[N x x x ,,,21?]表示的连续时间信号，其微分可以通过下式求得 1,,2,1,|)('1-?=?-≈ +=N k t x x t x k k t t k 其中t ?表示采样间隔。MA TLAB 中用diff 函数来计算差分 k k x x -+1。 连续时间信号的定积分可以由MATLAB 的qud 函数实现，调用格式为 quad ('function_name',a,b) 其中，function_name 为被积函数名，a 、b 为积分区间。 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．实验目的 ● 学会使用MATLAB 产生各种常见的连续时间信号与离散时间信号； ● 通过MATLAB 中的绘图工具对产生的信号进行观察，加深对常用信号的 理解； ● 熟悉MATLAB 的基本操作，以及一些基本函数的使用，为以后的实验奠 定基础。 2．实验内容 ⑴ 运行以上九个例子程序，掌握一些常用基本信号的特点及其MATLAB 实现方法；改变有关参数，进一步观察信号波形的变化。 ⑵ 在 k [10:10]=- 范围内产生并画出以下信号： a) 1f [k][k]δ=； b) 2f [k][k+2]δ=； c) 3f [k][k-4]δ=； d) 4f [k]2[k+2][k-4]δδ=-。 源程序： k=-10:10; f1k=[zeros(1,10),1,zeros(1,10)]; subplot(2,2,1) stem(k,f1k) title('f1[k]') f2k=[zeros(1,8),1,zeros(1,12)]; subplot(2,2,2) stem(k,f2k) title('f2[k]') f3k=[zeros(1,14),1,zeros(1,6)]; subplot(2,2,3) stem(k,f3k) title('f3[k]') f4k=2*f2k-f3k; subplot(2,2,4) stem(k,f4k) title('f4[k]') ⑶ 在 k [0:31]=范围内产生并画出以下信号： a) ()()k k 144f [k]sin cos π π=； b) ()2k 24f [k]cos π =； c) ()()k k 348f [k]sin cos π π=。 请问这三个信号的基波周期分别是多少？ 源程序： k=0:31; f1k=sin(pi/4*k).*cos(pi/4*k); subplot(3,1,1) stem(k,f1k) title('f1[k]') f2k=(cos(pi/4*k)).^2; subplot(3,1,2) stem(k,f2k) title('f2[k]') f3k=sin(pi/4*k).*cos(pi/8*k); subplot(3,1,3) stem(k,f3k) title('f3[k]') 其中f1[k]的基波周期是4, f2[k]的基波周期是4, f3[k]的基波周期是16。 实验一 基本信号的产生与运算 一、 实验目的 学习使用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('杨婕婕') 大连理工大学 本科实验报告 课程名称：___信号与系统实验学院：信息与通信工程学院专业：电子信息工程 班级： 学号： 学生姓名： 2012年12月11日 信号与系统实验 项目列表 信号的频谱图 Signals Frequency Spectrum 连续时间系统分析 Analysis for Continuous-time System 信号抽样 Signal Sampling 离散时间LTI系统分析 Analysis for Discrete-time LTI System 语音信号的调制解调 Modulation and Demodulation for Audio Signals Simulink?模拟信号的调制解调 Modulation and Demodulation for Analog Signals in Simulink ? 实验1信号的频谱图 一、实验目的 1. 掌握周期信号的傅里叶级数展开； 2. 掌握周期信号的有限项傅里叶级数逼近； 3. 掌握周期信号的频谱分析； 4. 掌握连续非周期信号的傅立叶变换； 5. 掌握傅立叶变换的性质。 二、实战演练（5道题） 1.已知周期三角信号如下图1-5所示，试求出该信号的傅里叶级数，利用MA TLAB编程 实现其各次谐波的叠加，并验证其收敛性。 解： 调试程序如下： clc clear t=-2:0.001:2; omega=pi; y=-(sawtooth(pi*t,0.5)/2+0.5)+1; plot(t,y),grid on; xlabel('t'),ylabel('周期三角波信号'); axis([-2 2 -0.5 1.5]) n_max=[1 3 5 11 47]; N=length(n_max); for k=1:N n=1:2: n_max(k); c=n.^2; b=4./(pi*pi*c); x=b*cos(omega*n'*t)+0.5; figure; plot(t,y,'b'); hold on; plot(t,x,'r'); hold off; xlabel('t'),ylabel('部分和的波形'); axis([-2 2 -0.5 1.5]);grid on; title(['最大谐波数=',num2str(n_max(k))]) end 运行结果如下： 实验一、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 编程实现对时域抽样定理的验证，加深抽样定理的理解。同时训练应用计算机分析问题的能力。 任务： 连续信号f(t)=cos(8*pi*t)+2*sin(40*pi*t)+cos(24*pi*t)，经过理想抽样后得到抽样信号fs(t)，通过理想低通滤波器后重构信号f(t)。 方法： 1、确定f(t)的最高频率fm 。对于无限带宽信号，确定最高频率fm 的方法：设其频谱的模降到10-5左右时的频率为fm 。 2、确定Nyquist 抽样间隔T N 。选定两个抽样时间：T S 实验一 基本信号的产生与运算 一、 实验目的 学习使用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('杨婕婕朱艺星'); 《信号与系统》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 （） plot(t,f) %用plot 函数绘制连续函数 title('f=exp(t)') %用title 函数设置图形的名称 axis([-5,5,-1,100]) %用axis 函数规定横纵坐标的范围 a=2时： t=-5:0.01:5 f=exp(2*t) %调用指数函数exp （） plot(t,f) title('f=exp(2*t)') axis([-5,5,-1,100]) a=-2时： t=-5:0.01:5 f=exp(-2*t) plot(t,f) title('f=exp(-2*t)') axis([-5,5,-1,100]) （4）)()(t R t f = t=-5:0.01:5 f=rectpuls(t,2) %用rectpuls(t,a)表示门函数，默认以零点为中心，宽度为a plot(t,f) title('f=R(t)') axis([-5 5 -0.5 1.5]) （5）)()(t Sa t f ω= ω=1时： 第二个信号实验题目 1（1）用数值法求门函数4()G t 的傅里叶变换，并给出门函数的幅频特性曲线和相频特性曲线。 （2）用符号法给出函数5()2()()3 t f t e u t -=的傅里叶变换。 （3）已知系统函数为34 2 1()3 s s H s s s ++=++，画出该系统的零极点图。 2 （1）用数值法给出函数5(2)2()(2)3 t f t e u t --=-幅频特性曲线和相频特性曲线。 （2）对函数5(2)2()(2)3 t f t e u t --=-进行采样，采样间隔为0.01。 （3）已知输入信号为()sin(100)f t t =，载波频率为1000Hz ，采样频率为5000 Hz ，试产生输入信号的调幅信号。 3（1）用符号法实现函数4()G t 的傅里叶变换，并给出门函数的幅频特性曲线和相频特性曲线。 （2）已知系统函数为34 2 1()3 s s H s s s ++= ++，输入信号为()sin(100)f t t =，求该系统的 稳态响应。 （3）已知输入信号为()sin(100)f t t =，载波频率为100Hz ，采样频率为400 Hz ，试产生输入信号的调频信号。 4（1）已知系统函数为23 1()3 s s H s s s ++= ++，画出该系统的零极点图。 （2）已知函数5()2()()3 t f t e u t -=用数值法给出函数(3)f t 的幅频特性曲线和相频特性曲线。 （3）实现系统函数3421 ()3 s s H s s s ++= ++的频率响应。 （4）已知输入信号为()cos(100)f t t =，载波频率为100Hz ，采样频率为400 Hz ， 试产生输入信号的调相信号。 5（1）用数值法给出函数5(2)2 ()(2)3 t f t e u t -+=+幅频特性曲线和相频特性曲线。 （2）用符号法实现函数 2 2i ω +的傅里叶逆变换。 （3）已知输入信号为()5sin(200)f t t =，载波频率为1000Hz ，采样频率为5000 Hz ， 试产生输入信号的调频信号。信号与系统实验题目及答案
信号与系统matlab实验及答案
信号与系统实验
信号与系统实验报告1
信号与系统综合实验项目doc信号与系统综合实验项目(竞
信号与系统 MATLAB实验报告
信号与系统实验报告_1(常用信号的分类与观察)
(完整word版)信号与系统matlab实验
信号与系统综合实验报告-带通滤波器的设计DOC
MATLAB实验报告
北京理工大学信号与系统实验实验报告
信号与系统MATLAB实验
信号与系统实验报告
matlab信号与系统实验报告
答案-信号与系统实验报告.
信号与系统MATLAB实验总汇
信号与系统综合实验项目信号与系统综合实验项目竞
matlab信号与系统实验报告
信号与系统 MATLAB实验报告
信号与系统实验二的题目及答案
- 信号与系统MATLAB实验
- 信号与系统 MATLAB实验报告
- 中南民大信号与系统MATLAB实验总汇
- (完整版)信号与系统Matlab实验作业
- 信号与系统 MATLAB实验报告
- 信号与系统matlab实验及答案
- 信号与系统MATLAB实验一报告
- matlab信号与系统实验报告
- MATLAB信号与系统实验报告
- 信号与系统MATLAB实验(教师版)
- 信号与系统matlab实验
- (完整word版)信号与系统matlab实验
- 信号与系统实验函数代码(MATLAB里的)
- 基于MATLAB的信号与系统实验教程
- 信号与系统与MATLAB实现
- MATLAB信号与系统实验报告
- 信号与系统Matlab实验作业
- 信号与系统MATLAB实验总汇
- 信号与系统matlab实验及答案
- 信号与系统Matlab实验2014版