1、实验一 离散信号的频谱分析报告班级 姓名 学号实验一 离散信号的频谱分析报告1 掌握采样频率的概念2 掌握信号频谱分析方法；3 掌握在计算机中绘制信号频谱图的方法。采样频率为 1000Hz,信号频率为 30Hz 的正弦信号y1（n） 对其进行 FFT 变换采样频率为 1000Hz,信号频率为 120Hz 的正弦信号y2（n）对其进行 FFT 变换采样频率为 1000Hz， 30Hz 的正弦信号和 120Hz 的混合信号 y3（n）。 对其进行 FFT 变换语音信号波形附录程序：fs=1000;%设定采样频率N=1024;n=0:N-1;t=n/fs;f0=30;%设定正弦信号频率%生成正弦信号
2、x=sin(2*pi*f0*t);figure(1);subplot(3,2,1);plot(t,x);%作正弦信号的时域波形xlabel(t);ylabel(y);title(正弦信号 30HZ 时域波形);grid;%进行 FFT 变换并做频谱图y=fft(x,N);%进行 fft 变换mag=abs(y);%求幅值f=(0:length(y)-1)*fs/length(y);%进行对应的频率转换figure(1);subplot(3,2,2);plot(f,mag);%做频谱图axis(0,100,0,500);xlabel(频率(Hz);ylabel(幅值);title(正弦信号 30
3、HZ 幅频谱图 N=1024);grid;%120HZf1=120;x=sin(2*pi*f1*t);figure(1);subplot(3,2,3);plot(t,x);%作正弦信号的时域波形xlabel(t);ylabel(y);title(正弦信号 120HZ 时域波形);grid;%进行 FFT 变换并做频谱图y=fft(x,N);%进行 fft 变换mag=abs(y);%求幅值f=(0:length(y)-1)*fs/length(y);%进行对应的频率转换figure(1);subplot(3,2,4);plot(f,mag);%做频谱图axis(0,200,0,600);xla
4、bel(频率(Hz);ylabel(幅值);title(正弦信号 120HZ 幅频谱图 N=1024);grid;%混合x=sin(2*pi*f0*t)+sin(2*pi*f1*t);figure(1);subplot(3,2,5);plot(t,x);%作正弦信号的时域波形xlabel(t);ylabel(y);title(正弦信号混合时域波形);grid;%进行 FFT 变换并做频谱图y=fft(x,N);%进行 fft 变换mag=abs(y);%求幅值f=(0:length(y)-1)*fs/length(y);%进行对应的频率转换figure(1);subplot(3,2,6);pl
5、ot(f,mag);%做频谱图axis(0,200,0,600);xlabel(频率(Hz);ylabel(幅值);title(正弦信号混合幅频谱图 N=1024);grid;fs=11025;y,fs,bits=wavread(C:UsersAdministratorDesktop20151030133647.WAV);sound(y,fs,bits);Y=fft(y); figure(2);subplot(2,3,1);plot(y);title(原始信号波形);subplot(2,3,2);plot(Y); title(原始信号频谱);subplot(2,3,3);plot(abs(Y);title(原始信号幅值);subplot(234);plot(angle(Y);title(原始信号相位); c=fft(y,44100);%进行 fft 变换mag=abs(c);%求幅值f=(0:length(c)-1)*fs/length(c);%进行对应的频率转换subplot(2,3,5);plot(f,mag);%做频谱图axis(0,40000,0,100);xlabel(频率(Hz);ylabel(幅值);title(频谱图 N=44100);grid;
