matlab - 周波数の大きさと位相をプロットする方法は?

Question

私は実用的な周波数ドメイン分析の紹介を読んでいました。次のコードを使用して、信号の周波数成分を観察する方法を示します。

Fs = 44100;
y = audioread('guitartune.wav');

NFFT = length(y);
Y = fft(y,NFFT);
F = ((0:1/NFFT:1-1/NFFT)*Fs).';

magnitudeY = abs(Y);        % Magnitude of the FFT
phaseY = unwrap(angle(Y));  % Phase of the FFT

helperFrequencyAnalysisPlot1(F,magnitudeY,phaseY,NFFT)

プロットする関数はhelperFrequencyAnalysisPlot1. 関数がどのようにプロットを行うかを知る必要がありますが、記事にはその方法が示されていません。それは実際にどのように行われますか？