控制系统的MATLAB仿真与设计课后答案

控制系统的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)

end

end

2M文件：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

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

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];

>>syms s

>>inv(s*eye(2)-a);

7>>f=solve('a*x^2+b*x+c')

8>>f=solve('x+y+z=1','x-y+z=2','2*x-y-z=1')

9>>y=dsolve('D2y+2*Dy+2*y=0','y(0)=1','Dy(0)=0') >>ezplot(y),grid on

10>>a=maple('simplify(sin(x)^2+cos(x)^2);')

11>>f=maple('laplace(exp(-3*t)*sin(t),t,s);')