實(shí)驗(yàn)五連續(xù)系統(tǒng)分析實(shí)驗(yàn)?zāi)康?/p>

深刻理解連續(xù)時(shí)間系統(tǒng)的系統(tǒng)函數(shù)在分析連續(xù)系統(tǒng)的時(shí)域特性、頻域特性及穩(wěn)定性中的重要作用及意義，掌握根據(jù)系統(tǒng)函數(shù)的零極點(diǎn)設(shè)計(jì)簡單的濾波器的方法。掌握利用MATLAB分析連續(xù)系統(tǒng)的時(shí)域響應(yīng)、頻響特性和零極點(diǎn)的基本方法。

實(shí)驗(yàn)原理

MATLAB提供了許多可用于分析線性時(shí)不變連續(xù)系統(tǒng)的函數(shù)，主要包含有系統(tǒng)函數(shù)、系統(tǒng)時(shí)域響應(yīng)、系統(tǒng)頻域響應(yīng)等分析函數(shù)。可見數(shù)值計(jì)算和理論計(jì)算曲線基本重合。誤差分析：可見誤差小于0.001，計(jì)算值與理論值契合度很高。3.

研究具有以下零極點(diǎn)的連續(xù)系統(tǒng)：

(a)1個(gè)極點(diǎn)s=—0.1，增益k=1。

(b)1個(gè)極點(diǎn)s=0，增益k=1。

(c)2個(gè)共軛極點(diǎn)，增益k=1。

(d)2個(gè)共軛極點(diǎn)，增益k=1。

(e)

零點(diǎn)在，極點(diǎn)在，增益k=1。

(f)

零點(diǎn)在，極點(diǎn)在，增益k=1。

完成下列任務(wù)：(1)

利用zpk和tf命令建立系統(tǒng)的系統(tǒng)函數(shù)，畫出系統(tǒng)的零極點(diǎn)圖。(2)

分析系統(tǒng)是否穩(wěn)定。若穩(wěn)定，畫出系統(tǒng)的幅頻特性曲線。(3)

畫出系統(tǒng)的沖激響應(yīng)波形。(4)

詳細(xì)列出根據(jù)零極點(diǎn)分析系統(tǒng)特性的過程。實(shí)驗(yàn)代碼：（a）%零極點(diǎn)圖subplot(3,1,1)b=[1];a=[1,0.1];z=roots(b);p=roots(a);sys=tf(b,a);pzmap(sys)%幅頻響應(yīng)subplot(3,1,2)b=[1];a=[1,0.1];[H,w]=freqs(b,a);plot(w,abs(H));xlabel('w');ylabel('幅頻響應(yīng)');%沖激響應(yīng)subplot(3,1,3)b=[1];a=[1,0.1];sys=tf(b,a);t=0:0.1:10;h=impulse(sys,t);plot(h);xlabel('t');ylabel('h(t)')（b）%零極點(diǎn)圖subplot(3,1,1)b=[1];a=[1,0];z=roots(b);p=roots(a);sys=tf(b,a);pzmap(sys)%幅頻響應(yīng)subplot(3,1,2)b=[1];a=[1,0];[H,w]=freqs(b,a);plot(w,abs(H));xlabel('w');ylabel('幅頻響應(yīng)');%沖激響應(yīng)subplot(3,1,3)b=[1];a=[1,0];sys=tf(b,a);t=0:0.1:10;h=impulse(sys,t);plot(h);xlabel('t');ylabel('h(t)')（c）%零極點(diǎn)圖subplot(3,1,1)b=[1];a=conv([1,5j],[1,-5j]);z=roots(b);p=roots(a);sys=tf(b,a);pzmap(sys)%幅頻響應(yīng)subplot(3,1,2)b=[1];a=conv([1,5j],[1,-5j]);[H,w]=freqs(b,a);plot(w,abs(H));xlabel('w');ylabel('幅頻響應(yīng)');%沖激響應(yīng)subplot(3,1,3)b=[1];a=conv([1,5j],[1,-5j]);sys=tf(b,a);t=0:0.1:10;h=impulse(sys,t);plot(h);xlabel('t');ylabel('h(t)')（d）%零極點(diǎn)圖subplot(3,1,1)b=[1];a=conv([1,0.5+5j],[1,0.5-5j]);z=roots(b);p=roots(a);sys=tf(b,a);pzmap(sys)%幅頻響應(yīng)subplot(3,1,2)b=[1];a=conv([1,0.5+5j],[1,0.5-5j]);[H,w]=freqs(b,a);plot(w,abs(H));xlabel('w');ylabel('幅頻響應(yīng)');%沖激響應(yīng)subplot(3,1,3)b=[1];a=conv([1,0.5+5j],[1,0.5-5j]);sys=tf(b,a);t=0:0.1:10;h=impulse(sys,t);plot(h);xlabel('t');ylabel('h(t)')（e）%零極點(diǎn)圖subplot(3,1,1)b=[1,-0.5];a=conv([1,0.1+5j],[1,0.1-5j]);z=roots(b);p=roots(a);sys=tf(b,a);pzmap(sys)%幅頻響應(yīng)subplot(3,1,2)b=[1,-0.5];a=conv([1,0.1+5j],[1,0.1-5j]);[H,w]=freqs(b,a);plot(w,abs(H));xlabel('w');ylabel('幅頻響應(yīng)');%沖激響應(yīng)subplot(3,1,3)b=[1,-0.5];a=conv([1,0.1+5j],[1,0.1-5j]);sys=tf(b,a);t=0:0.1:10;h=impulse(sys,t);plot(h);xlabel('t');ylabel('h(t)')（f）%零極點(diǎn)圖subplot(3,1,1)b=[1,-0.5];a=conv([1,-0.1+5j],[1,-0.1-5j]);z=roots(b);p=roots(a);sys=tf(b,a);pzmap(sys)%幅頻響應(yīng)subplot(3,1,2)b=[1,-0.5];a=conv([1,-0.1+5j],[1,-0.1-5j]);[H,w]=freqs(b,a);plot(w,abs(H));xlabel('w');ylabel('幅頻響應(yīng)');%沖激響應(yīng)subplot(3,1,3)b=[1,-0.5];a=conv([1,-0.1+5j],[1,-0.1-5j]);sys=tf(b,a);t=0:0.1:10;h=impulse(sys