内蒙古科技大学20082009学年第1学期答案

08-09 A答案
课程名称：《控制系统仿真》考试班级：自动化、电气2005
考试时间：2009年1月8日15时00分至17时00分标准制订人：李琦、梁丽
一、填空题（每空1分，共15分）
1、模型、实际系统、数学、物理、混合
2、rlocfind、nyquist
3、zp2ss、tf2ss
4、大、小、小、大
5、模块法、菜单法
二、选择题（每题1分，共25分）
BCADB ACDCA ABCAB DDCAA BCBCA
三、简答题（共28分）
1、r=10.2;
num0=[1.875e6 0.25e6];
den0=[1 24 102 22.4 20]; （0.5分）
numh=0.001;denh=1; （0.5分）
[num,den]=feedback(num0,den0,numh,denh); （1分）
[A，b，C，D]=tf2ss(num,den); （1分）
Tf=input('仿真时间Tf=');
h=input('计算步长h=');
x=[zeros(length(A),1)];y=0;t=0;
for i=1:Tf/h
k1=A*x+b*r;
k2=A*(x+h*k1/2)+b*r;
k3=A*(x+h*k2/2)+b*r;
k4=A*(x+h*k3)+b*r;
x=x+h*(k1+2*k2+2*k3+k4)/6;
y=[y;C*x];t=[t;t(i)+h]; （6分）
end
plot(t,y) （1分）
2、A=[-8 -16 -6;1 0 0;0 1 0];
B=[1;0;0];
C=[2 8 6];
D=0; （1分）
[num,den]=ss2tf(A,B,C,D); （3分）
[z,p,k]=ss2zp(A,B,C,D); （3分）
[R,P,H]=residue(num,den) （3分）
3、
x=0:0.01:2*pi; （1分）
y=sin(x)+tan(sin(x)); （1分）
plot(x,y) （2分）
xlabel('this is x axis') （1分）
ylabel('this is y axis') （1分）
title('this is example') （1分）
grid on （1分）
四、系统分析（共4题，每题8分，共32分）
1、num=500; （1分）
den=conv([1 10 50],[1 10]); （1分）
figure(1)
step(num,den) （3分）
figure(2)
impulse(num,den) （3分）
2、num0=conv([0.2],[1 0.5]); （1分）
den0=conv([1 0],conv([1 0.5],conv([1 0.7],[1 3]))); （1分）[num,den]=cloop(num0,den0); （1分）
[z,p,k]=tf2zp(num,den); （1分）
ii=find(real(p)>0);n1=length(ii); （1分）
if(n1>0)
disp('system is unstable')
else
disp('system is stable')
end （1分）判断语句1分（）pzmap(num,den) （2分）
3、num=[6 12 24 5]; （1分）
den=[1 2 3 4 1]; （1分）
[gm,pm,wcg,wcp]=margin(num,den) （3分）
bode(num,den) （3分）
4、for m=-1:1 （1分）
num=[1,m];
den=conv([1,2],conv([1,3],[1,4])); （1分）
[a,b,c,d]=tf2ss(num,den); （2分）
if(rank(obsv(a,c))==length(a)) （4分）
disp('系统是可观的')
else disp('系统是不可观的') end end。