- 1、本文档共8页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
MATLAB 上机测验题
姓名 学号________________班级:__________________
考试要求:
1、要求独立完成不得与他人共享,答卷雷同将做不及格处理。
2 、答卷用Word 文件递交,文件名为学号+姓名.doc,试卷写上姓名及学号。
3、答卷内容包括:
(1) 程序;
(2 ) 运行结果及其分析;
(3 ) 图也要粘贴在文档中。
上机考题:
1+z −1
一、系统传递函数为H (z ) ,按照以下要求求解:
1−z −1 +0.81z −2
1)求其极零点图,判断系统的稳定性,画出系统的频谱特性;
2 )当系统输入信号为 : x (n ) [5 =+cos(0.2πn) +2sin(0.7πn)] ,
0 ≤n ≤50 时,画出系统的输出。
答:1)
num=[1,1,0];
den=[1,-1,0.81];
[a,b,c,d]=tf2ss(num,den);
sys=ss(a,b,c,d,-1);
pzmap(sys);
分析:系统极零点仅有一个零点在单位圆上,其余均在单位圆内,所以系统稳定。
num=[1,1,0];
den=[1,-1,0.81];
[a,b,c,d]=tf2ss(num,den);
sys=ss(a,b,c,d,-1);
freqz(num,den,512)
2)
num=[1,1,0];
den=[1,-1,0.81];
[a,b,c,d]=tf2ss(num,den);
sys=ss(a,b,c,d,-1);
n=0:50;
e=5+cos(0.2*pi*n)+2*sin(0.7*pi*n);
lsim(sys,e,n)
分析:图中灰色波形为输入信号,蓝色为输出信号。
s3 +11s 2 +30s
二、系统传递函数为H (s ) ,
s4 +9s3 +45s2 +97s +50
1) 画出系统的零极点图,判断稳定性;
2) 给定频率范围为[0,10],步长为0.1,画出其频率响应;
3) 画出系统的单位脉冲响应。
答:1)
num=[1,11,30,0];
den=[1,9,45,97,50];
[a,b,c,d]=tf2ss(num,den);
sys=ss(a,b,c,d);
pzmap(sys)
分析:系统除一个零点在虚轴上,其余极零点均在虚轴左侧,所以系统稳定;
2 )
num=[1,11,30,0];
den=[1,9,45,97,50];
[a,b,c,d]=tf2ss(num,den);
sys=ss(a,b,c,d);
w=0:0.2*pi:20*pi;
%2=2*pi*f
freqs(num,den,w);
3 )
num=[1,11,30,0];
den=[1,9,45,97,50];
[a,b,c,d]=tf2ss(num,den);
sys=ss(a,b,c,d);
impulse(sys)
s+3
三、系统传递函数为H (s )= 2 ,初始状态为y (0) 2, y (0)=1 ,激励信号为
s +3s+2
− t
e(t) e 2 [ε(t) =−ε(t −6)] ,求系统的零输入、零状态和全响应。
答:
零输入响应:
num=[1,3];
den=[1,3,2];
x0=[4,-1];
[a,b,c,d]=tf2ss(num,den);
sys=ss(a,b,c,d);
initial(sys,x0)
title(零输入相应)
零状态响应:
%function epxl
function y=epxl(t)
if t=0
y=1;
els
文档评论(0)