资源描述:
《信号的频谱分析实验报告.doc》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库。
1、实验四信号的频谱分析一.实验目的1.掌握利用FFT分析连续周期,非周期信号的频谱,如周期,非周期方波,正弦信号等。理解CFS,CTFT与DFT(FFT)的关系。2.利用FFT分析离散周期,非周期信号的频谱,如周期,非周期方波,正弦信号等。理解DFS,DTFT与DFT(FFT)的关系,并讨论连续信号与离散信号频谱分析方法的异同。二.实验要求1.编写程序完成任意信号数字谱分析算法;2.编写实验报告。三.实验内容1.利用FFT,分析并画出频谱,改变采样间隔与截断长度,分析混叠与泄漏对单一频率成分信号频谱的影
2、响。(1)sin(100*pi*t)产生程序:closeall;clc;clear;t=0:0.0025:0.5-0.0025;f=400*t;w0=100*pi;y=sin(w0*t);a=fft(y);b=abs(a)/200;d=angle(a)*180/pi;subplot(311);plot(t,y);title('y=sin(wt)');xlabel('t');ylabel('y(t)');subplot(312);stem(f,b);title('振幅');xlabel('f');yla
3、bel('y(t)');subplot(313);stem(f,d);title('相位');xlabel('t');ylabel('y(t)');混叠closeall;clc;clear;t=0:0.0115:0.46-0.0115;f=(t/0.0115)*2;w0=100*pi;y=sin(w0*t);a=fft(y);b=abs(a)/40;d=angle(a)*180/pi;subplot(311);plot(t,y);title('y=sin(wt)');xlabel('t');ylabe
4、l('y(t)');subplot(312);stem(f,b);title('振幅');xlabel('f');ylabel('y(t)');subplot(313);stem(f,d);title('相位');xlabel('t');ylabel('y(t)');泄漏closeall;clc;clear;t=0:0.0025:0.5-0.0075;f=800*t;w0=100*pi;y=sin(w0*t);a=fft(y);b=abs(a)/198;d=angle(a)*180/pi;subplo
5、t(311);plot(t,y);title('y=sin(wt)');xlabel('t');ylabel('y(t)');subplot(312);stem(f,b);title('振幅');xlabel('f');ylabel('y(t)');subplot(313);stem(f,d);title('相位');xlabel('t');ylabel('y(t)');(2)cos(100*pi*t);closeall;clc;clear;t=0:0.0025:0.5-0.0025;f=800*t;
6、w0=100*pi;y=cos(w0*t);a=fft(y);b=abs(a)/200;d=angle(a)*180/pi;subplot(311);plot(t,y);title('y=cos(wt)');xlabel('t');ylabel('y(t)');gridon;holdon;subplot(312);stem(f,b);title('振幅');xlabel('f');ylabel('y(t)');gridon;holdon;subplot(313);stem(f,d);title('相位
7、');xlabel('f');ylabel('y(t)');混叠closeall;clc;clear;t=0:0.0115:0.46-0.0115;f=(t/0.0115)*2;w0=100*pi;y=cos(w0*t);a=fft(y);b=abs(a)/40;d=angle(a)*180/pi;subplot(311);plot(t,y);title('y=cos(wt)');xlabel('t');ylabel('y(t)');subplot(312);stem(f,b);title('振幅')
8、;xlabel('f');ylabel('y(t)');subplot(313);stem(f,d);title('相位');xlabel('t');ylabel('y(t)');泄漏closeall;clc;clear;t=0:0.0025:0.5-0.0075;f=800*t;w0=100*pi;y=cos(w0*t);a=fft(y);b=abs(a)/198;d=angle(a)*180/pi;subplot(311);plot(t,y);tit