计算机仿真课后答案.docx

1、计算机仿真课后答案第二章2.1x=15 22 33 94 85 77 60x(6)x(1 3 5)x(4:end)x(find(x70)2.2T=1 -2 3 -4 2 -3 ;n=length(T);TT=T;for k=n-1:-1:0B(:,n-k)=TT.k;endBtest=vander(T)2.3A=zeros(2,5);A(:)=-4:5L=abs(A)3islogical(L)X=A(L)2.4A=4,15,-45,10,6;56,0,17,-45,0find(A=10&A=20)2.5p1=conv(1,0,2,conv(1,4,1,1);p2=1 0 1 1;q,r=dec

2、onv(p1,p2);cq=商多项式为 ; cr= 余多项式为 ;disp(cq,poly2str(q,s),disp(cr,poly2str(r,s)2.6A=11 12 13;14 15 16;17 18 19;PA=poly(A)PPA=poly2str(PA,s)第三章3.1n=(-10:10);y=abs(n);plot(n,y,r.,MarkerSize,20)axis equalgrid onxlabel(n)3.2x=0:pi/100:2*pi;y=2*exp(-0.5*x).*sin(2*pi*x);plot(x,y),grid on;3.3t=0:pi/50:2*pi;x=

3、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.4theta=0:0.01:2*pi;rho=sin(2*theta).*cos(2*theta);polar(theta,rho,k);3.5x,y,z=sphere(20);z1=z;z1(:,1:4)=NaN;c1=ones(size(z1);surf(3*x,3*y,3*z1,c1);hold onz2

4、=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 onhold off第四章4.1for m=100:999m1=fix(m/100);m2=rem(fix(m/10),10);m3=rem(m,10);if m=m1*m1*m1+m2*m2*m2+m3*m3*m3disp(m)endend4.2functions,p=fcircle(r)s=pi*r*r;p=2*pi*r;4.3y=0;n=100;for i=

5、1:ny=y+1/i/i;endy4.4s=0;for i=1:5s=s+factor(i);ends4.5sum=0;i=1;while sum2000sum=sum+i;i=i+1;end;n=i-24.6function k=jcsum(n)k=0;for i=0:nk=k+2i;end或function k=jcsum1(n)k=0;i=0;while i0)disp( 系统是稳定的 );elsedisp( 系统是不稳定的 );end9.5num=7*1,5;den=conv(1,0,0,conv(1,10,1,1);gm,pm,wg,wc=margin(num,den)第十章10.1

6、ng0=1; dg0=10000*1 0 -1.1772;g0=tf(ng0,dg0); %满足开环增益的为校正系统的传递函数s=kw2s(0.7,0.5) %期望的闭环主导极点ngc=rg_lead(ng0,dg0,s);gc=tf(ngc,1)g0c=tf(g0*gc);rlocus(g0,g0c);b1=feedback(g0,1); b2=feedback(g0c,1); figure,step(b1,r-,b2,b);%未校正系统的闭环传递函数% 校正后系统的闭环传递函数grid on %绘制校正前后系统的单位阶跃响应曲线10.2KK=20;s1=-2+i*sqrt(6);a=1ng

7、0=10;dg0=conv(1,0,1,4);g0=tf(ng0,dg0);ngc,dgc,k=rg_lag(ng0,dg0,KK,s1,a);gc=tf(ngc,dgc)g0c=tf(KK*g0*gc);b1=feedback(k*g0,1);b2=feedback(g0c,1);step(b1,r-,b2,b);grid on10.3KK=128;s1=-2+i*2*sqrt(3);a=2ng0=10;dg0=conv(1,0,conv(1,2,1,8);g0=tf(ng0,dg0);ngc,dgc,k=rg_lag(ng0,dg0,KK,s1,a);gc=tf(ngc,dgc)g0c=t

8、f(KK*g0*gc);rlocus(g0,g0c);b1=feedback(k*g0,1);b2=feedback(g0c,1);figure,step(b1,r-,b2,b);grid on10.4ng0=1;dg0=conv(1,0,0,1,5);g0=tf(ng0,dg0);w=logspace(-3,3);KK=1;Pm=50;ngc,dgc=lead4(ng0,dg0,KK,Pm,w);gc=tf(ngc,dgc);g0c=tf(KK*g0*gc);bode(KK*g0,w);hold on,bode(g0c,w);grid on,hold offgm,pm,wcg,wcp=mar

9、gin(g0c)Kg=20*log10(gm)g1=feedback(g0c,1);bode(g1),grid on,mag,phase,w=bode(g1);a=find(mag=0.707*mag(1);wb=w(a(1)max(mag)b=find(mag=max(mag)wr=w(b)10.5KK=40; Pm=50;ng0= KK *1; dg0=conv(1,0,conv(1,1,1,4);g0=tf(ng0,dg0) ;w=logspace(-2,4);ngc,dgc=fg_lead_pm(ng0,dg0,Pm,w)gc=tf(ngc,dgc),g0c=tf(g0*gc);b1=

10、feedback(g0,1); b2=feedback(g0c,1);step(b1,r-, b2,b); grid onfigure, bode(g0,r-,g0c,b,w), grid on,gm,pm,wcg,wcp=margin(g0c), Km=20*log10(gm)10.6KK=200; bp=0.3; ts=0.7; delta=0.05;ng0=1;dg0=conv(1,0,conv(0.1,1,conv(0.02 1,conv(0.01,1,0.005 1);g0=tf(ng0,dg0);w=logspace(-4,3);t=0:0.1:3;mag,phase=bode(K

11、K*g0,w);gm0,pm0,wg0,wc0=margin(mag,phase,w),gm0=20*log10(gm0)%gm0 =-15.6769%2、确定期望的开环传递函数mr=0.6+2.5*bp;wc=ceil(2+1.5*(mr-1)+2.5*(mr-1)2)*pi/ts),h=(mr+1)/(mr-1)w1=2*wc/(h+1),w2=h*w1w1=wc/10; w2=25;ng1=1/w1,1; dg1=conv(1/w2,1,conv(1,0,1,0);g1=tf(ng1,dg1);g=polyval(ng1,j*wc)/polyval(dg1,j*wc); K=abs(1/g); % 剪切频率处幅值为 1，求 K 值g1=tf(K*g1)%3、确定反馈环节传递函数h=tf(dg1,ng1); Kh=1/K;h=tf(Kh*h) % 期望频率特性的倒特性%4、验算性能指标g2=feedback(KK*g0,h); %校正后，系统的开环传递函数 b1=feedback(KK*g0,1); b2=feedback(g2,1); bode(KK*g0,r-,g2,b,h,g,w);grid onfigure,step(b1, r-,b2, b,t); grid on,pos,tr,ts,tp=stepchar(b2,delta)