资源描述:
《matlab的信号与系统实验》由会员上传分享,免费在线阅读,更多相关内容在工程资料-天天文库。
1、理学院物理1203班Q1-1:修改程序Programi」,将dt改为0.2,再执行该程序,保存图形,看看所得图形的效果如何?dt=0.01时的信号波形dt=0.2时的信号波形这两幅图形冇什么区别,哪一幅图形看起來与实际信号波形更像?答:第一幅图看起来更加圆滑,第一幅的与实际信号更像Q1-2:以Q1_2为文件名存盘,产生实门信号g2(t)和信号/(f)=g2(r)COS10^0要求在图形中加上网格线,并使用函数axis。控制图形的时间范围在・2〜2秒之间。然后执行该程序,保存所的图形。t=-2:0.01:2;y=g(t
2、).*cos(10*pi*t);plot(t,y)gridon;axis([-2,2,-2,2])title(1y=g(t).*cos(10*pi*t)1)y=g(t).*cos(1O*pi*t)Ql-3:将实验原理中所给的单位冲激信号和单位阶跃信号的函数文件在MATLAB文件编辑器屮编写好,并分别以以文件名delta和Heaviside存入work文件夹中以便于使用。抄写函数文件delta如下:抄写函数文件u如下:n>0n<0仿照前面的示例程序的编写方法,编写一个MATLAB程序,以Ql_4为文件名存盘,使Z能够在
3、同一个图形窗口屮的两个子图屮分别绘制信号x[n]=0.5
4、n
5、和X(t)=COS(27Tt)[u(t)-U(t-3)]o要求选择的时间窗能够表现出信号的主要部分(或特征)。编写的程序Ql_4如下:-5:0.01:5;x=0.5.Aabs(n);t=-5:0.01:5;y=cos(2*pi*t).*(u(t)-u(t-3));subplot(121)plot(n,x)gridon,title('x(n)=0・5・八4匕s(n)')subplot(122)plot(n,y)gridon,title(1x(t)=cos(2
6、*pi*t)(u(t)-u(t-3))1)信号x
7、n]=0-5
8、n
9、的波形图和信号x(t)=cos(2^t)
10、u(t)-u(t-3)
11、的波形图x(n)=0.5abs(n)x(t)=cos(2*pi*t)(u(t)-u(t-3))Ql-5:根据示例程序的编程方法,编写一个MATLAB程序,以Ql_5为文件名存盘,由给定信号x(t)=e_0,5lu(t)求信号y(t)=x(1.5t+3),并绘制出x(t)和y(t)的图形。编写的程序Ql_5如下:t=-5:0.01:5;x=exp(-0.5壮).*u(t);y=exp(-
12、0.5*(1・5・*t+3))・*u(1・5・*t+3);subplot(121)plot(t,x)gridontitle(1x=x(t)1);subplot(122)plot(tzy)gridontitle(1y=x(1.5•壮+3)1);信号x(t)的波形图和信号y(t)=x(l・5t+3)的波形图x=x(t)y=x(1.5*t+3)Q1-6:给定两个离散时I'可序列x[n]=0.5n{u[n]-u[n-8]}h[n]=u[n]-u[n-8]编写程序Ql_6,计算它们的卷积,并分别绘制x[n]、h[n]和它们的卷
13、积y[n]的图形。编写的程序Ql_6如下:n=-10:1:10;x=0.5.An.*(u(n)-u(n-8));h=u(n)-u(n-8);y=conv(x,h);subplot(131)stem(n,x)gridontitle('x=0.5*n.*(u(n)-u(n-8))*);subplot(132)stem(n,h)gridontitle('h=u(n)-u(n-8)');subplot(133)n=-20:1:20;stem(n,y)gridontitle(1y=conv(x,h)1);信号x[n卜h[n]和
14、y[n]的波形图x=0.5*n*(u(n)-u(n-8))h=u(n)-u(n-8)尸conv(x,h)Q1-7编写程序Ql_7,分别用lsim和dsolve函数计算并绘制由如下微分方程表示的系统在输入信号为x(t)=(e'2t-e-3t)u(t)时的零状态响应和你手工计算得到的系统零状态响应曲线。缪+3響+2曲)=8咖drdt手工计算得到的系统零状态响应的数学表达式是:y=4*exp(-t)+8*exp(-2*t)-4*exp(-3*t)*(2*exp(t)-1)-exp(-2*t)*(8*t+8*exp(-t))
15、执行程序Ql_7得到的系统响应此处粘帖执行程序Ql_7得到的系统响应(2条曲线分为上下排列的2个子图来绘制)编写的程序Ql_7如下:t=0:0•01:8;x=8*(exp(-2*t)-exp(-3*t)(t);y=lsim(1,[132],x,t);subplot(211)plot(tzy)titie(16ZDD3iDdQl_7pAp^pAIp