离散傅里叶变换--频移特性

从wave文件读取一段音乐信号（此时是实信号），画出其幅频特性，然后频移6KHz，此时音乐信号变为复信号，再观察其幅频特性。

filename='C:\Windows\Media\Alarm01.wav';
[s1, fs, bits]=wavread(filename);%读取音乐信号，采样频率为fs
fs=fs/1000;%采样频率转换为kHz
s=s1';
x=s(20000:22047);%取音乐信号的一部分
xf=20*log10(abs(fft(x, 4096)));%求实信号的频谱
xf1(1:2048)=xf(2049:4096);
xf1(2049:4096)=xf(1:2048);
figure;
plot((-2048:2047) * fs/4096, xf1);
grid on;

y=x.*exp(2*j*pi*6/32*(0:2047));%频谱搬移
yf=20*log10(abs(fft(y, 4096)));%求复信号的频谱
yf1(1:2048)=yf(2049:4096);
yf1(2049:4096)=yf(1:2048);
figure;
plot((-2048:2047)* fs/4096, yf1);
grid on;