clearclcclose?all?%?清屏
N=176400;?%?采樣點數?一共4秒，每秒44100個點，采集176400個點
n=[0:N-1];
[x1fs1nbits]=wavread（‘D:\bin\au.wav‘）;%?采語音信號
%sound（x1fs1nbits）;?%?播放歌曲??用來試聽
x=x1（1:N）;
Ts=1/fs1;?%?抽樣時間間隔
X=fft（x）;?%?傅里葉變換
f=n/N*fs1;?%?把點數轉換成頻率
b=abs（X）;
subplot（211）;
plot（fb）;
title（‘?輸入語音的幅度譜?‘）;
xlabel（‘?頻率為?/Hz‘）;
%翻轉頻譜
f1=22050;???%中點就是44100的一半
y=x.*exp（1i*2*pi*f1.*n*Ts）;?%?進行翻轉
Y=fft（y）;
subplot（212）;
plot（fabs（Y））;?%?畫出翻轉后的幅度圖
title（‘?頻譜倒置語音的幅度譜?‘）;
xlabel（‘?