Неполучается построить 3D график

clear;
clc;
f0=30*10^6;%частота сигнала
f1=2*10^6;%частота следования бит
f=60*f0;%частота дискретизации
N = 1024; %колличество бит
tau=20*10^-6;
fd=5000;
w0 = 2*pi*f0;
In = randint(1,N+1); % информационный сигнал
T = 512*10^-6/N; % интервал времени, определяющий скачек фазы
 
t = (0:1/f:N/f1);
F = In( fix(t./T)+1); % расширяем инф. сигнал до кол-ва точек времени
S = sin(2*pi*f0*t + F*pi); % если F=1 - скачек фазы на пи.
plot(t,S,'r',t,F,':b','LineWidth',2);
grid on
axis ([min(t) max(t) -2 2])
xlabel('Время, с')
ylabel('Сигнал')
legend('ФМн сигнал','Кодируюущая последовательность')
t1 = (tau:1/f:N/f1);
F1 = In( fix(t1./T)+1); % расширяем инф. сигнал до кол-ва точек времени
 
S1 = sin(2*pi*(f0+fd)*t1 + F1*pi); % если F=1 - скачек фазы на пи.
 
S2=awgn(S1,100,'measured',100);
 
 
k=abs(xcorr(S,S2));
k1=k/max(k);
lags=-(N/f1):1/f:(N/f1);
lags1=0:10000/1843200:10000;
figure;
plot3(lags,lags1,k1);
grid on
q=[-20 -10 0 10 20 30 40 100];
to=[0.08e-6 0.06e-6 0.05e-6 0.05e-6 0.04e-6 0.04e-6 0.03e-6 0.01e-6];
fo=[196 196 196 196 196 196 196 196];
figure
plot(q,to);
axis ([-20 110 0 0.09e-6])
xlabel('Сигнал/шум, дБ')
ylabel('Ошибка, с')
legend('Ошибка по времени')
grid on
figure
plot(q,fo);
axis([-20 100 0 250])
xlabel('Сигнал/шум, дБ')
ylabel('Ошибка, Гц')
legend('Ошибка по частоте')
grid on