控制系统的matlab仿真与设计(王海英)前9章习题答案.doc

第二章：2.1 x=15 22 33 94 85 77 60x(6)x(1 3 5)x(4:end)x(find(x70)2.2 T=1 -2 3 -4 2 -3 ;n=length(T); TT=T; for k=n-1:-1:0 B(:,n-k)=TT.k; end B test=vander(T)2.3 A=zeros(2,5);A(:)=-4:5L=abs(A)3islogical(L)X=A(L)2.4 A=4,15,-45,10,6;56,0,17,-45,0find(A=10&A=20)2.5 p1=conv(1,0,2,conv(1,4,1,1);p2=1 0 1 1;q,r=deconv(p1,p2);cq=商多项式为 ; cr=余多项式为 ;disp(cq,poly2str(q,s),disp(cr,poly2str(r,s)2.6 A=11 12 13;14 15 16;17 18 19;PA=poly(A) PPA=poly2str(PA,s) 第三章：3.1 n=(-10:10);y=abs(n);plot(n,y,r.,MarkerSize,20)axis equalgrid onxlabel(n)3.2 x=0:pi/100:2*pi;y=2*exp(-0.5*x).*sin(2*pi*x);plot(x,y),grid on;3.3 t=0:pi/50:2*pi;x=8*cos(t);y=4*sqrt(2)*sin(t);z=-4*sqrt(2)*sin(t);plot3(x,y,z,p);title(Line in 3-D Space);text(0,0,0,origin);xlabel(X),ylable(Y),zlable(Z);grid;3.4 theta=0:0.01:2*pi;rho=sin(2*theta).*cos(2*theta);polar(theta,rho,k);3.5 x,y,z=sphere(20);z1=z;z1(:,1:4)=NaN;c1=ones(size(z1);surf(3*x,3*y,3*z1,c1);hold onz2=z;c2=2*ones(size(z2);c2(:,1:4)=3*ones(size(c2(:,1:4);surf(1.5*x,1.5*y,1.5*z2,c2);colormap(0,1,0;0.5,0,0;1,0,0);grid onhold off第四章：1.factor ：function f=factor(n)if n=1 f=1;else f=factor(n-1)*n;end2. fcircle ：functions,p=fcircle(r)s=pi*r*r;p=2*pi*r;3. jcsum1： function k=jcsum1(n)k=0;i=0;while i=n k=k+2i; i=i+1;end4. jcsum ：function k=jcsum(n)k=0;for i=0:n k=k+2i;end4.1 for m=100:999 m1=fix(m/100); m2=rem(fix(m/10),10); m3=rem(m,10); if m=m1*m1*m1+m2*m2*m2+m3*m3*m3 disp(m) endend4.2 s,p=fcircle(10)4.3 y=0;n=100;for i=1:n y=y+1/i/i;endy4.4 s=0;for i=1:5 s=s+factor(i);ends第五章：1. fxyz：function f=fxyz(u)x=u(1);y=u(2);z=u(3);f=x+y.2./x/4+z.2./y+2./z;5.1 A=2,1,-5,1;1,-5,0,7;0,2,1,-1;1,6,-1,-4;b=13,-9,6,0;x=Ab5.2 U,fmin=fminsearch(fxyz,0.5,0.5,0.5)5.3 X=linspace(0,2*pi,50);Y=sin(X);P=polyfit(X,Y,3)AX=linspace(0,2*pi,50);Y=sin(X);Y1=polyval(P,X)plot(X,Y,:O,X,Y1,-*)5.4 x=0:2.5:10;h=0:30:60;T=95,14,0,0,0;88,48,32,12,6;67,64,54,48,41;xi=0:0.5:10;hi=0:10:60;temps=interp2(x,h,T,xi,hi,cubic);mesh(xi,hi,temps);第六章：6.1 syms xy=finverse(1/tan(x)6.2 syms x yf=1/(1+x2); g=sin(y);fg=compose(f,g)6.3 syms xg=(exp(x)+x*sin(x)(1/2);dg=diff(g)6.4 F=int(int(x*exp(-x*y),x),y)6.7 f=solve(a*x2+b*x+c)6.8 f=solve(x+y+z=1,x-y+z=2,2*x-y-z=1)6.10 a=maple(simplify(sin(x)2+cos(x)2);)6.12 syms t x F=sin(x*t+2*t); L=laplace(F)第七章：1.ww:function sys,x0,str,ts=ww(t,x,u,flag)%定义连续系统的S函数A=0,1;-0.4,-0.2;B=0;0.2;C=1,0;D=0;switch flag,case 0, sys,x0,str,ts=mdlInitializeSizes(A,B,C,D);case 1, sys=mdlDerivatives(t,x,u,A,B,C,D);case 2, sys=mdlUpdate(t,x,u);case 3, sys=mdlOutputs(t,x,u,A,B,C,D);case 4, sys=mdlGetTimeOfNextVarHit(t,x,u);case 9, sys=mdlTerminate(t,x,u);otherwise error(Unhandled flag = ,num2str(flag);end%=function sys,x0,str,ts=mdlInitializeSizes(A,B,C,D)sizes = simsizes;sizes.NumContStates = 2;sizes.NumDiscStates = 0;sizes.NumOutputs = 1;sizes.NumInputs = 1;sizes.DirFeedthrough = 1;sizes.NumSampleTimes = 1; sys = simsizes(sizes);x0 = 0;0;str = ;ts = 0 0;%=function sys=mdlDerivatives(t,x,u,A,B,C,D)sys = A*x+B*u;%=function sys=mdlUpdate(t,x,u)sys = ;%=function sys=mdlOutputs(t,x,u,A,B,C,D)sys = C*x+D*u;%=function sys=mdlGetTimeOfNextVarHit(t,x,u)sampleTime = 1; sys = t + sampleTime;%=function sys=mdlTerminate(t,x,u)sys = ; 7.1 7.2 7.37.4 7.57.67.7第八章：8.1 （1） num=5; den=1,2,2;sys=tf(num,den)（2）s = tf(s);H = 5/(s2+2*s+2) ; H.inputdelay =2（3）h=tf(0.5,0,1,-0.5,0.5,0.1)8.2 num=2*1,0.5;den=1,0.2,1.01;sys=tf(num,den)z,p,k=tf2zp(num,den);zpk(z,p,k)A,B,C,D=tf2ss(num,den);ss(A,B,C,D)8.3 num=1,5; den=1,6,5,1; ts=0.1; sysc=tf(num,den); sysd=c2d(sysc,ts,tustin) 8.4 （1） %对系统方框图每个环节进行编号,有8个通道,列写每个通道传递函数r1=1;r2=2;c1=3;c2=4;G1=r1;G2=tf(1,c1,0);G3=1;%是分离点和汇合点的连线,不能合并,传函为1G4=-1;G5=1/r2;G6=tf(1,c2,0);G7=-1;G8=-1;%建立无连接的状态空间模型G=append(G1,G2,G3,G4,G5,G6,G7,G8)%写出系统的连接矩阵Q=1 4 0 %通道1的输入是通道4 2 1 7 %通道2的输入是通道1,7 3 2 0 4 2 0 5 3 8 6 5 0 7 5 0 8 6 0;%负号在传函中体现%列出系统的总的输入和输出端的编号inputs=1;outputs=6;%生成组合后系统的状态空间模型sys=connect(G,Q,inputs,outputs)（2）r1=1;r2=2;c1=3;c2=4;A,B,C,D=linmod(x84);num,den=ss2tf(A,B,C,D);sys=tf(num,den)8.5A=1,1,0;0,1,0;0,0,2;B=0,0;1,0;0,-2;n=size(A)Tc=ctrb(A,B);if n=rank(Tc) disp(系统完全能控);else disp(系统不完全能控);end第九章：1. routh：function rtab,info=routh(den)info=; vec1=den(1:2:length(den); nrT=length(vec1);vec2=den(2:2:length(den)-1); rtab=vec1; vec2, zeros(1,nrT-length(vec2);for k=1:length(den)-2, alpha(k)=vec1(1)/vec2(1); for i=1:length(vec2), a3(i)=rtab(k,i+1)-alpha(k)*rtab(k+1,i+1