控制系统的matlab仿真与设计

第二章：

2.1 x=[15 22 33 94 85 77 60]

x(6)

x([1 3 5])

x(4:end)

x(find(x>70))

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:5

L=abs(A)>3

islogical(L)

X=A(L)

2.4 A=[4,15,-45,10,6;56,0,17,-45,0]

find(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 equal

grid on

xlabel('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 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.factor ：

function f=factor(n)

if n<=1

f=1;

else

f=factor(n-1)*n;

end

2. fcircle ：

function[s,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+2^i;

i=i+1;

end

4. jcsum ：

function k=jcsum(n)

k=0;

for i=0:n

k=k+2^i;

end