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

成绩 实验报告 实验名称 连续时间系统的时域分析 实验班级 电子08-1 姓 名 *** 学号（后两位） ** 指导教师 张涛 实验日期 2010年月日学会用MATLAB求解系统的零状态响应冲激响应及阶跃响应与两函数的卷积运算。 其中：在一个图形窗口中，画出、以及卷积结果。要求每个坐标系有标题、坐标轴名称。 程序如下： 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); 编写程序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*