控制系统仿真与CAD实验报告专业版

2021-10-11 13:02:58| 浏览次数：

　控制系统仿真与CAD 实验报告专业版

　资料内容仅供参考，如有不当或者侵权，请联系本人改正或者删除。

　《

　控制系统仿真与 CAD 》

　实验课程汇报

　一、

　实验教学目的与基本要求

　资料内容仅供参考，如有不当或者侵权，请联系本人改正或者删除。

　上机实验是本课程重要的实践教学环节。

　实验的目的不仅仅是验证理论知识, 更重要的是通过上机加强学生的实验手段与实践技能, 掌握应用MATLAB/ Simulink 求解控制问题的方法, 培养学生分析问题、解决问题、应用知识的能力和创新精神, 全面提高学生的综合素质。

　通过对MATLAB/ Simulink进行求解, 基本掌握常见控制问题的求解方法与命令调用, 更深入地熟悉和了解MATLAB语言的强大的计算功能与其在控制领域的应用优势。

　上机实验最终以书面汇报的形式提交, 作为期末成绩的考核内容。

　二、

　题目及解答

　第一部分: MATLAB 必备基础知识、

　控制系统模型与转换、

　线性控制系统的计算机辅助分析

　>>f=inline("[- x(2)- x(3);x(1)+ a* x(2);b+ (x(1)- c)* x(3)]","t","x","flag","a","b","c");[t,x]=ode45(f,[0,100],[0;0;0],[],0.2,0.2,5.7);plot3(x(:,1),x(:,2),x(:,3)),grid,figure,plot(x(:,1),x(:,2)),grid

　资料内容仅供参考，如有不当或者侵权，请联系本人改正或者删除。

　 >>y=@(x)x(1)^2- 2* x(1)+ x(2);ff=optimset;ff.LargeScale="off";ff.TolFun=1e- 30;ff.TolX=1e- 15;ff.TolCon=1e- 20;x0=[1;1;1];xm=[0;0;0];xM=[];A=[];B=[];Aeq=[];Beq=[];[x,f,c,d]=fmincon(y,x0,A,B,Aeq,Beq,xm,xM,@wzhfc1,ff)

　Warning: Options LargeScale = "off" and Algorithm =

　"trust- region- reflective" conflict.

　资料内容仅供参考，如有不当或者侵权，请联系本人改正或者删除。

　 Ignoring Algorithm and running active- set algorithm. To run

　trust- region- reflective, set

　LargeScale = "on". To run active- set without this warning, use

　Algorithm = "active- set".

　> In fmincon at 456

　Local minimum possible. Constraints satisfied.

　fmincon stopped because the size of the current search direction is less than

　twice the selected value of the step size tolerance and constraints are

　satisfied to within the selected value of the constraint tolerance.

　 Active inequalities (to within options.TolCon = 1e- 20):

　lower

　upper

　 ineqlin

　 ineqnonlin

　x =

　 1.0000

　f =

　资料内容仅供参考，如有不当或者侵权，请联系本人改正或者删除。

　 - 1.0000

　c =

　d =

　 iterations: 5

　funcCount: 20

　lssteplength: 1

　stepsize: 3.9638e- 26

　algorithm: "medium- scale: SQP, Quasi- Newton, line- search"

　firstorderopt: 7.4506e- 09

　constrviolation: 0

　message: [1x766 char]

　 (a) >> s=tf("s");G=(s^3+ 4* s+ 2)/ (s^3* (s^2+ 2)* ((s^2+ 1)^3+ 2* s+ 5))

　G =

　s^3 +

　4 s +

　- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -

　资料内容仅供参考，如有不当或者侵权，请联系本人改正或者删除。

　 - - - -

　s^11 +

　5 s^9 +

　9 s^7 +

　2 s^6 +

　12 s^5 +

　4 s^4 +

　12 s^3

　 Continuous- time transfer function.

　(b)

　>> z=tf("z",0.1);H=(z^2+ 0.568)/ ((z- 1)* (z^2- 0.2* z+ 0.99))

　H =

　 z^2 +

　0.568

　- - - - - - - - - - - - - - - - - - - - - - - - - - - - -

　z^3 -

　1.2 z^2 +

　1.19 z -

　0.99

　Sample time: 0.1 seconds

　Discrete- time transfer function.

　 >> A=[0 1 0;0 0 1;- 15 - 4 - 13];B=[0 0 2]";C=[1 0 0];D=0;G=ss(A,B,C,D),Gs=tf(G),Gz=zpk(G)

　G =

　a =

　资料内容仅供参考，如有不当或者侵权，请联系本人改正或者删除。

　- 15

　 - 4

　- 13

　b =

　c =

　d =

　Continuous- time state- space model.

　Gs =

　- - - - - - - - - - - - - - - - - - - - - - -

　s^3 +

　13 s^2 +

　4 s +

　 Continuous- time transfer function.

　Gz =

　- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -

　资料内容仅供参考，如有不当或者侵权，请联系本人改正或者删除。

　 (s+ 12.78) (s^2 +

　0.2212s +

　1.174)

　 Continuous- time zero/ pole/ gain model.

　设采样周期为 0.01s

　>> z=tf("z",0.01);H=(z+ 2)/ (z^2+ z+ 0.16)

　H =

　z +

　- - - - - - - - - - - - - -

　z^2 +

　z +

　0.16

　 Sample time: 0.01 seconds

　Discrete- time transfer function.

　 >> syms J Kp Ki s;G=(s+ 1)/ (J* s^2+ 2* s+ 5);Gc=(Kp* s+ Ki)/ s;GG=feedback(G* Gc,1)

　GG =

　((Ki +

　Kp* s)* (s +

　1))/ (J* s^3 +

　(Kp +

　2)* s^2 +

　(Ki +

　Kp +

　5)* s

　资料内容仅供参考，如有不当或者侵权，请联系本人改正或者删除。

　Ki)

　 (a)>>s=tf("s");G=(211.87* s+ 317.64)/ ((s+ 20)* (s+ 94.34)* (s+ 0.1684));Gc=(169.6* s+ 400)/ (s* (s+ 4));H=1/ (0.01* s+ 1);GG=feedback(G* Gc,H),Gd=ss(GG),Gz=zpk(GG)

　GG =

　359.3 s^3 +

　3.732e04 s^2 +

　1.399e05 s +

　127056

　- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -

　 0.01 s^6 +

　2.185 s^5 +

　142.1 s^4 +

　2444 s^3 +

　4.389e04 s^2 +

　1.399e05 s +

　127056

　Continuous- time transfer function.

　 Gd =

　a =

　资料内容仅供参考，如有不当或者侵权，请联系本人改正或者删除。

　- 218.5

　- 111.1

　- 29.83

　- 16.74

　- 6.671

　- 3.029

　 128

　b =

　c =

　 1.097

　 3.559

　 1.668

　0.7573

　d =

　资料内容仅供参考，如有不当或者侵权，请联系本人改正或者删除。

　 Continuous- time state- space model.

　Gz =

　 35933.152 (s+ 100) (s+ 2.358) (s+ 1.499)

　- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -

　 (s^2 +

　3.667s +

　3.501) (s^2 +

　11.73s +

　339.1) (s^2 +

　203.1s +

　1.07e04)

　Continuous- time zero/ pole/ gain model.

　 (b)设采样周期为 0.1s

　>>z=tf("z",0.1);G=(35786.7* z^2+ 108444* z^3)/ ((1+ 4* z)* (1+ 20* z)* (1+ 74.04* z));Gc=z/ (1- z);H=z/ (0.5- z);GG=feedback(G* Gc,H),Gd=ss(GG),Gz=zpk(GG)

　GG =

　- 108444 z^5 +

　1.844e04 z^4 +

　1.789e04 z^3

　 - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -

　 1.144e05 z^5 +

　2.876e04 z^4 +

　274.2 z^3 +

　782.4 z^2 +

　47.52 z +

　0.5

　资料内容仅供参考，如有不当或者侵权，请联系本人改正或者删除。

　 Sample time: 0.1 seconds

　Discrete- time transfer function.

　 Gd =

　a =

　 - 0.2515

　- 0.00959

　 - 0.1095

　- 0.05318

　- 0.01791

　0.25

　 0.125

　 0.03125

　b =

　c =

　资料内容仅供参考，如有不当或者侵权，请联系本人改正或者删除。

　 0.3996

　 0.6349

　 0.1038

　0.05043

　0.01698

　d =

　- 0.9482

　Sample time: 0.1 seconds

　Discrete- time state- space model.

　 Gz =

　- 0.94821 z^3 (z- 0.5) (z+ 0.33)

　- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -

　(z+ 0.3035) (z+ 0.04438) (z+ 0.01355) (z^2 -

　0.11z +

　0.02396)

　 Sample time: 0.1 seconds

　Discrete- time zero/ pole/ gain model.

　 >>s=tf("s");g1=1/ (s+ 1);g2=s/ (s^2+ 2);g3=1/ s^2;g4=(4* s+ 2)/

　资料内容仅供参考，如有不当或者侵权，请联系本人改正或者删除。

　 (s+ 1)^2;g5=50;g6=(s^2+ 2)/ (s^3+ 14);G1=feedback(g1* g2,g4);G2=feedback(g3,g5);GG=3* feedback(G1* G2,g6)

　GG =

　 3 s^6 +

　6 s^5 +

　3 s^4 +

　42 s^3 +

　84 s^2 +

　42 s

　 - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -

　s^10 +

　3 s^9 +

　55 s^8 +

　175 s^7 +

　300 s^6 +

　1323 s^5 +

　2656 s^4 +

　3715 s^3 +

　7732 s^2 +

　5602 s +

　1400

　Continuous- time transfer function.

　 >>s=tf("s");T0=0.01;T1=0.1;T2=1;G=(s+ 1)^2* (s^2+ 2* s+ 400)/ ((s+ 5)^2* (s^2+ 3* s+ 100)* (s^2+ 3* s+ 2500));Gd1=c2d(G,T0),Gd2=c2d(G,T1),Gd3=c2d(G,T2),step(G),figure,step(Gd1),figure,step(Gd2),figure,step(Gd3)

　 Gd1 =

　 4.716e- 05 z^5 -

　0.0001396 z^4 +

　9.596e- 05 z^3 +

　8.18e- 05

　资料内容仅供参考，如有不当或者侵权，请联系本人改正或者删除。

　 z^2 -

　0.0001289 z +

　4.355e- 05

　 - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -

　 z^6 -

　5.592 z^5 +

　13.26 z^4 -

　17.06 z^3 +

　12.58 z^2 -

　5.032 z +

　0.8521

　 Sample time: 0.01 seconds

　Discrete- time transfer function.

　Gd2 =

　 0.0003982 z^5 -

　0.0003919 z^4 -

　0.000336 z^3 +

　0.0007842 z^2 -

　0.000766 z +

　0.0003214

　 - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -

　 z^6 -

　2.644 z^5 +

　4.044 z^4 -

　3.94 z^3 +

　2.549 z^2 -

　1.056 z +

　 Sample time: 0.1 seconds

　Discrete- time transfer function.

　Gd3 =

　资料内容仅供参考，如有不当或者侵权，请联系本人改正或者删除。

　 8.625e- 05 z^5 -

　4.48e- 05 z^4 +

　6.545e- 06 z^3 +

　1.211e - 05 z^2 -

　3.299e- 06 z +

　1.011e- 07

　- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -

　z^6 -

　0.0419 z^5 -

　0.07092 z^4 -

　0.0004549 z^3 +

　0.002495 z^2 -

　3.347e- 05 z +

　1.125e- 07

　Sample time: 1 seconds

　Discrete- time transfer function.

　资料内容仅供参考，如有不当或者侵权，请联系本人改正或者删除。

　10.

　 (a)>> G=tf(1,[1 2 1 2]);eig(G),pzmap(G)

　ans =

　 - 2.0000

　 - 0.0000 +

　1.0000i

　- 0.0000 -

　1.0000i

　资料内容仅供参考，如有不当或者侵权，请联系本人改正或者删除。

　系统为临界稳定。

　(b) >> G=tf(1,[6 3 2 1 1]);eig(G),pzmap(G)

　 ans =

　 - 0.4949 +

　0.4356i

　- 0.4949 -

　0.4356i

　 0.2449 +

　0.5688i

　 0.2449 -

　0.5688i

　有一对共轭复根在右半平面, 所以系统不稳定。

　资料内容仅供参考，如有不当或者侵权，请联系本人改正或者删除。

　 ans =

　- 2.0000

　 - 1.0000

　1.0000

　有两根在右半平面, 故系统不稳定。

　11.

　 (1) >> H=tf([- 3 2],[1 - 0.2 - 0.25 0.05]);pzmap(H),abs(eig(H"))

　 ans =

　 0.5000

　资料内容仅供参考，如有不当或者侵权，请联系本人改正或者删除。

　系统稳定。

　(2) >> H=tf([3 - 0.39 - 0.09],[1 - 1.7 1.04 0.268 0.024]);pzmap(H),abs(eig(H"))

　 ans =

　 1.1939

　0.1298

　资料内容仅供参考，如有不当或者侵权，请联系本人改正或者删除。

　系统不稳定。

　12.

　 (1)>> A=[- 0.2 0.5 0 0 0;0 - 0.5 1.6 0 0;0 0 - 14.3 85.8 0;0 0 0 - 33.3 100;0 0 0 0 - 10];B=[0 0 0 0 30]";C=zeros(1,5);D=0;G=ss(A,B,C,D),eig(G)

　 G =

　a =

　 - 0.2

　0.5

　 - 0.5

　1.6

　- 14.3

　 85.8

　- 33.3

　100

　- 10

　b =

　资料内容仅供参考，如有不当或者侵权，请联系本人改正或者删除。

　 ans =

　- 0.

　 - 0.5000

　- 14.3000

　- 33.3000

　- 10.0000

　c =

　d =

　Continuous- time state- space model.

　系统稳定。

　13.

　资料内容仅供参考，如有不当或者侵权，请联系本人改正或者删除。

　>> A=[- 5 2 0 0; 0 - 4 0 0; - 3 2 - 4 - 1; - 3 2 0 - 4]; A=sym(A);syms t;x=expm(A* t)* [1;2;0;1]

　x =

　4* exp(- 4* t) -

　3* exp(- 5* t)

　2* exp(- 4* t)

　12* exp(- 4* t) -

　18* exp(- 5* t) +

　3* t* exp(- 4* t) -

　4* t^2* (exp(- 4* t)/ (4* t) +

　exp(- 4* t)/ (2* t^2)) +

　8* t^2* (exp(- 4* t)/ 2 -

　exp(- 4* t)/ (2* t)) -

　16* t* (exp(- 4* t) -

　exp(- 4* t)/ (2* t))

　6* exp(- 4* t) -

　9* exp(- 5* t) -

　8* t* (exp(- 4* t) -

　exp(- 4* t)/ (2* t))

　 >> G=ss([- 5 2 0 0; 0 - 4 0 0; - 3 2 - 4 - 1; - 3 2 0 - 4],[1;2;0;1],eye(4),zeros(4,1));tt=0:0.01:2; xx=[];

　for i=1:length(tt)

　 t=tt(i); xx=[xx eval(x)];

　end

　y=impulse(G,tt); plot(tt,xx,tt,y,":")

　资料内容仅供参考，如有不当或者侵权，请联系本人改正或者删除。

　解析解和数值解的脉冲响应曲线如图所示, 可以看出她们完全一致。

　14.

　(a) >> s=tf("s");G=(s+ 6)* (s- 6)/ (s* (s+ 3)* (s+ 4- 4j)* (s+ 4+ 4j));rlocus(G),grid

　不存在 K 使得系统稳定。

　(b) >> G=tf([1,2,2],[1 1 14 8 0]);rlocus(G),grid

　资料内容仅供参考，如有不当或者侵权，请联系本人改正或者删除。

　放大根轨迹图像, 可以看到, 根轨迹与虚轴交点处, K 值为 5.53, 所以, 0<K<5.53 时, 系统稳定。

　15.

　 pade_app.m

　function Gr=pade_app(c,r,k)

　w=- c(r+ 2:r+ k+ 1)";vv=[c(r+ 1:- 1:1)";zeros(k- 1- r,1)];

　W=rot90(hankel(c(r+ k:- 1:r+ 1),vv));V=rot90(hankel(c(r:- 1:1)));

　资料内容仅供参考，如有不当或者侵权，请联系本人改正或者删除。

　 x=[1 (W\w)"];dred=x(k+ 1:- 1:1)/ x(k+ 1);

　y=[c(1) x(2:r+ 1)* V"+ c(2:r+ 1)];nred=y(r+ 1:- 1:1)/ x(k+ 1);

　Gr=tf(nred,dred);

　paderm.m

　function [n,d]=paderm(tau,r,k)

　c(1)=1;for i=2:r+ k+ 1,c(i)=- c(i- 1)* tau/ (i- 1);end

　Gr=pade_app(c,r,k);n=Gr.num{1}(k- r+ 1:end);d=Gr.den{1};

　 >> tau=2;[n,d]=paderm(tau,1,3);s=tf("s");G=tf(n,d)* (s- 1)/ (s+ 1)^5,rlocus(G)

　 G =

　- 1.5 s^2 +

　4.5 s -

　- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -

　s^8 +

　8 s^7 +

　29.5 s^6 +

　65.5 s^5 +

　95 s^4 +

　91 s^3 +

　55.5 s^2 +

　19.5 s +

　Continuous- time transfer function.

　资料内容仅供参考，如有不当或者侵权，请联系本人改正或者删除。

　由图得 0<K<3.68 可以使得闭环系统稳定。

　16.

　 (a)>>s=tf("s");G=8* (s+ 1)/ (s^2* (s+ 15)* (s^2+ 6* s+ 10));bode(G),figure,nyquist(G),figure,nichols(G),[Gm,y,wcg,wcp]=margin(G),figure,step(feedback(G,1))

　 Gm =

　30.4686

　资料内容仅供参考，如有不当或者侵权，请联系本人改正或者删除。

　 y =

　 4.2340

　wcg =

　 1.5811

　wcp =

　 0.2336

　资料内容仅供参考，如有不当或者侵权，请联系本人改正或者删除。

　系统稳定。

　(b)>>z=tf("z");G=0.45* (z+ 1.31)* (z+ 0.054)* (z- 0.957)/ (z* (z- 1)* (z- 0.368)* (z- 0.99));bode(G),figure,nyquist(G),figure,nichols(G),[Gm,y,wcg,wcp]=margin(G),figure,step(feedback(G,1))

　Warning: The closed- loop system is unstable.

　> In warning at 26

　In DynamicSystem.margin at 63

　Gm =

　资料内容仅供参考，如有不当或者侵权，请联系本人改正或者删除。

　0.9578

　y =

　- 1.7660

　wcg =

　 1.0464

　wcp =

　 1.0734

　资料内容仅供参考，如有不当或者侵权，请联系本人改正或者删除。

　系统不稳定。

　17.

　 >>s=tf("s");G=100* (1+ s/ 2.5)/ (s* (1+ s/ 0.5)* (1+ s/ 50));Gc=1000* (s+ 1)* (s+ 2.5)/ ((s+ 0.5)* (s+ 50));GG=G* Gc;nyquist(GG),grid,figure,bode(GG),figure,nichols(GG),grid,figure,step(feedback(GG,1))

　资料内容仅供参考，如有不当或者侵权，请联系本人改正或者删除。

　由奈氏图可得, 曲线不包围(- 1,j0)点, 而开环系统不含有不稳定极点, 所以根据奈氏稳定判据闭环系统是稳定的。

　用阶跃响应来验证, 可得系统是稳定的。

　第二部分: Simulink 在系统仿真中的应用、

　控制系统计算机辅助设计、

　控制工程中的仿真技术应用

　资料内容仅供参考，如有不当或者侵权，请联系本人改正或者删除。

　>> syms y t;

　y=dsolve("D4y+ 5* D3y+ 6* D2y+ 4* Dy+ 2* y=exp(- 3* t)+ exp(- 5* t)* sin(4* t+ pi/ 3)","y(0)=1","Dy(0)=1/ 2","D2y(0)=1/ 2","D3y(0)=1/ 5");tt=0:.05:10; yy=[];

　for k=1:length(tt)

　 ti=tt(k);

　 yy=[yy subs(y,"t",ti)];

　end

　plot(tout,yout,tt,yy,":")

　资料内容仅供参考，如有不当或者侵权，请联系本人改正或者删除。

　输出曲线及误差曲线

　资料内容仅供参考，如有不当或者侵权，请联系本人改正或者删除。

　442.3( )1 22 12.3( )3 24 3sin( )sin( )xxx t x ex xx x ex x  5.

　>> [A,B,C,D]=linmod("part2_5");G=ss(A,B,C,D)

　资料内容仅供参考，如有不当或者侵权，请联系本人改正或者删除。

　 Warning: Using a default value of 0.2 for maximum step size.

　The simulation step

　size will be equal to or less than this value.

　You can disable this diagnostic by

　setting "Automatic solver parameter selection" diagnostic to "none" in the

　Diagnostics page of the configuration parameters dialog.

　> In dlinmod at 172

　In linmod at 60

　a =

　- 100

　130

　- 100

　 200

　 - 0.88

　- 100

　 294.1

　- 29.41

　- 149.3

　 100

　 - 0.44

　 - 27.56

　1.045e+ 004

　 100

　 - 10

　 x10

　资料内容仅供参考，如有不当或者侵权，请联系本人改正或者删除。

　 x10

　 1.4

　11.76

　 1.4

　 19.61

　- 6.667

　 x10

　b =

　资料内容仅供参考，如有不当或者侵权，请联系本人改正或者删除。

　 x10

　c =

　x10

　130

　d =

　Continuous- time model.

　>>subplot(221),step(G),grid,subplot(222),bode(G),grid,subplot(223),nyquist(G),grid,subplot(224),nichols(G),grid

　阶跃响应和频率响应曲线

　资料内容仅供参考，如有不当或者侵权，请联系本人改正或者删除。

　>>s=tf("s");G=210* (s+ 1.5)/ ((s+ 1.75)* (s+ 16)* (s^2+ 3* s+ 11.25));Gc=52.5* (s+ 1.5)/ (s+ 14.86);GG=feedback(G* Gc,1);step(feedback(G,1)),figure,step(GG),xlim([85 95])

　>> [Gm,garma,wcg,wcp]=margin(G)

　Gm =

　 4.8921

　garma =

　资料内容仅供参考，如有不当或者侵权，请联系本人改正或者删除。

　60.0634

　wcg =

　 7.9490

　wcp =

　 3.9199

　>> [Gm,garma,wcg,wcp]=margin(G* Gc)

　Warning: The closed- loop system is unstable.

　> In warning at 26

　In DynamicSystem.margin at 63

　Gm =

　 0.8090

　garma =

　- 6.0615

　wcg =

　17.1659

　wcp =

　18.9029

　资料内容仅供参考，如有不当或者侵权，请联系本人改正或者删除。

　>>A=[0 1 0 0;0 0 1 0;- 3 1 2 3;2 1 0 0];B=[1 0;2 1;3 2;4 3];Q=diag([1 2 3 4]);R=eye(2);[K,P]=lqr(A,B,Q,R),eig(A- B* K)

　 K =

　- 0.0978

　1.2118

　1.8767

　0.7871

　 - 3.8819

　 - 0.4668

　2.6713

　1.0320

　P =

　 5.4400

　0.6152

　 - 2.3163

　0.0452

　0.6152

　1.8354

　 - 0.0138

　 - 0.7582

　 - 2.3163

　 - 0.0138

　1.9214

　 - 0.3859

　0.0452

　 - 0.7582

　 - 0.3859

　0.8540

　ans =

　- 12.2563

　 - 1.6786 +

　0.9981i

　- 1.6786 -

　0.9981i

　- 1.4627

　资料内容仅供参考，如有不当或者侵权，请联系本人改正或者删除。

　 >> A=[- 0.2 0.5 0 0 0;0 - 0.5 1.6 0 0;0 0 - 14.3 85.8 0;0 0 0 - 33.3 100;0 0 0 0 - 10];B=[0 0 0 0 30]";C=[1 0 0 0 0];P=[- 1 - 2 - 3 - 4 - 5];K=acker(A,B,P)

　 K =

　0.0004

　 - 0.0035

　0.3946

　 - 1.4433

　>> eig(A- B* K)

　ans =

　资料内容仅供参考，如有不当或者侵权，请联系本人改正或者删除。

　 - 5.0000

　 - 4.0000

　 - 3.0000

　 - 2.0000

　 - 1.0000

　>> A=[- 0.2 0.5 0 0 0;0 - 0.5 1.6 0 0;0 0 - 14.3 85.8 0;0 0 0 - 33.3 100;0 0 0 0 - 10];B=[0 0 0 0 30]";C=[1 0 0 0 0];D=0;P=[- 1 - 2 - 3 - 4 - 5]";L=place(A",C",P)",[xh,x,t]=simobsv(ss(A,B,C,D),L);

　plot(t,x,t,xh,":"),axis([0 15 - 8000 4000])

　 L =

　1.0e+ 04 *

　- 0.0043

　0.3216

　 - 6.9200

　1.5873

　 - 0.0002

　资料内容仅供参考，如有不当或者侵权，请联系本人改正或者删除。

　optpid.m:

　function y=optpidf(x)

　opt=simset("OutputVariables","y");

　assignin("base","Kp",x(1));

　assignin("base","Ki",x(2));

　assignin("base","Kd",x(3));

　try

　 [t_time,x_state,y_out]=sim("pidctrl_model",[0,30.000000],opt);

　资料内容仅供参考，如有不当或者侵权，请联系本人改正或者删除。

　 catch, y_out=10000; end

　y=y_out(end,1);

　设置相应参数后, 点击 Optimize 后得到如下结果

　ctrl_pars =

　 0.5326

　0.2294

　0.3561

　 yy =

　 1.2446

　即比例, 积分, 微分相应系数为: Kp =0.5326, Ki = 0.2294, Kd = 0.3561 再代入模型 xiti9.mdl 进行仿真运行, 得到曲线如下图

　>>plot(tout,yout)

　资料内容仅供参考，如有不足或者侵权，请联系本人删除。

　关于方案编制、合同签订、投标文件编制的经验总结

　一、编制施工方案前需要提供的资料 1、需要提供项目的地勘报告、水文地质及地下管线分布情况。

　2、项目的平面、断面等详细图纸。

　3、项目的现场作业带及两侧为何现状、地上建筑物等实际情况，进场道路是否满足运输、就近的预制场面积、水、电是否具备条件，位置与作业带的距离是否具备运输道路（现场实际考查或提供详细的现场平面布置图）

　4、项目的工作范围，总长度，规格型号、埋深、产品上的地面为何用途。

　5、项目的工期要求。

　6、项目所在地的特殊情况加以说明。

　7、针对不同项目，甲方的注重的测重点和本项目的特点、难点。

　8、施工技术方案、模具投入是否有要求，如有要求明确，未明确视为自主安排编制。

　9、管廊用途及管廊内包含哪些管线、设施等 10、项目所需材料哪些甲供、乙供，分别说明材料的名称、数量和详细技术参数。

　11、项目预制场地是甲方提供还是我们自行解决，费用由谁承担。

　12、前期各项手续办理和征占地由谁负责办理或协助。

　二、编制施工方案注意事项：

　1、现场地上及地下实际地质、管线、建筑物、地下水对开挖、是否

　资料内容仅供参考，如有不足或者侵权，请联系本人删除。

　设支护？ 2、进场前充分做好地下管线的详细堪测，已利于做好不同地质和地下管线交错分布的开挖施工方案。

　3、工作范围、工期要求、工程特点，施工难度了解清楚。

　4、现场预制厂地具备情况。

　5、有无异形段、曲线、预留件等其它要求。

　6、此管线用途，是否同仓敷设，包含哪些管线？ 7、结合不同项目特点，合理投入、安排切实可行的施工方法。

　8、写方案时要注意我们提供内容、工作范围，不清楚的不要乱写，直接关系到报价和成本。

　三、编制投标文件注意事项：

　1、重中之重逐条查看废标条件，按此条目检查不能出现废标。

　2、编写投标文件条理清晰、组织有序，方案有针对性、方法经济有效，技术先进、严格按招标文件的要求逐项对应编写。

　3、注意加减分项，对应着每条加减分项把所能加分项材料有序尽可能提供，所有减分项要避免出现，要认真核查。

　4、注意编写的投标文件不能出现致命错误或错字，如名称、地点、特征、长度、规格等等。

　5、把我们的亮点、技术先进等突出，招标方无论有无要求，公司有的对我们有利的证书证件、业绩、奖项等要尽可能提供到标书内。

　6、开标要求到场的人员和提供的原件要一件不少，缺一即废，没有强调必带原件尽可能带着，在开标后招标方有可能要看原件。

　资料内容仅供参考，如有不足或者侵权，请联系本人删除。

　 7、注意标书内按要求签字和盖章、不能出现遗漏和侥幸，不需要签字、盖章的也不要多此一举。

　8、封标要领悟招标方要求，决不允许出现不按要求封标，否则可能以不按要求即废。

　9、报价不能出现普通材料、机械台班费、单价远远高于市场单价。

　10、预算要全面、合理、不能出现漏项、少算量、漏取费。

　11、增加单价时注意，不能直接增加数量增加费用，人工及材料、机械单价调增不能明显太高。

　12、注招标文件里要求必须提供的，即使不全，没有，一项不提供此块就是零分。

　13、投标保证金提前提供、原件备查、按要求附标书中。

　14、图纸设计和结构配筋人员在设计图纸前要全面了解清楚，严紧、不能出现少、漏项、设计错误直接影响造价。

　四、谈判、签属合同前注意事项：

　1、甲乙双方的责任合同中一定要分清。

　2、是否有预付款、进度款、结算付款方式明确，结算条件、时间明确，尽量以对我们有利的方式签定。

　3、质量标准明确，是否有质保、质保期限、质保金比例，返还方式、期限合同中要明确。

　4、如出现违约，违约责任明确，出现争议的方法，最终仲裁或诉讼选利与我公司的地方的机构。

　5、明确不可抗力，不可预见、出现不可抗力或不可预见处理方法及

　资料内容仅供参考，如有不足或者侵权，请联系本人删除。

　双方履行的责任。

　6、双方的约定保密义务。

　7、合同的终止和中止条件及各方承担的责任。

　8、合同的变更要在合同中约定，合同量的增减和施工中的洽商、单价及取费如何调整在合同中明确。

　9、设计变更、地质条件变化、规格的变更等定价及取费，在合同中明确。

　10、分别因双方原因出现停、窝工等，如何承担责任及解决办法。

　11、合同中工程名称、量、地点、规格型号、价格、地质条件、工期、工程内容等不能出现错误，合同中明确列明。

　12、合同中明确相关交接人员。

　13、付款是否提供发票、什么样的发票、提供条件。

　14、工程竣工验收条件和合格标准，时间。

　15、其他未明确事项的协商解决。

　16、非法人签约合同，要附授权委托书，防止出现诈骗。

　项目实施方案编制要点

　项目简介包括；项目提出的主要理由，项目目标、建设内容、规模、方案等，项目总投资、资金构成及筹措方案，主要技术经济指标，承担单位基本情况。

　项目实施的意义和必要性产业国内外发展现状及存在的主要问题、拟解决的关键问题、拟达到的技术水平和市场前景。

　资料内容仅供参考，如有不足或者侵权，请联系本人删除。

　项目承担单位的基本情况注册资金及股东构成，企业组织架构近 3 年经营业绩和成长性（包括总产值、销售收入、利润上缴税金、研发投入、资产负债率、银行信用等级等、既要有历年数据，也要有增长分析）

　核心竞争力分析（企业的技术开发能力、经营管理水平、成本控制能力、资金筹措能力、市场开拓能力等）

　企业发展战略和规划项目的主要建设内容及预期目标产品技术研发及产业化等重点环节的项目内容。项目拟采用的工艺路线与技术特点、设备选型及主要技术经济指标等。项目在科技研发投入、技术水平、规模及创新团队建设等方面的预期目标，与方案的目标做好相应结合。

　项目实施条件

　环境保护、资源综合利用、节能措施、外部配套条件落实情况、项目实施基础等。

　项目投资及资金筹措项目总投资及测算依据，并按研发、产业化等环节进行分项说明。资金筹措方案、投资使用方案和年度投资计划。

　项目效益分析项目总投资经济社会效益分析（企业内生发展能力提升、对地方战略性新兴产业培育和产业结构调整升级的作用）。

　资料内容仅供参考，如有不足或者侵权，请联系本人删除。

　项目实施和考核目标

　项目组织实施方式、年度考核指标有关附件技术来源及技术先进性的有关证明文件项目备案、银行贷款、自筹资金证明、环境保护、节能审查、土地等项目实施条件支撑材料与项目相关的其他证明材料或文件等签订施工合同的十大注意事项

　建设工程施工合同是依法保护发承包双方权益的法律文件，是发承包双方在工程施工过程中的最高行为准则。为防范合同纠纷，在签订过程中，以下 10 个方面需要注意。

　一关于发包人与承包人 1、发包方需要满足两方面内容：①主体资格，即建设相关手续是否齐全。例：建设用地是否已经批准？是否列入投资计划？规划、设计是否得到批准？是否进行了招标等。②履约能力即资金问题。施工所需资金是否已经落实或可能落实等。

　2、承包方需要满足的内容有：①资质情况；②施工能力；③社会信誉；④财务情况。承包方的二级公司和工程处不能对外签订合同。上述内容是体现履约能力的指标，应认真的分析和判断。

　二合同价款应注意 1、“合同价款”的填写，应依据住建部规定，招标工程的合同价款由发包人、承包人依据中标通知书中的中标价格在协议书内约定。

　资料内容仅供参考，如有不足或者侵权，请联系本人删除。

　非招标工程合同价款由发包人承包人依据工程预算在协议书内约定。

　2、合同价款是双方共同约定的条款，要求第一要协议，第二要确定。暂定价、暂估价、概算价、都不能作为合同价款，约而不定的造价不能作为合同价款。

　三发包人工作与承包人工作条款应注意 1、双方各自工作的具体时间要填写准确。

　2、双方所做工作的具体内容和要求应填写详细。

　3、双方不按约定完成有关工作应赔偿对方损失的范围、具体责任和计算方法要填写清楚。

　四合同价款及调整价款应注意 1、填写合同价款及调整时应按《通用条款》所列的固定价格、可调价格、成本加酬金三种方式，约定一种写入本款。

　2、采用固定价格应注意明确包死价的种类。如：总价包死、单价包死，还是部分总价包死，以免履约过程中发生争议。

　3、采用固定价格必须把风险范围约定清楚。

　4、应当把风险费用的计算方法约定清楚。双方应约定一个百分比系数，也可采用绝对值法。

　5、对于风险范围以外的风险费用，应约定调整方法。

　五工程预付款条款应注意 1、填写约定工程预付款的额度应结合工程款、建设工期及包工包料情况来计算。

　资料内容仅供参考，如有不足或者侵权，请联系本人删除。

　 2、应准确填写发包人向承包人拨付款项的具体时间或相对时间。

　3、应填写约定扣回工程款的时间和比例。

　六工程进度款条款应注意 1、工程进度款的拨付应以发包方代表确认的已完工程量，相应的单价及有关计价依据计算。

　2、工程进度款的支付时间与支付方式以形象进度可选择：按月结算、分段结算、竣工后一次结算（小工程）及其它结算方式。

　七材料设备供应条款应注意 1、应详细填写材料设备供应的具体内容、品种、规格、数量、单价、质量等级、提供的时间和地点。

　2、应约定供应方承担的具体责任。

　3、双方应约定供应材料和设备的结算方法（可以选择预结法、现结法、后结法或其他方法）。

　八违约条款应注意 1、在合同中首先应约定发包人对《通用条款》中的预付款、工程进度款、竣工结算的违约应承担的具体违约责任。

　2、在合同中应约定承包人对《通用条款》违约应承担的具体违约责任。

　3、还应约定其它违约责任。

　4、违约金与赔偿金应约定具体数额和具体计算方法，要越具体越好，具有可操作性，以防止事后产生争议。

　九争议与工程分包条款应注意

　资料内容仅供参考，如有不足或者侵权，请联系本人删除。

　 1、争议的解决方式是选择仲裁方式，还是选择诉讼方式，双方应达成一致意见。

　2、如果选择仲裁方式，当事人可以自主选择仲裁机构。仲裁不受级别地域管辖限制。

　3、如果选择诉讼方式，应当选定有管辖权的人民法院（诉讼是地域管辖）。

　4、分包的工程项目须经发包人同意，禁止分包单位将其承包的工程再分包。

　十关于补充条款 1、需要补充新条款或哪条、哪款需要细化、补充或修改，可在《补充条款》内尽量补充，依次按顺序排列。

　2、补充条款必须符合国家、现行的法律、法规，另行签订的有关书面协议应与主体合同精神相一致。要杜绝“阴阳合同”。

　投标文件编制注意事项

　一、封面 1.封面格式是否与招标文件要求格式一致，文字打印是否有错字。

　2.封面标段、内容是否与所投标段、内容一致。

　3.企业法人或委托代理人是否按照规定签字或盖章，是否按规定加盖单位公章，投标单位名称是否与资格审查时的单位名称相符。

　4.投标日期是否正确。

　5.密封时检查好封套上格式是否与招标文件文件一致、招标人名称、项目名称、标段号、日期是否正确。

　6.封标时仔细核对招标文件看是否需要盖密封章或者法人章。看招

　资料内容仅供参考，如有不足或者侵权，请联系本人删除。

　标文件对封套是否有特殊要求如：正副本是要求写在什么位置。

　二、目录 1.目录编号、页码、标题是否与内容编号、页码（内容首页）、标题一致。

　三、投标书及投标书附录

　 1.投标书格式、标段、内容是否与招标文件规定相符，建设单位名称与招标单位名称是否正确。

　2.报价金额是否与“投标报价汇总表合计”、“投标报价汇总表”、“综合报价表”一致，大小写是否一致。

　3.投标书所示工期是否满足招标文件要求。

　4.投标书是否按已按要求盖公章。

　5.投标书日期是否正确，是否与封面所示吻合。

　6.投标保证金是否与招标文件要求一致。

　7.投标函附录中缺陷责任期是否满足招标文件要求。

　四、授权书、银行保函、信贷证明 1.授权书、银行保函、信贷证明是否按照招标文件要求格式填写。

　2.上述三项是否由法人正确签字或盖章。

　 3.委托代理人是否正确签字或盖章。

　 4.委托书日期是否正确。

　 5.委托权限是否满足招标文件要求，单位公章加盖完善。

　五、对招标文件及合同条款的确认和承诺

　 1.投标书承诺与招标文件要求是否吻合。

　2.承诺内容与投标书其他有关内容是否一致。

　3.对招标文件（含补遗书）及合同条款的确认和承诺，是否确认了全部内容和全部条款，不能只确认、承诺主要条款，用词要确切，不允许有保留或留有其他余地。

　资料内容仅供参考，如有不足或者侵权，请联系本人删除。

　 4.投标承诺检查招标人名称和投标日期是否准确。

　5.仔细阅读招标文件的要求看是否要求删除法定代表人或其委托代理人中的一个。

　6.检查好标书中是否有需要盖法人章或授权委托人私章的。

　六、施工组织及施工进度安排 1.工程概况是否准确描述。

　2.计划开竣工日期是否符合招标文件中工期安排与规定，分项工程的阶段工期、节点工期是否满足招标文件规定。工期提前要合理，要有相应措施，不能提前的决不提前，如铺架工程工期。

　3.工期的文字叙述、施工顺序安排与“形象进度图”、“横道图”、“网络图”是否一致，特别是铺架工程工期要针对具体情况仔细安排，以免造成与实际情况不符的现象。

　 4.总体部署：施工队伍及主要负责人与资审方案是否一致，文字叙述与“平面图”、“组织机构框图”、“人员简历”及职务等是否吻合。

　5.施工方案与施工方法、工艺是否匹配。

　6.检查项目清单中有的项是否在施工组织设计中有。

　7.施工横道图有工程名称的检查名称是否是该项目的名称。

　8.检查好施工组织设计中项目名称是否是本次投标工程的项目名称。

　七、项目管理机构 1、项目经理为建造师的，是否按要求将证明材料附齐全，特别是不要忘记附安考证。

　2、技术负责人根据招标文件要求附职称证 3、其他人员是否根据招标文件要求附岗位证或职称证八、企业有关资质、社会信誉

　资料内容仅供参考，如有不足或者侵权，请联系本人删除。

　 1.营业执照、资质证书、法人代表、安全资格是否齐全并满足招标文件要求。

　2.是否按要求放锁人表和入库证明 3.企业近年来从事过的类似工程主要业绩是否满足招标文件要求。

　4.财务状况表、近年财务决算表及审

　5.仔细阅读招标文件看是否所要求的资料都已经附齐。

推荐访问:专业版控制系统仿真

上一篇：金强驿小学禁毒教育自查报告(定)

下一篇：漂流项目规划设计书

控制系统仿真与CAD实验报告专业版

相关推荐

热门排行Top Ranking