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

下载本文档

30
0
约1.76万字
约 12页
2016-06-29 发布于贵州
举报
版权申诉

控制系统的Matlb仿真与设计课后答案.doc

1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
7、VIP文档为合作方或网友上传，每下载1次，网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

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

MATLAB课后习题答案 2.1 x=[15 22 33 94 85 77 60] x(6) x([1 3 5]) x(4:end) x(find(x70)) 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=10A=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 第四章 function f=factor(n) if n=1 f=1; else f=factor(n-1)*n; end function[s,p]=fcircle(r) s=pi*r*r; p=2*pi*r; function k=jcsum1(n) k=0;i=0; while i=n k=k+2^i; i=i+1; end function k=jcsum(n) k=0; for i=0:n k=k+2^i; end 4.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 4.2[s,p]=fcircle(10) 4.3y=0;n=100; for i=1:n y=y+1/i/i; end y 4.4s=0; for i=1:5 s=s+factor(i); end s 4.5sum=0;i=1; while sum2000 sum=sum+i; i=i+1; end; n=i-2 4.6jcsum(63) jcsum1(63) 4.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 4.3 y=0;n=100; for i=1:n y=y+1/i/i; end y 4.4 s=0; for i=1:5 s=s+factor(i); end s 4.5 sum=0;i=1; while sum2000 sum=sum+i; i=i+1; end; n=i-2 4.6 i=0;k=0; while i=63 k=k+2^i; i=i+1; end k i i=0;k=