实验五信号频谱分析II.docVIP

实验五 信号的频谱分析II 实验目的： 1、理解频谱分析时N的选取，熟练掌握横轴（频率）坐标的对应关系。 2、学会计算实际信号的幅度谱，并分析简单的特性。 一、基本练习 1、对余弦信号做傅里叶变换，改变频率，考察幅度谱。 M=100; Dt=1/M; t=0:Dt:1-Dt; %设1s时间长度，步长Dt=0.01s，t为M个元素的向量 x1=2*cos(2*pi*2*t); %频率为2Hz的余弦信号，M个点 x2=2*cos(2*pi*5*t); %频率为5Hz的余弦信号，M个点 x3=2*cos(2*pi*10*t); %频率为10Hz的余弦信号，M个点 x4=2*cos(2*pi*2.5*t); %频率为2.5Hz的余弦信号，M个点 x5=2*cos(2*pi*4.5*t); %频率为4.5Hz的余弦信号，M个点 subplot(5,2,1); plot(t,x1); title(频率2Hz); subplot(5,2,3); plot(t,x2); title(频率5Hz); subplot(5,2,5); plot(t,x3); title(频率10Hz); subplot(5,2,7); plot(t,x4); title(频率2.5Hz); subplot(5,2,9); plot(t,x5); title(频率4.5Hz); N=100; %设频谱长度1000点， f=(-N/2:N/2-1)/N*(1/Dt); %步长1/N*(1/Dt)=0.1Hz X1=fftshift(fft(x1,N)); AMP1=abs(X1); X2=fftshift(fft(x2,N)); AMP2=abs(X2); X3=fftshift(fft(x3,N)); AMP3=abs(X3); X4=fftshift(fft(x4,N)); AMP4=abs(X4); X5=fftshift(fft(x5,N)); AMP5=abs(X5); subplot(5,2,2); stem(f,AMP1); title(频率2Hz); axis([-50 50 0 100 ]); subplot(5,2,4); stem(f,AMP2); title(频率5Hz); axis([-50 50 0 100 ]); subplot(5,2,6); stem(f,AMP3); title(频率10Hz); axis([-50 50 0 100 ]); subplot(5,2,8); stem(f,AMP4); title(频率2.5Hz); axis([-50 50 0 100 ]); subplot(5,2,10); stem(f,AMP5); title(频率4.5Hz); axis([-50 50 0 100 ]); 将上面的程序运行后，发现频率2.5Hz和频率4.5Hz的频谱与其他相差很大。将stem改成plot，再将N改为1000，此时各余弦信号除频率不同外，频谱基本相似。因此在频谱分析时，如果允许，最好在频域的点数大一些，其中的道理随着后续课程会逐步理解。 2、频率轴的问题 前面程序中，f=(-N/2:N/2-1)/N*(1/Dt)，其中： （1）如果每两个点之间的频率差为1/N，(-N/2:N/2-1)/N是指横轴的频率范围是-0.5~0.5（一般范围的上限要减去一个间隔，即小学学过的植树问题）。这是归一化的频率。 （2）如果需要使频率范围在-0.5 fs~0.5 fs之间，则f=(-N/2:N/2-1)/N*fs。 （3）fs实际是时域波形两点间的间隔Dt之倒数，及fs=1/Dt。以后我们将学到fs称为采样频率，Dt称为采样间隔。 3、音频信号的读写和播放 [x b]=wavread(music1.wav); %读音频信号文件，x信号为，b为采样频率，b=11025。 wavplay(x ,b); %播放音频信号 x1=resample(x,a1,a2) ;%时域压缩或扩展，其中a1和a2为正整数。 wavwrite(x,b,ss.wav); %写音频信号文件 如果相对某个信号（或频谱）归一化，y=y/max(y)，最大值为1。下面求带宽用到。 二、实验报告内容（一个必做题，一个选做题） 1、对所给音频信号，进行时域压缩或扩展，画出时域波形与幅度谱（如下图），使其满足以下要求，并说明做了什么变换。 （1）将music1.wav的音调变低a倍（0.8a1），及变高a倍（1a1.3），比较变换前后的10%