《matlab作业.doc》由会员分享，可在线阅读，更多相关《matlab作业.doc（34页珍藏版）》请在金锄头文库上搜索。
1、实验一 产生信号波形的仿真实验一、实验目的：熟悉MATLAB软件的使用，并学会信号的表示和以及用MATLAB来产生信号并实现信号的可视化。 二、实验内容： 对信号进行时域分析，首先需要将信号随时间变化的规律用二维曲线表示出来。对于简单信号可以通过手工绘制其波形，但对于复杂的信号，手工绘制信号波形显得十分困难，且难以绘制精确的曲线。 用MATLAB软件的信号处理工具箱（Signal Processing Toolbox）来产生并表示信号。一种是用向量来表示信号，另一种则是用符合运算的方法来表示信号。用适当的MATLAB语句表示信号后，可以利用MATLAB的绘图命令绘制出直观的信号波形。产生以下信
2、号波形 3sin(x)、5exp(-x)、sin(x)/x、12abs(x)/a、sqrt(a*x) 1、如下图所示为3sin(x)的信号波形 ：2、如下图所示为5exp(-x)的信号波形 ：3、如下图所示为sin(x)/x的信号波形 ：4、如下图所示为12abs(x)/3的信号波形 ：5、如下图所示为sqrt(5*x)的信号波形 ： 实验二 连续时间信号卷积及MATLAB实现 一、实验目的：熟悉使用MATLAB软件来分析连续时间信号的卷积积分运算并用图形可视化相关结果。 二、实验内容： 1卷积积分 卷积积分在信号与线形系统分析中具有非常重要的意义，是信号与系统分析的基本方法之一。 连续时间信
3、号f1(t)和f2(t)的卷积积分（简称为卷积）f(t)定义为： 由此可得到两个与卷积相关的重要结论，即是： （1） ，即连续信号可分解为一系列幅度由 决定的冲激信号及其平移信号之和； （2）线形时不变连续系统，设其输入信号为 ，单位响应为 ，其零状态响应为 ，则有： 。 可见，连续信号卷积的计算对我们进行连续信号与系统的分析具有重要的意义。 用MATLAB实现连续信号 与 卷积的过程如下： （1）将连续信号 与 以时间间隔进行取样，得到离散序列 和 ； （2）构造 与 相对应的时间向量 和 ；（3）调用conv()函数计算卷积积分 的近似向量 ； （4）构造 对应的时间向量k。 下面即是利用
4、MATLAB实现连续时间卷积的通用函数sconv()，该程序在计算出卷积积分的数值近似的同时，还绘出 的时域波形图。需要注意的是，程序中是如何构造 的对应时间向量k的？另外，程序在绘制 波形图时采用的是plot命令而不是stem命令。 function f,k=sconv(f1,f2,k1,k2,p) %计算连续信号卷积积分f(t)=f1(t)*f2(t) % f: 卷积积分f(t)对应的非零样值向量 % k：f(t)的对应时间向量 % f1: f1(t)非零样值向量 % f2: f2(t)的非零样值向量 % k1: f1(t)的对应时间向量 % k2: f2(t)的对应时间向量 % p：取样
5、时间间隔 f=conv(f1,f2); %计算序列f1与f2的卷积和f f=f*p; k0=k1(1)+k2(1); %计算序列f非零样值的起点位置 k3=length(f1)+length(f2)-2; %计算卷积和f的非零样值的宽度 k=k0:p:k3*p; %确定卷积和f非零样值的时间向量 subplot(2,2,1) plot(k1,f1) %在子图1绘f1(t)时域波形图 title(f1(t) xlabel(t) ylabel(f1(t) subplot(2,2,2) plot(k2,f2) %在子图2绘f2(t)时波形图 title(f2(t) xlabel(t) ylabel(
6、f2(t) subplot(2,2,3) plot(k,f); %画卷积f(t)的时域波形 h=get(gca,position); h(3)=2.5*h(3); set(gca,position,h) %将第三个子图的横坐标范围扩为原来的2.5倍 title(f(t)=f1(t)*f2(t) xlabel(t) ylabel(f(t)三、实验部分: 1已知两连续时间信号如下图所示，试用MATLAB求 f(t)=f1(t)*f2(t)，并绘出 f(t) 的时域波形图。（设定取样时间间隔为p） 参考程序：p=0.01; k1=0:p:2; f1=0.5*k1; k2=k1; f2=f1; f=c
7、onv(f1,f2); %f1f2f f=f*p; k0=k1(1)+k2(1); %f k3=length(f1)+length(f2)-2; %f k=k0:p:k3*p; %f subplot(2,2,1) plot(k1,f1) %1f1(t) title(f1(t) xlabel(t) ylabel(f1(t) subplot(2,2,2) plot(k2,f2) %2f2(t) title(f2(t) xlabel(t) ylabel(f2(t) subplot(2,2,3) plot(k,f); %f(t) h=get(gca,position); h(3)=2.5*h(3);
8、set(gca,position,h) %2.5 title(f(t)=f1(t)*f2(t) xlabel(t) ylabel(f(t)A、如下图所示当P=0.5时的波形：B、如下图所示当P=0.1时的波形：C、如下图所示当P=0.05时的波形：D、如下图所示当P=0.01时的波形：2已知两连续时间信号如下图所示，试用MATLAB求f(t)=f1(t)*f2(t)，并绘出f(t)的时域波形图。（设定取样时间间隔为p） 参考程序： p=1.5; k1=0:p:2 f1=2*ones(1,length(k1) k2=0:p:4 f2=ones(1,length(k2) f=conv(f1,f2)
9、; %f1f2f f=f*p; k0=k1(1)+k2(1); %f k3=length(f1)+length(f2)-2; %f k=k0:p:k3*p; %f subplot(2,2,1) plot(k1,f1) %1f1(t) title(f1(t) xlabel(t) ylabel(f1(t) subplot(2,2,2) plot(k2,f2) %2f2(t) title(f2(t) xlabel(t) ylabel(f2(t) subplot(2,2,3) plot(k,f); %f(t) h=get(gca,position); h(3)=2.5*h(3); set(gca,po
10、sition,h) %2.5 title(f(t)=f1(t)*f2(t) xlabel(t) ylabel(f(t)A、如下图所示当P=0.5时的波形：B、如下图所示当P=0.1时的波形：C、如下图所示当P=0.05时的波形：D、如下图所示当P=0.8时的波形： E、如下图所示当P=0.9时的波形：F、如下图所示当P=1时的波形： 实验三 系统时域特性的仿真分析实验一、实验目的：通过使用MATLAB仿真软件对LTI系统的时域特性进行仿真分析对系统的冲激响应和零状态响应等有更深入的理解和掌握。 二、连续系统的冲激响应、阶跃响应及MATLAB实现 对于LTI连续系统，求解系统的冲激响应h(t)和
11、阶跃响应g(t)对我们进行连续系统的分析具有非常重要的意义。MATLAB为用户提供了专门用于求连续系统冲激响应和阶跃响应并绘制其时域波形的函数impulse（）和step（）。 在调用impulse（）和step（）函数时，我们需要用向量来对连续系统进行分析。 设描述连续系统的微分方程为： 则我们可用向量a和b来表示该系统，即： a=aN，aN-1，a1，a0 b=bN，bN-1，b1，b0 注意，向量a和b的元素一定要以微分方程中时间求导的降幂次序来排列，且缺项要用0来补齐。例如，对微分方程，则表示该系统的对应向量应为 a=1 3 2，b=1 0 1。 1 impulse（）函数 函数imp
12、ulse（）将绘出由向量a和b表示的连续系统在指定时间范围内的冲激响应h(t)的时域波形图，并能求出指定时间范围内冲激响应的数值解。impulse（）函数有如下几种调用格式： （1） impulse(b,a)：该调用格式以默认方式绘出向量a和b定义的连续系统的冲激响应的时 域波形。例如描述连续系统的微分方程为，运行如下MATLAB命令： a=1 5 6; b=3 2; impulse(b,a); 则绘出系统的冲激响应波形，如图1所示。 如图1 连续系统的冲激响应1 （2） impulse(b,a,t)：绘出系统在0t时间范围内冲激响应的时域波形。对上例，若运行命令impulse(b,a,10)
13、，则绘出系统在010秒范围内冲激响应的时域波形，如图2所示 如图2 连续系统的冲激响应2 （3） impulse(b,a,t1:p:t2)：绘出在t1t2时间范围内，且以时间间隔p均匀取样的冲激响应波形。对上例，若运行命令impulse(b,a,1:0.1:2)，则绘出12秒内，每隔0.1秒取样的冲激响应的时域波形，如图3所示 如图3 连续系统的冲激响应3 （4） y=impulse(b,a,t1:p:t2)：不绘出波形，而是求出系统冲激响应的数值解。对上例，若运行命令y=impulse(b,a,0:0.2:2)，则运行结果为： y = 3.0000 1.1604 0.3110 -0.0477 -0.1726 -0.1928 -0.1716
