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

1、實驗一 信號與系統(tǒng)的傅立葉分析一. 實驗目的 用傅立葉變換對信號和系統(tǒng)進行頻域分析。二.實驗儀器裝有matlab軟件的計算機三.實驗內(nèi)容及步驟（1）已知系統(tǒng)用下面差分方程描述：試在和兩種情況下用傅立葉變換分析系統(tǒng)的頻率特性。要求寫出系統(tǒng)的傳輸函數(shù)，并打印曲線。、當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.95y(n-1)傳輸函數(shù)零、極點分布');grid on；H,w=freqz(B,A,'

2、;whole'); subplot(1,3,2);plot(w/pi,abs(H),'linewidth',2);grid on; xlabel('omega/pi');ylabel('|H(ejomega)|');title('幅頻響應特性');axis(0,2,0,2.5);subplot(1,3,3);plot(w/pi,angle(H),'linewidth',2);grid on; xlabel('omega/pi');ylabel('phi(omega)');ti

3、tle('相頻響應特性');axis(-0.1,2.1,-1.5,1.5);a=0.5程序如上，圖如下（2）已知兩系統(tǒng)分別用下面差分方程描述：試分別寫出它們的傳輸函數(shù)，并分別打印曲線。當方程為的程序代碼：B=1,1;A=1; subplot(2,3,1);zplane(B,A);xlabel('實部Re');ylabel('虛部Im');title('y(n)=x(n)+x(n-1)傳輸函數(shù)零、極點分布');grid onH,w=freqz(B,A,'whole');subplot(2,3,2);plot(w/pi

4、,abs(H),'linewidth',2);grid on;xlabel('omega/pi');ylabel('|H(ejomega)|');title('幅頻響應特性');axis(0,2,0,2.2);subplot(2,3,3);plot(w/pi,angle(H),'linewidth',2);grid on;xlabel('omega/pi');ylabel('phi(omega)');title('相頻響應特性');axis(-0.1,2.1,-1.6

5、,1.6);當方程為的程序代碼：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ù)零、極點分布');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)|

6、9;);title('幅頻響應特性');axis(0,2,0,2.2);subplot(2,3,6);plot(w/pi,angle(H),'linewidth',2);grid on;xlabel('omega/pi');ylabel('phi(omega)');title('相頻響應特性'); axis(-0.1,2.1,-1.6,1.6);（3）已知信號，試分析它的頻域特性，要求打印曲線。B=1,0,0,-1;A=1,-1;subplot(1,3,1);zplane(B,A);xlabel('實部Re

7、');ylabel('虛部Im');title('x(n)=R3(n)傳輸函數(shù)零、極點分布');axis(-1.1,1.1,-1.5,1.5);grid onH,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('幅頻響應特性');axis(0,2,0,1.3);subplot(1

8、,3,3);plot(w/pi,angle(H),'linewidth',2);grid on;xlabel('omega/pi');ylabel('phi(omega)');title('相頻響應特性'); axis(-0.1,2.1,-2.1,2.1);（4）假設(shè)，將以2為周期進行延拓，得到，試分析它的頻率特性，并畫出它的幅頻特性。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,'.

9、9;)xlabel('n'),ylabel('xtide(n)');title('延拓后周期為2的周期序列')B=1;A=1;subplot(2,2,2);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('幅頻響應特性');axis(0,2,0,2);subplot(2,2