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

《实验一 离散信号的频谱分析报告》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

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;%设定正弦信号频率%生成正弦信号x=sin(2*pi*f0

2、*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('

3、正弦信号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

4、([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);

5、%进行对应的频率转换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('原始信号波形')

6、;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('幅值

7、');title('频谱图N=44100');grid;