实验一信号与系统的傅立叶分析报告

1、实用标准文案实验一信号与系统的傅立叶分析一 .实验目的用傅立叶变换对信号和系统进行频域分析。二 .实验仪器装有 matlab 软件的计算机三 .实验内容及步骤（1）已知系统用下面差分方程描述：y(n) x( n) ay (n 1)试在 a0.95 和 a0.5两种情况下用傅立叶变换分析系统的频率特性。要求写出系统的传输函数，并打印H (e jw ) w 曲线。、当 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) 传输函数零、极点分布');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(' 幅频响应特性');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(' 相频响应特性');axis(-0.1,2.1,-1.5,1.5);a=0.5 程序如上，图如下精彩文档实用标准文案（ 2）已知两系统分别用下面差分方程描述：y1 ( n)x(n)x(n1)y2 (n)x(n)x(n1)H (e jw ) w试分别写出它们的传输函数，并分别打印曲线。当方程为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)传输函数零、极点分布');grid onH,w=freqz(B,A,'whole');subplot(2,3,2);plot(w/pi,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,

5、angle(H),'linewidth',2);grid on;xlabel('omega/pi');ylabel('phi(omega)');title(' 相频响应特性');axis(-0.1,2.1,-1.6,1.6);当方程为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)传输函数零、极点分布&#

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(' 幅频响应特性');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(' 相频响应特性');axis(-0.1,2.1,-1.6,1.6);精彩文档实用标准文案（ 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)传输函数零、极点分布');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(' 幅频响应特性');axis(0,2,0,1.3);精彩文档实用标准文案subplot(1,3,3);plot(w/pi,angle(H),'linewidth',2);grid on;xlabel('omega/pi');ylabel(&

9、#39;phi(omega)');title(' 相频响应特性');axis(-0.1,2.1,-2.1,2.1);（ 4）假设 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);精彩文档实用标准文案zplane(B,A);xlabel('实部 Re');ylabel('虚部 Im');title('x(n)延拓后得到x (n)传输函数零、极点分布');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);subp