《实验一离散信号地频谱分析报告报告材料.doc》由会员分享,可在线阅读,更多相关《实验一离散信号地频谱分析报告报告材料.doc(7页珍藏版)》请在课桌文档上搜索。
1、实验一 离散信号的频谱分析报告1 掌握采样频率的概念2 掌握信号频谱分析方法;3 掌握在计算机中绘制信号频谱图的方法。采样频率为1000Hz,信号频率为30Hz的正弦信号y1n 对其进展FFT变换采样频率为1000Hz,信号频率为120Hz的正弦信号y2n对其进展FFT变换采样频率为1000Hz, 30Hz的正弦信号和120Hz的混合信号y3n。对其进展FFT变换语音信号波形附录程序:fs=1000;%设定采样频率N=1024;n=0:N-1;t=n/fs;f0=30;%设定正弦信号频率%生成正弦信号x=sin(2*pi*f0*t);figure(1);subplot(3,2,1);plot(
2、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(正弦信号30HZ幅频谱图N=1024);grid;%120HZf1=120;x=sin(2*pi*f1*t);figur
3、e(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);xlabel(频率(Hz);ylabel(幅值);title(正弦信号120HZ幅频谱图N=1024);grid;%混合x=sin
4、(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);plot(f,mag);%做频谱图axis(0,200,0,600);xlabel(频率(Hz);ylabel(幅值);title(正弦信号混
5、合幅频谱图N=1024);grid; fs=11025; 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;