控制系统的MATLAB仿真与设计课后答案>>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; 4>>theta=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第四章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) endend2M文件：function[s,p]=fcircle(r) s=pi*r*r; p=2*pi*r; 主程序：[s,p]=fcircle(10) 3>>y=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 5>>sum=0;i=1; while sum<2000 sum=sum+i; i=i+1; end; >>n=i-2 6for 循环 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) 第五章1>>A=[2,1,-5,1;1,-5,0,7;0,2,1,-1;1,6,-1,-4]; >>b=[13,-9,6,0]'; >>x=A\b 2M文件：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]) 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, '-*') 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); 第六章1>>syms x>>y=finverse(1/tan(x)) 2>>syms x y>>f=1/(1+x^2); g=sin(y); >>fg=compose(f,g) 3>>syms x>>g=(exp(x)+x*sin(x))^(1/2); >>dg=diff(g) 4>>F=int(int('x*exp(-x*y)', 'x' ), 'y' ) 5>>syms x>>F=ztrans(x*exp(-x*10)) 6>>a=[0 1;-2 -3...
