% Get the number of time instants
S=size(time); L=S(2);
% Compute the numerator
numerator = Qsignal(2:L).*Isignal(1:L-1) - Isignal(2:L).*Qsignal(1:L-1);
% Compute the denominator
denominator = Isignal(2:L).*Isignal(1:L-1) + Qsignal(2:L).*Qsignal(1:L-1);
% Take the ratio of the numerator and denominator
demodulation = numerator ./ denominator;
% Plot the signal
plot(time(2:L),8*demodulation);
% Annotations
title('Frequency Demodulation');
xlabel('Time (s)');
ylabel('Amplitude (V)');