资源描述：
《信号的频谱分析实验报告.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
