控制系統仿真與CAD 實驗報告專業版
資料內容僅供參考，如有不當或者侵權，請聯系本人改正或者刪除。
《
控制系統仿真與 CAD 》
實驗課程匯報
一、
實驗教學目的與基本要求
資料內容僅供參考，如有不當或者侵權，請聯系本人改正或者刪除。
上機實驗是本課程重要的實踐教學環節。
實驗的目的不僅僅是驗證理論知識, 更重要的是通過上機加強學生的實驗手段與實踐技能, 掌握應用MATLAB/ Simulink 求解控制問題的方法, 培養學生分析問題、 解決問題、 應用知識的能力和創新精神, 全面提高學生的綜合素質。
通過對MATLAB/ Simulink進行求解, 基本掌握常見控制問題的求解方法與命令調用, 更深入地熟悉和了解MATLAB語言的強大的計算功能與其在控制領域的應用優勢。
上機實驗最終以書面匯報的形式提交, 作為期末成績的考核內容。
二、
題目及解答
第一部分: MATLAB 必備基礎知識、
控制系統模型與轉換、
線性控制系統的計算機輔助分析
1.
>>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
資料內容僅供參考，如有不當或者侵權，請聯系本人改正或者刪除。
2.
>>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.
<stopping criteria details>
Active inequalities (to within options.TolCon = 1e- 20):
lower
upper
ineqlin
ineqnonlin
2
x =
1.0000
0
1.0000
f =
資料內容僅供參考，如有不當或者侵權，請聯系本人改正或者刪除。
- 1.0000
c =
4
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]
3.
(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 +
2
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
資料內容僅供參考，如有不當或者侵權，請聯系本人改正或者刪除。
- - - -
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.
4.
>> 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 =
x1
x2
x3
x1
0
1
0
x2
0
0
1
資料內容僅供參考，如有不當或者侵權，請聯系本人改正或者刪除。
x3
- 15
- 4
- 13
b =
u1
x1
0
x2
0
x3
2
c =
x1
x2
x3
y1
1
0
0
d =
u1
y1
0
Continuous- time state- space model.
Gs =
2
- - - - - - - - - - - - - - - - - - - - - - -
s^3 +
13 s^2 +
4 s +
15
Continuous- time transfer function.
Gz =
2
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
資料內容僅供參考，如有不當或者侵權，請聯系本人改正或者刪除。
(s+ 12.78) (s^2 +
0.2212s +
1.174)
Continuous- time zero/ pole/ gain model.
5.
設采樣周期為 0.01s
>> z=tf("z",0.01);H=(z+ 2)/ (z^2+ z+ 0.16)
H =
z +
2
- - - - - - - - - - - - - -
z^2 +
z +
0.16
Sample time: 0.01 seconds
Discrete- time transfer function.
6.
>> 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)
7.
(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 =
x1
x2
x3
x4
x5
x6
資料內容僅供參考，如有不當或者侵權，請聯系本人改正或者刪除。
x1
- 218.5
- 111.1
- 29.83
- 16.74
- 6.671
- 3.029
x2
128
0
0
0
0
0
x3
0
64
0
0
0
0
x4
0
0
32
0
0
0
x5
0
0
0
8
0
0
x6
0
0
0
0
2
0
b =
u1
x1
4
x2
0
x3
0
x4
0
x5
0
x6
0
c =
x1
x2
x3
x4
x5
x6
y1
0
0
1.097
3.559
1.668
0.7573
d =
u1
y1
0
資料內容僅供參考，如有不當或者侵權，請聯系本人改正或者刪除。
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 =
x1
x2
x3
x4
x5
x1
- 0.2515
- 0.00959
- 0.1095
- 0.05318
- 0.01791
x2
0.25
0
0
0
0
x3
0
0.25
0
0
0
x4
0
0
0.125
0
0
x5
0
0
0
0.03125
0
b =
u1
x1
1
x2
0
x3
0
x4
0
x5
0
c =
x1
x2
x3
x4
x5
資料內容僅供參考，如有不當或者侵權，請聯系本人改正或者刪除。
y1
0.3996
0.6349
0.1038
0.05043
0.01698
d =
u1
y1
- 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.
8.
>>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.
9.
>>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 +
0.
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
有一對共軛復根在右半平面, 所以系統不穩定。
(c) >> G=tf(1,[1 1 - 3 - 1 2]);eig(G),pzmap(G)
資料內容僅供參考，如有不當或者侵權，請聯系本人改正或者刪除。
ans =
- 2.0000
- 1.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
0.5000
0.
資料內容僅供參考，如有不當或者侵權，請聯系本人改正或者刪除。
系統穩定。
(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
1.1939
0.1298
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 =
x1
x2
x3
x4
x5
x1
- 0.2
0.5
0
0
0
x2
0
- 0.5
1.6
0
0
x3
0
0
- 14.3
85.8
0
x4
0
0
0
- 33.3
100
x5
0
0
0
0
- 10
b =
u1
資料內容僅供參考，如有不當或者侵權，請聯系本人改正或者刪除。
x1
0
ans =
- 0.
- 0.5000
- 14.3000
- 33.3000
- 10.0000
x2
0
x3
0
x4
0
x5
30
c =
x1
x2
x3
x4
x5
y1
0
0
0
0
0
d =
u1
y1
0
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 -
3
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
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 +
3
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 在系統仿真中的應用、
控制系統計算機輔助設計、
控制工程中的仿真技術應用
2.
資料內容僅供參考，如有不當或者侵權，請聯系本人改正或者刪除。
>> 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,":")
3.
資料內容僅供參考，如有不當或者侵權，請聯系本人改正或者刪除。
輸出曲線及誤差曲線
4.
資料內容僅供參考，如有不當或者侵權，請聯系本人改正或者刪除。
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 =
x1
x2
x3
x4
x5
x6
x1
0
0
0
0
0
0
x2
0
- 100
0
0
0
0
x3
130
0
- 100
0
0
0
x4
0
200
- 0.88
- 100
0
0
x5
0
0
0
0
- 100
0
x6
0
0
0
294.1
- 29.41
- 149.3
x7
0
100
- 0.44
0
0
0
x8
- 27.56
0
0
0
0
1.045e+ 004
x9
0
0
0
100
- 10
0
x10
0
0
0
0
0
0
資料內容僅供參考，如有不當或者侵權，請聯系本人改正或者刪除。
x7
x8
x9
x10
x1
0
1.4
0
0
x2
0
0
0
0
x3
0
0
0
0
x4
11.76
0
0
0
x5
0
1.4
0
0
x6
0
0
19.61
0
x7
0
0
0
0
x8
0
- 6.667
0
0
x9
0
0
0
0
x10
0
0
0
0
b =
u1
x1
0
x2
1
x3
0
x4
0
x5
0
x6
0
x7
0
x8
0
資料內容僅供參考，如有不當或者侵權，請聯系本人改正或者刪除。
x9
0
x10
0
c =
x1
x2
x3
x4
x5
x6
x7
x8
x9
x10
y1
130
0
0
0
0
0
0
0
0
0
d =
u1
y1
0
Continuous- time model.
>>subplot(221),step(G),grid,subplot(222),bode(G),grid,subplot(223),nyquist(G),grid,subplot(224),nichols(G),grid
階躍響應和頻率響應曲線
6.
資料內容僅供參考，如有不當或者侵權，請聯系本人改正或者刪除。
>>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
7.
資料內容僅供參考，如有不當或者侵權，請聯系本人改正或者刪除。
>>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
資料內容僅供參考，如有不當或者侵權，請聯系本人改正或者刪除。
8.
>> 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.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
資料內容僅供參考，如有不當或者侵權，請聯系本人改正或者刪除。
9.
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.封標時仔細核對招標文件看是否需要蓋密封章或者法人章?？凑?/p>
資料內容僅供參考，如有不足或者侵權，請聯系本人刪除。
標文件對封套是否有特殊要求如：正副本是要求寫在什么位置。
二、目錄 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.仔細閱讀招標文件看是否所要求的資料都已經附齊。
推薦訪問： 專業版 控制系統 仿真上一篇：操作系統文件管理實驗報告
下一篇：金強驛小學禁毒教育自查報告(定)
在偉大祖國73華誕之際，我參加了單位組織的“光影鑄魂”主題黨日活動，集中觀看了抗美援朝題材影片《長津湖》，再一次重溫這段悲壯歷史，再一次深刻感悟偉大抗美援朝精神。1950年10月，新中國剛剛成立一年，
根據省局黨組《關于舉辦習近平談治國理政（第四卷）讀書班的通知》要求，我中心通過專題學習、專題研討以及交流分享等形式，系統的對《習近平談治國理政》（第四卷)進行了深入的學習與交流，下面我就來談一談我個人
《習近平談治國理政》（第四卷）是在百年變局和世紀疫情相互疊加的大背景下，對以習近平同志為核心的黨中央治國理政重大戰略部署、重大理論創造、重大思想引領的系統呈現。它生動記錄了新一代黨中央領導集體統籌兩個
《真抓實干做好新發展階段“三農工作”》是《習近平談治國理政》第四卷中的文章，這是習近平總書記在2020年12月28日中央農村工作會議上的集體學習時的講話。文章指出，我常講，領導干部要胸懷黨和國家工作大
在《習近平談治國理政》第四卷中，習近平總書記強調，江山就是人民，人民就是江山，打江山、守江山，守的是人民的心。從嘉興南湖中駛出的小小紅船，到世界上最大的執政黨，在中國共產黨的字典里，“人民”一詞從來都
黨的十八大以來，習近平總書記以馬克思主義戰略家的博大胸襟和深謀遠慮，在治國理政和推動全球治理中牢固樹立戰略意識，在不同場合多次圍繞戰略策略的重要性，戰略和策略的關系，提高戰略思維、堅定戰略自信、強化戰
《習近平談治國理政》第四卷集中展示了以習近平同志為核心的黨中央在百年變局和世紀疫情相互疊加背景下，如何更好地堅持和發展中國特色社會主義而進行的生動實踐與理論探索；對于新時代堅持和發展什么樣的中國特色社
在黨組織的關懷下，我有幸參加了區委組織部組織的入黨積極分子培訓班。為期一周的學習，學習形式多樣，課程內容豐富，各位專家的講解細致精彩，對于我加深對黨的創新理論的認識、對黨的歷史的深入了解、對中共黨員的
《習近平談治國理政》第四卷《共建網上美好精神家園》一文中指出：網絡玩命是新形勢下社會文明的重要內容，是建設網絡強國的重要領域。截至2021年12月，我國網民規模達10 32億，較2020年12月增長4
剛剛召開的中國共產黨第十九屆中央委員會第七次全體會議上討論并通過了黨的十九屆中央委員會向中國共產黨第二十次全國代表大會的報告、黨的十九屆中央紀律檢查委員會向中國共產黨第二十次全國代表大會的工作報告和《