實驗一信號與系統(tǒng)的傅立葉分析報告

1、實用標(biāo)準(zhǔn)文案實驗一信號與系統(tǒng)的傅立葉分析一 .實驗?zāi)康挠酶盗⑷~變換對信號和系統(tǒng)進行頻域分析。二 .實驗儀器裝有 matlab 軟件的計算機三 .實驗內(nèi)容及步驟（1）已知系統(tǒng)用下面差分方程描述：y(n) x( n) ay (n 1)試在 a0.95 和 a0.5兩種情況下用傅立葉變換分析系統(tǒng)的頻率特性。要求寫出系統(tǒng)的傳輸函數(shù)，并打印H (e jw ) w 曲線。、當(dāng) a=0.95B=1;A=1,0.95;subplot(1,3,1);zplane(B,A);xlabel('實部 Re');ylabel('虛部 Im');title('y(n)=x(n)+0

2、.95y(n-1) 傳輸函數(shù)零、極點分布');grid on；H,w=freqz(B,A,'whole');subplot(1,3,2);plot(w/pi,abs(H),'linewidth',2);grid on;xlabel('omega/pi');ylabel('|H(ejomega)|');title(' 幅頻響應(yīng)特性');axis(0,2,0,2.5);subplot(1,3,3);plot(w/pi,angle(H),'linewidth',2);grid on;xlabel(

3、'omega/pi');ylabel('phi(omega)');title(' 相頻響應(yīng)特性');axis(-0.1,2.1,-1.5,1.5);a=0.5 程序如上，圖如下精彩文檔實用標(biāo)準(zhǔn)文案（ 2）已知兩系統(tǒng)分別用下面差分方程描述：y1 ( n)x(n)x(n1)y2 (n)x(n)x(n1)H (e jw ) w試分別寫出它們的傳輸函數(shù)，并分別打印曲線。當(dāng)方程為y1 (n) x(n) x(n 1) 的程序代碼：B=1,1;A=1; subplot(2,3,1);zplane(B,A);xlabel('實部 Re');yla

4、bel('虛部 Im');title('y(n)=x(n)+x(n-1)傳輸函數(shù)零、極點分布');grid onH,w=freqz(B,A,'whole');subplot(2,3,2);plot(w/pi,abs(H),'linewidth',2);精彩文檔實用標(biāo)準(zhǔn)文案grid on;xlabel('omega/pi');ylabel('|H(ejomega)|');title(' 幅頻響應(yīng)特性');axis(0,2,0,2.2);subplot(2,3,3);plot(w/pi,

5、angle(H),'linewidth',2);grid on;xlabel('omega/pi');ylabel('phi(omega)');title(' 相頻響應(yīng)特性');axis(-0.1,2.1,-1.6,1.6);當(dāng)方程為y2 (n) x( n) x(n 1) 的程序代碼：B=1,-1;A=1;subplot(2,3,4);zplane(B,A);xlabel('實部 Re');ylabel('虛部 Im');title('y(n)=x(n)-x(n-1)傳輸函數(shù)零、極點分布&#

6、39;);grid onH,w=freqz(B,A,'whole');subplot(2,3,5);plot(w/pi,abs(H),'linewidth',2);grid on;xlabel('omega/pi');ylabel('|H(ejomega)|');title(' 幅頻響應(yīng)特性');axis(0,2,0,2.2);subplot(2,3,6);plot(w/pi,angle(H),'linewidth',2);grid on;xlabel('omega/pi');yla

7、bel('phi(omega)');title(' 相頻響應(yīng)特性');axis(-0.1,2.1,-1.6,1.6);精彩文檔實用標(biāo)準(zhǔn)文案（ 3）已知信號 x(n)R3 (n) ，試分析它的頻域特性，要求打印 X ( ejw) w 曲線。B=1,0,0,-1;A=1,-1;subplot(1,3,1);zplane(B,A);xlabel('實部 Re');ylabel('虛部 Im');title('x(n)=R3(n)傳輸函數(shù)零、極點分布');axis(-1.1,1.1,-1.5,1.5);grid onH,w

8、=freqz(B,A,'whole');subplot(1,3,2);plot(w/pi,abs(H),'linewidth',2);grid on;xlabel('omega/pi');ylabel('|H(ejomega)|');title(' 幅頻響應(yīng)特性');axis(0,2,0,1.3);精彩文檔實用標(biāo)準(zhǔn)文案subplot(1,3,3);plot(w/pi,angle(H),'linewidth',2);grid on;xlabel('omega/pi');ylabel(&

9、#39;phi(omega)');title(' 相頻響應(yīng)特性');axis(-0.1,2.1,-2.1,2.1);（ 4）假設(shè) x( n)(n) ，將 x(n) 以 2 為周期進行延拓，得到x ( n ) ，試分析它的頻率特性，并畫出它的幅頻特性。n=-20:20;x=ones(1,1),zeros(1,1);xtide=x(mod(n,2)+1);grid on;subplot(2,2,1),stem(n,xtide,'.')xlabel('n'),ylabel('xtide(n)');title(' 延拓后周

10、期為2 的周期序列 ')B=1;A=1;subplot(2,2,2);精彩文檔實用標(biāo)準(zhǔn)文案zplane(B,A);xlabel('實部 Re');ylabel('虛部 Im');title('x(n)延拓后得到x (n)傳輸函數(shù)零、極點分布');axis(-1.1,1.1,-1.1,1.1);grid onH,w=freqz(B,A,'whole');subplot(2,2,3);plot(w/pi,abs(H),'linewidth',2);grid on;xlabel('omega/pi');ylabel('|H(ejomega)|');title(' 幅頻響應(yīng)特性');axis(0,2,0,2);subp