1、電 子 科 技 大 學實 驗 報 告學生姓名：彭淼 學 號：2903101008 指導教師：魏芳偉一、實驗室名稱：信號與系統實驗室二、實驗項目名稱：離散系統的轉移函數，零、極點分布和模擬三、實驗原理：離散系統的時域方程為其變換域分析方法如下：系統的頻率響應為 Z域 系統的轉移函數為 分解因式 ，其中和稱為零、極點。在MATLAB中，可以用函數z,p,K=tf2zp（num,den）求得有理分式形式的系統轉移函數的零、極點，用函數zplane（z，p）繪出零、極點分布圖；也可以用函數zplane（num，den）直接繪出有理分式形式的系統轉移函數的零、極點分布圖。四、實驗目的：1、加深對離散系統

2、轉移函數、零極點概念的理解；2、根據系統轉移函數求系統零極點分布。5、實驗內容：實驗內容（一）、使用實驗仿真系統(略實驗內容（二）、MATLAB仿真六、實驗器材（設備、元器件）：計算機、MATLAB軟件。七、實驗步驟：對系統系統1、 編程實現系統的參數輸入，繪出幅度頻率響應曲線和零、極點分布圖。2、 根據系統的零極點計算系統頻率響應的幅值和相位。定義omega=0:511*pi/256和unitcirc=exp(j*omega得到在單位圓上512個等分點，在這些點上將要對頻率響應求值。（a）定義polevectors1是一個2×512的矩陣，其中每一行包含這樣一些復數，這些復數是由u

3、nitcirc的相應列減去一個極點位置得到的。如果ps1是一個列向量，它包含了極點的位置，可以用下列命令來完成polevectors1 = ones(2, 1 * unitcirc ps1 * ones(1, 512用abs和astan2定義polelength1和poleangle1作為polevectors1中每一元素的幅值和相位。（b）類似與polevectors1定義zerovectors1，使得它是2×512的矩陣，其中包含從零點位置到unitcirc元素的向量。定義zerolength1和zeroangle1分別是這些向量的幅值和相位。（c）畫出polelength1和z

4、erolength1對于omega的圖。根據這些圖，預計在哪里有最大值和最小值？（d）利用MATLAB命令H = freqz (b,a,512,whole 導出幅值和相位，與（c）中的結果相比較。8、實驗數據及結果分析：1、幅度頻率響應曲線和零、極點分布圖Matlab 程序源代碼：a=1 -0.9 0.81;b=1;figure(1;zplane(b,a;title('零、極點分布圖'figure(2;Fs=1024;freqz(b,a,Fs; grid on;title('幅度頻率響應曲線'輸出圖像:2、Matlab 程序源代碼：w=0:511*pi/256;

5、H,w=freqz(b,a,w;unitcirc=exp(j*w;ps1=roots(a;polevectors1=repmat(unitcirc,2,1-repmat(ps1,1,512'polelength1=abs(polevectors1;poleangle1=angle(polevectors1;figure,plot(w',polelength1;title('極點-幅值'figure,plot(w',poleangle1;title('極點-相位'figuresubplot(211,plot(w',polelength1(:,1.*polelength1(:,2title('極點乘積-幅值' %極點乘積最小處就是幅頻響應的峰值subplot(212,plot(w',abs(H'title('幅值響應' %極點乘積最小處就是幅頻響應的峰值%上面兩圖說明，距極點乘積最小時，幅值響應最大輸出圖像:9、實驗結論：Matlab功能很強大，能快速方便地模擬出離散系統的轉移函數，零、