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

实验一 离散信号的频谱分析报告班级 姓名 学号实验一 离散信号的频谱分析报告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;%设定正弦信号频率%生成正弦信号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; 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);