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

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:\\Users\\Administrator\\Desktop\\20151030133647.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); plot(f,mag);%做频谱图 axis([0,40000,0,100]); xlabel('频率(Hz)'); ylabel('幅值');

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