function[]=OriginalSignal();
[x,fs,bits]=wavread('bluetooth_input.wav');%播放原始信號
wavplay(x,fs);
N=length(x);%返回采樣點數(shù)
t=(1:N)/fs;
df=fs/N;%采樣間隔
n1=1:N/2;
f=(n1-1)*df;%頻帶寬度
figure(2);
plot(x,'LineWidth',2),grid on;%信號的時域波形
title('原始信號的時域波形');
xlabel('時間/t');
ylabel('幅值/A');
y0=fft(x);%快速傅立葉變換
figure(3);
plot(f,20*log10(abs(y0(n1)))),grid on;%離散信號的頻譜圖
title('原始信號的頻譜圖');
xlabel('頻率/Hz');
ylabel('幅值/db');
grid;