连续时间系统的时域分析 实验报告.doc

成绩 实验报告 实验名称 连续时间系统的时域分析 实验班级 电子08-1 姓 名 * 学号（后两位） * 指导教师 张涛 实验日期 2010年5月6日 实验二 连续时间系统的时域分析一、实验目的：1、学习MATLAB语言的编程方法及熟悉MATLAB指令。2、深刻理解卷积运算，利用离散卷积实现连续卷积运算。3、学会用MATLAB求解系统的零状态响应、冲激响应及阶跃响应的方法。二、实验内容与步骤。1、 编写程序Q2_1，完成与两函数的卷积运算。其中：在一个图形窗口中，画出、以及卷积结果。要求每个坐标系有标题、坐标轴名称。程序如下：t0 = -3;t1 = 5;dt = 0.01;t = t0:dt:t1;x = exp(-2*t).*ut(t);h = ut(t)-ut(t-4);y = dt*conv(x,h); subplot(221)plot(t,x), grid on, title(Signal x(t), axis(t0,t1,-0.2,1.2)subplot(222)plot(t,h), grid on, title(Signal h(t), axis(t0,t1,-0.2,1.2)subplot(212)t = 2*t0:dt:2*t1; plot(t,y), grid on, title(The convolution of x(t) and h(t), axis(2*t0,2*t1,-0.1,0.6), xlabel(Time t sec);2、 编写程序Q2_2。若系统模型为： 其中 求零状态响应，画出波形。并与手工计算的零状态响应比较。程序如下：t0=0;t1=7;dt=0.01;a=1 0 3 4;b=1 3;sys=tf(b,a);t=t0:dt:t1;f=exp(-t).*ut(t);y=lsim(sys,f,t);plot(t,y);title(零状态响应);xlabel(Time(sec);ylabel(y(t); 3、编写程序Q2_3。已知某连续系统的微分方程为：，绘出该系统的冲激响应和阶跃响应的波形。程序如下：a=2 1 8;b=1;t=0:0.1:10;y=impulse(b,a,t);subplot(211);plot(t,y);title(单位冲激响应 );y=step(b,a,t);subplot(212);plot(t,y);title(单位阶跃响应 );xlabel(Time(sec) 4、编写程序Q2_4。验证卷积的相关性质。 验证性质：.选择信号x(t)的数学表达式为：x(t)、(t)和x(t)*(t)的波形.程序如下：t0=-2;t1=5;dt=0.01;t = t0:dt:t1;x = exp(-2*t).*ut(t);n=length(t);h=zeros(1,n);h(1,(0-t0)/dt+1)=1/dt;y = dt*conv(x,h); subplot(221)plot(t,x), grid on, title(Signal x(t), axis(t0,t1,-0.2,1.2)subplot(222)plot(t,h), grid on, title(Signal h(t), axis(t0,t1,-0.2,110)subplot(212)t = 2*t0:dt:2*t1; plot(t,y), grid on, title(The convolution of x(t) and h(t), axis(2*t0,2*t1,-0.1,1.5), xlabel(Time t sec);5、length()、ones()两条指令的作用是什么？说明其格式。答：length()指令是求向量的长度，length(x)返回向量x的长度,相当与求非空数组和空矩阵大小的最大值。Ones()指令全1矩阵，例如ones(n)表示一个n行n列的全1矩阵！6、（选作题）在知识扩展中，令输入信号为高频信号和低频信号的叠加（例如），通过滤波器电路，查看输出信号和输入信号的波形。初步了解滤波器的作用。程序如下：t0=0;t1=5;dt=0.01;a=1 70.7 2500;b=2500;sys=tf(b,a);t=t0:dt:t1;f=sin(2*pi*t)+0.1*sin(50*pi*t);y=lsim(sys,f,t);subplot(211)plot(t,f);title(f=f=sin(2*pi*t)+0.1*sin(50*pi*t);subplot(212);plot(t,y);title(滤波波形)xlabel(Time(sec);ylabel(y(t);三、实验心得：通过这次实验我掌握了MATLAB语言的编程方法及熟悉MATLAB指令，虽然在实验前已经做了充分的准备，在宿舍自己的