控制系统的MATLAB仿真和设计课后答案
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 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;