计算机仿真课后答案.doc_第1页
计算机仿真课后答案.doc_第2页
计算机仿真课后答案.doc_第3页
计算机仿真课后答案.doc_第4页
计算机仿真课后答案.doc_第5页
已阅读5页,还剩6页未读 继续免费阅读

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

第二章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:0 B(:,n-k)=TT.k; end B test=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=deconv(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=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=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: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) endend4.2functions,p=fcircle(r)s=pi*r*r;p=2*pi*r;4.3y=0;n=100;for i=1:n y=y+1/i/i;endy4.4s=0;for i=1:5 s=s+factor(i);ends4.5sum=0;i=1;while sum2000 sum=sum+i; i=i+1;end;n=i-24.6function k=jcsum(n)k=0;for i=0:n k=k+2i;end或function k=jcsum1(n)k=0;i=0;while i0) disp(系统是稳定的); else disp(系统是不稳定的); end9.5num=7*1,5;den=conv(1,0,0,conv(1,10,1,1);gm,pm,wg,wc=margin(num,den)第十章10.1ng0=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=1ng0=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=tf(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=margin(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=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(KK*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=feedbac

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论