资源描述:
《实验一信号与系统的傅立叶分析.pdf》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库。
1、⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯最新资料推荐⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯实验一信号与系统的傅立叶分析一.实验目的用傅立叶变换对信号和系统进行频域分析。二.实验仪器装有matlab软件的计算机三.实验内容及步骤(1)已知系统用下面差分方程描述:y(n)x(n)ay(n1)试在a0.95和a0.5两种情况下用傅立叶变换分析系统的频率特性。要求写出jwH(e)~w系统的传输函数,并打印曲线。、当a=0.95B=1;A=[1,0.95];subplot(1,3,1);zplane(B,A);x
2、label('实部Re');ylabel('虚部Im');title('y(n)=x(n)+0.95y(n-1)传输函数零、极点分布');gridon;[H,w]=freqz(B,A,'whole');subplot(1,3,2);plot(w/pi,abs(H),'linewidth',2);gridon;xlabel('omega/pi');ylabel('
3、H(e^j^omega)
4、');title('幅频响应特性');axis([0,2,0,2.5]);subplot(1,3,3);plot(w
5、/pi,angle(H),'linewidth',2);gridon;xlabel('omega/pi');ylabel('phi(omega)');title('相频响应特性');axis([-0.1,2.1,-1.5,1.5]);a=0.5程序如上,图如下1⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯最新资料推荐⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯(2)已知两系统分别用下面差分方程描述:y1(n)x(n)x(n1)y2(n)x(n)x(n1)jwH(e)~w试分别写出它们的传输函数,并分别打印
6、曲线。当方程为y1(n)x(n)x(n1)的程序代码:B=[1,1];A=1;subplot(2,3,1);zplane(B,A);xlabel('实部Re');ylabel('虚部Im');title('y(n)=x(n)+x(n-1)传输函数零、极点分布');gridon[H,w]=freqz(B,A,'whole');subplot(2,3,2);plot(w/pi,abs(H),'linewidth',2);gridon;2⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯最新资料推荐⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯
7、⋯⋯⋯⋯⋯⋯⋯⋯⋯xlabel('omega/pi');ylabel('
8、H(e^j^omega)
9、');title('幅频响应特性');axis([0,2,0,2.2]);subplot(2,3,3);plot(w/pi,angle(H),'linewidth',2);gridon;xlabel('omega/pi');ylabel('phi(omega)');title('相频响应特性');axis([-0.1,2.1,-1.6,1.6]);当方程为y2(n)x(n)x(n1)的程序代码:B
10、=[1,-1];A=1;subplot(2,3,4);zplane(B,A);xlabel('实部Re');ylabel('虚部Im');title('y(n)=x(n)-x(n-1)传输函数零、极点分布');gridon[H,w]=freqz(B,A,'whole');subplot(2,3,5);plot(w/pi,abs(H),'linewidth',2);gridon;xlabel('omega/pi');ylabel('
11、H(e^j^omega)
12、');title('幅频响应特性');axis
13、([0,2,0,2.2]);subplot(2,3,6);plot(w/pi,angle(H),'linewidth',2);gridon;xlabel('omega/pi');ylabel('phi(omega)');title('相频响应特性');axis([-0.1,2.1,-1.6,1.6]);3⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯最新资料推荐⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯jwx(n)R3(n)X(e)~w(3)已知信号,试分析它的频域特性,要求打印曲线。B=[1,0,0,-
14、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]);gridon[H,w]=freqz(B,A,'whole');subplot(1,3,2);plot(w/pi,abs(H),'linewidth',2