控制系统的MATAB仿真与设计课后答案.docVIP

第二章 1x=[15 22 33 94 85 77 60] x(6) x([1 3 5]) x(4:end) x(find(x70)) 2T=[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) 3A=zeros(2,5); A(:)=-4:5 L=abs(A)3 islogical(L) X=A(L) 4A=[4,15,-45,10,6;56,0,17,-45,0] find(A=10A=20) 5p1=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)]) 6A=[11 12 13;14 15 16;17 18 19]; PA=poly(A) PPA=poly2str(PA,s) 第三章 1n=(-10:10); y=abs(n); plot(n,y,r.,MarkerSize,20) axis equal grid on xlabel(n) 2x=0:pi/100:2*pi; y=2*exp(-0.5*x).*sin(2*pi*x); plot(x,y),grid on; 3t=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; 4theta=0:0.01:2*pi; rho=sin(2*theta).*cos(2*theta); polar(theta,rho,k); 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 on z2=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 on hold off 第四章 1for 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) end end M文件： function[s,p]=fcircle(r) s=pi*r*r; p=2*pi*r; 主程序： [s,p]=fcircle(10) 3y=0;n=100; for i=1:n y=y+1/i/i; end y 4 M文件： function f=factor(n) if n=1 f=1; else f=factor(n-1)*n; end 主程序： s=0; for i=1:5 s=s+factor(i); end s 5sum=0;i=1; while sum2000 sum=sum+i; i=i+1; end; n=i-2 for循环M文件： function k=jcsum(n) k=0; for i=0:n k=k+2^i; end 主程序： jcsum(63) While循环M文件： function k=jcsum1(n) k=0;i=0; while i=n k=k+2^i; i=i+1; end 主程序： jcsum1(63) 第五章 1A=[2,1,-5,1;1,-5,0,7;0,2,1,-1;1,6,-1,-4]; b=[13,-9,6,0]; x=A\b M文件： function f=fxyz(u) x=u(1);y=u(2);z=u(3); f=x+y.^2./x/4+z.^2./y+2./z; 主程序： [U,fmin]=fminsearch(fxyz,[0.5,0.5,0.5]) 3X=linspace(0,2*pi,50); Y=sin(X); P=p