实验一连续和离散系统分析一、实验目的学习连续系统和离散系统响应的matlab求解方法;二、实验主要仪器设备和材料计算机三、实验方法、步骤及结果测试实验方法：编程，上机调试，分析实验结果；步骤：编程实现上述各实验容四、实验结果1、某系统的传递函数为:Y(s)/X(x)1(s1)(s2)试求系统的冲激响应和阶跃响应2、编制程序求解下列两个系统的单位冲激响应和阶跃响应，并绘出其图形。要求分别用filter、conv、impz三种函数完成。给出理论计算结果和程序计算结果并讨论。°)y[n]0.75y[n1]0.125y[n2]x[n]x[n1]理论计算结果：单位冲激响应y(町=-5X(-0.25)%(托)4&X(-U.5)%(町In012345h(n)1-1.751.19-0.670.355-0.18单位阶跃响应y(斤)=-1x(-0.25)^(n)+2x(-U.5)nu(n)n012345y(n)1-0.750.44-0.2340.12-0.06程序计算结果:A:单位冲激响应(3)用impz函数bpnl■亠101114161020B:单位阶跃响应(1)用Filter函数(2)用Conv函数(3)用Impz函数(")y[n]0.25{x[n1]x[n2]x[n3]x[n4]}理论计算结果：单位冲激响应：y(n)=0.25[
(3)用Impz函数结论：y=filter(p,d,x)用来实现差分方程，d
表
关于同志近三年现实表现材料材料类招标技术评分表图表与交易pdf视力表打印pdf用图表说话 pdf
示差分方程输出y的系数，p表示输入x的系数，而x表示输入序列。输出结果长度数等于x的长度；而y=conv(x,h)是用来实现卷积的，对x序列和h序列进行卷积，输出的结果个数等于x的长度与h的长度之和减去1；y=impz(p,d,N)是用来实现冲击响应的，d和p的定义同filter，N表示冲击响应输出的序列个数。3已知描述某连续系统的微分方程为y"(t)5y'(t)6y(t)2x'(t)8x(t)求该系统的单位冲激响应程序：b=[28];a=[156];sys=tf(b,a);t=0:0.1:10;y=impulse(sys,t);plot(t,y);xlabel(时间(t)');ylabel(‘y(t单位冲激响应')4计算上述系统在输入为1)离散系统的单位冲激响应h(n);5已知系统函数H(z)爲尧占'求2)输入为x(n)u(n)，求系统的零状态响应程序：num=[12];den=[10.40.12];%systemmodely1=impz(num,den,20);subplot(121);stem(y1);title('impulseresponse');n=0:20;x2=ones(1,21);y1filter=filter(num,den,x2);subplot(122);stem(n,y1filter);title('filter_step');xlabel('x');ylabel('y');x(t)el(t)时的零状态响应程序：ts=0;te=10;dt=0.01;sys=tf([28],[156]);t=ts:dt:te;f=exp(-t);y=lsim(sys,f,t);Plot(t,y)xlabel('time(sec)')ylabel('y(t)')impulseresponse0.5'f1j1520filterfttep6已知某离散因果系统的系统函数H（z）击時，试分析该系统的幅频特性。程序：num=[11];den=[1-10.5];sys=tf(num,den);w=-10:0.001:10;H=freqs(num,den,w);subplot(2,1,1);plot(w,abs(H));title('幅频特性');Q—'IJ■LII1—-10-fl-6-J-2Q246B10■*L—|l*1M1INO出击』々口2d6B10subplot(2,1,2);plot(w,angle(H));title('相频特性’)；五、思考题1)连续系统响应的计算机求解可以分为哪些方法？各是什么原理?答：①连续时间系统零状态响应的求解yIsim(sys,f,t)t表示计算系统响应的抽样点向量，f是系统输入信号向量sys是LTI系统模型，借助tf函数获得。②连续系统冲激响应和阶跃响应的求解连续时间系统冲激响应可用impulse函数直接求出，其调用形式为:y=impulse(sys,t)连续时间系统阶跃响应可用step函数直接求出，其调用形式为：y=step(sys,t)t表示计算系统响应的抽样点向量，sys是LTI系统模型.2)matlab中用于离散系统求解的命令有哪些？各基于什么求解方法？答:①y=filter(p,d,x)用来实现差分方程，d表示差分方程输出y的系数，p表示输入x的系数，而x表示输入序列。输出结果长度数等于x的长度；y=conv(x,h)是用来实现卷积的，对x序列和h序列进行卷积，输出的结果个数等于x的长度与h的长度之和减去1；y=impz(p,d,N)是用来实现冲击响应的，d和p的定义同filter,N表示冲击响应输出的序列个数。
浙公网安备 33021202002483