实验一离散信号地频谱分析报告报告材料.doc
实验一 离散信号的频谱分析报告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(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);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);xlabel('频率(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);plot(f,mag);%做频谱图axis(0,200,0,600);xlabel('频率(Hz)');ylabel('幅值');title('正弦信号混合幅频谱图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;
