[信号2.docVIP

大连理工大学实验报告 学院（系）： 专业： 班级： 姓 名： 学号： 组： ___ 实验时间： 实验室： 实验台： 指导教师签字： 成绩： 实验2：连续时间系统分析 一、实验目的和要求 （同预习报告） 二、实验原理和内容 三、实验程序和结果 1.已知系统的微分方程为，计算该系统的单位冲激响应和单位阶跃响应。 解：代码如下： b=[1,4]; a=[1,3,2]; sys=tf(b,a); t=0:0.1:10; y=impulse(sys,t); subplot(2,1,1) plot(t,y); xlabel(时间) ylabel(单位冲激响应) b=[1 4]; a=[1 3 2]; sys=tf(b,a); t=0:0.1:10; y=step(sys,t); subplot(2,1,2) plot(t,y); xlabel(时间); ylabel(单位阶跃响应); 运行结果见上页右图： 2. 实现卷积，其中： 解：代码如下： p=0.01; % 取样时间间隔 nf=0:p:2; % f(t)对应的时间向量 f=2*((nf=0)-(nf=2)); % 序列f(n)的值 nh=0:p:5; % h(t)对应的时间向量 h=exp(-nh).*(nh=0); k=0:p:7 y=conv(f,h)*p; % 计算y(t)=f(t)*h(t) subplot(3,1,1),stairs(nf,f); % 绘制f(t)的波形 title(f(t));axis([0 3 0 2.1]); subplot(3,1,2),stairs(nh,h); % 绘制h(t)的波形 title(h(t));axis([0 3 0 1.1]); subplot(3,1,3),plot(k,y); % 绘制y(t)=f(t)*h(t)的波形 title(y(t)=f(t)*h(t)); axis([0 5 0 2.1]); 运行结果如下图： 3.已知二阶系统方程 对下列情况分别求单位冲激响应，并画出其波形。 a. b. c. d. 解：代码如下： b=[1/(1*1/3)]; a=[1,4/1,1/(1*1/3)]; sys=tf(b,a); t=0:0.1:10; y=impulse(sys,t); subplot(2,2,1) plot(t,y); xlabel(t); ylabel(h(t)); title(R=4W,L=1H,C=1/3F); axis([0 10 0 1]); b=[1/(1*1)]; a=[1,2/1,1/(1*1)]; sys=tf(b,a); t=0:0.1:10; y=impulse(sys,t); subplot(2,2,2) plot(t,y); xlabel(t); ylabel(h(t)); title(R=2Ω,L=1H,C=1F); axis([0 10 0 1]); b=[1/(1*1)]; a=[1,1/1,1/(1*1)]; sys=tf(b,a); t=0:0.1:10; y=impulse(sys,t); subplot(2,2,3) plot(t,y); xlabel(t); ylabel(h(t)); title(R=1Ω,L=1H,C=1F); axis([0 10 0 1]); b=[1/(1*1)]; a=[1,0/1,1/(1*1)]; sys=tf(b,a); t=0:0.1:10; y=impulse(sys,t); subplot(2,2,4) plot(t,y); xlabel(t); ylabel(h(t)); title(R=0Ω,L=1H,C=1F); axis([0 10 0 1]); 运行结果如下图： .4.求下列系统的零极点 （1） （2） 解：代码如