Need Help with Telephone Dialing Signal

So I'm suppose to take the file I have. Read it, play it and then plot the magnitude and Phase of the DFT signal versus the frequency in Hz. I think I got the first part right.

The problem I have is that I also need to include vertical dotted lines at each of the DTMF column and row frequencies. Then I need to decode the number dialed from the DFt. I did find this resource that seemed to help but I couldn't find the actual code. Any help would be greatly appreciated.

clear; clc;
[data,Fs] = audioread('tel.wav'); % Read the audio file and store the audio in data and sampling rate into fs
player = audioplayer(data,Fs);
play(player);

    NFFT = length(data);
  Y = fft(data,NFFT);
  F = ((0:1/NFFT:1-1/NFFT)*Fs).';
  magnitudeY = abs(Y);        % Magnitude of the FFT
  phaseY = unwrap(angle(Y));  % Phase of the FFT
  tone= abs(data)';

figure
subplot(2,1,1)       
plot(F,magnitudeY)
title('Magnitude')
line([1209 697], [0 500],'LineStyle',':'); % dotted line at DTMF?

subplot(2,1,2)       
plot(F,phaseY)       
title('Phase')

[Walter Roberson]

https://www.mathworks.com/matlabcentral/answers/?term=tag%3A%22dtmf%22

Answer 1

[Hannah Chamberlain]

I tried searching and couldn't figure it out. I'm sorry. I'm so lost.

Answer 2

[Star Strider]

Have fun with How can i detect frequency of dtmf tone (.wav) file.

[Hannah Chamberlain]

How would I load my sound file into it? I have it as a .wav file. How would I change the type? Would the code I haven't work?

[Star Strider]

Use the audioread function. It will give you the signal (in my code, ‘x’) and the sampling frequency (in my code ‘fs’). See the documentation on audioread for details on how to use it.