控制系统的matlab仿真与设计Word下载.docx

1、find（A=10&A=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,）2.6 A=11 12 13;14 15 16;17 18 19;PA=poly（A） PPA=poly2str（PA,） 第三章：3.1 n=（-10:10）y=abs（n）;plot（n,y,r.,MarkerSize,20）axis equalgrid onxlabel（n）3.2 x=0:pi/100:2*pi;y=2

2、*exp（-0.5*x）.*sin（2*pi*x）;plot（x,y）,grid on;3.3 t=0:pi/50: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 Spacetext（0,0,0,originX）,ylable（Y）,zlable（Zgrid;3.4 theta=0:0.01:rho=sin（2*theta）.*cos（2*theta）;polar（theta,rho,k3.5 x,y,z=sphere（20）;z1=z;z1（:,1:4）=NaN;c1=one

3、s（size（z1）;surf（3*x,3*y,3*z1,c1）;hold onz2=z;c2=2*ones（size（z2）;c2（:4）=3*ones（size（c2（:4）;surf（1.5*x,1.5*y,1.5*z2,c2）;colormap（0,1,0;0.5,0,0;1,0,0）;hold off第四章：1.factor ：function f=factor（n）if n=1 f=1;else f=factor（n-1）*n;end2. fcircle ：functions,p=fcircle（r）s=pi*r*r;p=2*pi*r;3. jcsum1： function k=j

4、csum1（n）k=0;i=0;while i=n k=k+2i; i=i+1;4. jcsum ：function k=jcsum（n）for i=0:n4.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） end4.2 s,p=fcircle（10）4.3 y=0;n=100;for i=1: y=y+1/i/i;y4.4 s=0;5 s=s+factor（i）;s第五章：1. fxyz：function f=fxyz（u）x=

5、u（1）;y=u（2）;z=u（3）;f=x+y.2./x/4+z.2./y+2./z;5.1 A=2,1,-5,1;1,-5,0,7;0,2,1,-1;1,6,-1,-4;b=13,-9,6,0x=Ab5.2 U,fmin=fminsearch（fxyz,0.5,0.5,0.5）5.3 X=linspace（0,2*pi,50）;Y=sin（X）;P=polyfit（X,Y,3）AX=linspace（0,2*pi,50）;Y1=polyval（P,X）plot（X,Y,:O,X,Y1,-*5.4 x=0:2.5:10;h=0:30:60T=95,14,0,0,0;88,48,32,12,6

6、;67,64,54,48,41;xi=0:0.5:10;hi=0:10:temps=interp2（x,h,T,xi,hi,cubicmesh（xi,hi,temps）;第六章：6.1 syms xy=finverse（1/tan（x）6.2 syms x yf=1/（1+x2）; g=sin（y）;fg=compose（f,g）6.3 syms xg=（exp（x）+x*sin（x）（1/2）;dg=diff（g）6.4 F=int（int（x*exp（-x*y）x）,y6.7 f=solve（a*x2+b*x+c6.8 f=solve（x+y+z=1x-y+z=22*x-y-z=16.10

7、 a=maple（simplify（sin（x）2+cos（x）2）;6.12 syms t x F=sin（x*t+2*t）; L=laplace（F）第七章：1.ww:function sys,x0,str,ts=ww（t,x,u,flag）%定义连续系统的S函数A=0,1;-0.4,-0.2;B=0;0.2;C=1,0;D=0;switch flag,case 0, sys,x0,str,ts=mdlInitializeSizes（A,B,C,D）;case 1, sys=mdlDerivatives（t,x,u,A,B,C,D）;case 2, sys=mdlUpdate（t,x,u）

8、;case 3, sys=mdlOutputs（t,x,u,A,B,C,D）;case 4, sys=mdlGetTimeOfNextVarHit（t,x,u）;case 9, sys=mdlTerminate（t,x,u）;otherwise error（Unhandled flag = ,num2str（flag）;%=function sys,x0,str,ts=mdlInitializeSizes（A,B,C,D）sizes = simsizes;sizes.NumContStates = 2;sizes.NumDiscStates = 0;sizes.NumOutputs = 1;s

9、izes.NumInputs = 1;sizes.DirFeedthrough = 1;sizes.NumSampleTimes = 1;sys = simsizes（sizes）;x0 = 0;0;str = ;ts = 0 0;function sys=mdlDerivatives（t,x,u,A,B,C,D）sys = A*x+B*u;function sys=mdlUpdate（t,x,u）sys = ;function sys=mdlOutputs（t,x,u,A,B,C,D）sys = C*x+D*u;function sys=mdlGetTimeOfNextVarHit（t,x,

10、u）sampleTime = 1;sys = t + sampleTime;function sys=mdlTerminate（t,x,u） 7.1 7.2 7.37.47.57.67.7第八章：8.1 （1） num=5; den=1,2,2;sys=tf（num,den） （2）s = tf（H = 5/（s2+2*s+2） ;H.inputdelay =2（3）h=tf（0.5,0,1,-0.5,0.5,0.1）8.2 num=2*1,0.5;den=1,0.2,1.01;sys=tf（num,den）z,p,k=tf2zp（num,den）;zpk（z,p,k）A,B,C,D=tf2s

11、s（num,den）;ss（A,B,C,D）8.3 num=1,5; den=1,6,5,1; ts=0.1; sysc=tf（num,den）; sysd=c2d（sysc,ts,tustin） 8.4 （1） %对系统方框图每个环节进行编号,有8个通道,列写每个通道传递函数r1=1;r2=2;c1=3;c2=4;G1=r1;G2=tf（1,c1,0）;G3=1;%是分离点和汇合点的连线,不能合并,传函为1G4=-1;G5=1/r2;G6=tf（1,c2,0）;G7=-1;G8=-1;%建立无连接的状态空间模型G=append（G1,G2,G3,G4,G5,G6,G7,G8）%写出系统的连接矩阵Q=1 4 0 %通道1的输入是通道4 2 1 7 %通道2的输入是通道1,7 3 2 0 4 2 0 5 3 8 6 5 0 7 5 0 8 6 0;%负号在传函中体现%列出系统的总的输入和输出端的编号inputs=1;outputs=6;%生成组合后系统的状态空间模型