控制系統(tǒng)仿真

1、精選優(yōu)質(zhì)文檔-傾情為你奉上5.2設(shè)，求函數(shù)在附近的最小值。解： fun=inline(x(1)+x(2)2/(4*x(1)+x(3)2/x(2)+2/x(3),x); x0=0.5,0.5,0.5; x fval=fminsearch(fun,x0)x = 0.5000 1.0000 1.0000fval =4.0000 函數(shù)在附近的最小值為：6.8求方程組的解。解： A=1 1 1;1 -1 1;2 -1 -1; b=1;2;1; B=A,b; rank(A),rank(B)ans = 3ans = 3 X=AbX = 0.6667 -0.50000.8333 方程組的解為：，6.11求函數(shù)

2、的拉普拉斯變換。解： syms t; ft=exp(-3*t)*sin(t); Fs=laplace(ft)Fs =1/(s + 3)2 + 1) 函數(shù)的拉普拉斯變換為：7.11單位負(fù)反饋系統(tǒng)的開環(huán)傳遞函數(shù)為應(yīng)用Simulink仿真系統(tǒng)構(gòu)建其階躍響應(yīng)曲線。解：模型仿真圖 1單位階躍響應(yīng)曲線圖 17.7用函數(shù)創(chuàng)建二階系統(tǒng)，為單位階躍信號，使用Simulink創(chuàng)建和仿真系統(tǒng)的模型。解：function sys,x0,str,ts = sfun1(t,x,u,flag)switch flag, case 0 sys,x0,str,ts=mdlInitializeSizes; case 3 sys=m

3、dlOutputs(t,x,u); case 1,2,4,9 sys=;endfunction sys,x0,str,ts=mdlInitializeSizes()sizes=simsizes;sizes.NumContStates=0;sizes.NumDiscStates=0;sizes.NumOutputs=1;sizes.NumInputs=-1;sizes.DirFeedthrough=1;sizes.NumSampleTimes=1;sys=simsizes(sizes);x0=;str=;ts=0 0;function sys=mdlOutputs(t,x,u)sys=u + (

4、exp(-u/10)*cos(39(1/2)*u)/10)/2 - (19*39(1/2)*exp(-u/10)*sin(39(1/2)*u)/10)/78 - 1/2;8.1創(chuàng)建連續(xù)二階系統(tǒng)和離散系統(tǒng)的傳遞函數(shù)模型。(1)解： num=5; den=1 2 2; sysc=tf(num,den)sysc = 5 - s2 + 2 s + 2Continuous-time transfer function.(2)解： s=tf(s); H=5/(s2+2*s+2); clear s=tf(s); sysc=5/(s2+2*s+2); sysc.inputdelay=2sysc = 5 exp

5、(-2*s) * - s2 + 2 s + 2Continuous-time transfer function.(3)解： G=tf(0.5 0,1 -1.5 0.5,-1)G = 0.5 z - z2 - 1.5 z + 0.5Sample time: unspecifiedDiscrete-time transfer function.8.2已知系統(tǒng)的傳遞函數(shù)為建立系統(tǒng)的傳遞函數(shù)模型，并轉(zhuǎn)換為零極點(diǎn)模型和狀態(tài)空間模型。解： num=2 1 den=1 0.2 1.01; G=tf(num,den)G = 2 s + 1 - s2 + 0.2 s + 1.01Continuous-time

6、 transfer function. G1=zpk(G)G1 = 2 (s+0.5) - (s2 + 0.2s + 1.01)Continuous-time zero/pole/gain model. G2=ss(G)G2 = a = x1 x2 x1 -0.2 -1.01 x2 1 0 b = u1 x1 2 x2 0 c = x1 x2 y1 1 0.5 d = u1 y1 0Continuous-time state-space model.8.4已知系統(tǒng)的方框圖如題圖所示。其中,計(jì)算系統(tǒng)的。題圖8.1系統(tǒng)的方框圖解： R1=1; R2=2; C1=3; C2=4; G1=tf(1,R

7、1); G2=tf(1,C1 0); G3=tf(1,R2); G4=tf(1,C2 0); sys1=series(G3,G4); sys2=feedback(sys1,1); sys3=series(G1,G2); sys4=feedback(sys3,1); SYS5=series(sys4,sys2); sys=feedback(sys5,tf(conv(R1,C1 0),1)sys = 1 - 24 s2 + 14 s + 1Continuous-time transfer function.9.1系統(tǒng)的傳遞函數(shù)為繪制出其根軌跡。伯德圖和奈奎斯圖。解： num=2 5 1; den=1 2 3; sys=tf(num,den); subplot(3,1,1) rlocus(sys) grid on subplot(3,1,2) bode(sys) grid on subplot(3,1,3) nyquist(sys) grid on9.5系統(tǒng)的開環(huán)傳遞函數(shù)為計(jì)算系統(tǒng)的幅值裕度和相角裕度。解： num=7 35; den=conv(conv(1 0 0,1 10),1 1); sys=tf(num,den)sys = 7 s + 35 - s4