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

合集下载
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

第二章

1>>x=[15 22 33 94 85 77 60]

>>x(6)

>>x([1 3 5])

>>x(4:end)

>>x(find(x>70))

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)

3>>A=zeros(2,5);

>>A(:)=-4:5

>>L=abs(A)>3

>>islogical(L)

>>X=A(L)

4>>A=[4,15,-45,10,6;56,0,17,-45,0] >>find(A>=10&A<=20)

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')]) 6>>A=[11 12 13;14 15 16;17 18 19];

>>PA=poly(A)

>>PPA=poly2str(PA,'s')

第三章

1>>n=(-10:10)';

>>y=abs(n);

>>plot(n,y,'r.','MarkerSize',20)

>>axis equal

>>grid on

>>xlabel('n')

2>>x=0:pi/100:2*pi;

>>y=2*exp(-0.5*x).*sin(2*pi*x);

>>plot(x,y),grid on;

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;

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

4 M文件:

function f=factor(n)

if n<=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;

相关文档
最新文档