时域转频域

[sample, fs] = audioread("C:\Users\dev\Documents\录音\录音.wav");
[sample_count, n] = size(sample);
fprintf('sample_count: %d, fs: %d',sample_count,fs);
if (n == 2)
	sample = sample(:, 1);      % 双通道变单通道
end
freq_signal = fft(sample, sample_count);
z1 = abs(freq_signal);
% 时间
f1 = (0:sample_count - 1) * fs / sample_count;
% 除2是因为fft之后得到的是对称的双边谱
z = z1(1 : sample_count / 2) * 2 / sample_count;
plot(f1(1 : sample_count / 2), z);
xlabel('频率（hz）');
ylabel('振幅');