Амплитудный и фазовый спектры аудио сигнала

[x,Fs,bits]=wavread('c:\1.wav');
fvtool(x)

N=2048
Y=wavread('C:\1.wav', N)
T=0:1/(N-1):1
M=length(T)
A=fftshift(fft(Y))
plot(T, abs(A)) %амплитудный спектр
plot(T, angle(A)) %фазовый спектр

T = 5;
tay=0.01;
t = 0 : tay : T;
f0=1/T;
x=wavread('C:\1.wav');
max=1/2*tay;
N=length(x);
fmax=-max:abs((2*max)/(N-1)):max;
A=fftshift(fft(x));
        
plot(fmax, abs(A)); %амплитудный спектр
plot(fmax, angle(A)); %фазовый спектр
grid on;

FftS=abs(fft(x,fs));% Амплитуды преобразования Фурье сигнала
FftS=2*FftS./fs;% Нормировка спектра по амплитуде
F=0:1:fs/2-1/fs;% Массив частот вычисляемого спектра Фурье
figure;
plot(F,FftS(1:length(F)));% Построение спектра Фурье сигнала
grid on;

fs = 1000;
t = 0:1/fs:1;
x = cos(2*pi*100.*t);
 
FftS=abs(fft(x,fs));% Амплитуды преобразования Фурье сигнала
FftS=2*FftS./fs;% Нормировка спектра по амплитуде
F=0:1:fs/2-1/fs;% Массив частот вычисляемого спектра Фурье
figure;
subplot(2,2,1)
plot(F,FftS(1:length(F)));% Построение спектра Фурье сигнала
grid on;
 
%% Другое основание
fs = 500;
FftS=abs(fft(x,fs));% Амплитуды преобразования Фурье сигнала
FftS=2*FftS./fs;% Нормировка спектра по амплитуде
F=0:1:fs/2-1/fs;% Массив частот вычисляемого спектра Фурье
subplot(2,2,2)
plot(F,FftS(1:length(F)));% Построение спектра Фурье сигнала
grid on;
 
%% Меньше длительность, старое основание
fs = 1000;
t = 0:1/fs:0.5;
x = cos(2*pi*100.*t);
 
FftS=abs(fft(x,fs));% Амплитуды преобразования Фурье сигнала
FftS=2*FftS./fs;% Нормировка спектра по амплитуде
F=0:1:fs/2-1/fs;% Массив частот вычисляемого спектра Фурье
subplot(2,2,3)
plot(F,FftS(1:length(F)));% Построение спектра Фурье сигнала
grid on;
 
%% Ещё меньше длительность, старое основание
t = 0:1/fs:0.1;
x = cos(2*pi*100.*t);
 
FftS=abs(fft(x,fs));% Амплитуды преобразования Фурье сигнала
FftS=2*FftS./fs;% Нормировка спектра по амплитуде
F=0:1:fs/2-1/fs;% Массив частот вычисляемого спектра Фурье
subplot(2,2,4)
plot(F,FftS(1:length(F)));% Построение спектра Фурье сигнала
grid on;

[x,Fs]=audioread('1.wav');
L = length(x);
 
n = 2^nextpow2(L);
Y = fft(x,n);
f = Fs*(0:(n/2))/n;
P = abs(Y/n);
 
plot(f,P(1:n/2+1))
xlabel('Frequency (f)')
ylabel('|P(f)|')

clear; clc;
load('train');
 
FftS=abs(fft(y,Fs));% Амплитуды преобразования Фурье сигнала
FftS=2*FftS./Fs;% Нормировка спектра по амплитуде
F=0:1:Fs/2-1/Fs;% Массив частот вычисляемого спектра Фурье
figure;
subplot(2,1,1)
plot(F,FftS(1:length(F)));% Построение спектра Фурье сигнала
grid on;
 
L = length(y); 
n = 2^nextpow2(L);
Y = fft(y,n);
f = Fs*(0:(n/2))/n;
P = abs(Y/n);
subplot(2,1,2)
plot(f,P(1:n/2+1))
xlabel('Frequency (f)')
ylabel('|P(f)|')
grid on;