- 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
- 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
PAGE \* MERGEFORMAT 40
一.实验性质
验证性实验
二.实验目的
掌握在时域中对连续和离散时间线性非时变系统响应进行分析的方法。
三.实验内容与步骤
(1)已知系统的微分方程为+3+2r(t)=3e(t),用MATLAB画出该系统的冲激响应以及系统在输入信号e(t)=eu(t)时的零状态的波形。(改变取样时间间隔p观察仿真效果)
a.解:由微分方程写出MATLAB冲击响应程序为:
a=[1 3 2];
b=[3];
impulse(b,a)
其输出波形为:
更改其时间间隔后有:
a=[1 3 2];
b=[3];
impulse(b,a,0:0.05:10)
图形如下:
(由图形变化可知,时间越长,图形的响应越接近零。)
b.解:由微分方程写出MATLAB阶跃响应程序为:
a=[1 3 2];
b=[3];
step(b,a)
c.连续系统对输入信号为e(t)=eu(t)的零状态响应r(t)为:
H(t)
e(t) r(t)*h(t)
程序为:
a=[1 3 2];
b=[3];
p=0.01;
t=0:p:10;
x=exp(-2*t);
lsim(b,a,x,t);
r(t)的波形
(2). 已知离散系统的差分方程为:y(n)+y(n-1)+0.25y(n-2)=x(n)用MATLAB画出该系统的单位函数响应。(写出相应得程序并画出波形。)
a.解:离散系统的单位函数响应h(n)
LTI离散时间系统
x[n]=δ[n] y[n]=h[n]
程序: a=[1 1 0.25];
b=[1];
impz(b,a)
四、实验分析
实验四 系统的零极点分析
一:实验性质
验证性实验
二:实验目的
1 掌握系统函数级零极点的概念;
2 掌握对连续和离散系统的稳定性进行分析的方法。
三:实验内容与步骤
已知一连续时间线性非时变系统的系统函数为H(s)=画出系统的零极点图并判断系统的稳定性。
已知一离散系统的时间线性非时变的系统函数为H(z)=,画出系统的零极点图并判断系统的稳定性。
(1)a.用MATLAB来绘制连续系统的零极点图。
解: ∵ H(s)=
∴ H(s)=
试绘出其零极点图。
在MATLAB中输入ljdt.m函数文件后保存文件,可直接调用函数来计算系统函数的零极点并绘制零极点图。
ljdt.m输入命令为:
function [p,z]=ljdt(D,N)
p=roots(D)
z=roots(N)
p=p;
z=z;
x=max(abs([p z]));
x=x+0.01;
y=x;
hold on
axis([-x x -y y]);
plot([-x x],[0 0])
plot([0 0],[-y y])
plot(real(p),imag(p),x)
plot(real(z),imag(z),o)
title(连续系统的零极点图)
text(0.2,x-0.2,虚轴)
text(y-0.2,0.2,实轴)
保存后对应的MATLAB命令如下:
a=[1 2 -3 2 1];
b=[1 0 -4];
ljdt(a,b);
其连续系统的零极点图为:
b.判断该系统是否定?
根据MATLAB绘出零极点图即可判断出是否稳定。
∵ H(s)=
写出命令:
clf
a=[1 2 -3 2 1];
b=[1 0 -4];
ljdt(a,b);
p =
-3.1300
0.7247 + 0.6890i
0.7247 - 0.6890i
-0.3195
z =
2.0000
-2.0000
H(s)有极点在右半平面,因此系统是一个不稳定系统。
(2).a. 用MATLAB来绘制离散系统的零极点图。
∵ H(z)=在分式的分子、分母在乘以得
∴ H(z)=
试绘出其零极点图。
在MATLAB中输入nljdt.m函数文件后保存文件,可直接调用函数来计算系统函数的零极点并绘制零极点图。
nljdt.m输入命令为:
function [p,z]=nljdt(D,N)
p=roots(D)
z=roots(N)
p=p
z=z
x=max(abs([p z 1]));
x=x+0.1;
y=x;
hold on
axis([-x x -y y]);
w=0:pi/300:2*pi;
ucircle=exp(j*w);
plot(ucircle)
plot([-x x],[0,0])
plot([0,0],[-y,y])
plot(real(p),imag(
文档评论(0)