摘要:,已達到對離散與連續系統的特性研究。
關鍵詞:連續系統,離散系統,衝激響應,matlab,零極點。
連續系統函式零極點與系統特性研究
連續時間系統的穩定性與系統零點無關,與系統的極點有關,而系統零點則影響系統單位衝激響應的幅度和相位。理解系統的零極點與系統的穩定性之間的關係有利於對系統的理解。
如果給定系統函式h(s),或給定系統微分方程(可以求出系統函式),通過系統函式可以零極點圖判斷系統的穩定性。
(1)可用matlab函式pzmap來畫出系統的零極點圖。函式pzmap的呼叫形式為
[p,z] = pzmap(sys)
其中呼叫變數sys為系統函式,而sys生成可以利用sys=tf(num,den),num表示n(s),den表示d(s)。返回變數p存放系統h(s)的極點,返回變數z存放系統h(s)的零點。
(2)可用matlab函式impulse來畫出系統的單位衝激響應h(t)。函式impulse的呼叫形式為
h=impulse(num,den,t);
其中呼叫變數num表示n(s),den表示d(s)。返回變數h存放系統的單位衝激響應h(t)。
從h(t)的圖形可以基本判斷系統穩定性和零極點的關係。
以系統為例進行研究:
畫出系統的零極點,並畫出系統單位衝激響應h(t)的波形圖。並與理論圖形相比較
理論分析:
由h(s)可知道原系統方程為y』』(t)-6y』(t)+5y(t)=x』(t)+x(t)則可求系統衝激響應如下:
h』』(t)-6h』(t)+5h(t)=0--------r^2-6r+5=0r1=5;r2=1;
h(t)=(ae^(5t)+be^(t))*u(t)
帶入h』』(t)-6h』(t)+5h(t)=&』(t)+&(t)中,有左右相等解得a=5/8,b=3/8.
所以解得h(t)=(5/8*e^(5t)+3/8*e^(t))u(t);
程式如下:
t=0:0.02:30;
a=0.625;b=0.375;c=5;d=1;
xt=a*exp(c*t)+b*exp(d*t);
plot(t,xt)
xlabel('time(s)');
title('impulse respone')
結果如下:
通過matlab 直接畫出衝擊響應:
程式如下:
num=[1 1];
den=[1 -6 5];
sys=tf(num,den);
figure(1);pzmap(sys);
t=0:0.02:30;
h=impulse(num,den,t);
figure(2);
plot(t,h);
xlabel('time(s)');
title('impulse respone')
系統零極點圖如下
系統衝級響應圖如下:
經與理論值比較,影象符合的很好。
注意觀察零極點和系統單位衝激響應h(t)的波形圖走向關係,大致判斷系統零極點和系統穩定性之間的關係。
注意到零點在系統虛軸左邊,極點在虛軸右邊系統是不穩定的。具體**在第三部如下。
3. 只改變零點或改變極點,觀察系統單位衝激響應h(t)的波形圖,得出你的結論。
1)先**極點:由上圖可知極點為t=1和t=5 圖形上公升,可知系統不穩定。
將極點改為t=-1和t=-5程式如下:
num=[1 1];
den=[1 6 5];
sys=tf(num,den);
figure(1);pzmap(sys);
t=0:0.02:30;
h=impulse(num,den,t);
figure(2);
plot(t,h);
xlabel('time(s)');
title('impulse respone')結果如下:
將極點改為t=-1,t=2程式如下:
num=[1 1];
den=[1 -1 -2];
sys=tf(num,den);
figure(1);pzmap(sys);
t=0:0.02:30;
h=impulse(num,den,t);
figure(2);
plot(t,h);
xlabel('time(s)');
title('impulse respone')
可以看到,當虛軸右邊有極點時系統不穩定,當極點都在虛軸左邊時系統是穩定的下面**一下當極點在虛軸上時的穩定性:
極點為t=0 ,t=-1;程式如下:
num=[1 1];
den=[1 1 0];
sys=tf(num,den);
figure(1);pzmap(sys);
t=0:0.02:30;
h=impulse(num,den,t);
figure(2);
plot(t,h);
xlabel('time(s)');
title('impulse respone')
可以看到系統是不穩定的。
結論:只有當最右邊的極點在虛軸的左邊時系統才是穩定的,否則系統是不穩定的。
2)再**零點:由以上知零點為t=-1時系統是不穩定的在**零點為t=1時能否改變系統的穩定性程式如下:
num=[1 -1];
den=[1 -6 5];
sys=tf(num,den);
figure(1);pzmap(sys);
t=0:0.02:30;
h=impulse(num,den,t);
figure(2);
plot(t,h);
xlabel('time(s)');
title('impulse respone')
結果如下:
零點位置改變並未改變系統的不穩定性,接下來**改變零點是否改變穩定系統的穩定性:
程式如下:
num=[1 1];
den=[1 6 5];
sys=tf(num,den);
figure(1);pzmap(sys);
t=0:0.02:30;
h=impulse(num,den,t);
figure(2);
plot(t,h);
xlabel('time(s)');
title('impulse respone')
結果如下:
當零點改為t=+1時程式如下:
num=[1 -1];
den=[1 6 5];
sys=tf(num,den);
figure(1);pzmap(sys);
t=0:0.02:30;
h=impulse(num,den,t);
figure(2);
plot(t,h);
xlabel('time(s)');
title('impulse respone')
結果如下:
可知系統依然穩定,由此可知零點位置的改變並不會影響系統的穩定性。
結論:零點位置不影響系統穩定性,極點位置影響系統的穩定性,當所有極點都在虛軸左邊時系統是穩定的,否則系統是不穩定的。
離散系統函式零極點與系統特性研究
離散時間系統的穩定性與系統零點無關,與系統的極點有關,而系統零點則影響系統單位脈衝響應的幅度和相位。理解系統的零極點與系統的穩定性之間的關係有利於對系統的理解。
如果給定系統函式h(s),或給定系統微分方程(可以求出系統函式),通過系統函式可以零極點圖判斷系統的穩定性。
(1)可用matlab函式zplane來畫出系統的零極點圖。函式zplane的呼叫形式為
[p,z] = pzmap(num,den)
其中,如果系統函式是,呼叫變數num表示n(s),den表示d(s)。返回變數p存放系統h(s)的極點,返回變數z存放系統h(s)的零點。
(2)可用matlab函式impz來畫出系統的單位脈衝響應h[k]。函式impz的呼叫形式為
h=impz(num,den,fs);
其中呼叫變數num表示n(s),den表示d(s),抽樣率為1/fs。返回變數h存放系統單位脈衝響應h[k]。從h[k]的圖形可以基本判斷系統穩定性和零極點的關係。
以系統為例研究:
畫出系統的零極點,並畫出系統單位脈衝響應h[k]的波形圖。並與理論圖形相比較。
理論分析:
y[k-3]+0.1628y[k-2]+0.3403y[k-1]+0.0149[k]=x[k-3]-3x[k-2]+3x[k-1]-x[k]
h[k-3]+0.1628h[k-2]+0.3403h[k-1]+0.0149h[k]=0---
--------r^(-3)+0.1628r^(-2)+0.3403r^(-1)+0.0149=0
直接求衝擊響應:
程式如下:
num=[1 -3 3 -1];
den=[1 0.1682 0.3403 0.0149 ];
[r,p,k]=residue(num,den)
執行結果如下:
r = 0.1359 - 2.6563i
0.1359 + 2.6563i
-3.4399
p = -0.0618 + 0.5753i
-0.0618 - 0.5753i
-0.0445
k = 1
可解得h[k]
程式如下:
num=[1 -3 3 -1];
den=[1 0.1628 0.3403 0.0149];
figure(1);
zplane(num,den);
h=impz(num,den,21);
figure(2);
stem(0:20,h);
xlabel('k');
零點圖如下:
衝激響應圖如下:
與理論圖符合的很好。
注意觀察零極點和系統單位脈衝響應h[k]的波形圖走向關係,大致判斷系統零極點和系統穩定性之間的關係。
由圖可以看出極點在單位圓內,零點在單位圓上,由衝擊響應知該系統是穩定的具體**如下。
只改變零點或改變極點,觀察系統單位脈衝響應h[k]的波形圖,得出結論。
先**極點:
程式如下:
num=[1 -3 3 -1];
den=[1 -3 -4 12];
figure(1);
zplane(num,den);
h=impz(num,den,21);
figure(2);
stem(0:20,h);
xlabel('k');
結果如下:
可以看出系統不穩定。
程式如下:
num=[1 -3 3 -1];
den=[1 -4 -0.5 2];
figure(1);
zplane(num,den);
方程的根與函式零點
方程的根與函式零點習題 含答案 一 單選題 1 已知函式與,則函式在區間上所有零點的和為 a b c d 2 已知函式,若函式恰好有兩個零點,則實數等於 為自然對數的底數 a b c d 3 已知定義在上的函式,若有兩個零點,則實數的取值範圍是 a b c d 4 已知點是曲線上任意一點,記直線 為...
方程的根與函式的零點
方程的根與函式的零點 導學案 班級小組姓名評價 學習目標 1 知識與技能 通過二次函式的影象與x軸的交點的個數判斷一元二次方程的根的個數,了解函式零點的概念,領會方程的根與函式零點之間的關係,掌握函式零點存在性判定定理,能結合影象求解零點問題 2 過程與方法 運用 由特殊到一般 的認識規律 3 情感...
方程的根與函式的零點
課題 方程的根與函式的零點 課型 新授課 複習引入,創設情境 我們學習了一元二次方程和二次函式,那麼方程的根與其對應的函式的圖象它們又有怎樣的關係呢?新知 1 函式零點的 1 請同學完成下面的 讓學生檢查是否畫正確了,並觀察可得出以下兩點 方程根 不相等 的個數就是函式圖象與軸交點的個數 方程的實數...