当前位置：文档之家 › 信号系统MATLAB实验报告

信号系统MATLAB实验报告

  • 格式：doc
  • 大小：290.50 KB
  • 文档页数：31

信号与系统实验报告桂林理工大学信息科学与工程学院 电子信息工程实验二 信号及其表示【实验目的】了解各种常用信号的表达方式 掌握部分绘图函数 【实验内容】一、绘出连续时间信号x(t)=te707.0 sin32t 关于t 的曲线，t 的范围为 0~30s ，并以0.1s 递增。

MATLAB 源程序为：t=0:0.1:30; %对时间变量赋值x=exp(-0.707*t).*sin(2/3.*t); %计算变量所对应得函数值 plot(t,x);grid; %绘制函数曲线 ylabel('x(t)');xlabel('Time(sec)')二、产生周期为0.02的方波。

MATLAB源程序为：Fs=100000;t=0:1/Fs:1;x1=square(2*pi*50*t,20);x2=square(2*pi*50*t,80);subplot(2,1,1),plot(t,x1),axis([0,0.2,-1.5,1.5]); subplot(2,1,2),plot(t,x2),axis([0,0.2,-1.5,1.5]);三、产生sinc(x)函数波形。

MATLAB源程序为：x=linspace(-4,4);y=sinc(x);plot(x,y)四、绘制离散时间信号的棒状图。

其中x(-1)=-1,x(0)=1,x(1)=2,x(2)=1,x(3)=0,x(4)=-1,其他时间x(n)=0。

MATLAB源程序为：n=-3:5; %定位时间变量x=[0,0,-1,1,2,1,-1,0,0];stem(n,x);grid; %绘制棒状图line([-3,5],[0,0]); %画X轴线xlabel('n');ylabel('x[n]')五、单位脉冲序列 δ（n-0n ）={00...1...0n n n n =≠直接实现：x=zeros(1,N);x(1,n0)=1;函数实现：利用单位脉冲序列)(0n n -δ的生成函数impseq,即 function[x,n]=impseq(n0,ns,nf) n=[ns:nf];x=[(n-n0)==0]; plot(n,x); stem(n,x);输入参数：impseq(0,0,9)——连续图形12345678900.10.20.30.40.50.60.70.80.91输入参数：impseq(0,0,9)——离散图形六、单位阶跃序列ε（n-0n ）={00...1...0n n n n ≥<直接实现：n=[ns:nf];x=[(n-n0)>=0];函数实现：利用单位阶跃序列)(0n n -ε的生成函数stepseq ，即 Function[x,n]=stepseq(n0,ns,nf) n=[ns:nf];x=[(n-n0)>=0]; plot(n,x);七、实指数序列 R a n a n x n∈∀=,,)( 直接实现：n=[ns:nf]:x=a.^n;函数实现：利用实指数序列na n x =)(的生成函数rexpseq,即 Function[x,n]=rexpseq(a,ns,nf) n=[ns:nf];x=a,^n:八、复指数序列 n e n x n j ∀=+,)()(ωδ直接实现：n=[ns:nf];x=exp((sigema+jw)*n); 函数实现：利用复指数序列nj en x )()(ωδ+=的生成函数cexpseq,即Function[x,n]=cexpseq(sigema,w,ns,nf) n=[ns:nf];x=exp((sigema+j*w)*n);0123456789-3000-2000-10000100020003000400050006000九、正（余）弦序列 n wn n x ∀+=),cos()(θ直接实现：n=[ns:nf];x=cos(w*n+sita);函数实现：利用正（余）弦序列x(n)=cos(wn+θ)的生成函数cosswq,即 Function[x,n]=cosseq(w,ns,nf,sita) n=[ns:nf];x=cos(w*n+sita);输入参数：cosseq(3.14,0,9,30)——连续信号0123456789-0.2-0.15-0.1-0.0500.050.10.150.2输入参数：cosseq(3.14,0,9,30)——离散信号0123456789实验三信号的运算【实验目的】了解信号处理的基本操作。

熟悉一些常用的序列及其应用。

掌握编写MATLAB函数进行一些常用运算。

【实验内容】一、信号的相加与相乘MATLAB源程序为：n1=[5:4]; %序列x1(n)的时间起始及终止位置n1s=-5;n1f=4;x1=[2,3,1,-1,3,4,2,1,-5,-3]; %序列x1(n)不同时间的幅度n2=[0:9]; %序列x2(n)的时间起始及终止位置n2s=0;n2f=9;x2=[1,1,1,1,1,1,1,1,1,1]; %序列x2(n)不同时间的幅度ns=min(n1s,n2s);nf=max(n1f,n2f); %求取新信号的时间起始及终止位置n=ns:nf;y1=zeros(1,length(n)); %延拓序列初始化y2=zeros(1,length(n));y1(find((n>=n1s)&(n<=n1f)==1))=x1; %给延拓序列y1赋值x1y2(find((n>=n2s)&(n<=n2f)==1))=x2; %给延拓序列y2赋值x2ya=y1+y2; %逐点相加yp=y1.*y2; %逐点相乘subplot(4,1,1),stem(n,y1,'.');line([n(1),n(end)],[0,0]);ylabel('x1(n)'); %画x轴subplot(4,1,2),stem(n,y2,'.');line([n(1),n(end)],[0,0]);ylabel('x2(n)');subplot(4,1,3),stem(n,ya,'.');line([n(1),n(end)],[0,0]);ylabel('x1(n)+x2(n)');subplot(4,1,4),stem(n,yp,'.');line([n(1),n(rnd)],[0,0]);ylabel('x1(n).x2(n)');二、序列移位与周期延拓序列移位的数学描述为：y(n)=x(n-m)。

序列移位的METLAB 实现：y=x;ny=nx-m 。

序列周期延拓的数学描述为))(()(n x n y M= ，其中M 表示延拓周期。

序列周期延拓的MA TLAB 实现：ny=nxs:nxf;y=x(mod(ny,M)+1)。

MATLAB 源程序为： N=24;M=8;m=3; n=0:N-1;x1=(0.8).^n; %生成指数序列x2=[(n>=0)&(n<M)]; %形成矩形序列)(n R M x=x1.*x2; %截取操作形成新序列x(n) xm=zeros(1,N); for k=m+1:m+Mxm(k)=x(k-m); %产生序列移位x(n-3)endxc=x(mod(n,M)+1); %产生x(n)的周期延拓8))((n x xcm=x(mod(n-m,M)+1); %产生移位序列x(n-3)的周期延拓8))3((-n x subplot(4,1,1),stem(n,x,'.');ylabel('x(n)'); subplot(4,1,2),stem(n,xm,'.');ylabel('x(n-3)'); subplot(4,1,3),stem(n,xc,'.');ylabel('x((n))_8'); subplot(4,1,4),stem(n,xcm,'.');ylabel('x((n-3))_8');三、序列反转与序列累加运算序列反转的数学描述为：y(n)=x(-n)。

序列反转的METLAB可由函数fliplr来实现。

该函数的功能将行向量左右翻转，其调用格式：y=fliplr(x)。

序列累加的数学描述为：y(n)=∑)(i x。

序列累加的MATLAB实现可由函数cumsum 来实现，调用格式为：y=cumsum(x)。

例、求序列x(n)=3e n2.0-的反转序列y(n)=x(-n),MATLAB源程序为：n=0:10; %x(n)序列的时间序列x=3*exp(-0.2*n); %x(n)序列大小y=fliplr(x); %x(n)序列反转n1=-fliplr(n); %时间序列的反转（反转点为原点）n2=fliplr(-(n-3)); %在指定位置m=3处的时间序列的反转subplot(2,1,1);stem(n,x); %绘制x(n)序列xlabel('n'),ylabel('x(n)');subplot(2,1,2);stem(n1,y); %绘制反转y(n)=x(-n)序列xlabel('n'),ylabel('y(n)=x(-n)');s=cumsum(x); %求累加序列s(n)figure(2);subplot(2,1,1);stem(n2,y); %绘制点为3处的反转序列y(n)=x(-n+3)xlabel('n'),ylabel('y(n)=x(-n+3)');subplot(2,1,2);stem(n,s); %绘制s(n)序列xlabel('n'),ylabel('s(n)');四、两序列的卷积运算两序列的卷积运算的数学描述为：y(n)=x 1(n)*x 2(n)=)()(21m n m x x m-∑两序列的卷积运算的MA TLAB 实现：y=conv(x1,x2)。

序列x1(n)和x2(n)必须长度有限。

例、计算下列卷积，并图示各序列及卷积结果。

)()(),()(),(*)()(1012011119.0n n n n n n n R h Rx h x y n===；).(),5()(),(*)()(10220522229.0n n n n n n R h Rx h x y n =-==-MATLAB 源程序如下：Nx=20;Nh=10;m=5; %设定Nx,Nh 和位移值m n=0:Nx-1;x1=(0.9).^n; %产生x1(n) x2=zeros(1,Nx+m);for k=m+1:m+Nx %产生x2(n)=x1(n-m) x2(k)=x1(k-m) %完成位移 endnh=0:Nh-1;h1=ones(1,Nh); %产生h1(n) h2=h1;y1=conv(x1,h1); %计算y1(n)=x1(n)*h1(n) y2=conv(x2,h2); %计算y2(n)=x2(n)*h2(n) x2 =0 0 0 0 0 1 0 0 0 0 0 0 0 0Columns 15 through 250 0 0 0 0 0 0 0 0 0 0x2 =Columns 1 through 80 0 0 0 0 1.0000 0.9000 0Columns 9 through 160 0 0 0 0 0 0 0Columns 17 through 240 0 0 0 0 0 0 0Column 25x2 =Columns 1 through 80 0 0 0 0 1.0000 0.9000 0.8100Columns 9 through 160 0 0 0 0 0 0 00 0 0 0 0 0 0 0Column 25x2 =Columns 1 through 80 0 0 0 0 1.0000 0.9000 0.8100Columns 9 through 160.7290 0 0 0 0 0 0 0Columns 17 through 240 0 0 0 0 0 0 0Column 25x2 =Columns 1 through 80 0 0 0 0 1.0000 0.9000 0.8100Columns 9 through 160.7290 0.6561 0 0 0 0 0 00 0 0 0 0 0 0 0Column 25x2 =Columns 1 through 80 0 0 0 0 1.0000 0.9000 0.8100Columns 9 through 160.7290 0.6561 0.5905 0 0 0 0 0Columns 17 through 240 0 0 0 0 0 0 0Column 25x2 =Columns 1 through 80 0 0 0 0 1.0000 0.9000 0.8100Columns 9 through 160.7290 0.6561 0.5905 0.5314 0 0 0 00 0 0 0 0 0 0 0Column 25x2 =Columns 1 through 80 0 0 0 0 1.0000 0.9000 0.8100Columns 9 through 160.7290 0.6561 0.5905 0.5314 0.4783 0 0 0Columns 17 through 240 0 0 0 0 0 0 0Column 25x2 =Columns 1 through 80 0 0 0 0 1.0000 0.9000 0.8100Columns 9 through 160.7290 0.6561 0.5905 0.5314 0.4783 0.4305 0 00 0 0 0 0 0 0 0Column 25x2 =Columns 1 through 80 0 0 0 0 1.0000 0.9000 0.8100Columns 9 through 160.7290 0.6561 0.5905 0.5314 0.4783 0.4305 0.3874 0Columns 17 through 240 0 0 0 0 0 0 0Column 25x2 =Columns 1 through 80 0 0 0 0 1.0000 0.9000 0.8100Columns 9 through 160.7290 0.6561 0.5905 0.5314 0.4783 0.4305 0.3874 0.34870 0 0 0 0 0 0 0Column 25x2 =Columns 1 through 80 0 0 0 0 1.0000 0.9000 0.8100Columns 9 through 160.7290 0.6561 0.5905 0.5314 0.4783 0.4305 0.3874 0.3487Columns 17 through 240.3138 0 0 0 0 0 0 0Column 25x2 =Columns 1 through 80 0 0 0 0 1.0000 0.9000 0.8100Columns 9 through 160.7290 0.6561 0.5905 0.5314 0.4783 0.4305 0.3874 0.34870.3138 0.2824 0 0 0 0 0 0Column 25x2 =Columns 1 through 80 0 0 0 0 1.0000 0.9000 0.8100Columns 9 through 160.7290 0.6561 0.5905 0.5314 0.4783 0.4305 0.3874 0.3487Columns 17 through 240.3138 0.2824 0.2542 0 0 0 0 0Column 25x2 =Columns 1 through 80 0 0 0 0 1.0000 0.9000 0.8100Columns 9 through 160.7290 0.6561 0.5905 0.5314 0.4783 0.4305 0.3874 0.34870.3138 0.2824 0.2542 0.2288 0 0 0 0Column 25x2 =Columns 1 through 80 0 0 0 0 1.0000 0.9000 0.8100Columns 9 through 160.7290 0.6561 0.5905 0.5314 0.4783 0.4305 0.3874 0.3487Columns 17 through 240.3138 0.2824 0.2542 0.2288 0.2059 0 0 0Column 25x2 =Columns 1 through 80 0 0 0 0 1.0000 0.9000 0.8100Columns 9 through 160.7290 0.6561 0.5905 0.5314 0.4783 0.4305 0.3874 0.34870.3138 0.2824 0.2542 0.2288 0.2059 0.1853 0 0Column 25x2 =Columns 1 through 80 0 0 0 0 1.0000 0.9000 0.8100Columns 9 through 160.7290 0.6561 0.5905 0.5314 0.4783 0.4305 0.3874 0.3487Columns 17 through 240.3138 0.2824 0.2542 0.2288 0.2059 0.1853 0.1668 0Column 25x2 =Columns 1 through 80 0 0 0 0 1.0000 0.9000 0.8100Columns 9 through 160.7290 0.6561 0.5905 0.5314 0.4783 0.4305 0.3874 0.34870.3138 0.2824 0.2542 0.2288 0.2059 0.1853 0.16680.1501Column 250 x2 =Columns 1 through 80 0 0 0 0 1.0000 0.90000.8100Columns 9 through 160.7290 0.6561 0.5905 0.5314 0.4783 0.4305 0.38740.3487Columns 17 through 240.3138 0.2824 0.2542 0.2288 0.2059 0.1853 0.16680.1501Column 250.1351实验四 连续系统的时域分析【实验目的】进一步熟悉系统在MATLAB 中的表示 掌握利用卷积求系统的响应【实验内容】一、连续LTI 系统的响应在连续时间情况下，系统对任意输入信号x(t)的响应（即系统的输出）为y(t)=T[x(t)],则y(t)为系统输入与系统单位冲激响应的卷积积分，即y(t)= T[x(t)]=x(t)*h(t)=τττd t h x )()(-⎰+∞∞-例：某LTI 系统的单位冲激响应tet h 1.0)(-=，输入x(t)={101....1.....0≤≤t 其他，初始条件为零，求系统的响应y(t)。

合集下载

信号与系统matlab实验报告

信号与系统matlab实验报告

信号与系统MATLAB实验报告实验目的本实验旨在通过MATLAB软件进行信号与系统的相关实验，探究信号与系统的特性与应用。

实验步骤1. 准备工作在正式进行实验之前，我们需要做一些准备工作。

首先，确保已经安装好MATLAB软件，并且熟悉基本的操作方法。

其次，准备好实验所需的信号与系统数据，可以是已知的标准信号，也可以是自己采集的实际信号。

2. 信号的生成与显示使用MATLAB编写代码，生成不同类型的信号。

例如，可以生成正弦信号、方波信号、三角波信号等。

通过绘制信号波形图，观察不同信号的特点和变化。

t = 0:0.1:10; % 时间范围f = 1; % 信号频率s = sin(2*pi*f*t); % 正弦信号plot(t, s); % 绘制信号波形图3. 系统的建模与分析根据实验需求，建立相应的系统模型。

可以是线性时不变系统，也可以是非线性时变系统。

通过MATLAB进行模型的建立和分析，包括系统的时域特性、频域特性、稳定性等。

sys = tf([1, 2], [1, 3, 2]); % 系统传递函数模型step(sys); % 绘制系统的阶跃响应图4. 信号与系统的运算对于给定的信号和系统，进行信号与系统的运算。

例如，进行信号的卷积运算、系统的响应计算等。

通过MATLAB实现运算，并分析结果的意义与应用。

x = [1, 2, 3]; % 输入信号h = [4, 5, 6]; % 系统响应y = conv(x, h); % 信号的卷积运算plot(y); % 绘制卷积结果的波形图5. 实验结果分析根据实验数据和分析结果，对实验进行结果总结与分析。

可以从信号的特性、系统的特性、运算结果等方面进行综合性的讨论和分析。

实验总结通过本次实验，我们学习了如何在MATLAB中进行信号与系统的实验。

通过生成信号、建立系统模型、进行运算分析等步骤，我们深入理解了信号与系统的基本原理和应用方法。

通过实验数据和结果分析，我们对信号与系统有了更深刻的认识，并掌握了MATLAB在信号与系统实验中的应用技巧。

MATLAB实验报告(1-4)

MATLAB实验报告(1-4)

信号与系统MATLAB第一次实验报告一、实验目的1.熟悉MATLAB软件并会简单的使用运算和简单二维图的绘制。

2.学会运用MATLAB表示常用连续时间信号的方法3.观察并熟悉一些信号的波形和特性。

4.学会运用MATLAB进行连续信号时移、反折和尺度变换。

5.学会运用MATLAB进行连续时间微分、积分运算。

6.学会运用MATLAB进行连续信号相加、相乘运算。

7.学会运用MATLAB进行连续信号的奇偶分解。

二、实验任务将实验书中的例题和解析看懂，并在MATLAB软件中练习例题，最终将作业完成。

三、实验内容1.MATLAB软件基本运算入门。

1). MATLAB软件的数值计算：算数运算向量运算：1.向量元素要用”[ ]”括起来，元素之间可用空格、逗号分隔生成行向量，用分号分隔生成列向量。

2.x=x0:step:xn.其中x0位初始值，step表示步长或者增量，xn 为结束值。

矩阵运算：1.矩阵”[ ]”括起来；矩阵每一行的各个元素必须用”,”或者空格分开；矩阵的不同行之间必须用分号”;”或者ENTER分开。

2.矩阵的加法或者减法运算是将矩阵的对应元素分别进行加法或者减法的运算。

3.常用的点运算包括”.*”、”./”、”.\”、”.^”等等。

举例：计算一个函数并绘制出在对应区间上对应的值。

2).MATLAB软件的符号运算：定义符号变量的语句格式为”syms 变量名”2.MATLAB软件简单二维图形绘制1).函数y=f(x)关于变量x的曲线绘制用语：>>plot(x,y)2).输出多个图像表顺序：例如m和n表示在一个窗口中显示m行n列个图像，p表示第p个区域，表达为subplot(mnp)或者subplot(m,n,p)3).表示输出表格横轴纵轴表达范围：axis([xmax,xmin,ymax,ymin])4).标上横轴纵轴的字母：xlabel(‘x’),ylabel(‘y’)5).命名图像就在subplot写在同一行或者在下一个subplot前：title(‘……’)6).输出：grid on举例1：举例2：3.matlab程序流程控制1).for循环：for循环变量=初值:增量:终值循环体End2).while循环结构：while 逻辑表达式循环体End3).If分支：(单分支表达式)if 逻辑表达式程序模块End(多分支结构的语法格式)if 逻辑表达式1程序模块1Else if 逻辑表达式2程序模块2…else 程序模块nEnd4).switch分支结构Switch 表达式Case 常量1程序模块1Case 常量2程序模块2……Otherwise 程序模块nEnd4.典型信号的MATLAB表示1).实指数信号：y=k*exp(a*t)举例：2).正弦信号：y=k*sin(w*t+phi)3).复指数信号：举例：4).抽样信号5).矩形脉冲信号:y=square(t,DUTY) (width默认为1)6).三角波脉冲信号：y=tripuls(t,width,skew)(skew的取值在-1~+1之间，若skew取值为0则对称)周期三角波信号或锯齿波:Y=sawtooth(t,width)5.单位阶跃信号的MATLAB表示6.信号的时移、反折和尺度变换:Xl=fliplr(x)实现信号的反折7.连续时间信号的微分和积分运算1).连续时间信号的微分运算：语句格式：d iff(function,’variable’,n)Function:需要进行求导运算的函数，variable:求导运算的独立变量，n：求导阶数2).连续时间信号的积分运算：语句格式：int(function,’variable’,a,b)Function:被积函数variable:积分变量a:积分下限b:积分上限(a&b默认是不定积分)8.信号的相加与相乘运算9.信号的奇偶分解四、小结这一次实验让我能够教熟悉的使用这个软件，并且能够输入简单的语句并输出相应的结果和波形图，也在一定程度上巩固了c语言的一些语法。

MATLAB信号与系统实验报告19472[五篇范文]

MATLAB信号与系统实验报告19472[五篇范文]

MATLAB信号与系统实验报告19472[五篇范文]第一篇：MATLAB信号与系统实验报告19472信号与系统实验陈诉（5）MATLAB 综合实验项目二连续系统的频域阐发目的：周期信号输入连续系统的响应可用傅里叶级数阐发。

由于盘算历程啰嗦，最适适用MATLAB 盘算。

通过编程实现对输入信号、输出信号的频谱和时域响应的盘算，认识盘算机在系统阐发中的作用。

任务：线性连续系统的系统函数为11)(+=ωωjj H，输入信号为周期矩形波如图 1 所示，用MATLAB 阐发系统的输入频谱、输出频谱以及系统的时域响应。

-3-2-1 0 1 2 300.511.52Time(sec)图 1要领：1、确定周期信号 f(t)的频谱nF&。

基波频率Ω。

2、确定系统函数 )(Ω jn H。

3、盘算输出信号的频谱n nF jn H Y&&)(Ω=4、系统的时域响应∑∞-∞=Ω=nt jnn eY t y&)（MATLAB 盘算为y=Y_n*exp(j*w0*n“*t);要求（画出 3 幅图）：1、在一幅图中画输入信号f(t)和输入信号幅度频谱|F(jω)|。

用两个子图画出。

2、画出系统函数的幅度频谱|H(jω)|。

3、在一幅图中画输出信号y(t)和输出信号幅度频谱|Y(jω)|。

用两个子图画出。

解:(1)阐发盘算：输入信号的频谱为(n)输入信号最小周期为=2，脉冲宽度，基波频率Ω=2π/ =π，所以(n)系统函数为因此输出信号的频谱为系统响应为(2)步伐：t=linspace(-3,3,300);tau_T=1/4;%n0=-20;n1=20;n=n0:n1;%盘算谐波次数20F_n=tau_T*Sa(tau_T*pi*n);f=2*(rectpuls(t+1.75,0.5)+rectpuls(t-0.25,0.5)+rectpuls(t-2.25,0.5));figure(1),subplot(2,1,1),line(t,f,”linewidth“,2);%输入信号的波形 axis([-3,3,-0.1,2.1]);grid onxlabel(”Time(sec)“,”fontsize“,8),title(”输入信号“,”fontweight“,”bold“)%设定字体巨细，文本字符的粗细text(-0.4,0.8,”f(t)“)subplot(2,1,2),stem(n,abs(F_n),”.“);%输入信号的幅度频谱xlabel(”n“,”fontsize“,8),title(”输入信号的幅度频谱“,”fontweight“,”bold“)text(-4.0,0.2,”|Fn|“)H_n=1./(i*n*pi+1);figure(2),stem(n,abs(H_n),”.“);%系统函数的幅度频谱xlabel(”n“,”fontsize“,8),title(”系统函数的幅度频谱“,”fontweight“,”bold“)text(-2.5,0.5,”|Hn|“)Y_n=H_n.*F_n;y=Y_n*exp(i*pi*n”*t);figure(3),subplot(2,1,1),line(t,y,“linewidth”,2);%输出信号的波形 axis([-3,3,0,0.5]);grid onxlabel(“Time(sec)”,“fontsize”,8),title(“输出信号”,“fontweight”,“bold”)text(-0.4,0.3,“y(t)”)subplot(2,1,2),stem(n,abs(Y_n),“.”);%输出信号的幅度频谱xlabel(“n”,“fontsize”,8),title(“输出信号的幅度频谱”,“fontweight”,“bold”)text(-4.0,0.2,“|Yn|”)(3)波形：-3-2-1 0 1 2 300.511.52Time(sec)输入信号f(t)-20-15-10-5 0 5 10 15 2000.10.20.30.4n输入信号的幅度频谱|Fn|-20-15-10-5 0 5 10 15 2000.10.20.30.40.50.60.70.80.91n系统函数的幅度频谱|Hn|-3-2-1 0 1 2 300.10.20.30.4Time(sec)输出信号y(t)-20-15-10-5 0 5 10 15 2000.10.20.30.4n输出信号的幅度频谱|Yn| 项目三连续系统的复频域阐发目的：周期信号输入连续系统的响应也可用拉氏变更阐发。

信号与系统-MATLAB实验报告

信号与系统-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 （）plot(t,f) %用plot 函数绘制连续函数 title('f=exp(t)') %用title 函数设置图形的名称 axis([-5,5,-1,100]) %用axis 函数规定横纵坐标的范围 a=2时： t=-5:0.01:5f=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:5f=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:20f=sin(t)./t %调用正弦函数sin （）,并用sin （t ）./t 实现抽样函数 plot(t,f)title('f(t)=Sa(t)') axis([-20,-20,-0.5,1.1])ω=5时： t=-20:0.01:20 f=sin(5*t)./(5*t) plot(t,f)title('f(t)=Sa(5*t)') axis([-20,-20,-0.5,1.1]) （6）)2()(ft Sin t f π=ω=1时： t=-10:0.01:10f=sin(t) %调用正弦函数sin （） plot(t,f); title('f=sin(t)') axis([-10,10,-2,2]) ω=5时： t=-10:0.01:10 f=sin(5*t) plot(t,f);title('f=sin(5*t)') axis([-10,10,-2,2])实验结果；（1）-1-0.500.51 1.52 2.53-0.500.511.5（2）-1-0.500.51 1.52 2.53-0.500.511.5f(t)=heaviside(t)（3） a=1时：-5-4-3-2-1012345 a=2时：f=exp(2*t)-5-4-3-2-1012345 a=-2时：-5-4-3-2-1012345（4）-5-4-3-2-1012345-0.500.511.5f=R(t)（5） ω=1时：-20-15-10-505101520-0.4-0.20.20.40.60.81ω=5时：-20-15-10-505101520-0.4-0.20.20.40.60.81f(t)=Sa(5*t)（6） ω=1时：-10-8-6-4-20246810-2-1.5-1-0.500.511.52ω=5时：-10-8-6-4-20246810-2-1.5-1-0.500.511.52f=sin(5*t)实验心得体会:(1) 在 MATLAB 中，是用连续信号在等时间间隔点的样值来近似地表示连续信号的，当取样时间间隔足够小时，这些离散的样值就能较好地近似出连续信号。