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

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

班级

姓名

学号

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

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;

%120HZ

f1=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:\Users\Administrator\Desktop\2015103 0133647.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;