离散信号的频谱分析实验报告.doc

东北大学实验题目：离散信号的频谱分析姓 名：班 级：自动化班学 号：日 期：2015.11.02离散信号的频谱分析实验报告一、实验目的1 掌握采样频率的概念2 掌握信号频谱分析方法3 掌握在计算机中绘制信号频谱图的方法二、实验内容 1、产生以下时间序列信号，并画出相应时域序列图: 采样频率为1000Hz,信号频率为30Hz的正弦信号y1（n） 。采样频率为1000Hz,信号频率为120Hz的正弦信号y2（n） 。采样频率为1000Hz， 30Hz的正弦信号和120Hz的混合信号y3（n） 。2、 分别对信号y1（n）， y2（n）和y3（n）进行FFT变换，画出其频谱图。3、自带耳麦，采用goldwave等软件录制一段语音，内容为“数字信号处理”，文件按*.wav格式存储，设置采样频率为11025Hz。4、对采集到的语音信号，进行FFT变换，画出其频谱图，并分析出自己语音的频谱范围。三、实验结果及分析1，采样频率为1000Hz,信号频率为30Hz的正弦信号频谱图2，采样频率为1000Hz,信号频率为120Hz的正弦信号频谱图3，采样频率为1000Hz， 30Hz的正弦信号和120Hz的混合信号图4，对采样频率为1000Hz,信号频率为30Hz的正弦信号进行FFT变换的频谱图5，对采样频率为1000Hz,信号频率为120Hz的正弦信号进行FFT变换的频谱图6，对采样频率为1000Hz， 30Hz的正弦信号和120Hz的混合信号进行FFT变换的频谱图7，对采集到的语音信号（录制的自己声音，内容为“数字信号处理”），进行FFT变换的频谱图。四、MATABLE程序代码 fs=1000;%采样频率为1000HzN=1024;n=0:N-1;t=n/fs;f1=30;f2=120;x1=sin(2*pi*f1*t);x2=sin(2*pi*f2*t);x3=sin(2*pi*f1*t)+sin(2*pi*f2*t)+2*randn(1,length(t);figure(1);plot(t,x1);title(origenal1);grid;figure(2);plot(t,x2);title(origenal2);grid;figure(3);plot(t,x3);title(origenal3);grid;y=fft(x1,N);%傅里叶变换mag=abs(y);f=(0:length(y)-1)*fs/length(y);figure(4);plot(f(1:N/2),mag(1:N/2);%绘制频谱图title(with noise1);grid;y=fft(x2,N);%傅里叶变换mag=abs(y);f=(0:length(y)-1)*fs/length(y);figure(5);plot(f(1:N/2),mag(1:N/2);%绘制频谱图title(with noise2);grid;y=fft(x3,N);%傅里叶变换mag=abs(y);f=(0:length(y)-1)*fs/length(y);figure(6);plot(f(1:N/2),mag(1:N/2);%绘制频谱图title(with noise3);grid;fs=11025;x1=audioread(D:new.wav);sound(x1,11025);y1=fft(x1,4096);figure(1)subplot(321);plot(x1);title(原始信号); xlabel(time n);ylabel(fuzhi n);subplot(322);plot(y1);title(原始信号频谱);